Dynamical processes in the stratosphere and upper troposphere and their influence on the distribution of trace gases in the polar atmosphere

Transport plays an important role in the distribution of long-lived gases such as ozone and water vapour in the atmosphere. Understanding of observed variability in these gases as well as prediction of the future changes depends therefore on our knowledge of the relevant atmospheric dynamics. This dissertation studies certain dynamical processes in the stratosphere and upper troposphere which influence the distribution of ozone and water vapour in the atmosphere. The planetary waves that originate in the troposphere drive the stratospheric circulation. They influence both the meridional transport of substances as well as parameters of the polar vortices. In turn, temperatures inside the polar vortices influence abundance of the Polar Stratospheric Clouds (PSC) and therefore the chemical ozone destruction. Wave forcing of the stratospheric circulation is not uniform during winter. The November-December averaged stratospheric eddy heat flux shows a significant anticorrelation with the January-February averaged eddy heat flux in the midlatitude stratosphere and troposphere. These intraseasonal variations are attributable to the internal stratospheric vacillations. In the period 1979-2002, the wave forcing exhibited a negative trend which was confined to the second half of winter only. In the period 1958-2002, area, strength and longevity of the Arctic polar vortices do not exhibit significant long-term changes while the area with temperatures lower than the threshold temperature for PSC formation shows statistically significant increase. However, the Arctic vortex parameters show significant decadal changes which are mirrored in the ozone variability. Monthly ozone tendencies in the Northern Hemisphere show significant correlations (|r|=0.7) with proxies of the stratospheric circulation. In the Antarctic, the springtime vortex in the lower stratosphere shows statistically significant trends in temperature, longevity and strength (but not in area) in the period 1979-2001. Analysis of the ozone and water vapour vertical distributions in the Arctic UTLS shows that layering below and above the tropopause is often associated with poleward Rossby wave-breaking. These observations together with calculations of cross-tropopause fluxes emphasize the importance of poleward Rossby wave breaking for the stratosphere-troposphere exchange in the Arctic.

Studies of Earth Dynamics with the Superconducting Gravimeter

The superconducting (or cryogenic) gravimeter (SG) is based on the levitation of a super­conducting sphere in a stable magnetic field created by current in superconducting coils. Depending on frequency, it is capable of detecting gravity variations as small as 10-11ms-2. For a single event, the detection threshold is higher, conservatively about 10-9 ms-2. Due to its high sensitivity and low drift rate, the SG is eminently suitable for the study of geodynamical phenomena through their gravity signatures. I present investigations of Earth dynamics with the superconducting gravimeter GWR T020 at Metsähovi from 1994 to 2005. The history and key technical details of the installation are given. The data processing methods and the development of the local tidal model at Metsähovi are presented. The T020 is a part of the worldwide GGP (Global Geodynamics Project) network, which consist of 20 working station. The data of the T020 and of other participating SGs are available to the scientific community. The SG T020 have used as a long-period seismometer to study microseismicity and the Earth s free oscillation. The annual variation, spectral distribution, amplitude and the sources of microseism at Metsähovi were presented. Free oscillations excited by three large earthquakes were analyzed: the spectra, attenuation and rotational splitting of the modes. The lowest modes of all different oscillation types are studied, i.e. the radial mode 0S0, the "football mode" 0S2, and the toroidal mode 0T2. The very low level (0.01 nms-1) incessant excitation of the Earth s free oscillation was detected with the T020. The recovery of global and regional variations in gravity with the SG requires the modelling of local gravity effects. The most important of them is hydrology. The variation in the groundwater level at Metsähovi as measured in a borehole in the fractured bedrock correlates significantly (0.79) with gravity. The influence of local precipitation, soil moisture and snow cover are detectable in the gravity record. The gravity effect of the variation in atmospheric mass and that of the non-tidal loading by the Baltic Sea were investigated together, as sea level and air pressure are correlated. Using Green s functions it was calculated that a 1 metre uniform layer of water in the Baltic Sea increases the gravity at Metsähovi by 31 nms-2 and the vertical deformation is -11 mm. The regression coefficient for sea level is 27 nms-2m-1, which is 87% of the uniform model. These studies are associated with temporal height variations using the GPS data of Metsähovi permanent station. Results of long time series at Metsähovi demonstrated high quality of data and correctly carried out offsets and drift corrections. The superconducting gravimeter T020 has been proved to be an eminent and versatile tool in studies of the Earth dynamics.

Computational fluid dynamics simulations in aerosol and nucleation studies

Nucleation is the first step in the formation of a new phase inside a mother phase. Two main forms of nucleation can be distinguished. In homogeneous nucleation, the new phase is formed in a uniform substance. In heterogeneous nucleation, on the other hand, the new phase emerges on a pre-existing surface (nucleation site). Nucleation is the source of about 30% of all atmospheric aerosol which in turn has noticeable health effects and a significant impact on climate. Nucleation can be observed in the atmosphere, studied experimentally in the laboratory and is the subject of ongoing theoretical research. This thesis attempts to be a link between experiment and theory. By comparing simulation results to experimental data, the aim is to (i) better understand the experiments and (ii) determine where the theory needs improvement. Computational fluid dynamics (CFD) tools were used to simulate homogeneous onecomponent nucleation of n-alcohols in argon and helium as carrier gases, homogeneous nucleation in the water-sulfuric acid-system, and heterogeneous nucleation of water vapor on silver particles. In the nucleation of n-alcohols, vapor depletion, carrier gas effect and carrier gas pressure effect were evaluated, with a special focus on the pressure effect whose dependence on vapor and carrier gas properties could be specified. The investigation of nucleation in the water-sulfuric acid-system included a thorough analysis of the experimental setup, determining flow conditions, vapor losses, and nucleation zone. Experimental nucleation rates were compared to various theoretical approaches. We found that none of the considered theoretical descriptions of nucleation captured the role of water in the process at all relative humidities. Heterogeneous nucleation was studied in the activation of silver particles in a TSI 3785 particle counter which uses water as its working fluid. The role of the contact angle was investigated and the influence of incoming particle concentrations and homogeneous nucleation on counting efficiency determined.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Dynamics of contour, object and face processing in the human visual cortex

The neural basis of visual perception can be understood only when the sequence of cortical activity underlying successful recognition is known. The early steps in this processing chain, from retina to the primary visual cortex, are highly local, and the perception of more complex shapes requires integration of the local information. In Study I of this thesis, the progression from local to global visual analysis was assessed by recording cortical magnetoencephalographic (MEG) responses to arrays of elements that either did or did not form global contours. The results demonstrated two spatially and temporally distinct stages of processing: The first, emerging 70 ms after stimulus onset around the calcarine sulcus, was sensitive to local features only, whereas the second, starting at 130 ms across the occipital and posterior parietal cortices, reflected the global configuration. To explore the links between cortical activity and visual recognition, Studies II III presented subjects with recognition tasks of varying levels of difficulty. The occipito-temporal responses from 150 ms onwards were closely linked to recognition performance, in contrast to the 100-ms mid-occipital responses. The averaged responses increased gradually as a function of recognition performance, and further analysis (Study III) showed the single response strengths to be graded as well. Study IV addressed the attention dependence of the different processing stages: Occipito-temporal responses peaking around 150 ms depended on the content of the visual field (faces vs. houses), whereas the later and more sustained activity was strongly modulated by the observers attention. Hemodynamic responses paralleled the pattern of the more sustained electrophysiological responses. Study V assessed the temporal processing capacity of the human object recognition system. Above sufficient luminance, contrast and size of the object, the processing speed was not limited by such low-level factors. Taken together, these studies demonstrate several distinct stages in the cortical activation sequence underlying the object recognition chain, reflecting the level of feature integration, difficulty of recognition, and direction of attention.

The ozone transfer between atmosphere and vegetation. A study on Scots pine in the field

Ozone (O3) is a reactive gas present in the troposphere in the range of parts per billion (ppb), i.e. molecules of O3 in 109 molecules of air. Its strong oxidative capacity makes it a key element in tropospheric chemistry and a threat to the integrity of materials, including living organisms. Knowledge and control of O3 levels are an issue in relation to indoor air quality, building material endurance, respiratory human disorders, and plant performance. Ozone is also a greenhouse gas and its abundance is relevant to global warming. The interaction of the lower troposphere with vegetated landscapes results in O3 being removed from the atmosphere by reactions that lead to the oxidation of plant-related components. Details on the rate and pattern of removal on different landscapes as well as the ultimate mechanisms by which this occurs are not fully resolved. This thesis analysed the controlling processes of the transfer of ozone at the air-plant interface. Improvement in the knowledge of these processes benefits the prediction of both atmospheric removal of O3 and its impact on vegetation. This study was based on the measurement and analysis of multi-year field measurements of O3 flux to Scots pine (Pinus sylvestris L.) foliage with a shoot-scale gas-exchange enclosure system. In addition, the analyses made use of simultaneous CO2 and H2O exchange, canopy-scale O3, CO2 and H2O exchange, foliage surface wetness, and environmental variables. All data was gathered at the SMEAR measuring station (southern Finland). Enclosure gas-exchange techniques such as those commonly used for the measure of CO2 and water vapour can be applied to the measure of ozone gas-exchange in the field. Through analysis of the system dynamics the occurring disturbances and noise can be identified. In the system used in this study, the possible artefacts arising from the ozone reactivity towards the system materials in combination with low background concentrations need to be taken into account. The main artefact was the loss of ozone towards the chamber walls, which was found to be very variable. The level of wall-loss was obtained from simultaneous and continuous measurements, and was included in the formulation of the mass balance of O3 concentration inside the chamber. The analysis of the field measurements in this study show that the flux of ozone to the Scots pine foliage is generated in about equal proportions by stomatal and non-stomatal controlled processes. Deposition towards foliage and forest is sustained also during night and winter when stomatal gas-exchange is low or absent. The non-stomatal portion of the flux was analysed further. The pattern of flux in time was found to be an overlap of the patterns of biological activity and presence of wetness in the environment. This was seen to occur both at the shoot and canopy scale. The presence of wetness enhanced the flux not only in the presence of liquid droplets but also during existence of a moisture film on the plant surfaces. The existence of these films and their relation to the ozone sinks was determined by simultaneous measurements of leaf surface wetness and ozone flux. The results seem to suggest ozone would be reacting at the foliage surface and the reaction rate would be mediated by the presence of surface wetness. Alternative mechanisms were discussed, including nocturnal stomatal aperture and emission of reactive volatile compounds. The prediction of the total flux could thus be based on a combination of a model of stomatal behaviour and a model of water absorption on the foliage surfaces. The concepts behind the division of stomatal and non-stomatal sinks were reconsidered. This study showed that it is theoretically possible that a sink located before or near the stomatal aperture prevents or diminishes the diffusion of ozone towards the intercellular air space of the mesophyll. This obstacle to stomatal diffusion happens only under certain conditions, which include a very low presence of reaction sites in the mesophyll, an extremely strong sink located on the outer surfaces or stomatal pore. The relevance, or existence, of this process in natural conditions would need to be assessed further. Potentially strong reactions were considered, including dissolved sulphate, volatile organic compounds, and apoplastic ascorbic acid. Information on the location and the relative abundance of these compounds would be valuable. The highest total flux towards the foliage and forest happens when both the plant activity and ambient moisture are high. The highest uptake into the interior of the foliage happens at large stomatal apertures, provided that scavenging reactions located near the stomatal pore are weak or non-existent. The discussion covers the methodological developments of this study, the relevance of the different controlling factors of ozone flux, the partition amongst its component, and the possible mechanisms of non-stomatal uptake.

Tree mortality and deadwood dynamics in late-successional boreal forests

Here I aimed at quantifying the main components of deadwood dynamics, i.e. tree mortality, deadwood pools, and their decomposition, in late-successional boreal forests. I focused on standing dead trees in three stand types dominated by Picea mariana and Abies balsamea in eastern Canada, and on standing and down dead trees in Picea abies-dominated stands in three areas in Northern Europe. Dead and living trees were measured on five sample plots of 1.6-ha size in each study area and stand type. Stem disks from dead trees were sampled to determine wood density and year of death, using dendrochronological methods. The results were applied to reconstruct past tree mortality and to model deadwood decay class dynamics. Site productivity, stand developmental stage, and the occurrence of episodic tree mortality influenced deadwood volume and quality. In all study areas tree mortality was continuous, leading to continuity in deadwood decay stage distribution. Episodic tree mortality due to either autogenic or allogenic causes influenced deadwood volume and quality in all but one study area. However, regardless of productivity and disturbance history deadwood was abundant, accounting for 20 53% of total wood volume in European study areas, and 15 27% of total standing volume in eastern Canada. Deadwood was a persistent structural component, since its expected residence time in early- and midstages of decay was 18 yr even in the area with the most rapid decomposition. The results indicated that in the absence of episodic tree mortality, stands may eventually develop to a steady state, in which deadwood volume fluctuates around an equilibrium state. However, in many forests deadwood is naturally variable, due to recurrent moderate-severity disturbances. This variability, the continuous tree mortality, and variation in rates of wood decomposition determine the dynamics and availability of deadwood as a habitat and carbon storage medium in boreal coniferous forest ecosystems.

Fen ecosystem carbon gas dynamics in changing hydrological conditions

Northern peatlands are thought to store one third of all soil carbon (C). Besides the C sink function, peatlands are one of the largest natural sources of methane (CH4) to the atmosphere. Climate change may affect the C gas dynamics as well as the labile C pool. Because the peatland C sequestration and CH4 emissions are governed by high water levels, changes in hydrology are seen as the driving factor in peatland ecosystem change. This study aimed to quantify the carbon dioxide (CO2) and CH4 dynamics of a fen ecosystem at different spatial scales: plant community components scale, plant community scale and ecosystem scale, under hydrologically normal and water level drawdown conditions. C gas exchange was measured in two fens in southern Finland applying static chamber and eddy covariance techniques. During hydrologically normal conditions, the ecosystem was a CO2 sink and CH4 source to the atmosphere. Sphagnum mosses and sedges were the most important contributors to the community photosynthesis. The presence of sedges had a major positive impact on CH4 emissions while dwarf shrubs had a slightly attenuating impact. C fluxes varied considerably between the plant communities. Therefore, their proportions determined the ecosystem scale fluxes. An experimental water level drawdown markedly reduced the photosynthesis and respiration of sedges and Sphagnum mosses and benefited shrubs. Consequently, changes were smaller at the ecosystem scale than at the plant group scale. The decrease in photosynthesis and the increase in respiration, mostly peat respiration, made the fen a smaller CO2 sink. CH4 fluxes were significantly lowered, close to zero. The impact of natural droughts was similar to, although more modest than, the impact of the experimental water level drawdown. The results are applicable to the short term impacts of the water level drawdown and to climatic conditions in which droughts become more frequent.

Stand structural dynamics on pristine and managed boreal peatlands

The objectives of this study were to investigate the stand structure and succession dynamics in Scots pine (Pinus sylvestris L.) stands on pristine peatlands and in Scots pine and Norway spruce (Picea abies (L.) Karst.) dominated stands on drained peatlands. Furthermore, my focus was on characterising how the inherent and environmental factors and the intermediate thinnings modify the stand structure and succession. For pristine peatlands, the study was based on inventorial stand data, while for drained peatlands, longitudinal data from repeatedly measured stands were utilised. The studied sites covered the most common peatland site types in Finland. They were classified into two categories according to the ecohydrological properties related to microsite variation and nutrient levels within sites. Tree DBH and age distributions in relation to climate and site type were used to study the stand dynamics on pristine sites. On drained sites, the Weibull function was used to parameterise the DBH distributions and mixed linear models were constructed to characterise the impacts of different ecological factors on stand dynamics. On pristine peatlands, both climate and the ecohydrology of the site proved to be crucial factors determining the stand structure and its dynamics. Irrespective of the vegetation succession, enhanced site productivity and increased stand stocking they significantly affected the stand dynamics also on drained sites. On the most stocked sites on pristine peatlands the inter-tree competition seemed to also be a significant factor modifying stand dynamics. Tree age and size diversity increased with stand age, but levelled out in the long term. After drainage, the stand structural unevenness increased due to the regeneration and/or ingrowth of the trees. This increase was more pronounced on sparsely forested composite sites than on more fully stocked genuine forested sites in Scots pine stands, which further undergo the formation of birch and spruce undergrowth beneath the overstory as succession proceeds. At 20-30 years after drainage the structural heterogeneity started to decrease, indicating increased inter-tree competition, which increased the mortality of suppressed trees within stand. Peatland stands are more dynamic than anticipated and are generally not characterized by a balanced, self-perpetuating structure. On pristine sites, various successional pathways are possible, whereas on drained sites the succession has more uniform trend. Typically, stand succession proceeds without any distinct developmental stages on pristine peatlands, whereas on drained peatlands, at least three distinct stages could be identified. Thinnings had only little impact on the stand succession. The new information on stand dynamics may be utilised, e.g. in forest management planning to facilitate the allocation of the growth resources to the desired crop component by appropriate silvicultural treatments, as well as assist in assessing the effects of the climate change on the forested boreal peatlands.