TIGGE: Comparison of the prediction of northern hemisphere extratropical cyclones by different ensemble prediction systems

The Observing System Research and Predictability Experiment (THORPEX) Interactive Grand Global Ensemble (TIGGE) is a World Weather Research Programme project. One of its main objectives is to enhance collaboration on the development of ensemble prediction between operational centers and universities by increasing the availability of ensemble prediction system (EPS) data for research. This study analyzes the prediction of Northern Hemisphere extratropical cyclones by nine different EPSs archived as part of the TIGGE project for the 6-month time period of 1 February 2008–31 July 2008, which included a sample of 774 cyclones. An objective feature tracking method has been used to identify and track the cyclones along the forecast trajectories. Forecast verification statistics have then been produced [using the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis as the truth] for cyclone position, intensity, and propagation speed, showing large differences between the different EPSs. The results show that the ECMWF ensemble mean and control have the highest level of skill for all cyclone properties. The Japanese Meteorological Administration (JMA), the National Centers for Environmental Prediction (NCEP), the Met Office (UKMO), and the Canadian Meteorological Centre (CMC) have 1 day less skill for the position of cyclones throughout the forecast range. The relative performance of the different EPSs remains the same for cyclone intensity except for NCEP, which has larger errors than for position. NCEP, the Centro de Previsão de Tempo e Estudos Climáticos (CPTEC), and the Australian Bureau of Meteorology (BoM) all have faster intensity error growth in the earlier part of the forecast. They are also very underdispersive and significantly underpredict intensities, perhaps due to the comparatively low spatial resolutions of these EPSs not being able to accurately model the tilted structure essential to cyclone growth and decay. There is very little difference between the levels of skill of the ensemble mean and control for cyclone position, but the ensemble mean provides an advantage over the control for all EPSs except CPTEC in cyclone intensity and there is an advantage for propagation speed for all EPSs. ECMWF and JMA have an excellent spread–skill relationship for cyclone position. The EPSs are all much more underdispersive for cyclone intensity and propagation speed than for position, with ECMWF and CMC performing best for intensity and CMC performing best for propagation speed. ECMWF is the only EPS to consistently overpredict cyclone intensity, although the bias is small. BoM, NCEP, UKMO, and CPTEC significantly underpredict intensity and, interestingly, all the EPSs underpredict the propagation speed, that is, the cyclones move too slowly on average in all EPSs.

Instability of surface-temperature filaments in strain and shear

The effects of uniform straining and shearing on the stability of a surface quasi-geostrophic temperature filament are investigated. Straining is shown to stabilize perturbations for wide filaments but only for a finite time until the filament thins to a critical width, after which some perturbations can grow. No filament can be stabilized in practice, since there are perturbations that can grow large for any strain rate. The optimally growing perturbations, defined as solutions that reach a certain threshold amplitude first, are found numerically for a wide range of parameter values. The radii of the vortices formed through nonlinear roll-up are found to be proportional to θ/s, where θ is the temperature anomaly of the filament and s the strain rate, and are not dependent on the initial size of the filament. Shearing is shown to reduce the normal-mode growth rates, but it cannot stabilize them completely when there are temperature discontinuities in the basic state; smooth filaments can be stabilized completely by shearing and a simple scaling argument provides the shear rate required. Copyright © 2010 Royal Meteorological Society

A shift to randomness of brain oscillations in people with autism

BACKGROUND: Resting-state functional magnetic resonance imaging (fMRI) enables investigation of the intrinsic functional organization of the brain. Fractal parameters such as the Hurst exponent, H, describe the complexity of endogenous low-frequency fMRI time series on a continuum from random (H = .5) to ordered (H = 1). Shifts in fractal scaling of physiological time series have been associated with neurological and cardiac conditions. METHODS: Resting-state fMRI time series were recorded in 30 male adults with an autism spectrum condition (ASC) and 33 age- and IQ-matched male volunteers. The Hurst exponent was estimated in the wavelet domain and between-group differences were investigated at global and voxel level and in regions known to be involved in autism. RESULTS: Complex fractal scaling of fMRI time series was found in both groups but globally there was a significant shift to randomness in the ASC (mean H = .758, SD = .045) compared with neurotypical volunteers (mean H = .788, SD = .047). Between-group differences in H, which was always reduced in the ASC group, were seen in most regions previously reported to be involved in autism, including cortical midline structures, medial temporal structures, lateral temporal and parietal structures, insula, amygdala, basal ganglia, thalamus, and inferior frontal gyrus. Severity of autistic symptoms was negatively correlated with H in retrosplenial and right anterior insular cortex. CONCLUSIONS: Autism is associated with a small but significant shift to randomness of endogenous brain oscillations. Complexity measures may provide physiological indicators for autism as they have done for other medical conditions.

Adaptive approximate Bayesian computation

Sequential techniques can enhance the efficiency of the approximate Bayesian computation algorithm, as in Sisson et al.'s (2007) partial rejection control version. While this method is based upon the theoretical works of Del Moral et al. (2006), the application to approximate Bayesian computation results in a bias in the approximation to the posterior. An alternative version based on genuine importance sampling arguments bypasses this difficulty, in connection with the population Monte Carlo method of Cappe et al. (2004), and it includes an automatic scaling of the forward kernel. When applied to a population genetics example, it compares favourably with two other versions of the approximate algorithm.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Life-history evolution under a production constraint

The recently formulated metabolic theory of ecology has profound implications for the evolution of life histories. Metabolic rate constrains the scaling of production with body mass, so that larger organisms have lower rates of production on a mass-specific basis than smaller ones. Here, we explore the implications of this constraint for life-history evolution. We show that for a range of very simple life histories, Darwinian fitness is equal to birth rate minus death rate. So, natural selection maximizes birth and production rates and minimizes death rates. This implies that decreased body size will generally be favored because it increases production, so long as mortality is unaffected. Alternatively, increased body size will be favored only if it decreases mortality or enhances reproductive success sufficiently to override the preexisting production constraint. Adaptations that may favor evolution of larger size include niche shifts that decrease mortality by escaping predation or that increase fecundity by exploiting new abundant food sources. These principles can be generalized to better understand the intimate relationship between the genetic currency of evolution and the metabolic currency of ecology.

The allometry of ornaments and weapons

Exaggerated male traits that have evolved under sexual selection include ornaments to attract mates and weapons to deter rivals. Data from studies of many such traits in diverse kinds of organisms show that they almost universally exhibit positive allometries. Both ornaments and weapons increase disproportionately with overall body size, resulting in scaling exponents within species that are consistently > 1.0 and usually in the range 1.5-2.5. We show how scaling exponents reflect the relative fitness advantages of ornaments vs. somatic growth by using a simple mathematical model of resource allocation during ontogeny. Because the scaling exponents are similar for the different taxonomic groups, it follows that the fitness advantages of investing in ornaments also are similar. The model also shows how selection for ornaments influences body size at first reproduction and explains why interspecific allometries have consistently lower exponents than intraspecific ones.

Ecology of sponge assemblages (Porifera) in the Wakatobi region, south-east Sulawesi, Indonesia: richness and abundance

The influence of sedimentation, depth and substratum angle on sponge assemblages in the Wakatobi region, south-eastern Sulawesi, Indonesia was considered. Sponge assemblages were sampled from two reef localities. The first reef (Sampela) was highly impacted by high sedimentation rates with fine sediment particles that settle slowly, while the second (Hoga) experienced only fast settling coarse sediment with lower overall sedimentation rates. Sponge assemblages were sampled (area occupied and numbers) on the reef fiat (0 m) and at 5 (reef crest), 10 and 15 m (15 m at Hoga only). Some significant (P < 0.001) differences were observed in the area occupied and the number of sponge patches between surface angles and sites. Significantly lower (t > 4.61, df = 9, P < 0.001) sponge numbers, percentage cover and richness were associated with the reef flat at both sites compared with all other depths at each site, with the exception of abundance of sponges on the reef flat at Sampela, which was much greater than at any other depth sampled. Species richness increased with depth at both sites but differences between surface angles were only recorded at Sampela, with higher species richness being found on vertical, inclined and horizontal surfaces respectively A total of 100 sponge species (total area sampled 52.5 m(2)) was reported from the two sites, with 58 species found at Sampela and 71 species at Hoga (41% of species shared). Multi-dimensional scaling (MDS) indicated differences in assemblage structure between sites and most depth intervals, but not substratum angles. A number of biological (e.g. competition and predation) and physical (e.g. sedimentation and aerial exposure) factors were considered to control sponge abundance and richness. Unexpectedly a significant (F-1,F-169 = 148.98, P < 0.001) positive linear relationship was found between sponge density and area occupied. In areas of high sponge coverage, the number of patches was also high, possibly due to fragmentation of large sponges produced as a result of predation and physical disturbance. The MDS results were also the same whether sponge numbers or percentage cover estimates were used, suggesting that although these different approaches yield different sorts of information, the same assemblage structure can be identified.

Are the spatial patterns of weeds scale-invariant?

In previous empirical and modelling studies of rare species and weeds, evidence of fractal behaviour has been found. We propose that weeds in modern agricultural systems may be managed close to critical population dynamic thresholds, below which their rates of increase will be negative and where scale-invariance may be expected as a consequence. We collected detailed spatial data on five contrasting species over a period of three years in a primarily arable field. Counts in 20×20 cm contiguous quadrats, 225,000 in 1998 and 84,375 thereafter, could be re-structured into a wide range of larger quadrat sizes. These were analysed using three methods based on correlation sum, incidence and conditional incidence. We found non-trivial scale invariance for species occurring at low mean densities and where they were strongly aggregated. The fact that the scale-invariance was not found for widespread species occurring at higher densities suggests that the scaling in agricultural weed populations may, indeed, be related to critical phenomena.

Temporal and spatial variability of mobile fauna on a submarine cliff and boulder scree complex: a community in flux

The mobile component of a community inhabiting a submarine boulder scree/cliff was investigated at Lough Hyne, Ireland at dawn, midday, dusk and night over a 1-week period. Line transects (50 m) were placed in the infralittoral (6 m) and circumlittoral (18 m) zones and also the interface between these two zones (12 m). The dominant mobile fauna of this cliff consisted of echinoderms (6 species), crustaceans (10 species) and fish (23 species). A different component community was identified at each time/depth interval using Multi-Dimensional Scaling (MDS) even though both species diversity (Shannon-Wiener indices) and richness (number of species) remained constant. These changes in community composition provided indirect evidence for migration by these mobile organisms. However, little evidence was found for migration between different zones with the exception of the several wrasse species. These species were observed to spend the daytime foraging in the deeper zone, but returned to the upper zone at night presumably for protection from predators. For the majority of species, migration was considered to occur to cryptic habitats such as holes and crevices. The number of organisms declined during the night, although crustacean numbers peaked, while fish and echinoderms were most abundant during day, possibly due to predator-prey interactions. This submarine community is in a state of flux, whereby, community characteristics, including trophic and energetic relationships, varied over small temporal (daily) and spatial (m) scales.

Multivariate relationships between groundwater chemistry and toxicity in an urban aquifer

Multivariate statistical methods were used to investigate file Causes of toxicity and controls on groundwater chemistry from 274 boreholes in an Urban area (London) of the United Kingdom. The groundwater was alkaline to neutral, and chemistry was dominated by calcium, sodium, and Sulfate. Contaminants included fuels, solvents, and organic compounds derived from landfill material. The presence of organic material in the aquifer caused decreases in dissolved oxygen, sulfate and nitrate concentrations. and increases in ferrous iron and ammoniacal nitrogen concentrations. Pearson correlations between toxicity results and the concentration of individual analytes indicated that concentrations of ammoinacal nitrogen, dissolved oxygen, ferrous iron, and hydrocarbons were important where present. However, principal component and regression analysis suggested no significant correlation between toxicity and chemistry over the whole area. Multidimensional Scaling was used to investigate differences in sites caused by historical use, landfill gas status, or position within the sample area. Significant differences were observed between sites with different historical land use and those with different gas status. Examination of the principal component matrix revealed that these differences are related to changes in the importance of reduced chemical species.

Simulating the formation of secondary organic aerosol from the photooxidation of aromatic hydrocarbons

The formation and composition of secondary organic aerosol (SOA) from the photooxidation of benzene, p-xylene, and 1,3,5-trimethylbenzene has been simulated using the Master Chemical Mechanism version 3.1 (MCM v3.1) coupled to a representation of the transfer of organic material from the gas to particle phase. The combined mechanism was tested against data obtained from a series of experiments conducted at the European Photoreactor (EUPHORE) outdoor smog chamber in Valencia, Spain. Simulated aerosol mass concentrations compared reasonably well with the measured SOA data only after absorptive partitioning coefficients were increased by a factor of between 5 and 30. The requirement of such scaling was interpreted in terms of the occurrence of unaccounted-for association reactions in the condensed organic phase leading to the production of relatively more nonvolatile species. Comparisons were made between the relative aerosol forming efficiencies of benzene, toluene, p-xylene, and 1,3,5-trimethylbenzene, and differences in the OH-initiated degradation mechanisms of these aromatic hydrocarbons. A strong, nonlinear relationship was observed between measured (reference) yields of SOA and (proportional) yields of unsaturated dicarbonyl aldehyde species resulting from ring-fragmenting pathways. This observation, and the results of the simulations, is strongly suggestive of the involvement of reactive aldehyde species in association reactions occurring in the aerosol phase, thus promoting SOA formation and growth. The effect of NO, concentrations on SOA formation efficiencies (and formation mechanisms) is discussed.

Simulating the detailed chemical composition of secondary organic aerosol formed on a regional scale during the TORCH 2003 campaign in the southern UK

Following on from the companion study (Johnson et al., 2006), a photochemical trajectory model (PTM) has been used to simulate the chemical composition of organic aerosol for selected events during the 2003 TORCH (Tropospheric Organic Chemistry Experiment) field campaign. The PTM incorporates the speciated emissions of 124 nonmethane anthropogenic volatile organic compounds (VOC) and three representative biogenic VOC, a highly-detailed representation of the atmospheric degradation of these VOC, the emission of primary organic aerosol (POA) material and the formation of secondary organic aerosol (SOA) material. SOA formation was represented by the transfer of semi and non-volatile oxidation products from the gas-phase to a condensed organic aerosol-phase, according to estimated thermodynamic equilibrium phase-partitioning characteristics for around 2000 reaction products. After significantly scaling all phase-partitioning coefficients, and assuming a persistent background organic aerosol (both required in order to match the observed organic aerosol loadings), the detailed chemical composition of the simulated SOA has been investigated in terms of intermediate oxygenated species in the Master Chemical Mechanism, version 3.1 ( MCM v3.1). For the various case studies considered, 90% of the simulated SOA mass comprises between ca. 70 and 100 multifunctional oxygenated species derived, in varying amounts, from the photooxidation of VOC of anthropogenic and biogenic origin. The anthropogenic contribution is dominated by aromatic hydrocarbons and the biogenic contribution by alpha-and beta-pinene (which also constitute surrogates for other emitted monoterpene species). Sensitivity in the simulated mass of SOA to changes in the emission rates of anthropogenic and biogenic VOC has also been investigated for 11 case study events, and the results have been compared to the detailed chemical composition data. The role of accretion chemistry in SOA formation, and its implications for the results of the present investigation, is discussed.

A novel experimental method for measuring vergence and accommodation responses to the main near visual cues in typical and atypical groups

Binocular disparity, blur, and proximal cues drive convergence and accommodation. Disparity is considered to be the main vergence cue and blur the main accommodation cue. We have developed a remote haploscopic photorefractor to measure simultaneous vergence and accommodation objectively in a wide range of participants of all ages while fixating targets at between 0.3 and 2 m. By separating the three main near cues, we can explore their relative weighting in three-, two-, one-, and zero-cue conditions. Disparity can be manipulated by remote occlusion; blur cues manipulated by using either a Gabor patch or a detailed picture target; looming cues by either scaling or not scaling target size with distance. In normal orthophoric, emmetropic, symptom-free, naive visually mature participants, disparity was by far the most significant cue to both vergence and accommodation. Accommodation responses dropped dramatically if disparity was not available. Blur only had a clinically significant effect when disparity was absent. Proximity had very little effect. There was considerable interparticipant variation. We predict that relative weighting of near cue use is likely to vary between clinical groups and present some individual cases as examples. We are using this naturalistic tool to research strabismus, vergence and accommodation development, and emmetropization.

Stereoscopic acuity and observation distance

The amount of depth perceived from a fixed pattern of horizontal disparities varies with viewing distance. We investigated whether thresholds for discriminating stereoscopic corrugations at a range of spatial frequencies were also affected by viewing distance or whether they were determined solely by the angular disparity in the stimulus prior to scaling. Although thresholds were found to be determined primarily by disparity over a broad range of viewing distances, they were on average a factor of two higher at the shortest viewing distance (28.5 cm) than at larger viewing distances (57 to 450 cm). We found the same pattern of results when subjects' accommodation was arranged to be the same at all viewing distances. The change in thresholds at close distances is in the direction expected if subjects' performance is limited by a minimum perceived depth.