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.

Metabolite mobilization in the stem galls of Parthenium hysterophorus induced by Epiblema strenuana inferred from the signatures of isotopic carbon and nitrogen and concentrations of total non-structural carbohydrates

Parthenium hysterophorus L. (Asteraceae) is a weed of national significance in Australia. Among the several arthropod agents introduced into Australia to control populations of P. hysterophorus biologically, Epiblema strenuana Walker (Lepidoptera: Tortricidae) is the most widespread and abundant agent. By intercepting the normal transport mechanisms of P. hysterophorus, the larvae of E. strenuana drain nutrients, other metabolic products, and energy, and place the host plant under intense metabolic stress. In this study, determinations of total non-structural carbohydrates (TNC) levels and carbon and nitrogen isotope ratios of fixed products in different parts of the plant tissue, including the gall, have been made to establish the function of gall as a sink for the nutrients. Values of δ13C and δ15N in galls were significantly different than those in proximal and distal stems, whereas the TNC levels were insignificant, when measured in the total population of P. hysterophorus, regardless of plant age. However, carbon, nitrogen, and TNC signatures presented significant results, when assayed in different developmental stages of P. hysterophorus. Carbon isotope ratios in galls were consistently more negative than those from the compared plant organs. Nitrogen isotope ratios in galls, on the contrary, were either similar to or less negative than the compared plant organs, especially within a single host-plant stage population (i.e., either rosette, preflowering, or flowering stage). TNC levels varied within compared plant populations. The stem distal to the gall functioned more efficiently as a nodal channel than the stem proximal to the gall, especially in the translocation of nitrogenous nutrients. Our findings indicate that the gall induced by E. strenuana functions as a sink for the assayed nutrients, although some variations have been observed in the patterns of nutrient mobilization. By creating a sink for the nutrients in the gall, E. strenuana is able to place the overall plant metabolism under stress, and this ability indicates E. strenuana has the necessary potential for use as a biological-control agent.

Effect of the Queensland Shark Control Program on non target species : whale, dugong, turtle and dolphin : a review

The Queensland Shark Control Program (QSCP) aims to protect swimmers at ten beach areas on the east coast of Queensland between Cairns (17°S) and the Gold coast (28°S). Since its inception in 1962 it has deployed shark nets and baited drumlines in a `mixed gear strategy' that adapts the type of gear to the characteristics of a site (e .g . extreme tidal range, high energy wave action, or proximity of turtle breeding areas) . The policy has provided swimmer protection, and the incidental capture of non-target species has been lower than that resulting from deployment of nets alone (Dudley 1997; Gribble et al. 1998b). The QSCP is the only major public-safety shark-control program to routinely use mixed gear. Both the New South Wales (Holt 1998) and KwaZulu-Natal (Dudley 1998) programs use nets exclusively, although the KwaZulu-Natal program has recently tested drumlines on an experimental basis (Dudley 1998; Dudley, personal communication).

Effect of irrigation amounts applied with subsurface drip irrigation on corn evapotranspiration, yield, water use efficiency, and dry matter production in a semiarid climate

Quantifying the local crop response to irrigation is important for establishing adequate irrigation management strategies. This study evaluated the effect of irrigation applied with subsurface drip irrigation on field corn (Zea mays L.) evapotranspiration (ETc), yield, water use efficiencies (WUE = yield/ETc, and IWUE = yield/irrigation), and dry matter production in the semiarid climate of west central Nebraska. Eight treatments were imposed with irrigation amounts ranging from 53 to 356 mm in 2005 and from 22 to 226 mm in 2006. A soil water balance approach (based on FAO-56) was used to estimate daily soil water and ETc. Treatments resulted in seasonal ETc of 580-663 mm and 466-656 mm in 2005 and 2006, respectively. Yields among treatments differed by as much as 22% in 2005 and 52% in 2006. In both seasons, irrigation significantly affected yields, which increased with irrigation up to a point where irrigation became excessive. Distinct relationships were obtained each season. Yields increased linearly with seasonal ETc (R 2 = 0.89) and ETc/ETp (R 2 = 0.87) (ETp = ETc with no water stress). The yield response factor (ky), which indicates the relative reduction in yield to relative reduction in ETc, averaged 1.58 over the two seasons. WUE increased non-linearly with seasonal ETc and with yield. WUE was more sensitive to irrigation during the drier 2006 season, compared with 2005. Both seasons, IWUE decreased sharply with irrigation. Irrigation significantly affected dry matter production and partitioning into the different plant components (grain, cob, and stover). On average, the grain accounted for the majority of the above-ground plant dry mass (≈59%), followed by the stover (≈33%) and the cob (≈8%). The dry mass of the plant and that of each plant component tended to increase with seasonal ETc. The good relationships obtained in the study between crop performance indicators and seasonal ETc demonstrate that accurate estimates of ETc on a daily and seasonal basis can be valuable for making tactical in-season irrigation management decisions and for strategic irrigation planning and management.

Identifying major contributing sources to odour annoyance using a non-specific gas sensor array

Identification of major contributors to odour annoyance in areas with multiple emission sources is necessary to address and resolve odour disputes. In an effort to develop an appropriate tool for this task, odour samples were collected on-site at a piggery and an abattoir (the major odour sources in the area) and at surrounding off-site areas, then analysed using a commercial non-specific chemical sensor array to develop an odour fingerprint database. The developed odour fingerprint database was analysed using two pattern recognition algorithms including a partial least squares-discriminant analysis (PLS-DA) and a Kohonen self-organising map (KSOM). The KSOM model could identify odour samples sourced from the piggery shed 15, piggery pond 8, piggery pond 9, abattoir, motel and others with mean percentage values of 77.5, 65.0, 90.2, 75.7, 44.8 and 64.6%, respectively.

Predicting the response of wheat (Triticum aestivum L.) to liquid and granular phosphorus fertilisers in Australian soils

Liquid forms of phosphorus (P) have been shown to be more effective than granular P for promoting cereal growth in alkaline soils with high levels of free calcium carbonate on Eyre Peninsula, South Australia. However, the advantage of liquid over granular P forms of fertiliser has not been fully investigated across the wide range of soils used for grain production in Australia. A glasshouse pot experiment tested if liquid P fertilisers were more effective for growing spring wheat (Triticum aestivum L.) than granular P (monoammonium phosphate) in 28 soils from all over Australia with soil pH (H2O) ranging from 5.2 to 8.9. Application of liquid P resulted in greater shoot biomass, as measured after 4 weeks' growth (mid to late tillering, Feeks growth stage 2-3), than granular P in 3 of the acidic to neutral soils and in 3 alkaline soils. Shoot dry matter responses of spring wheat to applied liquid or granular P were related to soil properties to determine if any of the properties predicted superior yield responses to liquid P. The calcium carbonate content of soil was the only soil property that significantly contributed to predicting when liquid P was more effective than granular P. Five soil P test procedures (Bray, Colwell, resin, isotopically exchangeable P, and diffusive gradients in thin films (DGT)) were assessed to determine their ability to measure soil test P on subsamples of soil collected before the experiment started. These soil test values were then related to the dry matter shoot yields to assess their ability to predict wheat yield responses to P applied as liquid or granular P. All 5 soil test procedures provided a reasonable prediction of dry matter responses to applied P as either liquid or granular P, with the resin P test having a slightly greater predictive capacity on the range of soils tested. The findings of this investigation suggest that liquid P fertilisers do have some potential applications in non-calcareous soils and confirm current recommendations for use of liquid P fertiliser to grow cereal crops in highly calcareous soils. Soil P testing procedures require local calibration for response to the P source that is going to be used to amend P deficiency.

The limit to wheat water-use efficiency in eastern Australia. I. Gradients in the radiation environment and atmospheric demand

In the wheatbelt of eastern Australia, rainfall shifts from winter dominated in the south (South Australia, Victoria) to summer dominated in the north (northern New South Wales, southern Queensland). The seasonality of rainfall, together with frost risk, drives the choice of cultivar and sowing date, resulting in a flowering time between October in the south and August in the north. In eastern Australia, crops are therefore exposed to contrasting climatic conditions during the critical period around flowering, which may affect yield potential, and the efficiency in the use of water (WUE) and radiation (RUE). In this work we analysed empirical and simulated data, to identify key climatic drivers of potential water- and radiation-use efficiency, derive a simple climatic index of environmental potentiality, and provide an example of how a simple climatic index could be used to quantify the spatial and temporal variability in resource-use efficiency and potential yield in eastern Australia. Around anthesis, from Horsham to Emerald, median vapour pressure deficit (VPD) increased from 0.92 to 1.28 kPa, average temperature increased from 12.9 to 15.2°C, and the fraction of diffuse radiation (FDR) decreased from 0.61 to 0.41. These spatial gradients in climatic drivers accounted for significant gradients in modelled efficiencies: median transpiration WUE (WUEB/T) increased southwards at a rate of 2.6% per degree latitude and median RUE increased southwards at a rate of 1.1% per degree latitude. Modelled and empirical data confirmed previously established relationships between WUEB/T and VPD, and between RUE and photosynthetically active radiation (PAR) and FDR. Our analysis also revealed a non-causal inverse relationship between VPD and radiation-use efficiency, and a previously unnoticed causal positive relationship between FDR and water-use efficiency. Grain yield (range 1-7 t/ha) measured in field experiments across South Australia, New South Wales, and Queensland (n = 55) was unrelated to the photothermal quotient (Pq = PAR/T) around anthesis, but was significantly associated (r2 = 0.41, P < 0.0001) with newly developed climatic index: a normalised photothermal quotient (NPq = Pq . FDR/VPD). This highlights the importance of diffuse radiation and vapour pressure deficit as sources of variation in yield in eastern Australia. Specific experiments designed to uncouple VPD and FDR and more mechanistic crop models might be required to further disentangle the relationships between efficiencies and climate drivers.

A systems evaluation of high-input management using fortified molasses for beef production in Australia's dry tropics

The potential of beef producers to profitably produce 500-kg steers at 2.5 years of age in northern Australia's dry tropics to meet specifications of high-value markets, using a high-input management (HIM) system was examined. HIM included targeted high levels of fortified molasses supplementation, short seasonal mating and the use of growth promotants. Using herds of 300-400 females plus steer progeny at three sites, HIM was compared at a business level to prevailing best-practice, strategic low-input management (SLIM) in which there is a relatively low usage of energy concentrates to supplement pasture intake. The data presented for each breeding-age cohort within management system at each site includes: annual pregnancy rates (range: 14-99%), time of conception, mortalities (range: 0-10%), progeny losses between confirmed pregnancy and weaning (range: 0-29%), and weaning rates (range: 14-92%) over the 2-year observation. Annual changes in weight and relative net worth were calculated for all breeding and non-breeding cohorts. Reasons for outcomes are discussed. Compared with SLIM herds, both weaning weights and annual growth were >= 30 kg higher, enabling 86-100% of HIM steers to exceed 500 kg at 2.5 years of age. Very few contemporary SLIM steers reached this target. HIM was most profitably applied to steers. Where HIM was able to achieve high pregnancy rates in yearlings, its application was recommended in females. Well managed, appropriate HIM systems increased profits by around $15/adult equivalent at prevailing beef and supplement prices. However, a 20% supplement price rise without a commensurate increase in values for young slaughter steers would generally eliminate this advantage. This study demonstrated the complexity of pro. table application of research outcomes to commercial business, even when component research suggests that specific strategies may increase growth and reproductive efficiency and/or be more pro. table. Because of the higher level of management required, higher costs and returns, and higher susceptibility to market changes and disease, HIM systems should only be applied after SLIM systems are well developed. To increase profitability, any strategy must ultimately either increase steer growth and sale values and/or enable a shift to high pregnancy rates in yearling heifers.

Non-specific conducting polymer-based array capable of monitoring odour emissions from a biofiltration system in a piggery building

A commercial non-specific gas sensor array system was evaluated in terms of its capability to monitor the odour abatement performance of a biofiltration system developed for treating emissions from a commercial piggery building. The biofiltration system was a modular system comprising an inlet ducting system, humidifier and closed-bed biofilter. It also included a gravimetric moisture monitoring and water application system for precise control of moisture content of an organic woodchip medium. Principal component analysis (PCA) of the sensor array measurements indicated that the biofilter outlet air was significantly different to both inlet air of the system and post-humidifier air. Data pre-processing techniques including normalising and outlier handling were applied to improve the odour discrimination performance of the non-specific gas sensor array. To develop an odour quantification model using the sensor array responses of the non-specific sensor array, PCA regression, artificial neural network (ANN) and partial least squares (PLS) modelling techniques were applied. The correlation coefficient (r(2)) values of the PCA, ANN, and PLS models were 0.44, 0.62 and 0.79, respectively.

Development of near infrared analysis of faeces to estimate non-grass proportions in diets selected by cattle grazing tropical pastures

Grass (monocots) and non-grass (dicots) proportions in ruminant diets are important nutritionally because the non-grasses are usually higher in nutritive value, particularly protein, than the grasses, especially in tropical pastures. For ruminants grazing tropical pastures where the grasses are C-4 species and most non-grasses are C-3 species, the ratio of C-13/C-12 in diet and faeces, measured as delta C-13 parts per thousand, is proportional to dietary non-grass%. This paper describes the development of a faecal near infrared (NIR) spectroscopy calibration equation for predicting faecal delta C-13 from which dietary grass and non-grass proportions can be calculated. Calibration development used cattle faeces derived from diets containing only C-3 non-grass and C-4 grass components, and a series of expansion and validation steps was employed to develop robustness and predictive reliability. The final calibration equation contained 1637 samples and faecal delta C-13 range (parts per thousand) of [12.27]-[27.65]. Calibration statistics were: standard error of calibration (SEC) of 0.78, standard error of cross-validation (SECV) of 0.80, standard deviation (SD) of reference values of 3.11 and R-2 of 0.94. Validation statistics for the final calibration equation applied to 60 samples were: standard error of prediction (SEP) of 0.87, bias of -0.15, R-2 of 0.92 and RPD of 3.16. The calibration equation was also tested on faeces from diets containing C-4 non-grass species or temperate C-3 grass species. Faecal delta C-13 predictions indicated that the spectral basis of the calibration was not related to C-13/C-12 ratios per se but to consistent differences between grasses and non-grasses in chemical composition and that the differences were modified by photosynthetic pathway. Thus, although the calibration equation could not be used to make valid faecal delta C-13 predictions when the diet contained either C-3 grass or C-4 non-grass, it could be used to make useful estimates of dietary non-grass proportions. It could also be ut :sed to make useful estimates of non-grass in mixed C-3 grass/non-grass diets by applying a modified formula to calculate non-grass from predicted faecal delta C-13. The development of a robust faecal-NIR calibration equation for estimating non-grass proportions in the diets of grazing cattle demonstrated a novel and useful application of NIR spectroscopy in agriculture.

Weeds of Australian rangelands

Despite recognition that non-native plant species represent a substantial risk to natural systems, there is currently no compilation of weeds that impact on the biodiversity of the rangelands within Australia. Using published and expert knowledge, this paper presents a list of 622 non-native naturalised species known to occur within the rangelands. Of these, 160 species (26%) are considered a current threat to rangeland biodiversity. Most of these plant species have been deliberately introduced for forage or other commercial use (e.g. nursery trade). Among growth forms, shrubs and perennial grasses comprise over 50% of species that pose the greatest risk to rangeland biodiversity. We identify regions within the rangelands containing both high biodiversity values and a high proportion of weeds and recommend these areas as priorities for weed management. Finally, we examine the resources available for weed detection and identification since detecting weeds in the early stages of invasion is the most cost effective method of reducing further impact.

Queensland Coastal Wetland Resources: Sand Bay to Keppel Bay.

Protection of coastal wetland environments is an important prerequisite to effective and sustainable inshore fisheries management and conservation of habitats for use by future generations. Mangroves, saltmarshes, seagrasses and non vegetated habitats directly support local and regional inshore and offshore fisheries through the provision of food, shelter, breeding and nursery grounds. As such, these wetland environments have significant economic value as well as their intrinsic aesthetic and ecological values. This report summarises the results of the mapping undertaken in the Central Queensland Coast from Sand Bay to Keppel Bay (hereafter referred to as the Study Area). The study was undertaken in order to: 1. document and map the coastal wetland communities along the Queensland coastline from Sand Bay (20.93°S, 149.04°E) to Keppel Bay (23.65°S, 151.07°E); 2. document levels of existing disturbance to and protection of the wetlands; 3. examine existing recreational and commercial fisheries in the region; and 4. evaluate the conservation values of the areas investigated from the viewpoint of fisheries productivity and as habitat for important and/or threatened species. Dataset URL Link: Queensland Coastal Wetlands Resources Mapping data. [Dataset]

Characterising insect plant host relationships facilitates understanding multiple host use

Weed biocontrol relies on host specificity testing, usually carried out under quarantine conditions to predict the future host range of candidate control agents. The predictive power of host testing can be scrutinised directly with Aconophora compressa, previously released against the weed Lantana camara L. (lantana) because its ecology in its new range (Australia) is known and includes the unanticipated use of several host species. Glasshouse based predictions of field host use from experiments designed a posteriori can therefore be compared against known field host use. Adult survival, reproductive output and egg maturation were quantified. Adult survival did not differ statistically across the four verbenaceous hosts used in Australia. Oviposition was significantly highest on fiddlewood (Citharexylum spinosum L.), followed by lantana, on which oviposition was significantly higher than on two varieties of Duranta erecta (‘‘geisha girl’’ and ‘‘Sheena’s gold’’; all Verbenaceae). Oviposition rates across Duranta varieties were not significantly different from each other but were significantly higher than on the two non-verbenaceous hosts (Jacaranda mimosifolia D. Don: Bignoneaceae (jacaranda) and Myoporum acuminatum R. Br.: Myoporaceae (Myoporum)). Production of adult A. compressa was modelled across the hosts tested. The only major discrepancy between model output and their relative abundance across hosts in the field was that densities on lantana in the field were much lower than predicted by the model. The adults may, therefore, not locate lantana under field conditions and/or adults may find lantana but leave after laying relatively few eggs. Fiddlewood is the only primary host plant of A. compressa in Australia, whereas lantana and the others are used secondarily or incidentally. The distinction between primary, secondary and incidental hosts of a herbivore species helps to predict the intensity and regularity of host use by that herbivore. Populations of the primary host plants of a released biological control agent are most likely to be consistently impacted by the herbivore, whereas secondary and incidental host plant species are unlikely to be impacted consistently. As a consequence, potential biocontrol agents should be released only against hosts to which they have been shown to be primarily adapted.

Leaf trait co-ordination in relation to construction cost, carbon gain and resource-use efficiency in exotic invasive and native woody vine species

Background and Aims: Success of invasive plant species is thought to be linked with their higher leaf carbon fixation strategy, enabling them to capture and utilize resources better than native species, and thus pre-empt and maintain space. However, these traits are not well-defined for invasive woody vines. Methods: In a glass house setting, experiments were conducted to examine how leaf carbon gain strategies differ between non-indigenous invasive and native woody vines of south-eastern Australia, by investigating their biomass gain, leaf structural, nutrient and physiological traits under changing light and moisture regimes. Key Results: Leaf construction cost (CC), calorific value and carbon : nitrogen (C : N) ratio were lower in the invasive group, while ash content, N, maximum photosynthesis, light-use efficiency, photosynthetic energyuse efficiency (PEUE) and specific leaf area (SLA) were higher in this group relative to the native group. Trait plasticity, relative growth rate (RGR), photosynthetic nitrogen-use efficiency and water-use efficiency did not differ significantly between the groups. However, across light resource, regression analyses indicated that at a common (same) leaf CC and PEUE, a higher biomass RGR resulted for the invasive group; also at a common SLA, a lower CC but higher N resulted for the invasive group. Overall, trait co-ordination (using pair-wise correlation analyses) was better in the invasive group. Ordination using 16 leaf traits indicated that the major axis of invasive-native dichotomy is primarily driven by SLA and CC (including its components and/or derivative of PEUE) and was significantly linked with RGR. Conclusions: These results demonstrated that while not all measures of leaf resource traits may differ between the two groups, the higher level of trait correlation and higher revenue returned (RGR) per unit of major resource need (CC) and use (PEUE) in the invasive group is in line with their rapid spread where introduced.

Net radiation dynamics: performance of 20 daily net radiation models as related to model structure and intricacy in two climates.

We compared daily net radiation (Rn) estimates from 19 methods with the ASCE-EWRI Rn estimates in two climates: Clay Center, Nebraska (sub-humid) and Davis, California (semi-arid) for the calendar year. The performances of all 20 methods, including the ASCE-EWRI Rn method, were then evaluated against Rn data measured over a non-stressed maize canopy during two growing seasons in 2005 and 2006 at Clay Center. Methods differ in terms of inputs, structure, and equation intricacy. Most methods differ in estimating the cloudiness factor, emissivity (e), and calculating net longwave radiation (Rnl). All methods use albedo (a) of 0.23 for a reference grass/alfalfa surface. When comparing the performance of all 20 Rn methods with measured Rn, we hypothesized that the a values for grass/alfalfa and non-stressed maize canopy were similar enough to only cause minor differences in Rn and grass- and alfalfa-reference evapotranspiration (ETo and ETr) estimates. The measured seasonal average a for the maize canopy was 0.19 in both years. Using a = 0.19 instead of a = 0.23 resulted in 6% overestimation of Rn. Using a = 0.19 instead of a = 0.23 for ETo and ETr estimations, the 6% difference in Rn translated to only 4% and 3% differences in ETo and ETr, respectively, supporting the validity of our hypothesis. Most methods had good correlations with the ASCE-EWRI Rn (r2 > 0.95). The root mean square difference (RMSD) was less than 2 MJ m-2 d-1 between 12 methods and the ASCE-EWRI Rn at Clay Center and between 14 methods and the ASCE-EWRI Rn at Davis. The performance of some methods showed variations between the two climates. In general, r2 values were higher for the semi-arid climate than for the sub-humid climate. Methods that use dynamic e as a function of mean air temperature performed better in both climates than those that calculate e using actual vapor pressure. The ASCE-EWRI-estimated Rn values had one of the best agreements with the measured Rn (r2 = 0.93, RMSD = 1.44 MJ m-2 d-1), and estimates were within 7% of the measured Rn. The Rn estimates from six methods, including the ASCE-EWRI, were not significantly different from measured Rn. Most methods underestimated measured Rn by 6% to 23%. Some of the differences between measured and estimated Rn were attributed to the poor estimation of Rnl. We conducted sensitivity analyses to evaluate the effect of Rnl on Rn, ETo, and ETr. The Rnl effect on Rn was linear and strong, but its effect on ETo and ETr was subsidiary. Results suggest that the Rn data measured over green vegetation (e.g., irrigated maize canopy) can be an alternative Rn data source for ET estimations when measured Rn data over the reference surface are not available. In the absence of measured Rn, another alternative would be using one of the Rn models that we analyzed when all the input variables are not available to solve the ASCE-EWRI Rn equation. Our results can be used to provide practical information on which method to select based on data availability for reliable estimates of daily Rn in climates similar to Clay Center and Davis.