Inferential, non-parametric statistics to assess quality of probabilistic forecast systems

Many statistical forecast systems are available to interested users. In order to be useful for decision-making, these systems must be based on evidence of underlying mechanisms. Once causal connections between the mechanism and their statistical manifestation have been firmly established, the forecasts must also provide some quantitative evidence of `quality’. However, the quality of statistical climate forecast systems (forecast quality) is an ill-defined and frequently misunderstood property. Often, providers and users of such forecast systems are unclear about what ‘quality’ entails and how to measure it, leading to confusion and misinformation. Here we present a generic framework to quantify aspects of forecast quality using an inferential approach to calculate nominal significance levels (p-values) that can be obtained either by directly applying non-parametric statistical tests such as Kruskal-Wallis (KW) or Kolmogorov-Smirnov (KS) or by using Monte-Carlo methods (in the case of forecast skill scores). Once converted to p-values, these forecast quality measures provide a means to objectively evaluate and compare temporal and spatial patterns of forecast quality across datasets and forecast systems. Our analysis demonstrates the importance of providing p-values rather than adopting some arbitrarily chosen significance levels such as p < 0.05 or p < 0.01, which is still common practice. This is illustrated by applying non-parametric tests (such as KW and KS) and skill scoring methods (LEPS and RPSS) to the 5-phase Southern Oscillation Index classification system using historical rainfall data from Australia, The Republic of South Africa and India. The selection of quality measures is solely based on their common use and does not constitute endorsement. We found that non-parametric statistical tests can be adequate proxies for skill measures such as LEPS or RPSS. The framework can be implemented anywhere, regardless of dataset, forecast system or quality measure. Eventually such inferential evidence should be complimented by descriptive statistical methods in order to fully assist in operational risk management.

The risk of pollen-mediated gene flow from exotic Corymbia plantations into native Corymbia populations in Australia

Sectors of the forest plantation industry in Australia are set to expand in the near future using species or hybrids of the spotted gums (Corymbia, Section Politaria). Plantations of these taxa have already been introduced across temperate and subtropical Australia, representing locally exotic introductions from native stands in Queensland and New South Wales. A literature review was undertaken to provide insights into the potential for pollen-mediated gene flow from these plantations into native populations. Three factors suggest that such gene flow is likely; (1) interspecific hybridisation within the genus has frequently been recorded, including between distantly related species from different sections, (2) apparent high levels of vertebrate pollinator activity may result in plantation pollen being moved over hundreds of kilometres, (3) much of the plantation estate is being established among closely related taxa and therefore few barriers to gene flow are expected. Across Australia, 20 of the 100 native Corymbia taxa were found to have regional level co-occurrence with plantations. These were located most notably within regions of north-east New South Wales and south-east Queensland, however, co-occurrence was also found in south-west Western Australia and eastern Victoria. The native species found to have co-occurrence were then assessed for the presence of reproductive barriers at each step in the process of gene flow that may reduce the number of species at risk even further. The available data suggest three risk categories exist for Corymbia. The highest risk was for gene flow from plantations of spotted gums to native populations of spotted gums. This was based on the expected limited existence of pre- and post-zygotic barriers, substantial long-distance pollen dispersal and an apparent broad period of flowering in Corymbia citriodora subsp. variegata plantations. The following risk category focussed on gene flow from Corymbia torelliana × C. c. variegata hybrid plantations into native C. c. variegata, as the barriers associated with the production and establishment of F1 hybrids have been circumvented. For the lowest risk category, Corymbia plantations may present a risk to other non-spotted gum species, however, further investigation of the particular cross-combinations is required. A list of research directions is provided to better quantify these risks. Empirical data will need to be combined within a risk assessment framework that will not only estimate the likelihood of exotic gene flow, but also consider the conservation status/value of the native populations. In addition, the potential impacts of pollen flow from plantations will need to be weighed up against their various economic and environmental benefits.

Variable growth rates of the tropical estuarine fish barramundi Lates calcarifer (Bloch) under different freshwater flow conditions

Relationships between freshwater flows and growth rates of the opportunistic predatory finfish barramundi Lates calcarifer in a dry tropical estuary were examined using data from a long-term tag-recapture programme. Lagged effects were not investigated. After accounting for length at release, time at liberty and seasonal variation (e.g. winter, spring, summer and autumn), growth rates were significantly and positively related to fresh water flowing to the estuary. Effects were present at relatively low levels of freshwater flow (i.e. 2.15 m3 s-1, the 5th percentile of the mean flow rate experienced by fish in the study during time at liberty). The analysis, although correlative, provides quantitative evidence to support the hypothesis that freshwater flows are important in driving the productivity of estuaries and can improve growth of species high in the trophic chain.

Optical peanut yield monitor: Development and testing

An optical peanut yield monitor was developed, fabricated, and field-tested. The overall system includes an optical mass-flow sensor, a GPS receiver, and a data acquisition system. The concept for the mass-flow sensor is based on that of the cotton yield-monitor sensor developed previously by Thomasson and Sui (2000). A modified version of the sensor was designed to be specific to peanut mass-flow measurement. Field testing of the peanut yield monitor was conducted in Australia during the May 2003 harvest. After subsequent minor modifications, the system was more extensively tested in Mississippi in October of 2003 and November of 2004. Test results showed that the output of the peanut mass-flow sensor was very strongly correlated with the harvested load weight, and the system's performance was stable and reliable during the tests.

Probabilistic forecasts of the onset of the north Australian wet season

The amount and timing of early wet-season rainfall are important for the management of many agricultural industries in north Australia. With this in mind, a wet-season onset date is defined based on the accumulation of rainfall to a predefined threshold, starting from 1 September, for each square of a 1° gridded analysis of daily rainfall across the region. Consistent with earlier studies, the interannual variability of the onset dates is shown to be well related to the immediately preceding July-August Southern Oscillation index (SOI). Based on this relationship, a forecast method using logistic regression is developed to predict the probability that onset will occur later than the climatological mean date. This method is expanded to also predict the probabilities that onset will be later than any of a range of threshold dates around the climatological mean. When assessed using cross-validated hindcasts, the skill of the predictions exceeds that of climatological forecasts in the majority of locations in north Australia, especially in the Top End region, Cape York, and central Queensland. At times of strong anomalies in the July-August SOI, the forecasts are reliably emphatic. Furthermore, predictions using tropical Pacific sea surface temperatures (SSTs) as the predictor are also tested. While short-lead (July-August predictor) forecasts are more skillful using the SOI, long-lead (May-June predictor) forecasts are more skillful using Pacific SSTs, indicative of the longer-term memory present in the ocean.

The spatiotemporal dynamics of Tribolium castaneum (Herbst): adult flight and gene flow.

Tribolium castaneum (Herbst) has been used as a model organism to develop and test important ecological and evolutionary concepts and is also a major pest of grain and grain products globally. This beetle species is assumed to be a good colonizer of grain storages through anthropogenic movement of grain, and active dispersal by flight is considered unlikely. Studies using T. castaneum have therefore used confined or walking insects. We combine an ecological study of dispersal with an analysis of gene flow using microsatellites to investigate the spatiotemporal dynamics and adult flight of T. castaneum in an ecological landscape in eastern Australia. Flying beetles were caught in traps at grain storages and in fields at least 1 km from the nearest stored grain at regular intervals for an entire year. Significantly more beetles were trapped at storages than in fields, and almost no beetles were caught in native vegetation reserves many kilometres from the nearest stored grain. Genetic differentiation between beetles caught at storages and in fields was low, indicating that although T. castaneum is predominantly aggregated around grain storages, active dispersal takes place to the extent that significant gene flow occurs between them, mitigating founder effects and genetic drift. By combining ecological and molecular techniques, we reveal much higher levels of active dispersal through adult flight in T. castaneum than previously thought. We conclude that the implications of adult flight to previous and future studies on this model organism warrant consideration.

Flow impacts on estuarine finfish fisheries of the Gulf of Carpentaria

The estuaries of Australia s tropical rivers support commercial fisheries for finfish and shellfish valued at over $220 million per annum. There are also significant tourism-related and local recreational and indigenous fisheries for icon species such as barramundi. Development of water resources in Australia's Tropical Rivers region is being considered for the Flinders, Mitchell, McArthur, Roper, Daly and Victoria catchments. Greater knowledge of the freshwater requirements of tropical aquatic ecosystems, including estuaries is crucial, so that the communities of catchments where water resource development occurs can be assured that the downstream effects of such development are considered and managed based on the best available knowledge.

Honeydew as a source of nutrition for Lysiphlebus testaceipes (Cresson) (Hymenoptera: Braconidae): Effect of adult diet on lifespan and egg load

Parasitoid survival and fecundity is generally enhanced with access to carbohydrate food sources. In many agricultural ecosystems, there is often a scarcity of suitable carbohydrates for parasitoids. This study compared the suitability of aphid honeydew and buckwheat nectar as diet for the aphid parasitoid Lysiphlebus testaceipes. Wasp lifespan and egg load were both increased with access to carbohydrates, but no differences were detected between the various carbohydrates diets tested. As aphid honeydew is a sufficient source of nutrition and L.testaceipes is a short-lived species, adding additional sources of carbohydrates like floral nectar strips to the agricultural landscape is unlikely to significantly increase the biological control exerted by L.testaceipes. © 2012 Australian Entomological Society.

Reciprocal and advanced generation hybrids between Corymbia citriodora and C-torelliana: forestry breeding and the risk of gene flow

Corymbia F1 hybrids have high potential for plantation forestry; however, little is known of their reproductive biology and potential for genetic pollution of native Corymbia populations. This study aims to quantify the influence of reproductive isolating barriers on the success of novel reciprocal and advanced generation Corymbia hybrids. Two maternal taxa, Corymbia citriodora subsp. citriodora and Corymbia torelliana, were pollinated using five paternal taxa, C. citriodora subsp. citriodora, C. torelliana, one C. torelliana x C. citriodora subsp. citriodora hybrid and two C. torelliana x C. citriodora subsp. variegata hybrids. Pollen tube, embryo and seed development were assessed. Reciprocal hybridisation between C. citriodora subsp. citriodora and C. torelliana was successful. Advanced generation hybrids were also created when C. citriodora subsp. citriodora or C. torelliana females were backcrossed with F1 hybrid taxa. Prezygotic reproductive isolation was identified via reduced pollen tube numbers in the style and reduced numbers of ovules penetrated by pollen tubes. Reproductive isolation was weakest within the C. citriodora subsp. citriodora maternal taxon, with two hybrid backcrosses producing equivalent capsule and seed yields to the intraspecific cross. High hybridising potential was identified between all Corymbia species and F1 taxa studied. This provides opportunities for advanced generation hybrid breeding, allowing desirable traits to be amplified. It also indicates risks of gene flow between plantation and native Corymbia populations.

Pleistocene isolation, secondary introgression and restricted contemporary gene flow in the pig-eye shark, Carcharhinus amboinensis across northern Australia

We examine the structure and phylogeography of the pig-eye shark (Carcharhinus amboinensis) common in shallow coastal environments in northern Australia using two types of genetic markers, two mitochondrial (control region and NADH hydrogenase 4) and two nuclear (microsatellite and Rag 1) DNA. Two populations were defined within northern Australia on the basis of mitochondrial DNA evidence, but this result was not supported by nuclear microsatellite or Rag 1 markers. One possibility for this structure might be sex-specific behaviours such as female philopatry, although we argue it is doubtful that sufficient time has elapsed for any potential signatures from this behaviour to be expressed in nuclear markers. It is more likely that the observed pattern represents ancient populations repeatedly isolated and connected during episodic sea level changes during the Pleistocene epoch, until current day with restricted contemporary gene flow maintaining population genetic structure. Our results show the need for an understanding of both the history and ecology of a species in order to interpret patterns in genetic structure.

Riverine flow and spawning requirements of Macquaria ambigua oriens: implications for conservation and management

Understanding the life-history attributes of aquatic species is integral to the development of environmental-flow strategies in regulated river systems. This is particularly important when species are under continual and increasing pressure from water-resource development. In this study, the water temperature and flow requirements for spawning of the Fitzroy River golden perch (Macquaria ambigua oriens) were investigated over 4 years at 22 sites in the Fitzroy River catchment. Eggs, larvae and young-of-year (YOY) M. ambigua oriens were sampled on a variety of flow events to determine the environmental requirements for spawning. Eggs and larvae of M. ambigua oriens were detected during natural flow events generally with a minimum of 1.5 m river rise and duration of 7 days. Spawning was associated with the peak and/or recession of the first or second post-winter flow event where water temperatures exceeded 248 degrees C. Our data suggests that it is important to protect a range of flows, not just flood flows, as previously documented for this species. The interaction of spawning flows with existing and future water-resource development should be considered to ensure maintenance of the population viability of M. ambigua oriens.

Association between feline immunodeficiency virus (FIV) plasma viral RNA load, concentration of acute phase proteins and disease severity

Veterinarians have few tools to predict the rate of disease progression in FIV-infected cats. In contrast, in HIV infection, plasma viral RNA load and acute phase protein concentrations are commonly used as predictors of disease progression. This study evaluated these predictors in cats naturally infected with FIV. In older cats (>5 years), log10 FIV RNA load was higher in the terminal stages of disease compared to the asymptomatic stage. There was a significant association between log10 FIV RNA load and both log10 serum amyloid A concentration and age in unwell FIV-infected cats. This study suggests that viral RNA load and serum amyloid A warrant further investigation as predictors of disease status and prognosis in FIV-infected cats.

Gene flow from Corymbia hybrids in northern New South Wales

Some of the most productive taxa for forestry are interspecific F1 hybrids grown as exotics in the tropics and subtropics. Attributes of resilience, adaptability and vigour which engender the hybrids for wood production, may also exacerbate the risk they present from gene flow to native species gene pools or to local ecologies as weeds. To determine the biological and genetic factors that influence the extent of hybridisation, we examine the distribution and genealogy of wildlings surrounding plantings of locally-exotic Corymbia torelliana (Section Cadageria) near native C. henryi (Section Maculatae) in northern New South Wales. Our study showed pre-mating and pre- and post-zygotic barriers were incomplete, with in situ generation and natural establishment of both F1 hybrids (n = 3) and advanced generation hybrids under the disturbed conditions bordering native forest. As hybrids were located on alluvial flats exposed to frost, they also likely have an extended ecological range relative to native C. henryi. Despite the likely generation of large viable seed crops on F1 trees at the site over many years, establishment success and survival of advanced generation hybrids may be low, as only 5 immature and no mature advanced generation hybrids were identified. Propagation and genetic analysis of a seed crop from one F1 wildling showed early survival and vigour of seedlings in cultivation was high, and that at least for some F1 in some seasons, backcrossing to the recurrent native C. henryi parent is favoured (60%), whereas selfing (10%) and crossing with other F1 (30%) was less frequent. Transport of seed by stingless bees probably accounted for long distance dispersal from C. torelliana, but this mechanism does not appear to supplement gravity-dispersal of seed from the F1. Coupled with other evidence from studies of bee behaviour, controlled pollination in Corymbia sp., and long-term fitness in second generation eucalypt hybrids, we anticipate gene flow via pollen rather than seed will be the greater challenge for managing the risk of introgression of C. torelliana ancestry into native species from the planted F1 hybrid. If large sources of F1 pollen become available to compete with native pollen, gene flow will probably be frequent and hybrids may establish in disturbed conditions and in habitats beyond the ecological range of their native parent. Further study is needed to determine the degree to which outbreeding depression and poor survival inhibits on-going gene flow.

Comparison of the impact of six heat-load management strategies on thermal responses and milk production of feed-pad and pasture fed dairy cows in a subtropical environment

Exposure to hot environments affects milk yield (MY) and milk composition of pasture and feed-pad fed dairy cows in subtropical regions. This study was undertaken during summer to compare MY and physiology of cows exposed to six heat-load management treatments. Seventy-eight Holstein-Friesian cows were blocked by season of calving, parity, milk yield, BW, and milk protein (%) and milk fat (%) measured in 2 weeks prior to the start of the study. Within blocks, cows were randomly allocated to one of the following treatments: open-sided iron roofed day pen adjacent to dairy (CID) + sprinklers (SP); CID only; non-shaded pen adjacent to dairy + SP (NSD + SP); open-sided shade cloth roofed day pen adjacent to dairy (SCD); NSD + sprinkler (sprinkler on for 45 min at 1100 h if mean respiration rate >80 breaths per minute (NSD + WSP)); open-sided shade cloth roofed structure over feed bunk in paddock + 1 km walk to and from the dairy (SCP + WLK). Sprinklers for CID + SP and NSD + SP cycled 2 min on, 12 min off when ambient temperature >26°C. The highest milk yields were in the CID + SP and CID treatments (23.9 L cow−1 day−1), intermediate for NSD + SP, SCD and SCP + WLK (22.4 L cow−1 day−1), and lowest for NSD + WSP (21.3 L cow−1 day−1) (P < 0.05). The highest (P < 0.05) feed intakes occurred in the CID + SP and CID treatments while intake was lowest (P < 0.05) for NSD + WSP and SCP + WLK. Weather data were collected on site at 10-min intervals, and from these, THI was calculated. Nonlinear regression modelling of MY × THI and heat-load management treatment demonstrated that cows in CID + SP showed no decline in MY out to a THI break point value of 83.2, whereas the pooled MY of the other treatments declined when THI >80.7. A combination of iron roof shade plus water sprinkling throughout the day provided the most effective control of heat load.