The importance of mesofauna and decomposition environment on leaf decomposition in three forests in southeastern Brazil

We examined the effects of soil mesofauna and the litter decomposition environment (above and belowground) on leaf decomposition rates in three forest types in southeastern Brazil. To estimate decomposition experimentally, we used litterbags with a standard substrate in a full-factorial experimental design. We used model selection to compare three decomposition models and also to infer the importance of forest type, decomposition environment, mesofauna, and their interactions on the decomposition process. Rather than the frequently used simple and double-exponential models, the best model to describe our dataset was the exponential deceleration model, which assumed a single organic compartment with an exponential decrease of the decomposition rate. Decomposition was higher in the wet than in the seasonal forest, and the differences between forest types were stronger aboveground. Regarding litter decomposition environment, decomposition was predominantly higher below than aboveground, but the magnitude of this effect was higher in the seasonal than in wet forests. Mesofauna exclusion treatments had slower decomposition, except aboveground into the Semi-deciduous Forest, where the mesofauna presence did not affect decomposition. Furthermore, the effect of mesofauna was stronger in the wet forests and belowground. Overall, our results suggest that, in a regional scale, both decomposers activity and the positive effect of soil mesofauna in decomposition are constrained by abiotic factors, such as moisture conditions.

Investigation of the pore blockage of a Brazilian dolomite during the sulfation reaction

The influence of the temperature and reaction time on the sulfation process of a dolomite is investigated in this paper. The sulfation effectiveness was evaluated and correlated with changes in the physical characteristics of a Brazilian dolomite during the reactive process. Calcination and sulfation experiments were performed under isothermal conditions for dolomite samples with average particle sizes of 545 mu m at temperatures of 750 degrees C, 850 degrees C and 950 degrees C at different times of sulfation. Thermogravimetric tests were applied to establish the reactivity variation of the dolomite in function of the time in the sulfation reaction and evaluate the methodology of the samples preparation. Porosimetry tests were performed to study the pore blockage of dolomite during the sulfation reaction. The highest values of BET surface area were 25.55 m(2)/g, 29.55 m(2)/g and 12.62 m(2)/g for calcined samples and after their sulfation processes, conversions of 51.5%, 61.9% and 42.8% were obtained at 750 degrees C, 850 degrees C and 950 degrees C, respectively. Considering the process as a whole, the best fit was provided by a first-order exponential decay equation. Moreover, the results have shown that it is possible to quantify the decreasing in the dolomite reactivity for sulfur dioxide sorption and understand the changes in the behavior of the sulfation process of limestones when applied to technologies, as fluidized bed combustor, in which sulfur dioxide is present. (C) 2011 Elsevier B. V. All rights reserved.

An ancient founder mutation in PROKR2 impairs human reproduction

Congenital gonadotropin-releasing hormone (GnRH) deficiency manifests as absent or incomplete sexual maturation and infertility. Although the disease exhibits marked locus and allelic heterogeneity, with the causal mutations being both rare and private, one causal mutation in the prokineticin receptor, PROKR2 L173R, appears unusually prevalent among GnRH-deficient patients of diverse geographic and ethnic origins. To track the genetic ancestry of PROKR2 L173R, haplotype mapping was performed in 22 unrelated patients with GnRH deficiency carrying L173R and their 30 first-degree relatives. The mutations age was estimated using a haplotype-decay model. Thirteen subjects were informative and in all of them the mutation was present on the same approximate to 123 kb haplotype whose population frequency is 10. Thus, PROKR2 L173R represents a founder mutation whose age is estimated at approximately 9000 years. Inheritance of PROKR2 L173R-associated GnRH deficiency was complex with highly variable penetrance among carriers, influenced by additional mutations in the other PROKR2 allele (recessive inheritance) or another gene (digenicity). The paradoxical identification of an ancient founder mutation that impairs reproduction has intriguing implications for the inheritance mechanisms of PROKR2 L173R-associated GnRH deficiency and for the relevant processes of evolutionary selection, including potential selective advantages of mutation carriers in genes affecting reproduction.

Quantification of cellular action fibers alignment from confocal microscopic images.

The cellular rheology has recently undergone a rapid development with particular attention to the cytoskeleton mechanical properties and its main components - actin filaments, intermediate filaments, microtubules and crosslinked proteins. However it is not clear what are the cellular structural changes that directly affect the cell mechanical properties. Thus, in this work, we aimed to quantify the structural rearrangement of these fibers that may emerge in changes in the cell mechanics. We created an image analysis platform to study smooth muscle cells from different arteries: aorta, mammary, renal, carotid and coronary and processed respectively 31, 29, 31, 30 and 35 cell image obtained by confocal microscopy. The platform was developed in Matlab (MathWorks) and it uses the Sobel operator to determine the actin fiber image orientation of the cell, labeled with phalloidin. The Sobel operator is used as a filter capable of calculating the pixel brightness gradient, point to point, in the image. The operator uses vertical and horizontal convolution kernels to calculate the magnitude and the angle of the pixel intensity gradient. The image analysis followed the sequence: (1) opens a given cells image set to be processed; (2) sets a fix threshold to eliminate noise, based on Otsu's method; (3) detect the fiber edges in the image using the Sobel operator; and (4) quantify the actin fiber orientation. Our first result is the probability distribution II(Δθ) to find a given fiber angle deviation (Δθ) from the main cell fiber orientation θ0. The II(Δθ) follows an exponential decay II(Δθ) = Aexp(-αΔθ) regarding to its θ0. We defined and determined a misalignment index α of the fibers of each artery kind: coronary αCo = (1.72 ‘+ or =’ 0.36)rad POT -1; renal αRe = (1.43 + or - 0.64)rad POT -1; aorta αAo = (1.42 + or - 0.43)rad POT -1; mammary αMa = (1.12 + or - 0.50)rad POT -1; and carotid αCa = (1.01 + or - 0.39)rad POT -1. The α of coronary and carotid are statistically different (p < 0.05) among all analyzed cells. We discussed our results correlating the misalignment index data with the experimental cell mechanical properties obtained by using Optical Magnetic Twisting Cytometry with the same group of cells.

Concentration quenching and thermal activation of the luminescence from terbium-doped a-SiC:H and c-AlN thin films

The effect of terbium (Tb) doping on the photoluminescence (PL) of crystalline aluminum nitride (c-AlN) and amorphous hydrogenated silicon carbide (a-SiC:H) thin films has been investigated for different Tb atomic concentrations. The samples were prepared by DC and RF magnetron reactive sputtering techniques covering the concentration range of Tb from 0.5 to 11 at.%. The Tb-related light emission versus the Tb concentration is reported for annealing temperatures of 450 °C, 750 °C and 1000 °C. In the low concentration region the intensity exhibits a linear increase and its slope is enhanced with the annealing temperature giving an activation energy of 0.106 eV in an Arrhenius plot. In the high concentration region an exponential decay is recorded which is almost independent on the host material, its structure and the annealing process.

Echtzeitsuche nach Neutrinoausbrüchen von Supernovae mit dem AMANDA-II-Detektor

Der AMANDA-II Detektor ist primär für den richtungsaufgelösten Nachweis hochenergetischer Neutrinos konzipiert. Trotzdem können auch niederenergetische Neutrinoausbrüche, wie sie von Supernovae erwartet werden, mit hoher Signifikanz nachgewiesen werden, sofern sie innerhalb der Milchstraße stattfinden. Die experimentelle Signatur im Detektor ist ein kollektiver Anstieg der Rauschraten aller optischen Module. Zur Abschätzung der Stärke des erwarteten Signals wurden theoretische Modelle und Simulationen zu Supernovae und experimentelle Daten der Supernova SN1987A studiert. Außerdem wurden die Sensitivitäten der optischen Module neu bestimmt. Dazu mussten für den Fall des südpolaren Eises die Energieverluste geladener Teilchen untersucht und eine Simulation der Propagation von Photonen entwickelt werden. Schließlich konnte das im Kamiokande-II Detektor gemessene Signal auf die Verhältnisse des AMANDA-II Detektors skaliert werden. Im Rahmen dieser Arbeit wurde ein Algorithmus zur Echtzeit-Suche nach Signalen von Supernovae als Teilmodul der Datennahme implementiert. Dieser beinhaltet diverse Verbesserungen gegenüber der zuvor von der AMANDA-Kollaboration verwendeten Version. Aufgrund einer Optimierung auf Rechengeschwindigkeit können nun mehrere Echtzeit-Suchen mit verschiedenen Analyse-Zeitbasen im Rahmen der Datennahme simultan laufen. Die Disqualifikation optischer Module mit ungeeignetem Verhalten geschieht in Echtzeit. Allerdings muss das Verhalten der Module zu diesem Zweck anhand von gepufferten Daten beurteilt werden. Dadurch kann die Analyse der Daten der qualifizierten Module nicht ohne eine Verzögerung von etwa 5 Minuten geschehen. Im Falle einer erkannten Supernova werden die Daten für die Zeitdauer mehrerer Minuten zur späteren Auswertung in 10 Millisekunden-Intervallen archiviert. Da die Daten des Rauschverhaltens der optischen Module ansonsten in Intervallen von 500 ms zur Verfgung stehen, ist die Zeitbasis der Analyse in Einheiten von 500 ms frei wählbar. Im Rahmen dieser Arbeit wurden drei Analysen dieser Art am Südpol aktiviert: Eine mit der Zeitbasis der Datennahme von 500 ms, eine mit der Zeitbasis 4 s und eine mit der Zeitbasis 10 s. Dadurch wird die Sensitivität für Signale maximiert, die eine charakteristische exponentielle Zerfallszeit von 3 s aufweisen und gleichzeitig eine gute Sensitivität über einen weiten Bereich exponentieller Zerfallszeiten gewahrt. Anhand von Daten der Jahre 2000 bis 2003 wurden diese Analysen ausführlich untersucht. Während die Ergebnisse der Analyse mit t = 500 ms nicht vollständig nachvollziehbare Ergebnisse produzierte, konnten die Resultate der beiden Analysen mit den längeren Zeitbasen durch Simulationen reproduziert und entsprechend gut verstanden werden. Auf der Grundlage der gemessenen Daten wurden die erwarteten Signale von Supernovae simuliert. Aus einem Vergleich zwischen dieser Simulation den gemessenen Daten der Jahre 2000 bis 2003 und der Simulation des erwarteten statistischen Untergrunds kann mit einem Konfidenz-Niveau von mindestens 90 % gefolgert werden, dass in der Milchstraße nicht mehr als 3.2 Supernovae pro Jahr stattfinden. Zur Identifikation einer Supernova wird ein Ratenanstieg mit einer Signifikanz von mindestens 7.4 Standardabweichungen verlangt. Die Anzahl erwarteter Ereignisse aus dem statistischen Untergrund beträgt auf diesem Niveau weniger als ein Millionstel. Dennoch wurde ein solches Ereignis gemessen. Mit der gewählten Signifikanzschwelle werden 74 % aller möglichen Vorläufer-Sterne von Supernovae in der Galaxis überwacht. In Kombination mit dem letzten von der AMANDA-Kollaboration veröffentlicheten Ergebnis ergibt sich sogar eine obere Grenze von nur 2.6 Supernovae pro Jahr. Im Rahmen der Echtzeit-Analyse wird für die kollektive Ratenüberhöhung eine Signifikanz von mindestens 5.5 Standardabweichungen verlangt, bevor eine Meldung über die Detektion eines Supernova-Kandidaten verschickt wird. Damit liegt der überwachte Anteil Sterne der Galaxis bei 81 %, aber auch die Frequenz falscher Alarme steigt auf bei etwa 2 Ereignissen pro Woche. Die Alarm-Meldungen werden über ein Iridium-Modem in die nördliche Hemisphäre übertragen, und sollen schon bald zu SNEWS beitragen, dem weltweiten Netzwerk zur Früherkennung von Supernovae.

Feasibility study of performing high precision gamma spectroscopy of lambda lambda hypernuclei in the PANDA experiment

Hypernuclear physics is currently attracting renewed interest, due tornthe important role of hypernuclei spectroscopy rn(hyperon-hyperon and hyperon-nucleon interactions) rnas a unique toolrnto describe the baryon-baryon interactions in a unified way and to rnunderstand the origin of their short-range.rnrnHypernuclear research will be one of the main topics addressed by the {sc PANDA} experimentrnat the planned Facility for Antiproton and Ion Research {sc FAIR}.rnThanks to the use of stored $overline{p}$ beams, copiousrnproduction of double $Lambda$ hypernuclei is expected at thern{sc PANDA} experiment, which will enable high precision $gamma$rnspectroscopy of such nuclei for the first time.rnAt {sc PANDA} excited states of $Xi^-$ hypernuclei will be usedrnas a basis for the formation of double $Lambda$ hypernuclei.rnFor their detection, a devoted hypernuclear detector setup is planned. This setup consists ofrna primary nuclear target for the production of $Xi^{-}+overline{Xi}$ pairs, a secondary active targetrnfor the hypernuclei formation and the identification of associated decay products and a germanium array detector to perform $gamma$ spectroscopy.rnrnIn the present work, the feasibility of performing high precision $gamma$rnspectroscopy of double $Lambda$ hypernuclei at the {sc PANDA} experiment has been studiedrnby means of a Monte Carlo simulation. For this issue, the designing and simulation of the devoted detector setup as well as of the mechanism to produce double $Lambda$ hypernuclei have been optimizedrntogether with the performance of the whole system. rnIn addition, the production yields of double hypernuclei in excitedrnparticle stable states have been evaluated within a statistical decay model.rnrnA strategy for the unique assignment of various newly observed $gamma$-transitions rnto specific double hypernuclei has been successfully implemented by combining the predicted energy spectra rnof each target with the measurement of two pion momenta from the subsequent weak decays of a double hypernucleus.rn% Indeed, based on these Monte Carlo simulation, the analysis of the statistical decay of $^{13}_{Lambda{}Lambda}$B has been performed. rn% As result, three $gamma$-transitions associated to the double hypernuclei $^{11}_{Lambda{}Lambda}$Bern% and to the single hyperfragments $^{4}_{Lambda}$H and $^{9}_{Lambda}$Be, have been well identified.rnrnFor the background handling a method based on time measurement has also been implemented.rnHowever, the percentage of tagged events related to the production of $Xi^{-}+overline{Xi}$ pairs, variesrnbetween 20% and 30% of the total number of produced events of this type. As a consequence, further considerations have to be made to increase the tagging efficiency by a factor of 2.rnrnThe contribution of the background reactions to the radiation damage on the germanium detectorsrnhas also been studied within the simulation. Additionally, a test to check the degradation of the energyrnresolution of the germanium detectors in the presence of a magnetic field has also been performed.rnNo significant degradation of the energy resolution or in the electronics was observed. A correlationrnbetween rise time and the pulse shape has been used to correct the measured energy. rnrnBased on the present results, one can say that the performance of $gamma$ spectroscopy of double $Lambda$ hypernuclei at the {sc PANDA} experiment seems feasible.rnA further improvement of the statistics is needed for the background rejection studies. Moreover, a more realistic layout of the hypernuclear detectors has been suggested using the results of these studies to accomplish a better balance between the physical and the technical requirements.rn

Prediction and measurement of forced response of a composite blade undergoing nonlinear vibration

The main objective of this project is to experimentally demonstrate geometrical nonlinear phenomena due to large displacements during resonant vibration of composite materials and to explain the problem associated with fatigue prediction at resonant conditions. Three different composite blades to be tested were designed and manufactured, being their difference in the composite layup (i.e. unidirectional, cross-ply, and angle-ply layups). Manual envelope bagging technique is explained as applied to the actual manufacturing of the components; problems encountered and their solutions are detailed. Forced response tests of the first flexural, first torsional, and second flexural modes were performed by means of a uniquely contactless excitation system which induced vibration by using a pulsed airflow. Vibration intensity was acquired by means of Polytec LDV system. The first flexural mode is found to be completely linear irrespective of the vibration amplitude. The first torsional mode exhibits a general nonlinear softening behaviour which is interestingly coupled with a hardening behaviour for the unidirectional layup. The second flexural mode has a hardening nonlinear behaviour for either the unidirectional and angle-ply blade, whereas it is slightly softening for the cross-ply layup. By using the same equipment as that used for forced response analyses, free decay tests were performed at different airflow intensities. Discrete Fourier Trasform over the entire decay and Sliding DFT were computed so as to visualise the presence of nonlinear superharmonics in the decay signal and when they were damped out from the vibration over the decay time. Linear modes exhibit an exponential decay, while nonlinearities are associated with a dry-friction damping phenomenon which tends to increase with increasing amplitude. Damping ratio is derived from logarithmic decrement for the exponential branch of the decay.

IN-VIVO CT DOSIMETRY DURING VIRTUAL COLONOSCOPY

Virtual colonoscopy (VC) is a minimally invasive means for identifying colorectal polyps and colorectal lesions by insufflating a patient’s bowel, applying contrast agent via rectal catheter, and performing multi-detector computed tomography (MDCT) scans. The technique is recommended for colonic health screening by the American Cancer Society but not funded by the Centers for Medicare and Medicaid Services (CMS) partially because of potential risks from radiation exposure. To date, no in‐vivo organ dose measurements have been performed for MDCT scans; thus, the accuracy of any current dose estimates is currently unknown. In this study, two TLDs were affixed to the inner lumen of standard rectal catheters used in VC, and in-vivo rectal dose measurements were obtained within 6 VC patients. In order to calculate rectal dose, TLD-100 powder response was characterized at diagnostic doses such that appropriate correction factors could be determined for VC. A third-order polynomial regression with a goodness of fit factor of R2=0.992 was constructed from this data. Rectal dose measurements were acquired with TLDs during simulated VC within a modified anthropomorphic phantom configured to represent three sizes of patients undergoing VC. The measured rectal doses decreased in an exponential manner with increasing phantom effective diameter, with R2=0.993 for the exponential regression model and a maximum percent coefficient of variation (%CoV) of 4.33%. In-vivo measurements yielded rectal doses ranged from that decreased exponentially with increasing patient effective diameter, in a manner that was also favorably predicted by the size specific dose estimate (SSDE) model for all VC patients that were of similar age, body composition, and TLD placement. The measured rectal dose within a younger patient was favorably predicted by the anthropomorphic phantom dose regression model due to similarities in the percentages of highly attenuating material at the respective measurement locations and in the placement of the TLDs. The in-vivo TLD response did not increase in %CoV with decreasing dose, and the largest %CoV was 10.0%.

Dynamics of the nasal microbiota in infancy: A prospective cohort study.

BACKGROUND Understanding the composition and dynamics of the upper respiratory tract microbiota in healthy infants is a prerequisite to investigate the role of the microbiota in patients with respiratory diseases. This is especially true in early life, when the immune system is in development. OBJECTIVE We sought to describe the dynamics of the upper respiratory tract microbiota in healthy infants within the first year of life. METHODS After exclusion of low-quality samples, microbiota characterization was performed by using 16S rDNA pyrosequencing of 872 nasal swabs collected biweekly from 47 unselected infants. RESULTS Bacterial density increased and diversity decreased within the first year of life (R(2) = 0.95 and 0.73, respectively). A distinct profile for the first 3 months of life was found with increased relative abundances of Staphlyococcaceae and Corynebacteriaceae (exponential decay: R(2) = 0.94 and 0.96, respectively). In addition, relative bacterial abundance and composition differed significantly from summer to winter months. The individual composition of the microbiota changed with increasing time intervals between samples and was best modeled by an exponential function (R(2) = 0.97). Within-subject dissimilarity in a 2-week time interval was consistently lower than that between subjects, indicating a personalized microbiota. CONCLUSION This study reveals age and seasonality as major factors driving the composition of the nasal microbiota within the first year of life. A subject's microbiota is personalized but dynamic throughout the first year. These data are indispensable to interpretation of cross-sectional studies and investigation of the role of the microbiota in both healthy subjects and patients with respiratory diseases. They might also serve as a baseline for future intervention studies.

Modeling the Volatility of Real GDP Growth: The Case of Japan Revisited

Previous studies (e.g., Hamori, 2000; Ho and Tsui, 2003; Fountas et al., 2004) find high volatility persistence of economic growth rates using generalized autoregressive conditional heteroskedasticity (GARCH) specifications. This paper reexamines the Japanese case, using the same approach and showing that this finding of high volatility persistence reflects the Great Moderation, which features a sharp decline in the variance as well as two falls in the mean of the growth rates identified by Bai and Perronâs (1998, 2003) multiple structural change test. Our empirical results provide new evidence. First, excess kurtosis drops substantially or disappears in the GARCH or exponential GARCH model that corrects for an additive outlier. Second, using the outlier-corrected data, the integrated GARCH effect or high volatility persistence remains in the specification once we introduce intercept-shift dummies into the mean equation. Third, the time-varying variance falls sharply, only when we incorporate the break in the variance equation. Fourth, the ARCH in mean model finds no effects of our more correct measure of output volatility on output growth or of output growth on its volatility.

Mean snow accumulation rate, and mean methane sulfonate, chloride and nitrate concentration in snow pits and firn cores of Dronning Maud Land

We quantified postdepositional losses of methane sulfonate (MSA-), nitrate, and chloride at the European Project for Ice Coring in Antarctica (EPICA) drilling site in Dronning Maud Land (DML) (75°S, 0°E). Analyses of four intermediate deep firn cores and 13 snow pits were considered. We found that about 26 ± 13% of the once deposited nitrate and typically 51 ± 20% of MSA- were lost, while for chloride, no significant depletion could be observed in firn older than one year. Assuming a first order exponential decay rate, the characteristic e-folding time for MSA- is 6.4 ± 3 years and 19 ± 6 years for nitrate. It turns out that for nitrate and MSA- the typical mean concentrations representative for the last 100 years were reached after 5.4 and 6.5 years, respectively, indicating that beneath a depth of around 1.2-1.4 m postdepositional losses can be neglected. In the area of investigation, only MSA- concentrations and postdepositional losses showed a distinct dependence on snow accumulation rate. Consequently, MSA- concentrations archived at this site should be significantly dependent on the variability of annual snow accumulation, and we recommend a corresponding correction. With a simple approach, we estimated the partial pressure of the free acids MSA, HNO3, and HCl on the basis of Henry's law assuming that ionic impurities of the bulk ice matrix are localized in a quasi-brine layer (QBL). In contrast to measurements, this approach predicts a nearly complete loss of MSA-, NO3 - , and Cl-.

Oxygen and silica flux from 9 multicorer profiles

Pore water profiles from 24 stations in the South Atlantic (located in the Guinea, Angola, Cape, Guyana, and Argentine basins) show good correlations of oxygen and silicon, suggesting microbially mediated dissolution of biogenic silica. We used simple analytical transport and reaction models to show the tight coupling of the reconstructed process kinetics of aerobic respiration and silicon regeneration. A generic transport and reaction model successfully reproduced the majority of Si pore water profiles from aerobic respiration rates, confirming that the dissolution of biogenic silica (BSi) occurs proportionally to O2 consumption. Possibly limited to well-oxygenated sediments poor in BSi, benthic Si fluxes can be inferred from O2 uptake with satisfactory accuracy. Compared to aerobic respiration kinetics, the solubility of BSi emerged as a less influential parameter for silicon regeneration. Understanding the role of bacteria for silicon regeneration requires further investigations, some of which are outlined. The proposed aerobic respiration control of benthic silicon cycling is suitable for benthic-pelagic models. The empirical relation of BSi dissolution to aerobic respiration can be used for regionalization assessments and estimates of the silicon budget to increase the understanding of global primary and export production patterns.

Swift heavy ion damage to sodium chloride: synergy between excitation and thermal spikes

Systematic data on the effect of irradiation with swift ions (Zn at 735 MeV and Xe at 929 MeV) on NaCl single crystals have been analysed in terms of a synergetic two-spike approach (thermal and excitation spikes). The coupling of the two spikes, simultaneously generated by the irradiation, contributes to the operation of a non-radiative exciton decay model as proposed for purely ionization damage. Using this scheme, we have accounted for the π-emission yield of self-trapped excitons and its temperature dependence under ion-beam irradiation. Moreover, the initial production rates of F-centre growth have also been reasonably simulated for irradiation at low temperatures ( < 100 K), where colour centre annealing and aggregation can be neglected.

Geometric diagnostics of complex patterns: Spiral defect chaos

Motivated by the observation of spiral patterns in a wide range of physical, chemical, and biological systems, we present an automated approach that aims at characterizing quantitatively spiral-like elements in complex stripelike patterns. The approach provides the location of the spiral tip and the size of the spiral arms in terms of their arc length and their winding number. In addition, it yields the number of pattern components (Betti number of order 1), as well as their size and certain aspects of their shape. We apply the method to spiral defect chaos in thermally driven Rayleigh- Bénard convection and find that the arc length of spirals decreases monotonically with decreasing Prandtl number of the fluid and increasing heating. By contrast, the winding number of the spirals is nonmonotonic in the heating. The distribution function for the number of spirals is significantly narrower than a Poisson distribution. The distribution function for the winding number shows approximately an exponential decay. It depends only weakly on the heating, but strongly on the Prandtl number. Large spirals arise only for larger Prandtl numbers. In this regime the joint distribution for the spiral length and the winding number exhibits a three-peak structure, indicating the dominance of Archimedean spirals of opposite sign and relatively straight sections. For small Prandtl numbers the distribution function reveals a large number of small compact pattern components.