Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

High-resolution simulation of stack pollutant emissions

[EN]A three-dimensional air pollution model for the short-term simulation of emission, transport and reaction of pollutants is presented. In the finite element simulation of these environmental processes over a complex terrain, a mesh generator capable of adapting itself to the topographic characteristics is essential, A local refinement of tetrahedra is used in order to capture the plume rise. Then a wind field is computed by using a mass-consistent model and perturbing its vertical component to introduce the plume rise effect. Finally, an Eulerian convection-diffusionreaction model is used to simulate the pollutant dispersion…

Simulación de la calidad del aire en la isla de Gran Canaria mediante el método de los elementos finitos y su validación con datos experimentales

[EN]Este estudio muestra una simulación 3D de la calidad del aire en la isla de Gran Canaria. Se simula el transporte difusión y reacción de los gases emitidos en dos centrales térmicas de la isla. Se presenta una estrategia de modelización y simulación mediante elementos finitos. Los resultados se comparan con datos experimentales. Se usa el método del Meccano para mallar la isla de Gran Canaria, a partir de una triangulación de la topografia de la isla. A continuación, con los datos de las estaciones de viento, se calcula un campo de viento para la zona de estudio…

Modelización de sólidos mediante estructura Octree

[EN]En este trabajo presentamos un procedimiento para la modelización de sólidos con geometríacompleja mediante una estructura octree [1]. Partiendo de una triangulación superficial del sólido, se construye una malla de hexaedros que aproxima la geometría con una precisión indicada. Para ello se emplea una división recursiva del espacio con una estructura octree. Inicialmente se lleva a cabo una traslación y escalado de la triangulación al cubo inicial [0,1]3 donde se construye el octree...

Non linear application of an in-house finite elements software

Questa tesi si pone come obiettivo l'analisi delle componenti di sollecitazione statica di un serbatoio, in acciaio API 5L X52, sottoposto a carichi di flessione e pressione interna attraverso il programma agli elementi finiti PLCd4, sviluppato presso l'International Center for Numerical Methods in Engineering (CIMNE - Barcelona). Questo tipo di analisi rientra nel progetto europeo ULCF, il cui traguardo è lo studio della fatica a bassissimo numero di cicli per strutture in acciaio. Prima di osservare la struttura completa del serbatoio è stato studiato il comportamento del materiale per implementare all'interno del programma una nuova tipologia di curva che rappresentasse al meglio l'andamento delle tensioni interne. Attraverso il lavoro di preparazione alla tesi è stato inserito all'interno del programma un algoritmo per la distribuzione delle pressioni superficiali sui corpi 3D, successivamente utilizzato per l'analisi della pressione interna nel serbatoio. Sono state effettuate analisi FEM del serbatoio in diverse configurazioni di carico ove si è cercato di modellare al meglio la struttura portante relativa al caso reale di "full scale test". Dal punto di vista analitico i risultati ottenuti sono soddisfacenti in quanto rispecchiano un corretto comportamento del serbatoio in condizioni di pressioni molto elevate e confermano la bontà del programma nell'analisi computazionale.

Resonances in the two centers Coulomb system

In this work we investigate the existence of resonances for two-centers Coulomb systems with arbitrary charges in two and three dimensions, defining them in terms of generalized complex eigenvalues of a non-selfadjoint deformation of the two-center Schrödinger operator. After giving a description of the bifurcation of the classical system for positive energies, we construct the resolvent kernel of the operators and we prove that they can be extended analytically to the second Riemann sheet. The resonances are then defined and studied with numerical methods and perturbation theory.

Automotive diesel engine transient operation: modeling, optimization and control

Traditionally, the study of internal combustion engines operation has focused on the steady-state performance. However, the daily driving schedule of automotive engines is inherently related to unsteady conditions. There are various operating conditions experienced by (diesel) engines that can be classified as transient. Besides the variation of the engine operating point, in terms of engine speed and torque, also the warm up phase can be considered as a transient condition. Chapter 2 has to do with this thermal transient condition; more precisely the main issue is the performance of a Selective Catalytic Reduction (SCR) system during cold start and warm up phases of the engine. The proposal of the underlying work is to investigate and identify optimal exhaust line heating strategies, to provide a fast activation of the catalytic reactions on SCR. Chapters 3 and 4 focus the attention on the dynamic behavior of the engine, when considering typical driving conditions. The common approach to dynamic optimization involves the solution of a single optimal-control problem. However, this approach requires the availability of models that are valid throughout the whole engine operating range and actuator ranges. In addition, the result of the optimization is meaningful only if the model is very accurate. Chapter 3 proposes a methodology to circumvent those demanding requirements: an iteration between transient measurements to refine a purpose-built model and a dynamic optimization which is constrained to the model validity region. Moreover all numerical methods required to implement this procedure are presented. Chapter 4 proposes an approach to derive a transient feedforward control system in an automated way. It relies on optimal control theory to solve a dynamic optimization problem for fast transients. From the optimal solutions, the relevant information is extracted and stored in maps spanned by the engine speed and the torque gradient.

Periodicities in the Hodgkin-Huxley model and versions of this model with stochastic input

A field of computational neuroscience develops mathematical models to describe neuronal systems. The aim is to better understand the nervous system. Historically, the integrate-and-fire model, developed by Lapique in 1907, was the first model describing a neuron. In 1952 Hodgkin and Huxley [8] described the so called Hodgkin-Huxley model in the article “A Quantitative Description of Membrane Current and Its Application to Conduction and Excitation in Nerve”. The Hodgkin-Huxley model is one of the most successful and widely-used biological neuron models. Based on experimental data from the squid giant axon, Hodgkin and Huxley developed their mathematical model as a four-dimensional system of first-order ordinary differential equations. One of these equations characterizes the membrane potential as a process in time, whereas the other three equations depict the opening and closing state of sodium and potassium ion channels. The membrane potential is proportional to the sum of ionic current flowing across the membrane and an externally applied current. For various types of external input the membrane potential behaves differently. This thesis considers the following three types of input: (i) Rinzel and Miller [15] calculated an interval of amplitudes for a constant applied current, where the membrane potential is repetitively spiking; (ii) Aihara, Matsumoto and Ikegaya [1] said that dependent on the amplitude and the frequency of a periodic applied current the membrane potential responds periodically; (iii) Izhikevich [12] stated that brief pulses of positive and negative current with different amplitudes and frequencies can lead to a periodic response of the membrane potential. In chapter 1 the Hodgkin-Huxley model is introduced according to Izhikevich [12]. Besides the definition of the model, several biological and physiological notes are made, and further concepts are described by examples. Moreover, the numerical methods to solve the equations of the Hodgkin-Huxley model are presented which were used for the computer simulations in chapter 2 and chapter 3. In chapter 2 the statements for the three different inputs (i), (ii) and (iii) will be verified, and periodic behavior for the inputs (ii) and (iii) will be investigated. In chapter 3 the inputs are embedded in an Ornstein-Uhlenbeck process to see the influence of noise on the results of chapter 2.

Polygonal cracks in the seasonal semi-translucent CO2 ice layer in Martian polar areas

In this paper, we use morphological and numerical methods to test the hypothesis that seasonally formed fracture patterns in the Martian polar regions result from the brittle failure of seasonal CO2 slab ice. The observations by the High Resolution Imaging Science Experiment (HiRISE) of polar regions of Mars show very narrow dark elongated linear patterns that are observed during some periods of time in spring, disappear in summer and re-appear again in the following spring. They are repeatedly formed in the same areas but they do not repeat the exact pattern from year to year. This leads to the conclusion that they are cracks formed in the seasonal ice layer. Some of models of seasonal surface processes rely on the existence of a transparent form of CO2 ice, so-called slab ice. For the creation of the observed cracks the ice is required to be a continuous media, not an agglomeration of relatively separate particles like a firn. The best explanation for our observations is a slab ice with relatively high transparency in the visible wavelength range. This transparency allows a solid state green-house effect to act underneath the ice sheet raising the pressure by sublimation from below. The trapped gas creates overpressure and the ice sheet breaks at some point creating the observed cracks. We show that the times when the cracks appear are in agreement with the model calculation, providing one more piece of evidence that CO2 slab ice covers polar areas in spring.

Design and implementation of a 3D computer game controller using inertial MEMS sensors

Though 3D computer graphics has seen tremendous advancement in the past two decades, most available mechanisms for computer interaction in 3D are high cost and targeted for industry and virtual reality applications. Recent advances in Micro-Electro-Mechanical-System (MEMS) devices have brought forth a variety of new low-cost, low-power, miniature sensors with high accuracy, which are well suited for hand-held devices. In this work a novel design for a 3D computer game controller using inertial sensors is proposed, and a prototype device based on this design is implemented. The design incorporates MEMS accelerometers and gyroscopes from Analog Devices to measure the three components of the acceleration and angular velocity. From these sensor readings, the position and orientation of the hand-held compartment can be calculated using numerical methods. The implemented prototype is utilizes a USB 2.0 compliant interface for power and communication with the host system. A Microchip dsPIC microcontroller is used in the design. This microcontroller integrates the analog to digital converters, the program memory flash, as well as the core processor, on a single integrated circuit. A PC running Microsoft Windows operating system is used as the host machine. Prototype firmware for the microcontroller is developed and tested to establish the communication between the design and the host, and perform the data acquisition and initial filtering of the sensor data. A PC front-end application with a graphical interface is developed to communicate with the device, and allow real-time visualization of the acquired data.

Focal laminar cortical infarcts following aneurysmal subarachnoid haemorrhage

INTRODUCTION: The aim of this prospective study was to analyse small band-like cortical infarcts after subarachnoid haemorrhage (SAH) using magnetic resonance imaging (MRI) with reference to additional digital subtraction angiography (DSA). METHODS: In a 5-year period between January 2002 and January 2007 10 out of 188 patients with aneurysmal SAH were evaluated (one patient Hunt and Hess grade I, one patient grade II, four patients grade III, two patients grade IV, and two patients grade V). The imaging protocol included serially performed MRI with diffusion- and perfusion-weighted images (DWI/PWI) at three time points after aneurysm treatment, and cerebral vasospasm (CVS) was analysed on follow-up DSA on day 7+/-3 after SAH. RESULTS: The lesions were located in the frontal lobe (n=10), in the insular cortex (n=3) and in the parietal lobe (n=1). The band-like infarcts occurred after a mean time interval of 5.8 days (range 3-10 days) and showed unexceptional adjacent thick sulcal clots. Seven out of ten patients with cortical infarcts had no or mild CVS, and in the remaining three patients DSA disclosed moderate (n=2) or severe (n=1) CVS. CONCLUSION: The infarct pattern after aneurysmal SAH includes cortical band-like lesions. In contrast to territorial infarcts or lacunar infarcts in the white matter which develop as a result of moderate or severe proximal and/or distal vasospasm visible on angiography, the cortical band-like lesions adjacent to sulcal clots may also develop without evidence of macroscopic vasospasm, implying a vasospastic reaction of the most distal superficial and intraparenchymal vessels.

Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Modeling asteroid collisions and impact processes

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.

Pollen records and lithology of profiles from Lobsigensee, Switzerland

Lobsigensee is a small kettle hole lake 15 km north-west of Bern on the Swiss Plateau, at an altitude of 514 m asl. Its surface is 2ha today, its maximum depth 2.7 m; it has no inlet and the overflow functions mainly during snow melting. The area was covered by Rhone ice during the Last Glaciation (map in Fig.2). Local geology, climate and vegetation are summarized in Figure 3A-C, the history of settlement in Figures 5-7. In order to reconstruct the vegetational and environmental history of the lake and its surroundings pollen analysis and other bio- and isotope stratigraphies were applied to twelve profiles cored across the basin with modified Livingstone corers (Fig.3 D). (1) The standard diagram: The central core LQ-90 is described as the standard pollen diagram (Chapter 3) with 10 local pollen assemblage zones of the Late-Glacial (local PAZ Ll to Ll0, from about 16'000(7) to 10'000 years BP) and 20 PAZ of the Holocene (local PAZ L11 to L30), see Figs. 8-10 and 20-24. Local PAZ L 1 to L3 are in the Late-Glacial clay and record the vegetational development after the ice retreat: L1 shows very low pollen concentration and high Pinus percentages due to long-distance transport and reworking; the latter mechanism is corroborated by the findings of thermophilous and pre-Quaternary taxa. Local PAZ L2 has a high di versi ty of non-arboreal pollen (NAP) and reflects the Late-Glacial steppe rich in heliophilous species. Local PAZ L3 is similar but additionally rich in Betula nana and Sal1x, thus reflecting a "shrub tundra". The PAZ L1 to L3 belong to the Oldest Dryas biozone. Local PAZ L4 to L 10 are found in the gyttja of the profundal or in the lake marl of the littoral and record the Late-Glacial forests. L4 is the shrub phase of reforestation with very high Junlperus and rapidly increasing Betula percentages. L5 is the PAZ with a first, L7 with a second dominance of tree-birches, separated by L6 showing a depression in the Betula curve. L4 to L7 can be assigned to the Balling biozone. Possible correlation of the Betula depression to the Older Dryas biozone is discussed. In local PAZ L8 Plnus immigrates and expands. L9 shows a facies difference in that Plnus dominates over Betula in littoral but not in profundal spectra. L8 and L9 belong to the Allerod biozone. In its youngest part the volcanic ash from Laach/Eifel is regularly found (11,000 BP). The local PAZ Ll0 corresponds to the Younger Dryas blozone. The merely slight increase of the NAP indicates that the pine forests of the lowland were not strongly affected by a cooler climate. In order to evaluate the significance of the littoral accumulation of coniferous pollen the littoral profile LQ-150 is compared to the profundal. Radiocarbon stratigraphies derived from different materials are presented in Figures 13 and 14 and in Tables 2 and 3. The hard-water errors in the gyttja samples and the carbonate samples are similar. The samples of terrestrial plant macrofossils are not affected by hard-water errors. Two plateaux of constant age appear in the age-depth relationship; their consequence for biostratigraphy as well as pollen concentration and influx diagrams are discussed. Radiocarbon ages of the Late-Glacial pollen zones are shown in Table 10. The Holocene vegetational history is recorded in the local PAZ L 11 to L30. After a Preboreal (PAZ L11) dominated by pine and birch the expansions of Corylus, Ulmus and Quercus are very rapid. Among these taxa Corylus dominates dur ing the Boreal (PAZ L 12 and L 1 3), whereas the components of the mixed oak forest dominate in the Older Atlantic (PAZ L14 to L16). In the Younger Atlantic (PAZ L 17 to L 19) Fagus and Alnus play an increasing, the mixed oak forest a decreasing role. During the period of local PAZ L19 Neolithic settlers lived on the shore of Lobsigensee. During the Subboreal (PAZ L20 and L21) and the Older Subatlantic (L22 to L25) strong fluctuations of Fagus and often antagonistic peaks of NAP, Alnus, Betula and Corylus can be interpreted as signs of human impact on vegetation. L23 is characterized not only by high values of NAP (especially apophytes and anthropochorous species) but also by the appearance of Juglans, Castanea and Secale which point to the Roman colonization of the area. For a certain period during the Younger Subatlantic (PAZ L26 to L30) the lake was used for retting hemp (Cannabis). Later the dominance of Quercus pollen indicates the importance of wood pastures. The youngest sediments reflect the wide-spread agricultural grass lands and the plantation of Pinus and Picea. Radiocarbon dates for the Holocene are given in Figure 23 and Table 4, the extrapolated ages of the Holocene pollen zones in Table 15. (2) The cross sections: Figures 25 and 26 give a summary of the litho- and palynostratigraphy of the two cross sections. Based on 11 Late-Glacial and 9 Holocene pollen diagrams (in addition to the standard ones), the consistency of the criteria for the definition of the pollen zones is examined in Tables 7 and 8 for the Late-Glacial and in Tables 11 to 14 for the Holocene. Sediment thicknesses across the basin for each pollen zone are presented in these tables as well as in Figures 43 to 45 for the Late-Glacial and in Figures 59 to 65 for the Holocene. Sediment focusing can explain differences between the gyttja cores of the profundal. Focusing is more than compensated for through "stretching" by carbonate precipitation on the littoral terrace. Pollen influx to the cross section are discussed (Chapters 4.1.5. and 4.2.3.). (3) The regional pollen zones: Based on some selected sites between Lake Geneva and Lake Constance regional pollen zones are proposed (Table 16, 17 and 19). (4) Paleoecology: Climatic change in the Late-Glacial can be inferred from Coleoptera, Trichoptera, Chironomidae and d18O of carbonates: a distinct warming is recorded around 12' 600 BP and around 10' 000 BP. The Younger Dryas biozone (10'700-10'000 BP) was the only cooling found in the Late-Glacial. The Betula depression often correlated wi th the Older Dryas biozone was possibl not colder but dryer than the previous period. During the Holocene the lowland site is not very sensitive to the minor climatic changes. Table 22 summarizes climatic and trophic changes before 8'000 BP as deduced from various biostratigraphies studied by a number of authors. Ostracods, Chironomids and fossil pigments indicate that anoxic conditions prevailed during the BoIling (possibly meromixis). Changes in the lake level are illustrated in Figure 74. A first lake-level lowering occurred in the early Holocene (10'000 to 9'000 BP), a second during the Atlantic (about 6'800 to 5'200 BP). The first "shrinking" of the lake volume resulted in a eutrophication recorded by laminations in the profundal and by pigments of Cyanophyceae. The second fall in water level corresponds to an increase of Nymphaeaceae. Human impact can be inferred in three ways: eutrophication of the lake (since the Neolithic), changes of terrestrial vegetation by deforestations (cyclicity of Fagus, see Figures 78 to 80), and enhanced erosion (increasing sedimentation rates by inwashed clay, particularly since the Roman Colonization, see Figures 49 and 81). Summary: This paper was planned as the final report on Lobsigensee. However, a number of issues are not answered but can only be asked more precisely, for example: (1) For the two periods with the highest rates of change, Le. the Bolling and the Preboreal biozones, pollen influx may reflect vegetation dynamics. Detailed investigations of these periods in annually laminated sediments are planned. (2) Biostratigraphies other than palynostratigraphy are needed to estimate the degree of linkage or independence in the development of terrestrial and lacustrine ecosystems. Often our sampling intervals were not identical, thus influencing our temporal resolution. (3) 6180- and 14C-stratigraPhies with high resolution will elucidate the leads and lags of these dynamic periods. Plateaux of constant age in the age-depth relationship have a strong bearing on both biological and geophysical understanding of Late-Glacial and early Holocene developments. (4) Numerical methods applied to the pollen diagrams of the cross section will help to quantify the significance of similari ties and dissimilarities across a single basin (with Prof. Birks). (5) Numerical methods applied to different sites on the Swiss Plateau and on the transect across the Alps will be helpful in evaluating the influence of different environmental factors (with Prof. Birks). (6) A new map 1: 1000 with 50cm-contour lines prov ided by Prof. Zurbuchen will be combined with a grid of cores sampling the transition from lake marl to peat enabling us to calculate paleo-volumes of the lake. This is interesting for the two "shrinking periods" (in Fig. 74A numbers 2-6 and 7-10), both accompanied by eutrophication. The pal eo-volume during the Neoli thic set tlement of the Cortaillod culture linked wi th an est l.mate of trophic change derived from diatoms (Prof. Smol in prep.) could possibly give an indication of the size of the human population of this period. (7) For the period with the antagonism between Fagus peaks and ABC-peaks close collaboration between palynologists, geochemists and archeologists should enable us to determine the influence of prehistoric and historic people on vegetation (collaboration with Prof. Stockli and Prof. Herzig). (8) The core LL-75 taken with a "cold letter box" will be analysed for major and trace elements by Dr. Sturm for 210pb and 137Cs by Prof.von Gunten and for pollen. We will see if our local PAZ L30 really corresponds to the surface sediment and if the small seepage lake reflects modern pollution.