Sampling density and proportion for the characterization of the variability of Oxisol attributes on different materials

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estimating the spatial scales of regionalized variables by nested sampling, hierarchical analysis of variance and residual maximum likelihood

The variogram is essential for local estimation and mapping of any variable by kriging. The variogram itself must usually be estimated from sample data. The sampling density is a compromise between precision and cost, but it must be sufficiently dense to encompass the principal spatial sources of variance. A nested, multi-stage, sampling with separating distances increasing in geometric progression from stage to stage will do that. The data may then be analyzed by a hierarchical analysis of variance to estimate the components of variance for every stage, and hence lag. By accumulating the components starting from the shortest lag one obtains a rough variogram for modest effort. For balanced designs the analysis of variance is optimal; for unbalanced ones, however, these estimators are not necessarily the best, and the analysis by residual maximum likelihood (REML) will usually be preferable. The paper summarizes the underlying theory and illustrates its application with data from three surveys, one in which the design had four stages and was balanced and two implemented with unbalanced designs to economize when there were more stages. A Fortran program is available for the analysis of variance, and code for the REML analysis is listed in the paper. (c) 2005 Elsevier Ltd. All rights reserved.

Estimating the spatial scales of regionalized variables by nested sampling, hierarchical analysis of variance and residual maximum likelihood

The variogram is essential for local estimation and mapping of any variable by kriging. The variogram itself must usually be estimated from sample data. The sampling density is a compromise between precision and cost, but it must be sufficiently dense to encompass the principal spatial sources of variance. A nested, multi-stage, sampling with separating distances increasing in geometric progression from stage to stage will do that. The data may then be analyzed by a hierarchical analysis of variance to estimate the components of variance for every stage, and hence lag. By accumulating the components starting from the shortest lag one obtains a rough variogram for modest effort. For balanced designs the analysis of variance is optimal; for unbalanced ones, however, these estimators are not necessarily the best, and the analysis by residual maximum likelihood (REML) will usually be preferable. The paper summarizes the underlying theory and illustrates its application with data from three surveys, one in which the design had four stages and was balanced and two implemented with unbalanced designs to economize when there were more stages. A Fortran program is available for the analysis of variance, and code for the REML analysis is listed in the paper. (c) 2005 Elsevier Ltd. All rights reserved.

The use of scaled semivariograms to plan soil sampling in sugarcane fields

Spatial sampling designs used to characterize the spatial variability of soil attributes are crucial for science studies. Sample planning for the interpolation of a regionalized variable may use several criteria, which could be best selected from an estimated semivariogram from a previously established grid. The objective of this study was to optimize the procedure for scaled semivariogram use to plan soil sampling in sugarcane fields in the Alfisol and Oxisol regions of Jaboticabal Town in So Paulo State, Brazil. A scaled semivariogram for several soil chemical attributes was estimated from the data obtained from two grids positioned on a sugarcane field area, sampled at a depth of 0.0-0.5 m. The research showed that regular grids with uniform intervals did not express the real spatial variability of the soil attributes of Oxisols and Alfisols in the study area. The calculated final sampling density based on the scaled parameters of the semivariogram was one sample for each 2 ha in Area 1 (convex landscape) and one sample for each 1 ha in Area 2 (linear landscape), as indicated by SANOS 0.1 software. The combined use of the simulation programs and scaled semivariograms can be used to define sampling points. These results may help in soil fertility mapping and thereby improve nutrient management in sugarcane crops.

Enhancing Sensitivity In Natural Product Nmr: Exploring Optimal Non-Uniform Sampling

Recent optimizations of NMR spectroscopy have focused their attention on innovations in new hardware, such as novel probes and higher field strengths. Only recently has the potential to enhance the sensitivity of NMR through data acquisition strategies been investigated. This thesis has focused on the practice of enhancing the signal-to-noise ratio (SNR) of NMR using non-uniform sampling (NUS). After first establishing the concept and exact theory of compounding sensitivity enhancements in multiple non-uniformly sampled indirect dimensions, a new result was derived that NUS enhances both SNR and resolution at any given signal evolution time. In contrast, uniform sampling alternately optimizes SNR (t < 1.26T2) or resolution (t~3T2), each at the expense of the other. Experiments were designed and conducted on a plant natural product to explore this behavior of NUS in which the SNR and resolution continue to improve as acquisition time increases. Possible absolute sensitivity improvements of 1.5 and 1.9 are possible in each indirect dimension for matched and 2x biased exponentially decaying sampling densities, respectively, at an acquisition time of ¿T2. Recommendations for breaking into the linear regime of maximum entropy (MaxEnt) are proposed. Furthermore, examination into a novel sinusoidal sampling density resulted in improved line shapes in MaxEnt reconstructions of NUS data and comparable enhancement to a matched exponential sampling density. The Absolute Sample Sensitivity derived and demonstrated here for NUS holds great promise in expanding the adoption of non-uniform sampling.

Performance Tuning Non-Uniform Sampling for Sensitivity Enhancement of Signal-Limited Biological NMR

Non-uniform sampling (NUS) has been established as a route to obtaining true sensitivity enhancements when recording indirect dimensions of decaying signals in the same total experimental time as traditional uniform incrementation of the indirect evolution period. Theory and experiments have shown that NUS can yield up to two-fold improvements in the intrinsic signal-to-noise ratio (SNR) of each dimension, while even conservative protocols can yield 20-40 % improvements in the intrinsic SNR of NMR data. Applications of biological NMR that can benefit from these improvements are emerging, and in this work we develop some practical aspects of applying NUS nD-NMR to studies that approach the traditional detection limit of nD-NMR spectroscopy. Conditions for obtaining high NUS sensitivity enhancements are considered here in the context of enabling H-1,N-15-HSQC experiments on natural abundance protein samples and H-1,C-13-HMBC experiments on a challenging natural product. Through systematic studies we arrive at more precise guidelines to contrast sensitivity enhancements with reduced line shape constraints, and report an alternative sampling density based on a quarter-wave sinusoidal distribution that returns the highest fidelity we have seen to date in line shapes obtained by maximum entropy processing of non-uniformly sampled data.

Recent Advances in Adaptive Sampling and Reconstruction for Monte Carlo Rendering

Monte Carlo integration is firmly established as the basis for most practical realistic image synthesis algorithms because of its flexibility and generality. However, the visual quality of rendered images often suffers from estimator variance, which appears as visually distracting noise. Adaptive sampling and reconstruction algorithms reduce variance by controlling the sampling density and aggregating samples in a reconstruction step, possibly over large image regions. In this paper we survey recent advances in this area. We distinguish between “a priori” methods that analyze the light transport equations and derive sampling rates and reconstruction filters from this analysis, and “a posteriori” methods that apply statistical techniques to sets of samples to drive the adaptive sampling and reconstruction process. They typically estimate the errors of several reconstruction filters, and select the best filter locally to minimize error. We discuss advantages and disadvantages of recent state-of-the-art techniques, and provide visual and quantitative comparisons. Some of these techniques are proving useful in real-world applications, and we aim to provide an overview for practitioners and researchers to assess these approaches. In addition, we discuss directions for potential further improvements.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

地球化学新方法在隐状矿床找矿预测中的应用研究

With the progress of prospecting, the need for the discovery of blind ore deposits become more and more urgent. To study and find out the method and technology for the discovery of blind and buried ores is now a priority task. New geochemical methods are key technology to discover blind ores. Information of mobile components related to blind ores were extracted using this new methods. These methods were tested and applied based on element' s mobile components migrating and enriched in geophysical-geochemical process. Several kinds of partial extraction techniques have tested based on element' s occurrence in hypergenic zone. Middle-large scale geochemical methods for exploration in forest and swamp have been tested. A serious of methods were tested and applied effetely about evaluation of regional geochemical anomaly, 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system instead of the normal net. 1. Element related with ores can be mobiled to migrate upwards and be absorpted by surface soil. These abnomal components can be concentrated by natural or artificial methods. These trace metalic ions partially exist in dissovlvable ion forms of active state, and partially have been absorbed by Fe-Mn oxide, soil and organic matter in the soil so that a series of reaction such as complex reaction have take place. Employing various partial extraction techniques, metallic ions related with the phase of the blind ores can be extracted, such as the technique of organic complex extraction, Fe-Mn oxide extraction and the extraction technique of metallic ions of various absorption phases. 2.1:200000 regional geochemical evaluation anomaly methods: Advantageous ore-forming areas were selected firstly. Center, concentration, morphological feature, belt of anomaly were choosed then. Geological and geochemical anomalies were combined. And geological and geochemical background information were restrained. Xilekuduke area in Fuyun sheet , Zhaheba area in Qiakuerte sheet, the west-north part in Ertai sheet and Hongshanzui anomaly in Daqiao sheet were selected as target areas, in Alertai, in the north of Xinjiang. in Xilekuduke area, 1:25000 soil geochemical methods sampling based on the net in dendritic water system was carried out. Cu anomaly and copper mineralization were determined in the center area. Au , Cu anomalies and high polarization anomaly were determined in the south part. Prospecting by primary halo and organic complex extraction were used to prognosis blind ore in widely rang outcrop of bedrock. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system were used in transported overburden outside of mining area. Shallow seismic method and primary halo found a new blind orebody in mining area. A mineralization site was fou and outside of Puziwan gold mine, in the north of Shanxi province. Developing middle-large scale geochemical exploration method is a key technique based 1:200000 regional geochemical exploration. Some conditions were tested as Sampling density , distribution sites of sample, grain size of sample and occurrence of element for exploration. 1:50000 exploration method was advanced to sample clast sediment supplement clast sediment in valley. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system was applied to sample residual material in A or C horizon. 1:2000 primary or soil halo methods used to check anomalies and determine mineralization. Daliang gold mineralization in the northern Moerdaoga was found appling these methods. Thermomagnetic method was tested in miniqi copper-polymetallic ore. Process methods such as grain size of sample, heated temperature, magnetic separating technique were tested. A suite of Thermomagnetic geochemical method was formed. This method was applied in Xiangshan Cu~Ni deposit which is cover by clast or Gobi in the eastern Xinjiang. Element's content and contrast of anomaly with Thermomagnetic geochemical method were higher than soil anomaly. Susceptibility after samples were heated could be as a assessment conference for anomaly. In some sectors thermo-magnetic Cu, Ni, Ti anomalious were found outside deposits area. There were strong anomal ies response up ore tested by several kind of partial extraction methods include Thermomagnetic, enzyme leach and other partial extractions in Kalatongke Cu-Ni deposit in hungriness area in the northern of Xinjiang. Element's anomalies of meobile were mainly in Fe-Mn oxide and salt. A Copper mineralization site in Xilekuduke anomaly area had been determined. A blind ore was foung by shallow seismic and geochemical method and a mineralization site was found outside this mining area in Puziwan gold deposit in shanxi province. A Gold mineralization site was found by 1:50000 geochemical exploration in Daliang, Inner Mongolia.

Using geostatistical Bayesian Updating to integrate airborne radiometrics and soil geochemistry to improve mapping for mineral exploration

Mineral exploration programmes around the world use data from remote sensing, geophysics and direct sampling. On a regional scale, the combination of airborne geophysics and ground-based geochemical sampling can aid geological mapping and economic minerals exploration. The fact that airborne geophysical and traditional soil-sampling data are generated at different spatial resolutions means that they are not immediately comparable due to their different sampling density. Several geostatistical techniques, including indicator cokriging and collocated cokriging, can be used to integrate different types of data into a geostatistical model. With increasing numbers of variables the inference of the cross-covariance model required for cokriging can be demanding in terms of effort and computational time. In this paper a Gaussian-based Bayesian updating approach is applied to integrate airborne radiometric data and ground-sampled geochemical soil data to maximise information generated from the soil survey, to enable more accurate geological interpretation for the exploration and development of natural resources. The Bayesian updating technique decomposes the collocated estimate into a production of two models: prior and likelihood models. The prior model is built from primary information and the likelihood model is built from secondary information. The prior model is then updated with the likelihood model to build the final model. The approach allows multiple secondary variables to be simultaneously integrated into the mapping of the primary variable. The Bayesian updating approach is demonstrated using a case study from Northern Ireland where the history of mineral prospecting for precious and base metals dates from the 18th century. Vein-hosted, strata-bound and volcanogenic occurrences of mineralisation are found. The geostatistical technique was used to improve the resolution of soil geochemistry, collected one sample per 2 km2, by integrating more closely measured airborne geophysical data from the GSNI Tellus Survey, measured over a footprint of 65 x 200 m. The directly measured geochemistry data were considered as primary data in the Bayesian approach and the airborne radiometric data were used as secondary data. The approach produced more detailed updated maps and in particular maximized information on mapped estimates of zinc, copper and lead. Greater delineation of an elongated northwest/southeast trending zone in the updated maps strengthened the potential to investigate stratabound base metal deposits.

Cartografia geoquímica da ilha de Santiago com uma densidade de amostragem média/baixa

A ilha de Santiago é a maior do arquipélago de Cabo Verde, apresenta uma área de 991 km2, com um comprimento e largura máximos de 54,9 km e 29 km, respectivamente, e 1392 m de maior altitude. As condições climáticas e a erosão são alguns dos problemas naturais do arquipélago de Cabo Verde. Além disso, a intervenção humana no ambiente superficial revela-se muitas vezes inadequada e com forte impacto. A construção em solos com aptidão agrícola ou florestal, as actividades industriais, a deposição de materiais sólidos ou líquidos de forma inadequada, as práticas agrícolas incorrectas e intensivas, o uso abusivo de pesticidas e fertilizantes, a rega com água contaminada, a sobreexploração de aquíferos que levam muitas vezes à salinização, etc. têm enorme impacto em termos de contaminação de solos, águas superficiais e subterrâneas nas áreas envolventes. Por conseguinte, o conhecimento da variabilidade geoquímica natural é fundamental para a resolução de questões de índole económica, ambiental e de ordenamento do território, médica, e jurídica. A necessidade de construir uma base de dados de geoquímica georeferenciada que caracterize o ambiente superficial da ilha de Santiago foi a principal motivação para a realização deste estudo. Realizou-se um levantamento geoquímico de 337 amostras de sedimentos de corrente e 249 amostras de solos na ilha de Santiago, tendo sido seguidas as recomendações do Projecto Internacional IGCP 259 não só na fase de amostragem, mas também nas fases seguintes de preparação, análise, tratamento dos dados e elaboração de mapas. Determinaram-se os teores, na fracção < 2 mm, para 36 elementos – 9 elementos maiores (Al, Ca, Fe, K, Mg, Mn, Na, P, Ti) e 27 elementos vestigiais (Ag, As, Au, B, Ba, Bi, Cd, Co, Cr, Cu, Ga, Hg, La, Mo, Ni, Pb, S, Sb, Sc, Se, Sr, Th, Tl, U, V, W, Zn). Efectuou-se ainda a análise textural e estudou-se a composição mineralógica de cerca de 25% das amostras. Analisaram-se, também, 103 amostras de rochas, colhidas nas várias formações da ilha de Santiago, tendo sido determinados os teores de K2O, Na2O, Fe2O3(T), MnO, Sc, Cr, Co, Zn, Ga, As, Br, Rb, Zr, Sb, Cs, Ba, Hf, Ta, W, Th e U e REE, a fim de se fazer uma comparação com os teores destes elementos encontrados nos solos e sedimentos de corrente, averiguando se a sua variação é ou não essencialmente condicionada pela geoquímica da rocha-mãe. Os padrões geoquímicos obtidos através dos mapas de distribuição espacial foram correlacionados com a natureza da rocha mãe, o tipo de solo, e ainda com algumas fontes de contaminação. A interpretação dos resultados foi realizada não só pela observação dos mapas geoquímicos, mas também após análise estatística dos conjuntos de dados obtidos, e apoiada em informação diversa disponível. A utilização da Análise de Componentes Principais permitiu distinguir associações entre elementos químicos, quer de origem geogénica quer antropogénica. Foram ainda elaborados mapas de distribuição espacial de vários índices multielementares de importância ambiental, como o Índice de Acidificação Al/(Ca+Mg+K), Índice de Combi, Índice de Avaliação de Risco Ambiental e Índices de Enriquecimento/Contaminação para vários grupos de elementos considerados como “primary pollutant metals”.

Particulate matter flux deposition in the vicinity of a coal-fired power plant

This paper describes the methodology adopted to assess local air quality impact in the vicinity of a coal power plant located in the south of Portugal. Two sampling areas were selected to assess the deposition flux of dust fallout and its potential spatial heterogeneity. The sampling area was divided into two subareas: the inner, with higher sampling density and urban and suburban characteristics, inside a 6-km circle centered on the stacks, and an outer subarea, mainly rural, with lower sampling density within a radius of 20 km. Particulate matter deposition was studied in the vicinity of the coal fired power plant during three seasonal sampling campaigns. For the first one, the average annual flux of dust fallout was 22.51 g/(m2 yr), ranging from 4.20 to 65.94 g/(m2 yr); for the second one was 9.47 g/(m2 yr), ranging from 0.78 to 32.72 g/(m2 yr) and for the last one was 38.42 g/(m2 yr), ranging from 1.41 to 117.48 g/(m2 yr). The fallout during the second campaign turned out to be much lower than for others. This was in part due to meteorological local patterns but mostly due to the fact that the power plant was not working at full power during the second sampling campaign.155

Systematics and evolution of the genus Deuterocohnia Mez (Bromeliaceae)

The present study investigates the systematics and evolution of the Neotropical genus Deuterocohnia Mez (Bromeliaceae). It provides a comprehensive taxonomic revision as well as phylogenetic analyses based on chloroplast and nuclear DNA sequences and presents a hypothesis on the evolution of the genus. A broad morphological, anatomical, biogeographical and ecological overview of the genus is given in the first part of the study. For morphological character assessment more than 700 herbarium specimens from 39 herbaria as well as living plant material in the field and in the living collections of botanical gardens were carefully examined. The arid habitats, in which the species of Deuterocohnia grow, are reflected by the morphological and anatomical characters of the species. Important characters for species delimitation were identified, like the length of the inflorescence, the branching order, the density of flowers on partial inflorescences, the relation of the length of the primary bracts to that of the partial inflorescence, the sizes of floral bracts, sepals and petals, flower colour, the presence or absence of a pedicel, the curvature of the stamina and the petals during anthesis. After scrutinizing the nomenclatural history of the taxa belonging to Deuterocohnia – including the 1992 syonymized genus Abromeitiella – 17 species, 4 subspecies and 4 varieties are accepted in the present revision. Taxonomic changes were made in the following cases: (I) New combinations: A. abstrusa (A. Cast.) N. Schütz is re-established – as defined by Castellanos (1931) – and transfered to D. abstrusa; D. brevifolia (Griseb.) M.A. Spencer & L.B. Sm. includes accessions of the former D. lorentziana (Mez) M.A. Spencer & L.B. Sm., which are not assigned to D. abstrusa; D. bracteosa W. Till is synonymized to D. strobilifera Mez; D. meziana Kuntze ex Mez var. carmineo-viridiflora Rauh is classified as a subspecies of D. meziana (ssp. carmineo-viridiflora (Rauh) N. Schütz); D. pedicellata W. Till is classified as a subspecies of D. meziana (ssp. pedicellata (W. Till) N. Schütz); D. scapigera (Rauh & L. Hrom.) M.A. Spencer & L.B. Sm ssp. sanctae-crucis R. Vásquez & Ibisch is classified as a species (D. sanctae-crucis (R. Vásquez & Ibisch) N. Schütz); (II) New taxa: a new subspecies of D. meziana Kuntze ex Mez is established; a new variety of D. scapigera is established; (the new taxa will be validly published elsewhere); (III) New type: an epitype for D. longipetala was chosen. All other species were kept according to Spencer and Smith (1992) or – in the case of more recently described species – according to the protologue. Beside the nomenclatural notes and the detailed descriptions, information on distribution, habitat and ecology, etymology and taxonomic delimitation is provided for the genus and for each of its species. An key was constructed for the identification of currently accepted species, subspecies and varieties. The key is based on easily detectable morphological characters. The former synonymization of the genus Abromeitiella into Deuterocohnia (Spencer and Smith 1992) is re-evalutated in the present study. Morphological as well as molecular investigations revealed Deuterocohnia incl. Abromeitiella as being monophyletic, with some indications that a monophyletic Abromeitiella lineage arose from within Deuterocohnia. Thus the union of both genera is confirmed. The second part of the present thesis describes and discusses the molecular phylogenies and networks. Molecular analyses of three chloroplast intergenic spacers (rpl32-trnL, rps16-trnK, trnS-ycf3) were conducted with a sample set of 119 taxa. This set included 103 Deuterocohnia accessions from all 17 described species of the genus and 16 outgroup taxa from the remainder of Pitcairnioideae s.str. (Dyckia (8 sp.), Encholirium (2 sp.), Fosterella (4 sp.) and Pitcairnia (2 sp.)). With its high sampling density, the present investigation by far represents the most comprehensive molecular study of Deuterocohnia up till now. All data sets were analyzed separately as well as in combination, and various optimality criteria for phylogenetic tree construction were applied (Maximum Parsimony, Maximum Likelihood, Bayesian inferences and the distance method Neighbour Joining). Congruent topologies were generally obtained with different algorithms and optimality criteria, but individual clades received different degrees of statistical support in some analyses. The rps16-trnK locus was the most informative among the three spacer regions examined. The results of the chloroplast DNA analyses revealed a highly supported paraphyly of Deuterocohnia. Thus, the cpDNA trees divide the genus into two subclades (A and B), of which Deuterocohnia subclade B is sister to the included Dyckia and Encholirium accessions, and both together are sister to Deuterocohnia subclade A. To further examine the relationship between Deuterocohnia and Dyckia/Encholirium at the generic level, two nuclear low copy markers (PRK exon2-5 and PHYC exon1) were analysed with a reduced taxon set. This set included 22 Deuterocohnia accessions (including members of both cpDNA subclades), 2 Dyckia, 2 Encholirium and 2 Fosterella species. Phylogenetic trees were constructed as described above, and for comparison the same reduced taxon set was also analysed at the three cpDNA data loci. In contrast to the cpDNA results, the nuclear DNA data strongly supported the monophyly of Deuterocohnia, which takes a sister position to a clade of Dyckia and Encholirium samples. As morphology as well as nuclear DNA data generated in the present study and in a former AFLP analysis (Horres 2003) all corroborate the monophyly of Deuterocohnia, the apparent paraphyly displayed in cpDNA analyses is interpreted to be the consequence of a chloroplast capture event. This involves the introgression of the chloroplast genome from the common ancestor of the Dyckia/ Encholirium lineage into the ancestor of Deuterocohnia subclade B species. The chloroplast haplotypes are not species-specific in Deuterocohnia. Thus, one haplotype was sometimes shared by several species, where the same species may harbour different haplotypes. The arrangement of haplotypes followed geographical patterns rather than taxonomic boundaries, which may indicate some residual gene flow among populations from different Deuteroccohnia species. Phenotypic species coherence on the background of ongoing gene flow may then be maintained by sets of co-adapted alleles, as was suggested by the porous genome concept (Wu 2001, Palma-Silva et al. 2011). The results of the present study suggest the following scenario for the evolution of Deuterocohnia and its species. Deuterocohnia longipetala may be envisaged as a representative of the ancestral state within the genus. This is supported by (1) the wide distribution of this species; (2) the overlap in distribution area with species of Dyckia; (3) the laxly flowered inflorescences, which are also typical for Dyckia; (4) the yellow petals with a greenish tip, present in most other Deuterocohnia species. The following six extant lineages within Deuterocohnia might have independently been derived from this ancestral state with a few changes each: (I) D. meziana, D. brevispicata and D. seramisiana (Bolivia, lowland to montane areas, mostly reddish-greenish coloured, very laxly to very densely flowered); (II) D. strobilifera (Bolivia, high Andean mountains, yellow flowers, densely flowered); (III) D. glandulosa (Bolivia, montane areas, yellow-greenish flowers, densely flowered); (IV) D. haumanii, D. schreiteri, D. digitata, and D. chrysantha (Argentina, Chile, E Andean mountains and Atacama desert, yellow-greenish flowers, densely flowered); (V) D. recurvipetala (Argentina, foothills of the Andes, recurved yellow flowers, laxly flowered); (VI) D. gableana, D. scapigera, D. sanctae-crucis, D. abstrusa, D. brevifolia, D. lotteae (former Abromeitiella species, Bolivia, Argentina, higher Andean mountains, greenish-yellow flowers, inflorescence usually simple). Originating from the lower montane Andean regions, at least four lineages of the genus (I, II, IV, VI) adapted in part to higher altitudes by developing densely flowered partial inflorescences, shorter flowers and – in at least three lineages (II, IV, VI) – smaller rosettes, whereas species spreading into the lowlands (I, V) developed larger plants, laxly flowered, amply branched inflorescences and in part larger flowers (I).

Are Cape floral clades the same age? Contemporaneous origins of two lineages in the genistoids s.l. (Fabaceae)

The hypothesis that the elements of the modern species-rich flora of the Cape Floristic Region (CFR), South Africa, originated more or less simultaneously at the Miocene/Pliocene boundary, in response to the development of a mediterranean climate, has been challenged by numerous molecular dating estimates of Cape floral clades. These studies reveal a more gradual emergence, with the oldest clades originating in the Eocene, but others appearing later, some as recently as the Pliocene. That there are factors which might affect the dates recovered, such as choice of calibration point, analysis method, sampling density and the delimitation of Cape floral clades, suggests a need for further critical evaluation of the age estimates presented to date. In this study, the dates of origin of two Cape floral clades (the legume Crotalarieae p.p. and Podalyrieae) are estimated, constrained by a shared calibration point in a single analysis using an rDNA ITS phylogeny in which 633 taxa are sampled. The results indicate that these two clades arose contemporaneously 44-46 mya, not at the Miocene/Pliocene boundary as had been previously supposed. The contemporaneous origin of these Cape floral clades suggests that additional more inclusive analyses are needed before rejecting the hypothesis that a. single environmental trigger explains the establishment of Cape floral clades. (c) 2007 Elsevier Inc. All rights reserved.

Undersampled Critical Branching Processes on Small-World and Random Networks Fail to Reproduce the Statistics of Spike Avalanches

The power-law size distributions obtained experimentally for neuronal avalanches are an important evidence of criticality in the brain. This evidence is supported by the fact that a critical branching process exhibits the same exponent t~3=2. Models at criticality have been employed to mimic avalanche propagation and explain the statistics observed experimentally. However, a crucial aspect of neuronal recordings has been almost completely neglected in the models: undersampling. While in a typical multielectrode array hundreds of neurons are recorded, in the same area of neuronal tissue tens of thousands of neurons can be found. Here we investigate the consequences of undersampling in models with three different topologies (two-dimensional, small-world and random network) and three different dynamical regimes (subcritical, critical and supercritical). We found that undersampling modifies avalanche size distributions, extinguishing the power laws observed in critical systems. Distributions from subcritical systems are also modified, but the shape of the undersampled distributions is more similar to that of a fully sampled system. Undersampled supercritical systems can recover the general characteristics of the fully sampled version, provided that enough neurons are measured. Undersampling in two-dimensional and small-world networks leads to similar effects, while the random network is insensitive to sampling density due to the lack of a well-defined neighborhood. We conjecture that neuronal avalanches recorded from local field potentials avoid undersampling effects due to the nature of this signal, but the same does not hold for spike avalanches. We conclude that undersampled branching-process-like models in these topologies fail to reproduce the statistics of spike avalanches.