Tailoring polyolefin properties by metallocene catalysis

This doctoral thesis deals with the syntheses of olefin homo- and copolymers using different kind of metallocene catalyst. Ethene, propene, 1-hexene, 1-hexadecene, vinylcyclohexane and phenylnorbornene were homo- or copolymerized with the catalysts. The unbridged benzyl substituted zirconium dichloride catalysts (1-4), ansa- bridged acenaphtyl substituted zirconium dichloride catalysts, ( 5, 6), rac- and meso-ethylene-bis(1-indenyl)zirconium dichlorides, (rac- and meso-8), rac-ethylene-bis(1-indenyl)hafnium dichloride, ( 12), bis(9-fluorenyl)hafnium dichloride (14 ) enantiomerically pure (R)- phenylethyl[(9-fluorenyl-1-indenyl)]ZrCl2, (11), 14 and asymmetric dimethylsilyl[(3-benzylindenyl-(2-methylbenzen[e]indenyl)] zirconium dichloride, (13), were prepared in our laboratory. Dimethylsilyl-bis(1-indenyl)zirconium dichloride, (9), isopropylidene(9-fluorenyl-cyclopentadienyl)zirconium dichloride, (10), and were obtained commercially. The solid-state structures of the catalysts rac- and meso-1 were determined by X-ray crystallography. Computational methods were used for the structure optimization of the catalyst rac- and meso-1 in order to compare the theoretical calculations with the experimental results. Polymerization experiments were conducted in a highly purified autoclave system using low pressures (< 5 bar) of gaseous monomers. The experiments were designed to attain the optimal catalytic activity and a uniform copolymer composition. The prepared homo- and copolymers were characterized by the gel permeation chromatography, GPC, differential scanning calorimetry, DSC, nuclear magnetic resonance, NMR, and Fourier transform infrared spectrometry, FTIR . Molar mass (Mw, Mn), molar mass distribution (Mw/Mn), tacticity, comonomer content, melting temperature, glass transition temperature, and end group structures and content were determined. A special attention was paid on the correlation of the polymer properties with the catalyst structures and polymerization conditions. An intramolecular phenyl coordination was found in phenyl substituted benzyl zirconocenes 1-3 explaining the decreased activity of the catalysts. Novel copolymers poly(propene-co-phenylnorbornene) and poly(propene co-vinylcyclohexane), were synthesized and high molar mass poly(ethene-co-1-hexene) and poly(ethene-co-1-hexadecene) copolymers with elastic properties were prepared. Activation of a hafnocene catalyst was studied with UV-Vis spectrometry and activation process for the synthesis of ultra high molar mass poly(1-hexene) was found out.

Antarctic ice shelf melting and its impact on the global sea ice-ocean system

Earth s ice shelves are mainly located in Antarctica. They cover about 44% of the Antarctic coastline and are a salient feature of the continent. Antarctic ice shelf melting (AISM) removes heat from and inputs freshwater into the adjacent Southern Ocean. Although playing an important role in the global climate, AISM is one of the most important components currently absent in the IPCC climate model. In this study, AISM is introduced into a global sea ice-ocean climate model ORCA2-LIM, following the approach of Beckmann and Goosse (2003; BG03) for the thermodynamic interaction between the ice shelf and ocean. This forms the model ORCA2-LIM-ISP (ISP: ice shelf parameterization), in which not only all the major Antarctic ice shelves but also a number of minor ice shelves are included. Using these two models, ORCA2-LIM and ORCA2-LIM-ISP, the impact of addition of AISM and increasing AISM have been investigated. Using the ORCA2-LIM model, numerical experiments are performed to investigate the sensitivity of the polar sea ice cover and the Antarctic Circumpolar Current (ACC) transport through Drake Passage (DP) to the variations of three sea ice parameters, namely the thickness of newly formed ice in leads (h0), the compressive strength of ice (P*), and the turning angle in the oceanic boundary layer beneath sea ice (θ). It is found that the magnitudes of h0 and P* have little impact on the seasonal sea ice extent, but lead to large changes in the seasonal sea ice volume. The variation in turning angle has little impact on the sea ice extent and volume in the Arctic but tends to reduce them in the Antarctica when ignored. The magnitude of P* has the least impact on the DP transport, while the other two parameters have much larger influences. Numerical results from ORCA2-LIM and ORCA2-LIM-ISP are analyzed to investigate how the inclusion of AISM affects the representation of the Southern Ocean hydrography. Comparisons with data from the World Ocean Circulation Experiment (WOCE) show that the addition of AISM significantly improves the simulated hydrography. It not only warms and freshens the originally too cold and too saline bottom water (AABW), but also warms and enriches the salinity of the originally too cold and too fresh warm deep water (WDW). Addition of AISM also improves the simulated stratification. The close agreement between the simulation with AISM and the observations suggests that the applied parameterization is an adequate way to include the effect of AISM in a global sea ice-ocean climate model. We also investigate the models capability to represent the sea ice-ocean system in the North Atlantic Ocean and the Arctic regions. Our study shows both models (with and without AISM) can successfully reproduce the main features of the sea ice-ocean system. However, both tend to overestimate the ice flux through the Nares Strait, produce a lower temperature and salinity in the Hudson Bay, Baffin Bay and Davis Strait, and miss the deep convection in the Labrador Sea. These deficiencies are mainly attributed to the artificial enlargement of the Nares Strait in the model. In this study, the impact of increasing AISM on the global sea ice-ocean system is thoroughly investigated. This provides a first idea regarding changes induced by increasing AISM. It is shown that the impact of increasing AISM is global and most significant in the Southern Ocean. There, increasing AISM tends to freshen the surface water, to warm the intermediate and deep waters, and to freshen and warm the bottom water. In addition, increasing AISM also leads to changes in the mixed layer depths (MLD) in the deep convection sites in the Southern Ocean, deepening in the Antarctic continental shelf while shoaling in the ACC region. Furthermore, increasing AISM influences the current system in the Southern Ocean. It tends to weaken the ACC, and strengthen the Antarctic coastal current (ACoC) as well as the Weddell Gyre and the Ross Gyre. In addition to the ocean system, increasing AISM also has a notable impact on the Antarctic sea ice cover. Due to the cooling of seawater, sea ice concentration and thickness generally become higher. In austral winter, noticeable increases in sea ice concentration mainly take place near the ice edge. In regards with sea ice thickness, large increases are mainly found along the coast of the Weddell Sea, the Bellingshausen and Amundsen Seas, and the Ross Sea. The overall thickening of sea ice leads to a larger volume of sea ice in Antarctica. In the North Atlantic, increasing AISM leads to remarkable changes in temperature, salinity and density. The water generally becomes warmer, more saline and denser. The most significant warming occurs in the subsurface layer. In contrast, the maximum salinity increase is found at the surface. In addition, the MLD becomes larger along the Greenland-Scotland-Iceland ridge. Global teleconnections due to AISM are studied. The AISM signal is transported with the surface current: the additional freshwater from AISM tends to enhance the northward spreading of the surface water. As a result, more warm and saline water is transported from the tropical region to the North Atlantic Ocean, resulting in warming and salt enrichment there. It would take about 30 40 years to establish a systematic noticeable change in temperature, salinity and MLD in the North Atlantic Ocean according to this study. The changes in hydrography due to increasing AISM are compared with observations. Consistency suggests that increasing AISM is highly likely a major contributor to the recent observed changes in the Southern Ocean. In addition, the AISM might contribute to the salinity contrast between the North Atlantic and North Pacific, which is important for the global thermohaline circulation.

Characterisation and processing of amorphous binary mixtures with low glass transition temperature

The number of drug substances in formulation development in the pharmaceutical industry is increasing. Some of these are amorphous drugs and have glass transition below ambient temperature, and thus they are usually difficult to formulate and handle. One reason for this is the reduced viscosity, related to the stickiness of the drug, that makes them complicated to handle in unit operations. Thus, the aim in this thesis was to develop a new processing method for a sticky amorphous model material. Furthermore, model materials were characterised before and after formulation, using several characterisation methods, to understand more precisely the prerequisites for physical stability of amorphous state against crystallisation. The model materials used were monoclinic paracetamol and citric acid anhydrate. Amorphous materials were prepared by melt quenching or by ethanol evaporation methods. The melt blends were found to have slightly higher viscosity than the ethanol evaporated materials. However, melt produced materials crystallised more easily upon consecutive shearing than ethanol evaporated materials. The only material that did not crystallise during shearing was a 50/50 (w/w, %) blend regardless of the preparation method and it was physically stable at least two years in dry conditions. Shearing at varying temperatures was established to measure the physical stability of amorphous materials in processing and storage conditions. The actual physical stability of the blends was better than the pure amorphous materials at ambient temperature. Molecular mobility was not related to the physical stability of the amorphous blends, observed as crystallisation. Molecular mobility of the 50/50 blend derived from a spectral linewidth as a function of temperature using solid state NMR correlated better with the molecular mobility derived from a rheometer than that of differential scanning calorimetry data. Based on the results obtained, the effect of molecular interactions, thermodynamic driving force and miscibility of the blends are discussed as the key factors to stabilise the blends. The stickiness was found to be affected glass transition and viscosity. Ultrasound extrusion and cutting were successfully tested to increase the processability of sticky material. Furthermore, it was found to be possible to process the physically stable 50/50 blend in a supercooled liquid state instead of a glassy state. The method was not found to accelerate the crystallisation. This may open up new possibilities to process amorphous materials that are otherwise impossible to manufacture into solid dosage forms.

The effect of temperature on height growth of Scots pine in northern Finland

The effect of temperature on height growth of Scots pine in the northern boreal zone in Lapland was studied in two different time scales. Intra-annual growth was monitored in four stands in up to four growing seasons using an approximately biweekly measurement interval. Inter-annual growth was studied using growth records representing seven stands and five geographical locations. All the stands were growing on a dry to semi-dry heath that is a typical site type for pine stands in Finland. The applied methodology is based on applied time-series analysis and multilevel modelling. Intra-annual elongation of the leader shoot correlated with temperature sum accumulation. Height growth ceased when, on average, 41% of the relative temperature sum of the site was achieved (observed minimum and maximum were 38% and 43%). The relative temperature sum was calculated by dividing the actual temperature sum by the long-term mean of the total annual temperature sum for the site. Our results suggest that annual height growth ceases when a location-specific temperature sum threshold is attained. The positive effect of the mean July temperature of the previous year on annual height increment proved to be very strong at high latitudes. The mean November temperature of the year before the previous had a statistically significantly effect on height increment in the three northernmost stands. The effect of mean monthly precipitation on annual height growth was statistically insignificant. There was a non-linear dependence between length and needle density of annual shoots. Exceptionally low height growth results in high needle-density, but the effect is weaker in years of average or good height growth. Radial growth and next year s height growth are both largely controlled by current July temperature. Nevertheless, their growth variation in terms of minimum and maximum is not necessarily strongly correlated. This is partly because height growth is more sensitive to changes in temperature. In addition, the actual effective temperature period is not exactly the same for these two growth components. Yet, there is a long-term balance that was also statistically distinguishable; radial growth correlated significantly with height growth with a lag of 2 years. Temperature periods shorter than a month are more effective variables than mean monthly values, but the improvement is on the scale of modest to good when applying Julian days or growing-degree-days as pointers.

Temperature sensitivity of soil organic matter decomposition in boreal soils

The temperature sensitivity of decomposition of different soil organic matter (SOM) fractions was studied with laboratory incubations using 13C and 14C isotopes to differentiate between SOM of different age. The quality of SOM and the functionality and composition of microbial communities in soils formed under different climatic conditions were also studied. Transferring of organic layers from a colder to a warmer climate was used to assess how changing climate, litter input and soil biology will affect soil respiration and its temperature sensitivity. Together, these studies gave a consistent picture on how warming climate will affect the decomposition of different SOM fractions in Finnish forest soils: the most labile C was least temperature sensitive, indicating that it is utilized irrespective of temperature. The decomposition of intermediate C, with mean residence times from some years to decades, was found to be highly temperature sensitive. Even older, centennially cycling C was again less temperature sensitive, indicating that different stabilizing mechanisms were limiting its decomposition even at higher temperatures. Because the highly temperature sensitive, decadally cycling C, forms a major part of SOM stock in the organic layers of the studied forest soils, these results mean that these soils could lose more carbon during the coming years and decades than estimated earlier. SOM decomposition in boreal forest soils is likely to increase more in response to climate warming, compared to temperate or tropical soils, also because the Q10 is temperature dependent. In the northern soils the warming will occur at a lower temperature range, where Q10 is higher, and a similar increase in temperature causes a higher relative increase in respiration rates. The Q10 at low temperatures was found to be inversely related to SOM quality. At higher temperatures respiration was increasingly limited by low substrate availability.

Preparation and characterization of potato starch films plasticized with polyols

The present study investigated the potato starches and polyols which were used to prepare edible films. The amylose content and the gelatinization properties of various potato starches extracted from different potato cultivars were determined. The amylose content of potato starches varied between 11.9 and 20.1%. Onset temperatures of gelatinization of potato starches in excess water varied independently of the amylose content from 58 to 61°C determined using differential scanning calorimetry (DSC). The crystallinity of selected native starches with low, medium and high amylose content was determined by X-ray diffraction. The relative crystallinity was found to be around 10 13% in selected native potato starches containing 13 17% water. The glass transition temperature, crystallization melting behavior and relaxations of polyols, erythritol, sorbitol and xylitol, were determined using (DSC), dielectric analysis (DEA) and dynamic mechanical analysis (DMA). The glass transition temperatures of xylitol and sorbitol decreased as a result of water plasticization. Anhydrous amorphous erythritol crystallized rapidly. Edible films were obtained from solutions containing gelatinized starch, plasticizer (polyol or binary polyol mixture) and water by casting and evaporating water at 35°C. The present study investigated effects of plasticizer type and content on physical and mechanical properties of edible films stored at various relative water vapor pressures (RVP). The crystallinity of edible films with low, medium and high amylose content was determined by X-ray diffraction and they were found to be practically amorphous. Water sorption and water vapor permeability (WVP) of films was affected by the type and content of plasticizer. Water vapor permeability of films increased with increasing plasticizer content and storage RVP. Generally, Young's modulus and tensile strength decreased with increasing plasticizer and water content with a concurrent increase in elongation at break of films. High contents of xylitol and sorbitol resulted in changes in physical and mechanical properties of films probably due to phase separation and crystallization of xylitol and sorbitol which was not observed when binary polyol mixtures were used as plasticizers. The mechanical properties and the water vapor permeability (WVP) of the films were found to be independent of the amylose content.

Stretching single DNA-molecules with temperature-stabilized optical tweezers

This thesis consists of two parts; in the first part we performed a single-molecule force extension measurement with 10kb long DNA-molecules from phage-λ to validate the calibration and single-molecule capability of our optical tweezers instrument. Fitting the worm-like chain interpolation formula to the data revealed that ca. 71% of the DNA tethers featured a contour length within ±15% of the expected value (3.38 µm). Only 25% of the found DNA had a persistence length between 30 and 60 nm. The correct value should be within 40 to 60 nm. In the second part we designed and built a precise temperature controller to remove thermal fluctuations that cause drifting of the optical trap. The controller uses feed-forward and PID (proportional-integral-derivative) feedback to achieve 1.58 mK precision and 0.3 K absolute accuracy. During a 5 min test run it reduced drifting of the trap from 1.4 nm/min in open-loop to 0.6 nm/min in closed-loop.

Parameters Affecting the Extraction of Polycyclic Aromatic Hydrocarbons with Pressurised Hot Water

Pressurised hot water extraction (PHWE) exploits the unique temperature-dependent solvent properties of water minimising the use of harmful organic solvents. Water is environmentally friendly, cheap and easily available extraction medium. The effects of temperature, pressure and extraction time in PHWE have often been studied, but here the emphasis was on other parameters important for the extraction, most notably the dimensions of the extraction vessel and the stability and solubility of the analytes to be extracted. Non-linear data analysis and self-organising maps were employed in the data analysis to obtain correlations between the parameters studied, recoveries and relative errors. First, pressurised hot water extraction (PHWE) was combined on-line with liquid chromatography-gas chromatography (LC-GC), and the system was applied to the extraction and analysis of polycyclic aromatic hydrocarbons (PAHs) in sediment. The method is of superior sensitivity compared with the traditional methods, and only a small 10 mg sample was required for analysis. The commercial extraction vessels were replaced by laboratory-made stainless steel vessels because of some problems that arose. The performance of the laboratory-made vessels was comparable to that of the commercial ones. In an investigation of the effect of thermal desorption in PHWE, it was found that at lower temperatures (200ºC and 250ºC) the effect of thermal desorption is smaller than the effect of the solvating property of hot water. At 300ºC, however, thermal desorption is the main mechanism. The effect of the geometry of the extraction vessel on recoveries was studied with five specially constructed extraction vessels. In addition to the extraction vessel geometry, the sediment packing style and the direction of water flow through the vessel were investigated. The geometry of the vessel was found to have only minor effect on the recoveries, and the same was true of the sediment packing style and the direction of water flow through the vessel. These are good results because these parameters do not have to be carefully optimised before the start of extractions. Liquid-liquid extraction (LLE) and solid-phase extraction (SPE) were compared as trapping techniques for PHWE. LLE was more robust than SPE and it provided better recoveries and repeatabilities than did SPE. Problems related to blocking of the Tenax trap and unrepeatable trapping of the analytes were encountered in SPE. Thus, although LLE is more labour intensive, it can be recommended over SPE. The stabilities of the PAHs in aqueous solutions were measured using a batch-type reaction vessel. Degradation was observed at 300ºC even with the shortest heating time. Ketones and quinones and other oxidation products were observed. Although the conditions of the stability studies differed considerably from the extraction conditions in PHWE, the results indicate that the risk of analyte degradation must be taken into account in PHWE. The aqueous solubilities of acenaphthene, anthracene and pyrene were measured, first below and then above the melting point of the analytes. Measurements below the melting point were made to check that the equipment was working, and the results were compared with those obtained earlier. Good agreement was found between the measured and literature values. A new saturation cell was constructed for the solubility measurements above the melting point of the analytes because the flow-through saturation cell could not be used above the melting point. An exponential relationship was found between the solubilities measured for pyrene and anthracene and temperature.

Role of Temperature in the Biological Activity of a Boreal Forest

In northern latitudes, temperature is the key factor driving the temporal scales of biological activity, namely the length of the growing season and the seasonal efficiency of photosynthesis. The formation of atmospheric concentrations of biogenic volatile organic compounds (BVOCs) are linked to the intensity of biological activity. However, interdisciplinary knowledge of the role of temperature in the biological processes related to the annual cycle and photosynthesis and atmospheric chemistry is not fully understood. The aim of this study was to improve understanding of the role of temperature in these three interlinked areas: 1) onset of growing season, 2) photosynthetic efficiency and 3) BVOC air concentrations in a boreal forest. The results present a cross-section of the role of temperature on different spatial (southern northern boreal), structural (tree forest stand - forest) and temporal (day-season- year) scales. The fundamental status of the Thermal Time model in predicting the onset of spring recovery was confirmed. However, it was recommended that sequential models would be more appropriate tools when the onset of the growing season is estimated under a warmer climate. A similar type of relationship between photosynthetic efficiency and temperature history was found in both southern and northern boreal forest stands. This result draws attention to the critical question of the seasonal efficiency of coniferous species to emit organic compounds under a warmer climate. New knowledge about the temperature dependence of the concentrations of biogenic volatile organic compounds in a boreal forest stand was obtained. The seasonal progress and the inter-correlation of BVOC concentrations in ambient air indicated a link to biological activity. Temperature was found to be the main driving factor for the concentrations. However, in addition to temperature, other factors may play a significant role here, especially when the peak concentrations are studied. There is strong evidence that the spring recovery and phenological events of many plant species have already advanced in Europe. This study does not fully support this observation. In a boreal forest, changes in the annual cycle, especially the temperature requirement in winter, would have an impact on the atmospheric BVOC composition. According to this study, more joint phenological and BVOC field observations and laboratory experiments are still needed to improve these scenarios.

Microbial activities in boreal soils: Biodegradation of organic contaminants at low temperature and ammonia oxidation

This thesis deals with the response of biodegradation of selected anthropogenic organic contaminants and natural autochthonous organic matter to low temperature in boreal surface soils. Furthermore, the thesis describes activity, diversity and population size of autotrophic ammonia-oxidizing bacteria (AOB) in a boreal soil used for landfarming of oil-refinery wastes, and presents a new approach, in which the particular AOB were enriched and cultivated in situ from the landfarming soil onto cation exchange membranes. This thesis demonstrates that rhizosphere fraction of natural forest humus soil and agricultural clay loam soil from Helsinki Metropolitan area were capable of degrading of low to moderate concentrations (0.2 50 µg cm-3) of PCP, phenanthrene and 2,4,5-TCP at temperatures realistic to boreal climate (-2.5 to +15 °C). At the low temperatures, the biodegradation of PCP, phenanthrene and 2,4,5-TCP was more effective (Q10-values from 1.6 to 7.6) in the rhizosphere fraction of the forest soil than in the agricultural soil. Q10-values of endogenous soil respiration (carbon dioxide evolution) and selected hydrolytic enzyme activities (acetate-esterase, butyrate-esterase and β-glucosidase) in acid coniferous forest soil were 1.6 to 2.8 at temperatures from -3 to +30 °C. The results indicated that the temperature dependence of decomposition of natural autochthonous soil organic matter in the studied coniferous forest was only moderate. The numbers of AOB in the landfarming (sandy clay loam) soil were determined with quantitative polymerase chain reaction (real-time PCR) and with Most Probable Number (MPN) methods, and potential ammonium oxidation activity was measured with the chlorate inhibition technique. The results indicated presence of large and active AOB populations in the heavily oil-contaminated and urea-fertilised landfarming soil. Assessment of the populations of AOB with denaturing gradient gel electrophoresis (DGGE) profiling and sequence analysis of PCR-amplified 16S rRNA genes showed that Nitrosospira-like AOB in clusters 2 and 3 were predominant in the oily landfarming soil. This observation was supported by fluorescence in situ hybridization (FISH) analysis of the AOB grown on the soil-incubated cation-exchange membranes. The results of this thesis expand the suggested importance of Nitrosospira-like AOB in terrestrial environments to include chronically oil-contaminated soils.

Optical Studies of the Antarctic Glacio-Oceanic System

The Antarctic system comprises of the continent itself, Antarctica, and the ocean surrounding it, the Southern Ocean. The system has an important part in the global climate due to its size, its high latitude location and the negative radiation balance of its large ice sheets. Antarctica has also been in focus for several decades due to increased ultraviolet (UV) levels caused by stratospheric ozone depletion, and the disintegration of its ice shelves. In this study, measurements were made during three Austral summers to study the optical properties of the Antarctic system and to produce radiation information for additional modeling studies. These are related to specific phenomena found in the system. During the summer of 1997-1998, measurements of beam absorption and beam attenuation coefficients, and downwelling and upwelling irradiance were made in the Southern Ocean along a S-N transect at 6°E. The attenuation of photosynthetically active radiation (PAR) was calculated and used together with hydrographic measurements to judge whether the phytoplankton in the investigated areas of the Southern Ocean are light limited. By using the Kirk formula the diffuse attenuation coefficient was linked to the absorption and scattering coefficients. The diffuse attenuation coefficients (Kpar) for PAR were found to vary between 0.03 and 0.09 1/m. Using the values for KPAR and the definition of the Sverdrup critical depth, the studied Southern Ocean plankton systems were found not to be light limited. Variabilities in the spectral and total albedo of snow were studied in the Queen Maud Land region of Antarctica during the summers of 1999-2000 and 2000-2001. The measurement areas were the vicinity of the South African Antarctic research station SANAE 4, and a traverse near the Finnish Antarctic research station Aboa. The midday mean total albedos for snow were between 0.83, for clear skies, and 0.86, for overcast skies, at Aboa and between 0.81 and 0.83 for SANAE 4. The mean spectral albedo levels at Aboa and SANAE 4 were very close to each other. The variations in the spectral albedos were due more to differences in ambient conditions than variations in snow properties. A Monte-Carlo model was developed to study the spectral albedo and to develop a novel nondestructive method to measure the diffuse attenuation coefficient of snow. The method was based on the decay of upwelling radiation moving horizontally away from a source of downwelling light. This was assumed to have a relation to the diffuse attenuation coefficient. In the model, the attenuation coefficient obtained from the upwelling irradiance was higher than that obtained using vertical profiles of downwelling irradiance. The model results were compared to field measurements made on dry snow in Finnish Lapland and they correlated reasonably well. Low-elevation (below 1000 m) blue-ice areas may experience substantial melt-freeze cycles due to absorbed solar radiation and the small heat conductivity in the ice. A two-dimensional (x-z) model has been developed to simulate the formation and water circulation in the subsurface ponds. The model results show that for a physically reasonable parameter set the formation of liquid water within the ice can be reproduced. The results however are sensitive to the chosen parameter values, and their exact values are not well known. Vertical convection and a weak overturning circulation is generated stratifying the fluid and transporting warmer water downward, thereby causing additional melting at the base of the pond. In a 50-year integration, a global warming scenario mimicked by a decadal scale increase of 3 degrees per 100 years in air temperature, leads to a general increase in subsurface water volume. The ice did not disintegrate due to the air temperature increase after the 50 year integration.

Effects of climate and morphology on temperature conditions of lakes

To a large extent, lakes can be described with a one-dimensional approach, as their main features can be characterized by the vertical temperature profile of the water. The development of the profiles during the year follows the seasonal climate variations. Depending on conditions, lakes become stratified during the warm summer. After cooling, overturn occurs, water cools and an ice cover forms. Typically, water is inversely stratified under the ice, and another overturn occurs in spring after the ice has melted. Features of this circulation have been used in studies to distinguish between lakes in different areas, as basis for observation systems and even as climate indicators. Numerical models can be used to calculate temperature in the lake, on the basis of the meteorological input at the surface. The simple form is to solve the surface temperature. The depth of the lake affects heat transfer, together with other morphological features, the shape and size of the lake. Also the surrounding landscape affects the formation of the meteorological fields over the lake and the energy input. For small lakes the shading by the shores affects both over the lake and inside the water body bringing limitations for the one-dimensional approach. A two-layer model gives an approximation for the basic stratification in the lake. A turbulence model can simulate vertical temperature profile in a more detailed way. If the shape of the temperature profile is very abrupt, vertical transfer is hindered, having many important consequences for lake biology. One-dimensional modelling approach was successfully studied comparing a one-layer model, a two-layer model and a turbulence model. The turbulence model was applied to lakes with different sizes, shapes and locations. Lake models need data from the lakes for model adjustment. The use of the meteorological input data on different scales was analysed, ranging from momentary turbulent changes over the lake to the use of the synoptical data with three hour intervals. Data over about 100 past years were used on the mesoscale at the range of about 100 km and climate change scenarios for future changes. Increasing air temperature typically increases water temperature in epilimnion and decreases ice cover. Lake ice data were used for modelling different kinds of lakes. They were also analyzed statistically in global context. The results were also compared with results of a hydrological watershed model and data from very small lakes for seasonal development.