Bayesian model comparison in cosmology with Population Monte Carlo

We use Bayesian model selection techniques to test extensions of the standard flat LambdaCDM paradigm. Dark-energy and curvature scenarios, and primordial perturbation models are considered. To that end, we calculate the Bayesian evidence in favour of each model using Population Monte Carlo (PMC), a new adaptive sampling technique which was recently applied in a cosmological context. The Bayesian evidence is immediately available from the PMC sample used for parameter estimation without further computational effort, and it comes with an associated error evaluation. Besides, it provides an unbiased estimator of the evidence after any fixed number of iterations and it is naturally parallelizable, in contrast with MCMC and nested sampling methods. By comparison with analytical predictions for simulated data, we show that our results obtained with PMC are reliable and robust. The variability in the evidence evaluation and the stability for various cases are estimated both from simulations and from data. For the cases we consider, the log-evidence is calculated with a precision of better than 0.08. Using a combined set of recent CMB, SNIa and BAO data, we find inconclusive evidence between flat LambdaCDM and simple dark-energy models. A curved Universe is moderately to strongly disfavoured with respect to a flat cosmology. Using physically well-motivated priors within the slow-roll approximation of inflation, we find a weak preference for a running spectral index. A Harrison-Zel'dovich spectrum is weakly disfavoured. With the current data, tensor modes are not detected; the large prior volume on the tensor-to-scalar ratio r results in moderate evidence in favour of r=0.

Impacts of level of utilisation by grazing on an Astrebla (Mitchell grass) grassland in north-western Queensland between 1984 and 2010. 1. Herbage mass and population dynamics of Astrebla spp

Managing large variations in herbage production, resulting from highly variable seasonal rainfall, provides a major challenge for the sustainable management of Astrebla (Mitchell grass) grasslands in Australia. A grazing study with sheep was conducted between 1984 and 2010 on an Astrebla grassland in northern Queensland to describe the effects of a range of levels of utilisation of the herbage at the end of the summer growing season (April–May in northern Australia) on the sustainability of these grasslands. In unreplicated paddocks, sheep numbers were adjusted annually to achieve 0, 10, 20, 30, 50 and 80% utilisation of the herbage mass at the end of the summer over the ensuing 12 months. Higher levels of utilisation reduced both total and Astrebla spp. herbage mass because of the effects of higher utilisation on Astrebla spp. and this effect was accentuated by drought. The tussock density of Astrebla spp. varied widely among years but with few treatment differences until 2005 when density was reduced at the 50% level of utilisation. A major change in density resulted from a large recruitment of Astrebla spp. in 1989 that influenced its density for the remainder of the study. Basal area of the tussocks fluctuated among years, with increases due to rainfall and decreases during droughts. Seasonal rainfall was more influential than level of utilisation in changes to the basal area of perennial grasses. Drought resulted in the death of Astrebla spp. tussocks and this effect was accentuated at higher levels of utilisation. A series of three grazing exclosures were used to examine the recovery of the density and basal area of Astrebla spp. after it had been reduced by 80% utilisation over the preceding 9 years. This recovery study indicated that, although grazing exclusion was useful in the recovery of Astrebla spp., above-average rainfall was the major factor driving increases in the basal area of perennial grasses. Spring values of the Southern Oscillation Index and associated rainfall probabilities were considered to have potential for understanding the dynamics of Astrebla spp. It was concluded that Astrebla grassland remained sustainable after 26 years when grazed at up to 30% utilisation, while, at 50% utilisation, they became unsustainable after 20 years. Results from this study emphasised the need to maintain the population of Astrebla spp. tussocks.

The (Non)Effects of Lethal Population Control on the Diet of Australian Dingoes

Top-predators contribute to ecosystem resilience, yet individuals or populations are often subject to lethal control to protect livestock, managed game or humans from predation. Such management actions sometimes attract concern that lethal control might affect top-predator function in ways ultimately detrimental to biodiversity conservation. The primary function of a predator is predation, which is often investigated by assessing their diet. We therefore use data on prey remains found in 4,298 Australian dingo scats systematically collected from three arid sites over a four year period to experimentally assess the effects of repeated broad-scale poison-baiting programs on dingo diet. Indices of dingo dietary diversity and similarity were either identical or near-identical in baited and adjacent unbaited treatment areas in each case, demonstrating no control-induced change to dingo diets. Associated studies on dingoes' movement behaviour and interactions with sympatric mesopredators were similarly unaffected by poison-baiting. These results indicate that mid-sized top-predators with flexible and generalist diets (such as dingoes) may be resilient to ongoing and moderate levels of population control without substantial alteration of their diets and other related aspects of their ecological function.

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.

Proportion of Jews to general population of the European provinces of Russia

Digital image

Population density of the European provinces of Russia

Digital image

Livelihood, land use and environment interactions in the highlands of East Africa

This study aims at improving understanding of the interactions of livelihoods and the environment focusing on both socio-economic and biodiversity implications of land use change in the context of population pressure, global and local markets, climate change, cultural and regional historical factors in the highlands of East Africa. The study is based on three components (1) two extensive livelihood surveys, one on Mt. Kilimanjaro in Tanzania and the other in the Taita Hills of Kenya, (2) a land use change study of the southern slopes of Mt. Kilimanjaro focusing on land use trends between 1960s and 1980s and 1980s and 2000 and (3) a bird diversity study focusing on the potential impacts of the future land use change on birds in the main land use types on the slopes and the adjacent plains of Mt. Kilimanjaro. In addition, information on the highlands in Embu and the adjacent lowlands in Mbeere of Kenya are added to the discussion. Some general patterns of livelihood, land use and environment interactions can be found in the three sites. However, the linkages are very complex. Various external factors at different times in history have influenced most of the major turning points. Farmers continually make small adaptations to their farming practices, but the locally conceived alternatives are too few. Farmers lack specific information and knowledge on the most suitable crops, market opportunities and the quality requirements for growing the crops for markets. Population growth emerges as the most forceful driver of land use and environmental change. The higher altitudes have become extremely crowded with population densities in some areas higher than typical urban population densities. Natural vegetation has almost totally been replaced by farmland. Decreasing farm size due to population pressure is currently threatening the viability of whole farming systems. In addition, capital-poor intensification has lead to soil fertility depletion. Agricultural expansion to the agriculturally marginal lowlands has created a new and distinct group of farmers struggling constantly with climate variability causing frequent crop failures. Extensification to the fragile drylands is the major cause of fragmentation and loss of wildlife habitat. The linkages between livelihoods, land use and the environment generally point to degradation of the environment leading to reduced environmental services and ecosystem functions. There is no indication that the system is self-regulating in this respect. Positive interventions will be needed to maintain ecosystem integrity.

Population dynamics of Thaumastocoris peregrinus in Eucalyptus plantations of South Africa

Thaumastocoris peregrinus is a sap-sucking insect that infests non-native Eucalyptus plantations in Africa, New Zealand, South America and parts of Southern Europe, in addition to street trees in parts of its native range of Australia. In South Africa, pronounced fluctuations in the population densities have been observed. To characterise spatiotemporal variability in T. peregrinus abundance and the factors that might influence it, we monitored adult population densities at six sites in the main eucalypt growing regions of South Africa. At each site, twenty yellow sticky traps were monitored weekly for 30 months, together with climatic data. We also characterised the influence of temperature on growth and survival experimentally and used this to model how temperature may influence population dynamics. T. peregrinus was present throughout the year at all sites, with annual site-specific peaks in abundance. Peaks occurred during autumn (February–April) for the Pretoria site, summer (November–January) for the Zululand site and spring (August–October) for the Tzaneen, Sabie and Piet Retief monitoring sites. Temperature (both experimental and field-collected), humidity and rainfall were mostly weakly, or not at all, associated with population fluctuations. It is clear that a complex interaction of these and other factors (e.g. host quality) influence population fluctuations in an annual, site specific cycle. The results obtained not only provide insights into the biology of T. peregrinus, but will also be important for future planning of monitoring and control programs using semiochemicals, chemical insecticides or biological control agents.

Effect of water quality and season on the population dynamics of Cabomba caroliniana in subtropical Queensland, Australia

Cabomba caroliniana is a submersed aquatic macrophyte that originates from the Americas and is currently invading temperate, subtropical, and tropical freshwater habitats around the world. Despite being a nuisance in many countries, little is known about its ecology. We monitored C. caroliniana populations in three reservoirs in subtropical Queensland, Australia, over 5.5 years. Although biomass, stem length, and plant density of the C. caroliniana stands fluctuated over time, they did not exhibit clear seasonal patterns. Water depth was the most important environmental factor explaining C. caroliniana abundance. Plant biomass was greatest at depths from 2–4 m and rooted plants were not found beyond 5 m. Plant density was greatest in shallow water and decreased with depth, most likely as a function of decreasing light and increasing physical stress. We tested the effect of a range of water physico-chemical parameters. The concentration of phosphorus in the water column was the variable that explained most of the variation in C. caroliniana population parameters. We found that in subtropical Australia, C. caroliniana abundance does not appear to be affected by seasonal conditions but is influenced by other environmental variables such as water depth and nutrient loading. Therefore, further spread will more likely be governed by local habitat rather than climatic conditions.

Environmental change and anthropogenic impact on lake sediments during the Holocene in the Finnish - Karelian inland area

This thesis discusses the prehistoric human disturbance during the Holocene by means of case studies using detailed high-resolution pollen analysis from lake sediment. The four lakes studied are situated between 61o 40' and 61o 50' latitudes in the Finnish Karelian inland area and vary between 2.4 and 28.8 ha in size. The existence of Early Metal Age population was one important question. Another study question concerned the development of grazing, and the relationship between slash-and-burn cultivation and permanent field cultivation. The results were presented as pollen percentages and pollen concentrations (grains cm 3). Accumulation values (grains cm 2 yr 1) were calculated for Lake Nautajärvi and Lake Orijärvi sediment, where the sediment accumulation rate was precisely determined. Sediment properties were determined using loss-on-ignition (LOI) and magnetic susceptibility (k). Dating methods used include both conventional and AMS 14C determinations, paleomagnetic dating and varve choronology. The isolation of Lake Kirjavalampi on the northern shore of Lake Ladoga took place ca. 1460 1300 BC. The long sediment cores from Finland, Lake Kirkkolampi and Lake Orijärvi in southeastern Finland and Lake Nautajärvi in south central Finland all extended back to the Early Holocene and were isolated from the Baltic basin ca. 9600 BC, 8600 BC and 7675 BC, respectively. In the long sediment cores, the expansion of Alnus was visible between 7200 - 6840 BC. The spread of Tilia was dated in Lake Kirkkolampi to 6600 BC, in Lake Orijärvi to 5000 BC and at Lake Nautajärvi to 4600 BC. Picea is present locally in Lake Kirkkolampi from 4340 BC, in Lake Orijärvi from 6520 BC and in Lake Nautajärvi from 3500 BC onwards. The first modifications in the pollen data, apparently connected to anthropogenic impacts, were dated to the beginning of the Early Metal Period, 1880 1600 BC. Anthropogenic activity became clear in all the study sites by the end of the Early Metal Period, between 500 BC AD 300. According to Secale pollen, slash-and-burn cultivation was practised around the eastern study lakes from AD 300 600 onwards, and at the study site in central Finland from AD 880 onwards. The overall human impact, however, remained low in the studied sites until the Late Iron Age. Increasing human activity, including an increase in fire frequency was detected from AD 800 900 onwards in the study sites in eastern Finland. In Lake Kirkkolampi, this included cultivation on permanent fields, but in Lake Orijärvi, permanent field cultivation became visible as late as AD 1220, even when the macrofossil data demonstrated the onset of cultivation on permanent fields as early as the 7th century AD. On the northern shore of Lake Ladoga, local activity became visible from ca. AD 1260 onwards and at Lake Nautajärvi, sediment the local occupation was traceable from 1420 AD onwards. The highest values of Secale pollen were recorded both in Lake Orijärvi and Lake Kirjavalampi between ca. AD 1700 1900, and could be associated with the most intensive period of slash-and-burn from AD 1750 to 1850 in eastern Finland.

Effect of water quality and season on the population dynamics of Cabomba caroliniana in subtropical Queensland, Australia

Cabomba caroliniana is a submersed aquatic macrophyte that originates from the Americas and is currently invading temperate, subtropical, and tropical freshwater habitats around the world. Despite being a nuisance in many countries, little is known about its ecology. We monitored C. caroliniana populations in three reservoirs in subtropical Queensland, Australia, over 5.5 years. Although biomass, stem length, and plant density of the C. caroliniana stands fluctuated over time, they did not exhibit clear seasonal patterns. Water depth was the most important environmental factor explaining C. caroliniana abundance. Plant biomass was greatest at depths from 2–4 m and rooted plants were not found beyond 5 m. Plant density was greatest in shallow water and decreased with depth, most likely as a function of decreasing light and increasing physical stress. We tested the effect of a range of water physico-chemical parameters. The concentration of phosphorus in the water column was the variable that explained most of the variation in C. caroliniana population parameters. We found that in subtropical Australia, C. caroliniana abundance does not appear to be affected by seasonal conditions but is influenced by other environmental variables such as water depth and nutrient loading. Therefore, further spread will more likely be governed by local habitat rather than climatic conditions.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Optimising initial population density, growth time and nitrogen nutrition for assessing resistance of wheat cultivars to root-lesion nematode (Pratylenchus thornei)

Pratylenchus thornei is a major pathogen of wheat in Australia. Two glasshouse experiments with four wheat cultivars that had different final populations (Pf) of P. thornei in the field were used to optimise conditions for assessing resistance. With different initial populations (Pi) ranging up to 5250 P. thornei/kg soil, Pf of P. thornei increased to 16 weeks after sowing, and then decreased at 20 weeks in some cultivar x Pi combinations. The population dynamics of P. thornei up to 16 weeks were best described by a modified exponential equation P f (t) = aP i e kt where P f (t) is the final population density at time t, P i is the initial population density, a is the proportion of P i that initiates population development, and k is the intrinsic rate of increase of the population. The cultivar GS50a had very low k values at Pi of 5250 and 1050 indicating its resistance, Suneca and Potam had high k values indicating susceptibility, whereas intolerant Gatcher had a low value at the higher Pi and a high value at the lower Pi. Nitrate fertiliser increased plant growth and Pf values of susceptible cultivars, but in unplanted soil it decreased Pf. Nematicide (aldicarb 5 mg/kg soil) killed P. thornei more effectively in planted than in unplanted soil and increased plant growth particularly in the presence of N fertiliser. In both experiments, the wheat cultivars Suneca and Potam were more susceptible than the cultivar GS50a reflecting field results. The method chosen to discriminate wheat cultivars was to assess Pf after growth for 16 weeks in soil with Pi ~1050–5250 P. thornei/kg soil and fertilised with 200 mg NO3–N/kg soil.