Phosphorus retention in forest soils and the functioning of buffer zones used in forestry

Phosphorus (P) retention properties of soils typical for boreal forest, i.e. podzolic soil and peat soils, vary significantly, but the range of this variation has not been sufficiently documented. To assess the usefulness of buffer zones used in forestry in removing P from the discharge by chemical sorption in soil, and to estimate the risk of P leaching after forestry operations, more data is needed on soil P retention properties. P retention properties of soils were studied at clear-cut areas, unharvested buffer zones adjoining the clear-cut and at peatland buffer zone areas. Desorption-sorption isotherms were determined for the humus layer, the mineral soil horizons E, B and C of the Podzol profile and for the surface layer peat (0-15 cm) and the subsurface layer peat (15-30 cm). The efficiency of buffer zones in retaining P was studied at six peatland buffer zone areas by adding P-containing solute in the inflow. A tracer study was conducted at one of the buffer zone areas to determine the allocation of the added P in soil and vegetation. Measured sorption or desorption rather than parameter values of fitted sorption equations described P desorption and sorption behaviour in soil. The highest P retention efficiency was in the B horizon and consequently, if contact occurred or was established between the soluble P in the water and the soil B horizon, the risk of P leaching was low. Humus layer was completely incapable of retaining P after clear-cutting. In the buffer zones, the decrease in P retention properties in the humus layer and the low amount of P sorbed by it indicated that the importance of the layer in the functioning of buffer zones is low. The peatland buffer zone areas were efficient in retaining soluble P from inflow. P sorption properties of the peat soil at the buffer zone areas varied largely but the contribution of P sorption in the peat was particularly important during high flow in spring, when the vegetation was not fully developed. Factors contributing to efficient P retention were large buffer size and low hydrological load whereas high hydrological load combined with the formation of preferential flow paths, especially during early spring or late autumn was disadvantageous. However, small buffer zone areas, too, may be efficient in reducing P load.

Phosphorus in agricultural soils of Finland - characterization of reserves and retention in mineral soil profiles

An overwhelming majority of all the research on soil phosphorus (P) has been carried out with soil samples taken from the surface soils only, and our understanding of the forms and the reactions of P at a soil profile scale is based on few observations. In Finland, the interest in studying the P in complete soil profiles has been particularly small because of the lack of tradition in studying soil genesis, morphology, or classification. In this thesis, the P reserves and the retention of orthophosphate phosphorus (PO4-P) were examined in four cultivated mineral soil profiles in Finland (three Inceptisols and one Spodosol). The soils were classified according to the U.S. Soil Taxonomy and soil samples were taken from the genetic horizons in the profiles. The samples were analyzed for total P concentration, Chang and Jackson P fractions, P sorption properties, concentrations of water-extractable P, and for concentrations of oxalate-extractable Al and Fe. Theoretical P sorption capacities and degrees of P saturation were calculated with the data from the oxalate-extractions and the P fractionations. The studied profiles can be divided into sections with clearly differing P characteristics by their master horizons Ap, B and C. The C (or transitional BC) horizons below an approximate depth of 70 cm were dominated by, assumingly apatitic, H2SO4-soluble P. The concentration of total P in the C horizons ranged from 729 to 810 mg kg-1. In the B horizons between the depths of 30 and 70 cm, a significant part of the primary acid-soluble P has been weathered and transformed to secondary P forms. A mean weathering rate of the primary P in the soils was estimated to vary between 230 and 290 g ha-1 year-1. The degrees of P saturation in the B and C horizons were smaller than 7%, and the solubility of PO4-P was negligible. The P conditions in the Ap horizons differed drastically from those in the subsurface horizons. The high concentrations of total P (689-1870 mg kg-1) in the Ap horizons are most likely attributable to long-term cultivation with positive P balances. A significant proportion of the P in the Ap horizons occurred in the NH4F- and NaOH-extractable forms and as organic P. These three P pools, together with the concentrations of oxalate-extractable Al and Fe, seem to control the dynamics of PO4-P in the soils. The degrees of P saturation in the Ap horizons were greater (8-36%) than in the subsurface horizons. This was also reflected in the sorption experiments: Only the Ap horizons were able to maintain elevated PO4-P concentrations in the solution phase − all the subsoil horizons acted as sinks for PO4-P. Most of the available sorption capacity in the soils is located in the B horizons. The results suggest that this capacity could be utilized in reducing the losses of soluble P from excessively fertilized soils by mixing highly sorptive material from the B horizons with the P-enriched surface soil. The drastic differences in the P characteristics observed between adjoining horizons have to be taken into consideration when conducting soil sampling. Sampling of subsoils has to be made according to the genetic horizons or at small depth increments. Otherwise, contrasting materials are likely to be mixed in the same sample; and the results of such samples are not representative of any material present in the studied profile. Air-drying of soil samples was found to alter the results of the sorption experiments and the water extractions. This indicates that the studies on the most labile P forms in soil should be carried out with moist samples.

Modeling Asymmetries in Financial Data with Multiplicative Error Models

This thesis addresses modeling of financial time series, especially stock market returns and daily price ranges. Modeling data of this kind can be approached with so-called multiplicative error models (MEM). These models nest several well known time series models such as GARCH, ACD and CARR models. They are able to capture many well established features of financial time series including volatility clustering and leptokurtosis. In contrast to these phenomena, different kinds of asymmetries have received relatively little attention in the existing literature. In this thesis asymmetries arise from various sources. They are observed in both conditional and unconditional distributions, for variables with non-negative values and for variables that have values on the real line. In the multivariate context asymmetries can be observed in the marginal distributions as well as in the relationships of the variables modeled. New methods for all these cases are proposed. Chapter 2 considers GARCH models and modeling of returns of two stock market indices. The chapter introduces the so-called generalized hyperbolic (GH) GARCH model to account for asymmetries in both conditional and unconditional distribution. In particular, two special cases of the GARCH-GH model which describe the data most accurately are proposed. They are found to improve the fit of the model when compared to symmetric GARCH models. The advantages of accounting for asymmetries are also observed through Value-at-Risk applications. Both theoretical and empirical contributions are provided in Chapter 3 of the thesis. In this chapter the so-called mixture conditional autoregressive range (MCARR) model is introduced, examined and applied to daily price ranges of the Hang Seng Index. The conditions for the strict and weak stationarity of the model as well as an expression for the autocorrelation function are obtained by writing the MCARR model as a first order autoregressive process with random coefficients. The chapter also introduces inverse gamma (IG) distribution to CARR models. The advantages of CARR-IG and MCARR-IG specifications over conventional CARR models are found in the empirical application both in- and out-of-sample. Chapter 4 discusses the simultaneous modeling of absolute returns and daily price ranges. In this part of the thesis a vector multiplicative error model (VMEM) with asymmetric Gumbel copula is found to provide substantial benefits over the existing VMEM models based on elliptical copulas. The proposed specification is able to capture the highly asymmetric dependence of the modeled variables thereby improving the performance of the model considerably. The economic significance of the results obtained is established when the information content of the volatility forecasts derived is examined.

Modeling Efficient Classification as a Process of Confidence Assessment and Delegation

In visual object detection and recognition, classifiers have two interesting characteristics: accuracy and speed. Accuracy depends on the complexity of the image features and classifier decision surfaces. Speed depends on the hardware and the computational effort required to use the features and decision surfaces. When attempts to increase accuracy lead to increases in complexity and effort, it is necessary to ask how much are we willing to pay for increased accuracy. For example, if increased computational effort implies quickly diminishing returns in accuracy, then those designing inexpensive surveillance applications cannot aim for maximum accuracy at any cost. It becomes necessary to find trade-offs between accuracy and effort. We study efficient classification of images depicting real-world objects and scenes. Classification is efficient when a classifier can be controlled so that the desired trade-off between accuracy and effort (speed) is achieved and unnecessary computations are avoided on a per input basis. A framework is proposed for understanding and modeling efficient classification of images. Classification is modeled as a tree-like process. In designing the framework, it is important to recognize what is essential and to avoid structures that are narrow in applicability. Earlier frameworks are lacking in this regard. The overall contribution is two-fold. First, the framework is presented, subjected to experiments, and shown to be satisfactory. Second, certain unconventional approaches are experimented with. This allows the separation of the essential from the conventional. To determine if the framework is satisfactory, three categories of questions are identified: trade-off optimization, classifier tree organization, and rules for delegation and confidence modeling. Questions and problems related to each category are addressed and empirical results are presented. For example, related to trade-off optimization, we address the problem of computational bottlenecks that limit the range of trade-offs. We also ask if accuracy versus effort trade-offs can be controlled after training. For another example, regarding classifier tree organization, we first consider the task of organizing a tree in a problem-specific manner. We then ask if problem-specific organization is necessary.

Modeling ecological and evolutionary processes in spatially structured populations

Many species inhabit fragmented landscapes, resulting either from anthropogenic or from natural processes. The ecological and evolutionary dynamics of spatially structured populations are affected by a complex interplay between endogenous and exogenous factors. The metapopulation approach, simplifying the landscape to a discrete set of patches of breeding habitat surrounded by unsuitable matrix, has become a widely applied paradigm for the study of species inhabiting highly fragmented landscapes. In this thesis, I focus on the construction of biologically realistic models and their parameterization with empirical data, with the general objective of understanding how the interactions between individuals and their spatially structured environment affect ecological and evolutionary processes in fragmented landscapes. I study two hierarchically structured model systems, which are the Glanville fritillary butterfly in the Åland Islands, and a system of two interacting aphid species in the Tvärminne archipelago, both being located in South-Western Finland. The interesting and challenging feature of both study systems is that the population dynamics occur over multiple spatial scales that are linked by various processes. My main emphasis is in the development of mathematical and statistical methodologies. For the Glanville fritillary case study, I first build a Bayesian framework for the estimation of death rates and capture probabilities from mark-recapture data, with the novelty of accounting for variation among individuals in capture probabilities and survival. I then characterize the dispersal phase of the butterflies by deriving a mathematical approximation of a diffusion-based movement model applied to a network of patches. I use the movement model as a building block to construct an individual-based evolutionary model for the Glanville fritillary butterfly metapopulation. I parameterize the evolutionary model using a pattern-oriented approach, and use it to study how the landscape structure affects the evolution of dispersal. For the aphid case study, I develop a Bayesian model of hierarchical multi-scale metapopulation dynamics, where the observed extinction and colonization rates are decomposed into intrinsic rates operating specifically at each spatial scale. In summary, I show how analytical approaches, hierarchical Bayesian methods and individual-based simulations can be used individually or in combination to tackle complex problems from many different viewpoints. In particular, hierarchical Bayesian methods provide a useful tool for decomposing ecological complexity into more tractable components.

Short-term effects of variable retention on epigaeic spiders and carabid beetles in Finland

Forestry has influenced forest dwelling organisms for centuries in Fennoscandia. For example, in Finland ca. 30% of the threatened species are threatened because of forestry. Nowadays forest management recommendations include practices aimed at maintaining biodiversity in harvesting, such as green-tree retention. However, the effects of these practices have been little studied. In variable retention, different numbers of trees are retained, varying from green-tree retention (at least a few live standing trees in clear-cuts) to thinning (only individual trees removed). I examined the responses of ground-dwelling spiders and carabid beetles to green-tree retention (with small and large tree groups), gap felling and thinning aimed at an uneven age structure of trees. The impacts of these harvesting methods were compared to those of clear-cutting and uncut controls. I aimed to test the hypothesis that retaining more trees positively affects populations of those species of spiders and carabids that were present before harvesting. The data come from two studies. First, spiders were collected with pitfall traps in south-central Finland in 1995 (pre-treatment) and 1998 (after-treatment) in order to examine the effects of clear-cutting, green-tree retention (with 0.01-0.02-ha sized tree groups), gap felling (with three 0.16-ha sized openings in a 1-ha stand), thinning aiming at an uneven age structure of trees and uncut control. Second, spiders and carabids were caught with pitfall traps in eastern Finland in 1998-2001 (pre-treatment and three post-treatment years) in eleven 0.09-0.55-ha sized retention-tree groups and clear-cuts adjacent to them. Original spider and carabid assemblages were better maintained after harvests that retained more trees. Thinning maintained forest spiders well. However, gap felling and large retention-tree groups maintained some forest spider and carabid species in the short-term, but negatively affected some species over time. However, use of small retention-tree groups was associated with negative effects on forest spider populations. Studies are needed on the long-term effects of variable retention on terrestrial invertebrates; especially those directed at defining appropriate retention patch size and on the importance of structural diversity provided by variable retention for invertebrate populations. However, the aims of variable retention should be specified first. For example, are retention-tree groups planned to constitute life-boats , stepping-stones or to create structural diversity? Does it suffice that some species are maintained, or do we want to preserve the most sensitive ones, and how are these best defined? Moreover, the ecological benefits and economic costs of modified logging methods should be compared to other approaches aimed at maintaining biodiversity.

Theoretical modeling of ionospheric electrodynamics including induction effects

This thesis deals with theoretical modeling of the electrodynamics of auroral ionospheres. In the five research articles forming the main part of the thesis we have concentrated on two main themes: Development of new data-analysis techniques and study of inductive phenomena in the ionospheric electrodynamics. The introductory part of the thesis provides a background for these new results and places them in the wider context of ionospheric research. In this thesis we have developed a new tool (called 1D SECS) for analysing ground based magnetic measurements from a 1-dimensional magnetometer chain (usually aligned in the North-South direction) and a new method for obtaining ionospheric electric field from combined ground based magnetic measurements and estimated ionospheric electric conductance. Both these methods are based on earlier work, but contain important new features: 1D SECS respects the spherical geometry of large scale ionospheric electrojet systems and due to an innovative way of implementing boundary conditions the new method for obtaining electric fields can be applied also at local scale studies. These new calculation methods have been tested using both simulated and real data. The tests indicate that the new methods are more reliable than the previous techniques. Inductive phenomena are intimately related to temporal changes in electric currents. As the large scale ionospheric current systems change relatively slowly, in time scales of several minutes or hours, inductive effects are usually assumed to be negligible. However, during the past ten years, it has been realised that induction can play an important part in some ionospheric phenomena. In this thesis we have studied the role of inductive electric fields and currents in ionospheric electrodynamics. We have formulated the induction problem so that only ionospheric electric parameters are used in the calculations. This is in contrast to previous studies, which require knowledge of the magnetospheric-ionosphere coupling. We have applied our technique to several realistic models of typical auroral phenomena. The results indicate that inductive electric fields and currents are locally important during the most dynamical phenomena (like the westward travelling surge, WTS). In these situations induction may locally contribute up to 20-30% of the total ionospheric electric field and currents. Inductive phenomena do also change the field-aligned currents flowing between the ionosphere and magnetosphere, thus modifying the coupling between the two regions.