Space partitioning for scalable K-means

K-Means is a popular clustering algorithm which adopts an iterative refinement procedure to determine data partitions and to compute their associated centres of mass, called centroids. The straightforward implementation of the algorithm is often referred to as `brute force' since it computes a proximity measure from each data point to each centroid at every iteration of the K-Means process. Efficient implementations of the K-Means algorithm have been predominantly based on multi-dimensional binary search trees (KD-Trees). A combination of an efficient data structure and geometrical constraints allow to reduce the number of distance computations required at each iteration. In this work we present a general space partitioning approach for improving the efficiency and the scalability of the K-Means algorithm. We propose to adopt approximate hierarchical clustering methods to generate binary space partitioning trees in contrast to KD-Trees. In the experimental analysis, we have tested the performance of the proposed Binary Space Partitioning K-Means (BSP-KM) when a divisive clustering algorithm is used. We have carried out extensive experimental tests to compare the proposed approach to the one based on KD-Trees (KD-KM) in a wide range of the parameters space. BSP-KM is more scalable than KDKM, while keeping the deterministic nature of the `brute force' algorithm. In particular, the proposed space partitioning approach has shown to overcome the well-known limitation of KD-Trees in high-dimensional spaces and can also be adopted to improve the efficiency of other algorithms in which KD-Trees have been used.

Peroxy radical partitioning during the AMMA radical intercomparison exercise

Peroxy radicals were measured onboard two scientific aircrafts during the AMMA (African Monsoon Multidisciplinary Analysis) campaign in summer 2006. This paper reports results from the flight on 16 August 2006 during which measurements of HO2 by laser induced fluorescence spectroscopy at low pressure (LIF-FAGE) and total peroxy radicals (RO2* = HO2+ΣRO2, R = organic chain) by two similar instruments based on the peroxy radical chemical amplification (PeRCA) technique were subject of a blind intercomparison. The German DLR-Falcon and the British FAAM-BAe-146 flew wing tip to wing tip for about 30 min making concurrent measurements on 2 horizontal level runs at 697 and 485 hPa over the same geographical area in Burkina Faso. A full set of supporting measurements comprising photolysis frequencies, and relevant trace gases like CO, NO, NO2, NOy, O3 and a wider range of VOCs were collected simultaneously. Results are discussed on the basis of the characteristics and limitations of the different instruments used. Generally, no data bias are identified and the RO2* data available agree quite reasonably within the instrumental errors. The [RO2*]/[HO2] ratios, which vary between 1:1 and 3:1, as well as the peroxy radical variability, concur with variations in photolysis rates and in other potential radical precursors. Model results provide additional information about dominant radical formation and loss processes.

Free atmospheric CO2 enrichment increased above ground biomass but did not affect symbiotic N2-fixation and soil carbon dynamics in a mixed deciduous stand in Wales

Through increases in net primary production (NPP), elevated CO2 is hypothesizes to increase the amount of plant litter entering the soil. The fate of this extra carbon on the forest floor or in mineral soil is currently not clear. Moreover, increased rates of NPP can be maintained only if forests can escape nitrogen limitation. In a Free atmospheric CO2 Enrichment (FACE) experiment near Bangor, Wales, 4 ambient CO2 and 4 FACE plots were planted with patches of Betula pendula, Alnus glutinosa and Fagus sylvatica on a former arable field. Four years after establishment, only a shallow L forest floor litter layer had formed due to intensive bioturbation. Total soil C and N contents increased irrespective of treatment and species as a result of afforestation. We could not detect an additional C sink in the soil, nor were soil C stabilization processes affected by FACE. We observed a decrease of leaf N content in Betula and Alnus under FACE, while the soil C/N ratio decreased regardless of CO2 treatment. The ratio of N taken up from the soil and by N2-fixation in Alnus was not affected by FACE. We infer that increased nitrogen use efficiency is the mechanism by which increased NPP is sustained under elevated CO2 at this site.

Individual biomass factors for beech, oak and pine in Slovakia: a comparative study in young naturally regenerated stands

Biomass conversion and expansion factors (BCEF) which convert tree stem volume to whole tree biomass and biomass allocation patterns in young trees were studied in order to estimate tree and stand biomass in naturally regenerated forests. European beech (Fagus sylvatica L.), Sessile oak (Quercus petraea (Mattuschka) Liebl.) and Scots pine (Pinus sylvestris L.) stands were compared. Seven forest stands of each species were chosen to cover their natural distribution in Slovakia. Species specific BCEF are presented, generally showing a steep decrease in all species in the smallest trees, with the only exception in the case of branch BCEF in beech which grows with increasing tree size. The values of BCEF for all tree compartments stabilise in all species once trees reach about 60-70mm diameter at base. As they grow larger, all species increase their allocation to stem and branches, while decreasing the relative growth of roots and foliage. There are, however, clear differences between species and also between broadleaves and conifers in biomass allocation. This research shows that species specific coefficients must be used if we are to reduce uncertainties in estimates of carbon stock changes by afforestation and reforestation activities.

Effect of soil waterlogging on below-ground biomass allometric relations in Norway spruce

An increasing importance is assigned to the estimation and verification of carbon stocks in forests. Forestry practice has several long-established and reliable methods for the assessment of aboveground biomass; however we still miss accurate predictors of belowground biomass. A major windthrow event exposing the coarse root systems of Norway spruce trees allowed us to assess the effects of contrasting soil stone and water content on belowground allocation. Increasing stone content decreases root/shoot ratio, while soil waterlogging leads to an increase in this ratio. We constructed allometric relationships for belowground biomass prediction and were able to show that only soil waterlogging significantly impacts model parameters. We showed that diameter at breast height is a reliable predictor of belowground biomass and, once site-specific parameters have been developed, it is possible to accurately estimate belowground biomass in Norway spruce.

Fine root biomass and turnover in southern taiga estimated by root inclusion nets

Fine roots play an important part in forest carbon, nutrient and water cycles. The turnover of fine roots constitutes a major carbon input to soils. Estimation of fine root turnover is difficult, labour intensive and is often compounded by artefacts created by soil disturbance. In this work, an alternative approach of using inclusion nets installed in an undisturbed soil profile was used to measure fine root production and was compared to the in-growth core method. There was no difference between fine root production estimated by the two methods in three southern taiga sites with contrasting soil conditions and tree species composition in the Central Forest State Biosphere Reserve, Russia. Expressed as annual production over standing biomass, Norway spruce fine root turnover was in the region of 0.10 to 0.24 y-1. The inclusion net technique is suitable for field based assessment of fine root production. There are several advantages over the in-growth core method, due to non-disturbance of the soil profile and its potential for very high rate of replication.

A simple, coherent framework for partitioning uncertainty in climate predictions

A simple and coherent framework for partitioning uncertainty in multi-model climate ensembles is presented. The analysis of variance (ANOVA) is used to decompose a measure of total variation additively into scenario uncertainty, model uncertainty and internal variability. This approach requires fewer assumptions than existing methods and can be easily used to quantify uncertainty related to model-scenario interaction - the contribution to model uncertainty arising from the variation across scenarios of model deviations from the ensemble mean. Uncertainty in global mean surface air temperature is quantified as a function of lead time for a subset of the Coupled Model Intercomparison Project phase 3 ensemble and results largely agree with those published by other authors: scenario uncertainty dominates beyond 2050 and internal variability remains approximately constant over the 21st century. Both elements of model uncertainty, due to scenario-independent and scenario-dependent deviations from the ensemble mean, are found to increase with time. Estimates of model deviations that arise as by-products of the framework reveal significant differences between models that could lead to a deeper understanding of the sources of uncertainty in multi-model ensembles. For example, three models are shown diverging pattern over the 21st century, while another model exhibits an unusually large variation among its scenario-dependent deviations.

Belowground biomass functions and expansion factors in high elevation Norway spruce

Biomass allocation to above- and belowground compartments in trees is thought to be affected by growth conditions. To assess the strength of such influences, we sampled six Norway spruce forest stands growing at higher altitudes. Within these stands, we randomly selected a total of 77 Norway spruce trees and measured volume and biomass of stem, above- and belowground stump and all roots over 0.5 cm diameter. A comparison of our observations with models parameterised for lower altitudes shows that models developed for specific conditions may be applicable to other locations. Using our observations, we developed biomass functions (BF) and biomass conversion and expansion factors (BCEF) linking belowground biomass to stem parameters. While both BF and BCEF are accurate in belowground biomass predictions, using BCEF appears more promising as such factors can be readily used with existing forest inventory data to obtain estimates of belowground biomass stock. As an example, we show how BF and BCEF developed for individual trees can be used to estimate belowground biomass at the stand level. In combination with existing aboveground models, our observations can be used to quantify total standing biomass of high altitude Norway spruce stands.

Adaptation of tree growth to elevated CO2: quantitative trait loci for biomass in Populus

Information on the genetic variation of plant response to elevated CO2 (e[CO2]) is needed to understand plant adaptation and to pinpoint likely evolutionary response to future high atmospheric CO2 concentrations.• Here, quantitative trait loci (QTL) for above- and below-ground tree growth were determined in a pedigree – an F2 hybrid of poplar (Populus trichocarpa and Populus deltoides), following season-long exposure to either current day ambient CO2 (a[CO2]) or e[CO2] at 600 µl l−1, and genotype by environment interactions investigated.• In the F2 generation, both above- and below-ground growth showed a significant increase in e[CO2]. Three areas of the genome on linkage groups I, IX and XII were identified as important in determining above-ground growth response to e[CO2], while an additional three areas of the genome on linkage groups IV, XVI and XIX appeared important in determining root growth response to e[CO2].• These results quantify and identify genetic variation in response to e[CO2] and provide an insight into genomic response to the changing environment

An isotope dilution model for partitioning phenylalanine uptake by the liver of lactating dairy cows

An isotope dilution model for partitioning phenylalanine uptake by the liver of the lactating dairy cow was constructed and solved in the steady state. If assumptions are made, model solution permits calculation of the rate of phenylalanine uptake from portal vein and hepatic arterial blood supply, phenylalanine release into the hepatic vein, phenylalanine oxidation and synthesis, and degradation of hepatic constitutive and export proteins. The model requires the measurement of plasma fow rate through the liver in combination with phenylalanine concentrations and plateau isotopic enrichments in arterial, portal and hepatic plasma during a constant infusion of [1-13C]phenylalanine tracer. The model can be applied to other amino acids with similar metabolic fates and will provide a means for assessing the impact of hepatic metabolism on amino acid availability to peripheral tissues. This is of particular importance for the dairy cow when considering the requirements for milk protein synthesis and the negative environmental impact of excessive nitrogen excretion.

Mineral partitioning in milk and milk permeates at high temperature

The soluble phase of milk was separated at 20 and 80°C using ultrafiltration. The resulting permeates were then subjected to further ultrafiltration and dialysis at close to these two temperatures. It was found that pH, Ca2+ and soluble Ca decreased as the separation temperature increased both in original UF permeates and in dialysates obtained from these permeates, but P decreased only slightly. The major reason for these changes was due to the precipitation of calcium phosphate/citrate complexes onto the casein micelle with concomitant release of H+. The pH of both permeates and dialysates from milk at 20°C were slightly higher than for milk. When UF permeates collected at 20 and 80°C, were each dialysed at both these temperatures, the dialysate collected at 80°C showed much less temperature dependence for pH and ionic calcium compared with that collected at 20°C. This is in contrast to milk, which shows considerable temperature dependence for pH and ionic calcium. Further experiments revealed that the pH and Ca2+ concentration of permeates showed high temperature dependence above the temperature at which they were separated, but a much lower temperature dependence below that temperature. These findings suggest that dialysis and UF of milk at high temperature provide the best means yet for estimating the pH and ionic calcium of milk at that temperature.

Purified phenolics from hydrothermal treatments of biomass: ability to protect sunflower bulk oil and model food emulsions from oxidation

The phenolic fractions released during hydrothermal treatment of selected feedstocks (corn cobs, eucalypt wood chips, almond shells, chestnut burs, and white grape pomace) were selectively recovered by extraction with ethyl acetate and washed with ethanol/water solutions. The crude extracts were purified by a relatively simple adsorption technique using a commercial polymeric, nonionic resin. Utilization of 96% ethanol as eluting agent resulted in 47.0-72.6% phenolic desorption, yielding refined products containing 49-60% w/w phenolics (corresponding to 30-58% enrichment with respect to the crude extracts). The refined extracts produced from grape pomace and from chestnut burs were suitable for protecting bulk oil and oil-in-water and water-in-oil emulsions. A synergistic action with bovine serum albumin in the emulsions was observed.

Biomass functions and expansion factors in young Norway spruce (Picea abies [L.] Karst) trees

Roots, stems, branches and needles of 160 Norway spruce trees younger than 10 years were sampled in seven forest stands in central Slovakia in order to establish their biomassfunctions (BFs) and biomassexpansionfactors (BEFs). We tested three models for each biomass pool based on the stem base diameter, tree height and the two parameters combined. BEF values decreased for all spruce components with increasing height and diameter, which was most evident in very young trees under 1 m in height. In older trees, the values of BEFs did tend to stabilise at the height of 3–4 m. We subsequently used the BEFs to calculate dry biomass of the stands based on average stem base diameter and tree height. Total stand biomass grew with increasing age of the stands from about 1.0 Mg ha−1 at 1.5 years to 44.3 Mg ha−1 at 9.5 years. The proportion of stem and branch biomass was found to increase with age, while that of needles was fairly constant and the proportion of root biomass did decrease as the stands grew older.

RAMI4PILPS: An intercomparison of formulations for the partitioning of solar radiation in land surface models

[1] Remotely sensed, multiannual data sets of shortwave radiative surface fluxes are now available for assimilation into land surface schemes (LSSs) of climate and/or numerical weather prediction models. The RAMI4PILPS suite of virtual experiments assesses the accuracy and consistency of the radiative transfer formulations that provide the magnitudes of absorbed, reflected, and transmitted shortwave radiative fluxes in LSSs. RAMI4PILPS evaluates models under perfectly controlled experimental conditions in order to eliminate uncertainties arising from an incomplete or erroneous knowledge of the structural, spectral and illumination related canopy characteristics typical for model comparison with in situ observations. More specifically, the shortwave radiation is separated into a visible and near-infrared spectral region, and the quality of the simulated radiative fluxes is evaluated by direct comparison with a 3-D Monte Carlo reference model identified during the third phase of the Radiation transfer Model Intercomparison (RAMI) exercise. The RAMI4PILPS setup thus allows to focus in particular on the numerical accuracy of shortwave radiative transfer formulations and to pinpoint to areas where future model improvements should concentrate. The impact of increasing degrees of structural and spectral subgrid variability on the simulated fluxes is documented and the relevance of any thus emerging biases with respect to gross primary production estimates and shortwave radiative forcings due to snow and fire events are investigated.