Location and scale mixtures of Gaussians with flexible tail behaviour: Properties, inference and application to multivariate clustering

The family of location and scale mixtures of Gaussians has the ability to generate a number of flexible distributional forms. The family nests as particular cases several important asymmetric distributions like the Generalized Hyperbolic distribution. The Generalized Hyperbolic distribution in turn nests many other well known distributions such as the Normal Inverse Gaussian. In a multivariate setting, an extension of the standard location and scale mixture concept is proposed into a so called multiple scaled framework which has the advantage of allowing different tail and skewness behaviours in each dimension with arbitrary correlation between dimensions. Estimation of the parameters is provided via an EM algorithm and extended to cover the case of mixtures of such multiple scaled distributions for application to clustering. Assessments on simulated and real data confirm the gain in degrees of freedom and flexibility in modelling data of varying tail behaviour and directional shape.

Large-scale characterization of drought pattern: a continent-wide modelling approach applied to the Australian wheatbelt – spatial and temporal trends

* Plant response to drought is complex, so that traits adapted to a specific drought type can confer disadvantage in another drought type. Understanding which type(s) of drought to target is of prime importance for crop improvement. * Modelling was used to quantify seasonal drought patterns for a check variety across the Australian wheatbelt, using 123 yr of weather data for representative locations and managements. Two other genotypes were used to simulate the impact of maturity on drought pattern. * Four major environment types summarized the variability in drought pattern over time and space. Severe stress beginning before flowering was common (44% of occurrences), with (24%) or without (20%) relief during grain filling. High variability occurred from year to year, differing with geographical region. With few exceptions, all four environment types occurred in most seasons, for each location, management system and genotype. * Applications of such environment characterization are proposed to assist breeding and research to focus on germplasm, traits and genes of interest for target environments. The method was applied at a continental scale to highly variable environments and could be extended to other crops, to other drought-prone regions around the world, and to quantify potential changes in drought patterns under future climates.

Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes

Despite international protection of white sharks Carcharodon carcharias, important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with 6 nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (F-ST = 0.142, p < 0.0001), implying female reproductive philopatry. This concurs with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited micro satellite markers (F-ST = 0.009, p < 0.05), suggesting that males may also exhibit some degree of reproductive philopatry; 5 sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal, or migration resulting in breeding, may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and close to thresholds at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least 1, possibly 2, orders of magnitude below our historical effective size estimates. Population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail.

Predictive modelling to aid the regional-scale management of a vertebrate pest

Extensive resources are allocated to managing vertebrate pests, yet spatial understanding of pest threats, and how they respond to management, is limited at the regional scale where much decision-making is undertaken. We provide regional-scale spatial models and management guidance for European rabbits (Oryctolagus cuniculus) in a 260,791 km(2) region in Australia by determining habitat suitability, habitat susceptibility and the effects of the primary rabbit management options (barrier fence, shooting and baiting and warren ripping) or changing predation or disease control levels. A participatory modelling approach was used to develop a Bayesian network which captured the main drivers of suitability and spread, which in turn was linked spatially to develop high resolution risk maps. Policy-makers, rabbit managers and technical experts were responsible for defining the questions the model needed to address, and for subsequently developing and parameterising the model. Habitat suitability was determined by conditions required for warren-building and by above-ground requirements, such as food and harbour, and habitat susceptibility by the distance from current distributions, habitat suitability, and the costs of traversing habitats of different quality. At least one-third of the region had a high probability of being highly suitable (support high rabbit densities), with the model supported by validation. Habitat susceptibility was largely restricted by the current known rabbit distribution. Warren ripping was the most effective control option as warrens were considered essential for rabbit persistence. The anticipated increase in disease resistance was predicted to increase the probability of moderately suitable habitat becoming highly suitable, but not increase the at-risk area. We demonstrate that it is possible to build spatial models to guide regional-level management of vertebrate pests which use the best available knowledge and capture fine spatial-scale processes.

The sensitivity of the northern hemisphere summer mean meridional circulation to changes in the large scale eddy forcing

A zonally averaged version of the Goddard Laboratory for Atmospheric Sciences (GLAS) climate model is used to study the sensitivity of the northern hemisphere (NH) summer mean meridional circulation to changes in the large scale eddy forcing. A standard solution is obtained by prescribing the latent heating field and climatological horizontal transports of heat and momentum by the eddies. The radiative heating and surface fluxes are calculated by model parameterizations. This standard solution is compared with the results of several sensitivity studies. When the eddy forcing is reduced to 0.5 times or increased to 1.5 times the climatological values, the strength of the Ferrel cells decrease or increase proportionally. It is also seen that such changes in the eddy forcing can influence the strength of theNH Hadley cell significantly. Possible impact of such changes in the large scale eddy forcing on the monsoon circulation via changes in the Hadley circulation is discussed. Sensitivity experiments including only one component of eddy forcing at a time show that the eddy momentum fluxes seem to be more important in maintaining the Ferrel cells than the eddy heat fluxes. In the absence of the eddy heat fluxes, the observed eddy momentum fluxes alone produce subtropical westerly jets which are weaker than those in the standard solution. On the other hand, the observed eddy heat fluxes alone produce subtropical westerly jets which are stronger than those in the standard solution.

Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification

In recent years, concern has arisen over the effects of increasing carbon dioxide (CO2) in the earth's atmosphere due to the burning of fossil fuels. One way to mitigate increase in atmospheric CO2 concentration and climate change is carbon sequestration to forest vegeta-tion through photosynthesis. Comparable regional scale estimates for the carbon balance of forests are therefore needed for scientific and political purposes. The aim of the present dissertation was to improve methods for quantifying and verifying inventory-based carbon pool estimates of the boreal forests in the mineral soils. Ongoing forest inventories provide a data based on statistically sounded sampling for estimating the level of carbon stocks and stock changes, but improved modelling tools and comparison of methods are still needed. In this dissertation, the entire inventory-based large-scale forest carbon stock assessment method was presented together with some separate methods for enhancing and comparing it. The enhancement methods presented here include ways to quantify the biomass of understorey vegetation as well as to estimate the litter production of needles and branches. In addition, the optical remote sensing method illustrated in this dis-sertation can be used to compare with independent data. The forest inventory-based large-scale carbon stock assessment method demonstrated here provided reliable carbon estimates when compared with independent data. Future ac-tivity to improve the accuracy of this method could consist of reducing the uncertainties regarding belowground biomass and litter production as well as the soil compartment. The methods developed will serve the needs for UNFCCC reporting and the reporting under the Kyoto Protocol. This method is principally intended for analysts or planners interested in quantifying carbon over extensive forest areas.

Ecological processes and large-scale climate relationships in northern coniferous forests

Ilmasto vaikuttaa ekologisiin prosesseihin eri tasoilla. Suuren mittakaavan ilmastoprosessit, yhdessä ilmakehän ja valtamerien kanssa, säätelevät paikallisia sääilmiöitä suurilla alueilla (mantereista pallopuoliskoihin). Tämä väistöskirja pyrkii selittämään kuinka suuren mittakaavan ilmasto on vaikuttanut tiettyihin ekologisiin prosesseihin pohjoisella havumetsäalueella. Valitut prosessit olivat puiden vuosilustojen kasvu, metsäpalojen esiintyminen ja vuoristomäntykovakuoriaisen aiheuttamat puukuolemat. Suuren mittakaavan ilmaston löydettiin vaikuttaneen näiden prosessien esiintymistiheyteen, kestoon ja levinneisyyteen keskeisten sään muuttujien välityksellä hyvin laajoilla alueilla. Tutkituilla prosesseilla oli vahva yhteys laajan mittakaavan ilmastoon. Yhteys on kuitenkin ollut hyvin dynaaminen ja muuttunut 1900-luvulla ilmastonmuutoksen aiheuttaessa muutoksia suuren mittakaavan ja alueellisten ilmastoprosessien välisiin sisäisiin suhteisiin.

Commercial-scale Alternaria brown spot resistance screening as the first step in breeding new mandarins for Australia

Rapid screening tests and an appreciation of the simple genetic control of Alternaria brown spot (ABS) susceptibility have existed for many years, and yet the application of this knowledge to commercial-scale breeding programs has been limited. Detached leaf assays were first demonstrated more than 40 years ago and reliable data suggesting a single gene determining susceptibility has been emerging for at least 20 years. However it is only recently that the requirement for genetic resistance in new hybrids has become a priority, following increased disease prevalence in Australian mandarin production areas previously considered too dry for the pathogen. Almost all of the high-fruit-quality parents developed so far by the Queensland-based breeding program are susceptible to ABS necessitating the screening of their progeny to avoid commercialisation of susceptible hybrids. This is done effectively and efficiently by spraying 3-6 month old hybrid seedlings with a spore suspension derived from a toxin-producing field isolate of Alternaria alternate, then incubating these seedlings in a cool room at 25°C and high humidity for 5 days. Susceptible seedlings show clear disease symptoms and are discarded. Analysis of observed and expected segregation ratios loosely support the hypothesis for a single dominant gene for susceptibility, but do not rule out the possibility of alternative genetic models. After implementing the routine screening for ABS resistance for three seasons we now have more than 20,000 hybrids growing in field progeny blocks that have been screened for resistance to the ABS disease.

Large-scale expansion of no-take closures within the Great Barrier Reef has not enhanced fishery production

A rare opportunity to test hypotheses about potential fishery benefits of large-scale closures was initiated in July 2004 when an additional 28.4% of the 348 000 km2 Great Barrier Reef (GBR) region of Queensland, Australia was closed to all fishing. Advice to the Australian and Queensland governments that supported this initiative predicted these additional closures would generate minimal (10%) initial reductions in both catch and landed value within the GBR area, with recovery of catches becoming apparent after three years. To test these predictions, commercial fisheries data from the GBR area and from the two adjacent (non-GBR) areas of Queensland were compared for the periods immediately before and after the closures were implemented. The observed means for total annual catch and value within the GBR declined from pre-closure (2000–2003) levels of 12 780 Mg and Australian $160 million, to initial post-closure (2005–2008) levels of 8143 Mg and $102 million; decreases of 35% and 36% respectively. Because the reference areas in the non-GBR had minimal changes in catch and value, the beyond-BACI (before, after, control, impact) analyses estimated initial net reductions within the GBR of 35% for both total catch and value. There was no evidence of recovery in total catch levels or any comparative improvement in catch rates within the GBR nine years after implementation. These results are not consistent with the advice to governments that the closures would have minimal initial impacts and rapidly generate benefits to fisheries in the GBR through increased juvenile recruitment and adult spillovers. Instead, the absence of evidence of recovery in catches to date currently supports an alternative hypothesis that where there is already effective fisheries management, the closing of areas to all fishing will generate reductions in overall catches similar to the percentage of the fished area that is closed.

Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

Small-scale capture, transport and tank adaptation of live, medium-sized Scombrids using "Tuna Tubes"

The transport of live fish is a crucial step to establish fish culture in captivity, and is especially challenging for species that have not been commonly cultured before, therefore transport and handling methods need to be optimized and tailored. This study describes the use of tuna tubes for small-scale transport of medium-sized pelagic fish from the Scombridae family. Tuna tubes are an array of vertical tubes that hold the fish, while fresh seawater is pumped up the tubes and through the fish mouth and gills, providing oxygen and removing wastes. In this study, 19 fish were captured using rod and line and 42% of the captured fish were transported alive in the custom-designed tuna tubes to an on-shore holding tank: five mackerel tuna (Euthynnus affinis) and three leaping bonito (Cybiosarda elegans). Out of these, just three (15.8% of total fish) acclimatized to the tank's condition. Based on these results, we discuss an improved design of the tuna tubes that has the potential to increase survival rates and enable a simple and low cost method of transporting of live pelagic fish.

Achieving medium-high density in low scale development: The Queensland experience in innovative ‘fine-grained’ urbanism

Worldwide population growth and economic agglomeration is driving increasing urban density within larger metropolitan conurbations. Population growth and housing diversity and affordability issues in Queensland have seen an increasing demand for more diverse and higher density development. Under Queensland’s flexible planning regulatory provisions, a level of ‘medium’ to ‘high density’ is being achieved by a focus on fine-grained urban design, low scale development, lot diversity, and delivery of single dwelling products. This for Queensland (and Australia) has been an unprecedented innovation in urban and dwelling design. Dwellings are being delivered on lots with zero regulatory minimum sizes providing for a range of new products including ‘apartments on the ground’. This paper reviews recent and nascent demonstrations of EDQ’s fine-grained urbanism principles, identifiable with historical ‘vernacular suburbanism’. The paper introduces and defines a concept of a ‘natural density’ linking human scale built form with walkability. The paper challenges the notion that (sub)urban development, outside major city centres, needs to be of a higher scale to achieve density and diversity aspirations. ‘Natural density’ provides a means of achieving the increasing demand for more diverse and higher density development.

A simple method for large-scale preparation of chick embryo tRNA

A simple method for preparing bulk quantities of tRNA from chick embryo has been developed. In this method chick embryos were homogenized in a buffer of pH 4.5, followed by deproteinization with phenol. The aqueous layer was allowed to separate under gravity. The resulting aqueous layer, after two more phenol treatments, was directly passed through a DEAE-cellulose column and the tRNA eluted therefrom with 1 Image NaCl. The tRNA prepared by this method was as active as the one prepared at neutral pH.