Nitrous oxide emissions from selected natural and managed northern ecosystems

Microbial activity in soils is the main source of nitrous oxide (N2O) to the atmosphere. Nitrous oxide is a strong greenhouse gas in the troposphere and participates in ozone destructive reactions in the stratosphere. The constant increase in the atmospheric concentration, as well as uncertainties in the known sources and sinks of N2O underline the need to better understand the processes and pathways of N2O in terrestrial ecosystems. This study aimed at quantifying N2O emissions from soils in northern Europe and at investigating the processes and pathways of N2O from agricultural and forest ecosystems. Emissions were measured in forest ecosystems, agricultural soils and a landfill, using the soil gradient, chamber and eddy covariance methods. Processes responsible for N2O production, and the pathways of N2O from the soil to the atmosphere, were studied in the laboratory and in the field. These ecosystems were chosen for their potential importance to the national and global budget of N2O. Laboratory experiments with boreal agricultural soils revealed that N2O production increases drastically with soil moisture content, and that the contribution of the nitrification and denitrification processes to N2O emissions depends on soil type. Laboratory study with beech (Fagus sylvatica) seedlings demonstrated that trees can serve as conduits for N2O from the soil to the atmosphere. If this mechanism is important in forest ecosystems, the current emission estimates from forest soils may underestimate the total N2O emissions from forest ecosystems. Further field and laboratory studies are needed to evaluate the importance of this mechanism in forest ecosystems. The emissions of N2O from northern forest ecosystems and a municipal landfill were highly variable in time and space. The emissions of N2O from boreal upland forest soil were among the smallest reported in the world. Despite the low emission rates, the soil gradient method revealed a clear seasonal variation in N2O production. The organic topsoil was responsible for most of the N2O production and consumption in this forest soil. Emissions from the municipal landfill were one to two orders of magnitude higher than those from agricultural soils, which are the most important source of N2O to the atmosphere. Due to their small areal coverage, landfills only contribute minimally to national N2O emissions in Finland. The eddy covariance technique was demonstrated to be useful for measuring ecosystem-scale emissions of N2O in forest and landfill ecosystems. Overall, more measurements and integration between different measurement techniques are needed to capture the large variability in N2O emissions from natural and managed northern ecosystems.

Substance flow analysis in Finland : Four case studies on N and P flows

Nitrogen (N) and phosphorus (P) are essential elements for all living organisms. However, in excess, they contribute to several environmental problems such as aquatic and terrestrial eutrophication. Globally, human action has multiplied the volume of N and P cycling since the onset of industrialization. The multiplication is a result of intensified agriculture, increased energy consumption and population growth. Industrial ecology (IE) is a discipline, in which human interaction with the ecosystems is investigated using a systems analytical approach. The main idea behind IE is that industrial systems resemble ecosystems, and, like them, industrial systems can then be described using material, energy and information flows and stocks. Industrial systems are dependent on the resources provided by the biosphere, and these two cannot be separated from each other. When studying substance flows, the aims of the research from the viewpoint of IE can be, for instance, to elucidate the ways how the cycles of a certain substance could be more closed and how the flows of a certain substance could be decreased per unit of production (= dematerialization). In Finland, N and P are studied widely in different ecosystems and environmental emissions. A holistic picture comparing different societal systems is, however, lacking. In this thesis, flows of N and P were examined in Finland using substance flow analysis (SFA) in the following four subsystems: I) forest industry and use of wood fuels, II) food production and consumption, III) energy, and IV) municipal waste. A detailed analysis at the end of the 1990s was performed. Furthermore, historical development of the N and P flows was investigated in the energy system (III) and the municipal waste system (IV). The main research sources were official statistics, literature, monitoring data, and expert knowledge. The aim was to identify and quantify the main flows of N and P in Finland in the four subsystems studied. Furthermore, the aim was to elucidate whether the nutrient systems are cyclic or linear, and to identify how these systems could be more efficient in the use and cycling of N and P. A final aim was to discuss how this type of an analysis can be used to support decision-making on environmental problems and solutions. Of the four subsystems, the food production and consumption system and the energy system created the largest N flows in Finland. For the creation of P flows, the food production and consumption system (Paper II) was clearly the largest, followed by the forest industry and use of wood fuels and the energy system. The contribution of Finland to N and P flows on a global scale is low, but when compared on a per capita basis, we are one of the largest producers of these flows, with relatively high energy and meat consumption being the main reasons. Analysis revealed the openness of all four systems. The openness is due to the high degree of internationality of the Finnish markets, the large-scale use of synthetic fertilizers and energy resources and the low recycling rate of many waste fractions. Reduction in the use of fuels and synthetic fertilizers, reorganization of the structure of energy production, reduced human intake of nutrients and technological development are crucial in diminishing the N and P flows. To enhance nutrient recycling and replace inorganic fertilizers, recycling of such wastes as wood ash and sludge could be promoted. SFA is not usually sufficiently detailed to allow specific recommendations for decision-making to be made, but it does yield useful information about the relative magnitude of the flows and may reveal unexpected losses. Sustainable development is a widely accepted target for all human action. SFA is one method that can help to analyse how effective different efforts are in leading to a more sustainable society. SFA's strength is that it allows a holistic picture of different natural and societal systems to be drawn. Furthermore, when the environmental impact of a certain flow is known, the method can be used to prioritize environmental policy efforts.

Export of organic matter, sulphate and base cations from boreal headwater catchments downstream to the coast: impacts of land use and climate

The terrestrial export of dissolved organic matter (DOM) is associated with climate, vegetation and land use, and thus is under the influence of climatic variability and human interference with terrestrial ecosystems, their soils and hydrological cycles. The present study provides an assessment of spatial variation of DOM concentrations and export, and interactions between DOM, catchment characteristics, land use and climatic factors in boreal catchments. The influence of catchment characteristics, land use and climatic drivers on the concentrations and export of total organic carbon (TOC), total organic nitrogen (TON) and dissolved organic phosphorus (DOP) was estimated using stream water quality, forest inventory and climatic data from 42 Finnish pristine forested headwater catchments, and water quality monitoring, GIS land use, forest inventory and climatic data from the 36 main Finnish rivers (and their sub-catchments) flowing to the Baltic Sea. Moreover, the export of DOM in relation to land use along a European climatic gradient was studied using river water quality and land use data from four European areas. Additionally, the role of organic and minerogenic acidity in controlling pH levels in Finnish rivers and pristine streams was studied by measuring organic anion, sulphate (SO4) and base cation (Ca, Mg, K and Na) concentrations. In all study catchments, TOC was a major fraction of DOM, with much lower proportions of TON and DOP. Moreover, most of TOC and TON was in a dissolved form. The correlation between TOC and TON concentrations was strong and TOC concentrations explained 78% of the variation in TON concentrations in pristine headwater streams. In a subgroup of 20 headwater catchments with similar climatic conditions and low N deposition in eastern Finland, the proportion of peatlands in the catchment and the proportion of Norway spruce (Picea abies Karsten) of the tree stand had the strongest correlation with the TOC and TON concentrations and export. In Finnish river basins, TOC export increased with the increasing proportion of peatland in the catchment, whereas TON export increased with increasing extent of agricultural land. The highest DOP concentrations and export were recorded in river basins with a high extent of agricultural land and urban areas, reflecting the influence of human impact on DOP loads. However, the most important predictor for TOC, TON and DOP export in Finnish rivers was the proportion of upstream lakes in the catchment. The higher the upstream lake percentage, the lower the export indicating organic matter retention in lakes. Molar TOC:TON ratio decreased from headwater catchments covered by forests and peatlands to the large river basins with mixed land use, emphasising the effect of the land use gradient on the stoichiometry of rivers. This study also demonstrated that the land use of the catchments is related to both organic and minerogenic acidity in rivers and pristine headwater streams. Organic anion dominated in rivers and streams situated in northern Finland, reflecting the higher extent of peatlands in these areas, whereas SO4 dominated in southern Finland and on western coastal areas, where the extent of fertile areas, agricultural land, urban areas, acid sulphate soils, and sulphate deposition is highest. High TOC concentrations decreased pH values in the stream and river water, whereas no correlation between SO4 concentrations and pH was observed. This underlines the importance of organic acids in controlling pH levels in Finnish pristine headwater streams and main rivers. High SO4 concentrations were associated with high base cation concentrations and fertile areas, which buffered the effects of SO4 on pH.

Distribution of non-biting midges (Diptera, Chironomidae) in subarctic lakes of Finnish Lapland applications in lake classification and palaeolimnology

Climate change contributes directly or indirectly to changes in species distributions, and there is very high confidence that recent climate warming is already affecting ecosystems. The Arctic has already experienced the greatest regional warming in recent decades, and the trend is continuing. However, studies on the northern ecosystems are scarce compared to more southerly regions. Better understanding of the past and present environmental change is needed to be able to forecast the future. Multivariate methods were used to explore the distributional patterns of chironomids in 50 shallow (≤ 10m) lakes in relation to 24 variables determined in northern Fennoscandia at the ecotonal area from the boreal forest in the south to the orohemiarctic zone in the north. Highest taxon richness was noted at middle elevations around 400 m a.s.l. Significantly lower values were observed from cold lakes situated in the tundra zone. Lake water alkalinity had the strongest positive correlation with the taxon richness. Many taxa had preference for lakes either on tundra area or forested area. The variation in the chironomid abundance data was best correlated with sediment organic content (LOI), lake water total organic carbon content, pH and air temperature, with LOI being the strongest variable. Three major lake groups were separated on the basis of their chironomid assemblages: (i) small and shallow organic-rich lakes, (ii) large and base-rich lakes, and (iii) cold and clear oligotrophic tundra lakes. Environmental variables best discriminating the lake groups were LOI, taxon richness, and Mg. When repeated, this kind of an approach could be useful and efficient in monitoring the effects of global change on species ranges. Many species of fast spreading insects, including chironomids, show a remarkable ability to track environmental changes. Based on this ability, past environmental conditions have been reconstructed using their chitinous remains in the lake sediment profiles. In order to study the Holocene environmental history of subarctic aquatic systems, and quantitatively reconstruct the past temperatures at or near the treeline, long sediment cores covering the last 10000 years (the Holocene) were collected from three lakes. Lower temperature values than expected based on the presence of pine in the catchment during the mid-Holocene were reconstructed from a lake with great water volume and depth. The lake provided thermal refuge for profundal, cold adapted taxa during the warm period. In a shallow lake, the decrease in the reconstructed temperatures during the late Holocene may reflect the indirect response of the midges to climate change through, e.g., pH change. The results from three lakes indicated that the response of chironomids to climate have been more or less indirect. However, concurrent shifts in assemblages of chironomids and vegetation in two lakes during the Holocene time period indicated that the midges together with the terrestrial vegetation had responded to the same ultimate cause, which most likely was the Holocene climate change. This was also supported by the similarity in the long-term trends in faunal succession for the chironomid assemblages in several lakes in the area. In northern Finnish Lapland the distribution of chironomids were significantly correlated with physical and limnological factors that are most likely to change as a result of future climate change. The indirect and individualistic response of aquatic systems, as reconstructed using the chironomid assemblages, to the climate change in the past suggests that in the future, the lake ecosystems in the north do not respond in one predictable way to the global climate change. Lakes in the north may respond to global climate change in various ways that are dependent on the initial characters of the catchment area and the lake.

Effect of pulse profile variations on measurement of eccentricity in orbits of Cen X-3 and SMC X-1

It has long been argued that better timing precision allowed by satellites like Rossi X-ray Timing Explorer (RXTE) will allow us to measure the orbital eccentricity and the angle of periastron of some of the bright persistent high-mass X-ray binaries (HMXBs) and hence a possible measurement of apsidal motion in these system. Measuring the rate of apsidal motion allows one to estimate the apsidal motion constant of the mass losing companion star and hence allows for the direct testing of the stellar structure models for these giant stars present in the HMXBs. In the present paper, we use the archival RXTE data of two bright persistent sources, namely Cen X-3 and SMC X-1, to measure the very small orbital eccentricity and the angle of periastron. We find that the small variations in the pulse profiles of these sources, rather than the intrinsic time resolution provided by RXTE, limit the accuracy with which we can measure arrival time of the pulses from these sources. This influences the accuracy with which one can measure the orbital parameters, especially the very small eccentricity and the angle of periastron in these sources. The observations of SMC X-1 in the year 2000 were taken during the high-flux state of the source and we could determine the orbital eccentricity and omega using this data set.

Multi-GNSS precise point positioning with raw single-frequency and dual-frequency measurement models

The emergence of multiple satellite navigation systems, including BDS, Galileo, modernized GPS, and GLONASS, brings great opportunities and challenges for precise point positioning (PPP). We study the contributions of various GNSS combinations to PPP performance based on undifferenced or raw observations, in which the signal delays and ionospheric delays must be considered. A priori ionospheric knowledge, such as regional or global corrections, strengthens the estimation of ionospheric delay parameters. The undifferenced models are generally more suitable for single-, dual-, or multi-frequency data processing for single or combined GNSS constellations. Another advantage over ionospheric-free PPP models is that undifferenced models avoid noise amplification by linear combinations. Extensive performance evaluations are conducted with multi-GNSS data sets collected from 105 MGEX stations in July 2014. Dual-frequency PPP results from each single constellation show that the convergence time of undifferenced PPP solution is usually shorter than that of ionospheric-free PPP solutions, while the positioning accuracy of undifferenced PPP shows more improvement for the GLONASS system. In addition, the GLONASS undifferenced PPP results demonstrate performance advantages in high latitude areas, while this impact is less obvious in the GPS/GLONASS combined configuration. The results have also indicated that the BDS GEO satellites have negative impacts on the undifferenced PPP performance given the current “poor” orbit and clock knowledge of GEO satellites. More generally, the multi-GNSS undifferenced PPP results have shown improvements in the convergence time by more than 60 % in both the single- and dual-frequency PPP results, while the positioning accuracy after convergence indicates no significant improvements for the dual-frequency PPP solutions, but an improvement of about 25 % on average for the single-frequency PPP solutions.

Assessment of code bias variations of BDS triple-frequency signals and their impacts on ambiguity resolution for long baselines

Carrier phase ambiguity resolution over long baselines is challenging in BDS data processing. This is partially due to the variations of the hardware biases in BDS code signals and its dependence on elevation angles. We present an assessment of satellite-induced code bias variations in BDS triple-frequency signals and the ambiguity resolutions procedures involving both geometry-free and geometry-based models. First, since the elevation of a GEO satellite remains unchanged, we propose to model the single-differenced fractional cycle bias with widespread ground stations. Second, the effects of code bias variations induced by GEO, IGSO and MEO satellites on ambiguity resolution of extra-wide-lane, wide-lane and narrow-lane combinations are analyzed. Third, together with the IGSO and MEO code bias variations models, the effects of code bias variations on ambiguity resolution are examined using 30-day data collected over the baselines ranging from 500 to 2600 km in 2014. The results suggest that although the effect of code bias variations on the extra-wide-lane integer solution is almost ignorable due to its long wavelength, the wide-lane integer solutions are rather sensitive to the code bias variations. Wide-lane ambiguity resolution success rates are evidently improved when code bias variations are corrected. However, the improvement of narrow-lane ambiguity resolution is not obvious since it is based on geometry-based model and there is only an indirect impact on the narrow-lane ambiguity solutions.

Physical properties of meteorites and their role in planetology

Together with cosmic spherules, interplanetary dust particles and lunar samples returned by Apollo and Luna missions, meteorites are the only source of extraterrestrial material on Earth. The physical properties of meteorites, especially their magnetic susceptibility, bulk and grain density, porosity and paleomagnetic information, have wide applications in planetary research and can reveal information about origin and internal structure of asteroids. Thus, an expanded database of meteorite physical properties was compiled with new measurements done in meteorite collections across Europe using a mobile laboratory facility. However, the scale problem may bring discrepancies in the comparison of asteroid and meteorite properties. Due to inhomogenity, the physical properties of meteorites studied on a centimeter or millimeter scale may differ from those of asteroids determined on kilometer scales. Further difference may arise from shock effects, space and terrestrial weathering and from difference in material properties at various temperatures. Close attention was given to the reliability of the paleomagnetic and paleointensity information in meteorites and the methodology to test for magnetic overprints was prepared and verified.

Investigations of planetary boundary layer processes and particle formation in the atmosphere of planet Mars

The planet Mars is the Earth's neighbour in the Solar System. Planetary research stems from a fundamental need to explore our surroundings, typical for mankind. Manned missions to Mars are already being planned, and understanding the environment to which the astronauts would be exposed is of utmost importance for a successful mission. Information of the Martian environment given by models is already now used in designing the landers and orbiters sent to the red planet. In particular, studies of the Martian atmosphere are crucial for instrument design, entry, descent and landing system design, landing site selection, and aerobraking calculations. Research of planetary atmospheres can also contribute to atmospheric studies of the Earth via model testing and development of parameterizations: even after decades of modeling the Earth's atmosphere, we are still far from perfect weather predictions. On a global level, Mars has also been experiencing climate change. The aerosol effect is one of the largest unknowns in the present terrestrial climate change studies, and the role of aerosol particles in any climate is fundamental: studies of climate variations on another planet can help us better understand our own global change. In this thesis I have used an atmospheric column model for Mars to study the behaviour of the lowest layer of the atmosphere, the planetary boundary layer (PBL), and I have developed nucleation (particle formation) models for Martian conditions. The models were also coupled to study, for example, fog formation in the PBL. The PBL is perhaps the most significant part of the atmosphere for landers and humans, since we live in it and experience its state, for example, as gusty winds, nightfrost, and fogs. However, PBL modelling in weather prediction models is still a difficult task. Mars hosts a variety of cloud types, mainly composed of water ice particles, but also CO2 ice clouds form in the very cold polar night and at high altitudes elsewhere. Nucleation is the first step in particle formation, and always includes a phase transition. Cloud crystals on Mars form from vapour to ice on ubiquitous, suspended dust particles. Clouds on Mars have a small radiative effect in the present climate, but it may have been more important in the past. This thesis represents an attempt to model the Martian atmosphere at the smallest scales with high resolution. The models used and developed during the course of the research are useful tools for developing and testing parameterizations for larger-scale models all the way up to global climate models, since the small-scale models can describe processes that in the large-scale models are reduced to subgrid (not explicitly resolved) scale.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Orbital evolution and orbital phase resolved spectroscopy of the HMXB pulsar 4U 1538-52 with RXTE-PCA and BeppoSAX

We report here results from detailed timing and spectral studies of the high mass X-ray binary pulsar 4U 1538-52 over several binary periods using observations made with the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX satellites. Pulse timing analysis with the 2003 RXTE data over two binary orbits confirms an eccentric orbit of the system. Combining the orbitial parameters determined from this observation with the earlier measurements we did not find any evidence of orbital decay in this X-ray binary. We have carried out orbital phase resolved spectroscopy to measure changes in the spectral parameters with orbital phase, particularly the absorption column density and the iron line flux. The RXTE-PCA spectra in the 3-20 keV energy range were fitted with a power law and a high energy cut-off along with a Gaussian line at similar to 6.4 keV, whereas the BeppoSAX spectra needed only a power law and Gaussian emission line at similar to 6.4keV in the restricted energy range of 0.3-10.0 keV. An absorption along the line of sight was included for both the RXTE and BeppoSAX data. The variation of the free spectral parameters over the binary orbit was investigated and we found that the variation of the column density of absorbing material in the line of sight with orbital phase is in reasonable agreement with a simple model of a spherically symmetric stellar wind from the companion star.

Enantiodiscrimination and extraction of short and long range homo- and hetero-nuclear residual dipolar couplings by a spin selective correlation experiment

A two dimensional correlation experiment for the measurement of short and long range homo- and hetero- nuclear residual dipolar couplings (RDCs) from the broad and featureless proton NMR spectra including C-13 satellites is proposed. The method employs a single natural abundant C-13 spin as a spy nucleus to probe all the coupled protons and permits the determination of RDCs of negligible strengths. The technique has been demonstrated for the study of organic chiral molecules aligned in chiral liquid crystal, where additional challenge is to unravel the overlapped spectrum of enantiomers. The significant advantage of the method is demonstrated in better chiral discrimination using homonuclear RDCs as additional parameters. (C) 2010 Elsevier B.V. All rights reserved.

Carbon dioxide and methane exchange between a boreal pristine lake and the atmosphere

Lakes serve as sites for terrestrially fixed carbon to be remineralized and transferred back to the atmosphere. Their role in regional carbon cycling is especially important in the Boreal Zone, where lakes can cover up to 20% of the land area. Boreal lakes are often characterized by the presence of a brown water colour, which implies high levels of dissolved organic carbon from the surrounding terrestrial ecosystem, but the load of inorganic carbon from the catchment is largely unknown. Organic carbon is transformed to methane (CH4) and carbon dioxide (CO2) in biological processes that result in lake water gas concentrations that increase above atmospheric equilibrium, thus making boreal lakes as sources of these important greenhouse gases. However, flux estimates are often based on sporadic sampling and modelling and actual flux measurements are scarce. Thus, the detailed temporal flux dynamics of greenhouse gases are still largely unknown. ----- One aim here was to reveal the natural dynamics of CH4 and CO2 concentrations and fluxes in a small boreal lake. The other aim was to test the applicability of a measuring technique for CO2 flux, i.e. the eddy covariance (EC) technique, and a computational method for estimation of primary production and community respiration, both commonly used in terrestrial research, in this lake. Continuous surface water CO2 concentration measurements, also needed in free-water applications to estimate primary production and community respiration, were used over two open water periods in a study of CO2 concentration dynamics. Traditional methods were also used to measure gas concentration and fluxes. The study lake, Valkea-Kotinen, is a small, humic, headwater lake within an old-growth forest catchment with no local anthropogenic disturbance and thus possible changes in gas dynamics reflect the natural variability in lake ecosystems. CH4 accumulated under the ice and in the hypolimnion during summer stratification. The surface water CH4 concentration was always above atmospheric equilibrium and thus the lake was a continuous source of CH4 to the atmosphere. However, the annual CH4 fluxes were small, i.e. 0.11 mol m-2 yr-1, and the timing of fluxes differed from that of other published estimates. The highest fluxes are usually measured in spring after ice melt but in Lake Valkea-Kotinen CH4 was effectively oxidised in spring and highest effluxes occurred in autumn after summer stratification period. CO2 also accumulated under the ice and the hypolimnetic CO2 concentration increased steadily during stratification period. The surface water CO2 concentration was highest in spring and in autumn, whereas during the stable stratification it was sometimes under atmospheric equilibrium. It showed diel, daily and seasonal variation; the diel cycle was clearly driven by light and thus reflected the metabolism of the lacustrine ecosystem. However, the diel cycle was sometimes blurred by injection of hypolimnetic water rich in CO2 and the surface water CO2 concentration was thus controlled by stratification dynamics. The highest CO2 fluxes were measured in spring, autumn and during those hypolimnetic injections causing bursts of CO2 comparable with the spring and autumn fluxes. The annual fluxes averaged 77 (±11 SD) g C m-2 yr-1. In estimating the importance of the lake in recycling terrestrial carbon, the flux was normalized to the catchment area and this normalized flux was compared with net ecosystem production estimates of -50 to 200 g C m-2 yr-1 from unmanaged forests in corresponding temperature and precipitation regimes in the literature. Within this range the flux of Lake Valkea-Kotinen yielded from the increase in source of the surrounding forest by 20% to decrease in sink by 5%. The free water approach gave primary production and community respiration estimates of 5- and 16-fold, respectively, compared with traditional bottle incubations during a 5-day testing period in autumn. The results are in parallel with findings in the literature. Both methods adopted from the terrestrial community also proved useful in lake studies. A large percentage of the EC data was rejected, due to the unfulfilled prerequisites of the method. However, the amount of data accepted remained large compared with what would be feasible with traditional methods. Use of the EC method revealed underestimation of the widely used gas exchange model and suggests simultaneous measurements of actual turbulence at the water surface with comparison of the different gas flux methods to revise the parameterization of the gas transfer velocity used in the models.