A hidden Markov model approach for determining vessel activity from vessel monitoring system data.

Many fisheries worldwide have adopted vessel monitoring systems (VMS) for compliance purposes. An added benefit of these systems is that they collect a large amount of data on vessel locations at very fine spatial and temporal scales. This data can provide a wealth of information for stock assessment, research, and management. However, since most VMS implementations record vessel location at set time intervals with no regard to vessel activity, some methodology is required to determine which data records correspond to fishing activity. This paper describes a probabilistic approach, based on hidden Markov models (HMMs), to determine vessel activity. A HMM provides a natural framework for the problem and, by definition, models the intrinsic temporal correlation of the data. The paper describes the general approach that was developed and presents an example of this approach applied to the Queensland trawl fishery off the coast of eastern Australia. Finally, a simulation experiment is presented that compares the misallocation rates of the HMM approach with other approaches.

Statistical methods for analysis of multi-harvest data from perennial pasture variety selection trials

Variety selection in perennial pasture crops involves identifying best varieties from data collected from multiple harvest times in field trials. For accurate selection, the statistical methods for analysing such data need to account for the spatial and temporal correlation typically present. This paper provides an approach for analysing multi-harvest data from variety selection trials in which there may be a large number of harvest times. Methods are presented for modelling the variety by harvest effects while accounting for the spatial and temporal correlation between observations. These methods provide an improvement in model fit compared to separate analyses for each harvest, and provide insight into variety by harvest interactions. The approach is illustrated using two traits from a lucerne variety selection trial. The proposed method provides variety predictions allowing for the natural sources of variation and correlation in multi-harvest data.

Temporal and Spatial Patterns in Recruitment of Three Penaeid Prawns in Moreton Bay, Queensland, Australia

The size at recruitment, temporal and spatial distribution, and abiotic factors influencing abundance of three commercially important species of penaeid prawns in the sublittoral trawl grounds of Moreton Bay (Queensland, Australia) were compared. Metapenaeus bennettae and Penaeus plebejus recruit to the trawl grounds at sizes which are relatively small (14-15 mm carapace length, CL) and below that at which prawns are selected for, and retained, in the fleet's cod-ends. In contrast, Penaeus esculenlus recruit at the relatively large size of 27 mm CL from February to May, well above the size ranges selected for. Recruitment of M. bennettae extends over several months, September-October and February March, and was thus likely to be bi-annual, while the recruitment period of P. plebejus was distinct, peaking in October-November each year. Size classes of M . bennettae were the most spatially stratified of the three species. Catch rates of recruits were negatively correlated with depth for all three species, and were also negatively correlated with salinity for M. bennettae.

Correlation of plasma IGF-I concentrations and growth rate in aquacultured finfish: a tool for assessing the potential of new diets

A recently developed radioimmunoassay (RIA) for measuring insulin-like growth factor (IGF-I) in a variety of fish species was used to investigate the correlation between growth rate and circulating IGF-I concentrations of barramundi (Lates calcarifer), Atlantic salmon (Salmo salar) and Southern Bluefin tuna (Thunnus maccoyii). Plasma IGF-I concentration significantly increased with increasing ration size in barramundi and IGF-I concentration was positively correlated to growth rates obtained in Atlantic salmon (r2=0.67) and barramundi (r2=0.65) when fed a variety of diet formulations. IGF-I was also positively correlated to protein concentration (r2=0.59). This evidence suggested that measuring IGF-I concentration may provide a useful tool for monitoring fish growth rate and also as a method to rapidly assess different aquaculture diets. However, no such correlation was demonstrated in the tuna study probably due to seasonal cooling of sea surface temperature shortly before blood was sampled. Thus, some recommendations for the design and sampling strategy of nutritional trials where IGF-I concentrations are measured are discussed

Temporal expression of G-protein-coupled receptor 54 (GPR54), gonadotropin-releasing hormones (GnRH), and dopamine receptor D2 (drd2) in pubertal female grey mullet, Mugil cephalus.

The G-protein-coupled receptor 54 (muGPR54) cDNA was cloned from the brain of the grey mullet, and its expression level, as well as those of the gonadotropin-releasing hormones (GnRH1, GnRH2, GnRH3) and dopamine receptor D2 (drd2), in the brain, pituitary and ovary of pubertal fish (early, intermediate, advanced) were determined by real-time quantitative RT-PCR (QPCR). The muGPR54 cDNA has an open reading frame of 1140 bp with a predicted 380 amino acid peptide, containing seven putative transmembrane domains and putative N-glycosylation and protein kinase C phosphorylation sites. QPCR results showed that the early stage of puberty in grey mullet is characterized by significantly high levels of expression of GPR54, GnRH and drd2 in the brain relative to the intermediate and advanced stages, except for GnRH1 that increased at the advanced stage of puberty. In the pituitary, drd2 expression declined significantly at the advanced stage relative to levels at the intermediate stage. Ovarian expression of GPR54 significantly increased from the intermediate stage of puberty relative to the early stage while that of GnRH1 acutely increased at the advanced stage of puberty. The ovarian expression of drd2 decreased as puberty progressed, but the changes were not significant. The results suggest the possible role of GPR54 and GnRH in positively regulating pubertal development in grey mullet and the dopaminergic inhibition of reproductive function mediated by drd2.

Spatio-temporal patterns in maturation of the chokka squid (Loligo vulgaris reynaudii) off the coast of South Africa

Knowledge of the temporal and spatial characteristics of chokka squid (Loligo vulgaris reynaudii) biology in South African waters is limited, so the possibility of there being a geographically fragmented stock was examined by investigating the distribution of maturity patterns for the species, covering all known spawning areas and using both historical and recent data. Gonadosomatic indices (GSI) varied between year-round consistency and apparent seasonal peaks in both summer and winter; there was no clear spatial pattern. Monthly percentage maturity provided further evidence for two peak reproductive periods each year, although mature squid were present throughout. Sex ratios demonstrated great variability between different areas and life history stages. Male-biased sex ratios were only apparent on the inshore spawning grounds and ranged between 1.118:1 and 4.267:1. Size at sexual maturity was also seasonal, squid maturing smaller in winter/spring than in summer/autumn. Also, squid in the east matured smaller than squid in the west. Although the results from the present study do not provide conclusive evidence of distinct geographic populations, squid likely spawn over a significantly larger area of the Agulhas Bank than previously estimated, and squid on the west coast of South Africa may return to spawn on the western portion of the Agulhas Bank. It remains likely, however, that the east and west coast populations are a single stock and that migration of juveniles to the west coast and their subsequent return as sub-adults is an integral but non-essential and variable part of the life history.

Temporal variation in arthropod sampling effectiveness: The case for using the beat sheet method in cotton

Predatory insects and spiders are key elements of integrated pest management (IPM) programmes in agricultural crops such as cotton. Management decisions in IPM programmes should to be based on a reliable and efficient method for counting both predators and pests. Knowledge of the temporal constraints that influence sampling is required because arthropod abundance estimates are likely to vary over a growing season and within a day. Few studies have adequately quantified this effect using the beat sheet, a potentially important sampling method. We compared the commonly used methods of suction and visual sampling to the beat sheet, with reference to an absolute cage clamp method for determining the abundance of various arthropod taxa over 5 weeks. There were significantly more entomophagous arthropods recorded using the beat sheet and cage clamp methods than by using suction or visual sampling, and these differences were more pronounced as the plants grew. In a second trial, relative estimates of entomophagous and phytophagous arthropod abundance were made using beat sheet samples collected over a day. Beat sheet estimates of the abundance of only eight of the 43 taxa examined were found to vary significantly over a day. Beat sheet sampling is recommended in further studies of arthropod abundance in cotton, but researchers and pest management advisors should bear in mind the time of season and time of day effects.

Spatial and temporal variation and effects of changes in management in discard rates from the commercial reef line fishery of the Great Barrier Reef, Australia

Discarding in commercially exploited fisheries has received considerable attention in the last decade, though only more recently in Australia. The Reef Line fishery (RLF) of the Great Barrier Reef (GBR) in Australia is a large-scale multi-sector, multi-species, highly regulated hook and line fishery with the potential for high levels of discarding. We used a range of data sources to estimate discard rates and discard quantities for the two main target groups of the RLF, the coral trout, Plectropomus spp, and the red throat emperor, Lethrinus miniatus, and investigated possible effects on discarding of recent changes in management of the fishery. Fleet-wide estimates of total annual quantities discarded from 1989 to 2003 were 292-622 t and 33-95 t for coral trout and red throat emperor, respectively. Hypothetical scenarios of high-grading after the introduction of a total allowable commercial catch for coral trout resulted in increases in discard quantities up to 3895 t, while no high-grading still meant 421 t were discarded. Increasing the minimum size limit of red throat emperor from 35 to 38 cm also increased discards to an estimated 103 t. We provide spatially and temporally explicit estimates of discarding for the two most important species in the GBR RLF of Australia to demonstrate the importance of accounting for regional variation in quantification of discarding. Effects of management changes on discarding are also highlighted. This study provides a template for exploring discarding levels for other species in the RLF and elsewhere.

Characterization of disease resistance gene candidates of the nucleotide binding site (NBS) type from banana and correlation of a transcriptional polymorphism with resistance to Fusarium oxysporum f.sp. cubense race 4

Most plant disease resistance (R) genes encode proteins with a nucleotide binding site and leucine-rich repeat structure (NBS-LRR). In this study, degenerate primers were used to amplify genomic NBS-type sequences from wild banana (Musa acuminata ssp. malaccensis) plants resistant to the fungal pathogen Fusarium oxysporum formae specialis (f. sp.) cubense (FOC) race 4. Five different classes of NBS-type sequences were identified and designated as resistance gene candidates (RGCs). The deduced amino acid sequences of the RGCs revealed the presence of motifs characteristic of the majority of known plant NBS-LRR resistance genes. Structural and phylogenetic analyses grouped the banana RGCs within the non-TIR (homology to Toll/interleukin-1 receptors) subclass of NBS sequences. Southern hybridization showed that each banana RGC is present in low copy number. The expression of the RGCs was assessed by RT-PCR in leaf and root tissues of plants resistant or susceptible to FOC race 4. RGC1, 3 and 5 showed a constitutive expression profile in both resistant and susceptible plants whereas no expression was detected for RGC4. Interestingly, RGC2 expression was found to be associated only to FOC race 4 resistant lines. This finding could assist in the identification of a FOC race 4 resistance gene.

Spatial and temporal patterns of feral pig diggings in rainforests of north Queensland

Feral pigs (Sus scrofa) are believed to have a severe negative impact on the ecological values of tropical rainforests in north Queensland, Australia. Most perceptions of the environmental impacts of feral pigs focus on their disturbance of the soil or surface material (diggings). Spatial and temporal patterns of feral pig diggings were identified in this study: most diggings occurred in the early dry season and predominantly in moist soil (swamp and creek) microhabitats, with only minimal pig diggings found elsewhere through the general forest floor. The overall mean daily pig diggings were 0.09% of the rainforest floor. Most diggings occurred 3-4 months after the month of maximum rainfall. Most pig diggings were recorded in highland swamps, with over 80% of the swamp areas dug by pigs at some time during the 18-month study period. These results suggest that management of feral pig impacts should focus on protecting swamp and creek microhabitats in the rainforest, which are preferred by pigs for digging and which have a high environmental significance.

Temporal patterns in incidence and abundance of Aconophora compressa (Hemiptera: Membracidae), a biological control agent for Lantana camara, on target and nontarget plants

The membracid Aconophora compressa Walker, a biological control agent released in 1995 to control Lantana camara (Verbenaceae) in Australia, has since been collected on several nontarget plant species. Our survey suggests that sustained populations of A. compressa are found only on the introduced nontarget ornamental Citharexylum spinosum (Verbenaceae) and the target weed L. camara. It is found on other nontarget plant species only when populations on C. spinosum and L. camara are high, suggesting that the presence of populations on nontarget species may be a spill-over effect. Some of the incidence and abundance on nontarget plants could have been anticipated from host specificity studies done on this agent before release, whereas others could not. This raises important issues about predicting risks posed by weed biological control agents and the need for long-term postintroduction monitoring on nontarget species.

Using harvest statistics to monitor temporal variation in kangaroo density and harvest rate

Management of the commercial harvest of kangaroos relies on quotas set annually as a proportion of regular estimates of population size. Surveys to generate these estimates are expensive and, in the larger states, logistically difficult; a cheaper alternative is desirable. Rainfall is a disappointingly poor predictor of kangaroo rate of increase in many areas, but harvest statistics (sex ratio, carcass weight, skin size and animals shot per unit time) potentially offer cost-effective indirect monitoring of population abundance (and therefore trend) and status (i.e. under-or overharvest). Furthermore, because harvest data are collected continuously and throughout the harvested areas, they offer the promise of more intensive and more representative coverage of harvest areas than aerial surveys do. To be useful, harvest statistics would need to have a close and known relationship with either population size or harvest rate. We assessed this using longterm (11-22 years) data for three kangaroo species (Macropus rufus, M. giganteus and M. fuliginosus) and common wallaroos (M. robustus) across South Australia, New South Wales and Queensland. Regional variation in kangaroo body size, population composition, shooter efficiency and selectivity required separate analyses in different regions. Two approaches were taken. First, monthly harvest statistics were modelled as a function of a number of explanatory variables, including kangaroo density, harvest rate and rainfall. Second, density and harvest rate were modelled as a function of harvest statistics. Both approaches incorporated a correlated error structure. Many but not all regions had relationships with sufficient precision to be useful for indirect monitoring. However, there was no single relationship that could be applied across an entire state or across species. Combined with rainfall-driven population models and applied at a regional level, these relationships could be used to reduce the frequency of aerial surveys without compromising decisions about harvest management.

Reassessing spatial and temporal dynamics of kangaroo populations

Long-running datasets from aerial surveys of kangaroos (Macropus giganteus, Macropus [uliginosus, Macropus robustus and Macropus rufus) across Queensland, New South Wales and South Australia have been analysed, seeking better predictors of rates of increase which would allow aerial surveys to be undertaken less frequently than annually. Early models of changes in kangaroo numbers in response to rainfall had shown great promise, but much variability. We used normalised difference vegetation index (NDVI) instead, reasoning that changes in pasture condition would provide a better predictor than rainfall. However, except at a fine scale, NDVI proved no better; although two linked periods of rainfall proved useful predictors of rates of increase, this was only in some areas for some species. The good correlations reported in earlier studies were a consequence of data dominated by large droughtinduced adult mortality, whereas over a longer time frame and where changes between years are less dramatic, juvenile survival has the strongest influence on dynamics. Further, harvesting, density dependence and competition with domestic stock are additional and important influences and it is now clear that kangaroo movement has a greater influence on population dynamics than had been assumed. Accordingly, previous conclusions about kangaroo populations as simple systems driven by rainfall need to be reassessed. Examination of this large dataset has permitted descriptions of shifts in distribution of three species across eastern Australia, changes in dispersion in response to rainfall, and an evaluation of using harvest statistics as an index of density and harvest rate. These results have been combined into a risk assessment and decision theory framework to identify optimal monitoring strategies.

Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection.

In order to investigate the effect of long term recurrent selection on the pattern of gene diversity, thirty randomly-selected individuals from the progenitors (p) and four selection cycles (C0, C3, C6 and C11) were sampled for DNA analysis from the tropical maize (Zea mays L.) breeding populations, Atherton 1 (AT1) and Atherton 2 (AT2). Fifteen polymorphic Simple Sequence Repeat markers amplified a total of 284 and 257 alleles in AT1 and AT2 populations, respectively. Reductions in the number of alleles were observed at advanced selection cycles. About 11 and 12% of the alleles in AT1 and AT2 populations respectively, were near to fixation. However, a higher number of alleles (37% in AT1 and 33% in AT2) were close to extinction. Fisher's exact test and analysis of molecular variance (AMOVA) showed significant population differentiations. Gene diversity estimates and AMOVA revealed increased genetic differentiations at the expense of loss of heterozygosity. Population differentiations were mainly due to fixation of complementary alleles at a locus in the two breeding populations. The estimates of effective population at an advanced selection cycle were close to the population size predicted by the breeding method.

Water and thermal regimes for field pea in Australia and their implications for breeding.

There is a large gap between the refined approaches to characterise genotypes and the common use of location and season as a coarse surrogate for environmental characterisation of breeding trials. As a framework for breeding, the aim of this paper is quantifying the spatial and temporal patterns of thermal and water stress for field pea in Australia. We compiled a dataset for yield of the cv. Kaspa measured in 185 environments, and investigated the associations between yield and seasonal patterns of actual temperature and modelled water stress. Correlations between yield and temperature indicated two distinct stages. In the first stage, during crop establishment and canopy expansion before flowering, yield was positively associated with minimum temperature. Mean minimum temperature below similar to 7 degrees C suggests that crops were under suboptimal temperature for both canopy expansion and radiation-use efficiency during a significant part of this early growth period. In the second stage, during critical reproductive phases, grain yield was negatively associated with maximum temperature over 25 degrees C. Correlations between yield and modelled water supply/demand ratio showed a consistent pattern with three phases: no correlation at early stages of the growth cycle, a progressive increase in the association that peaked as the crop approached the flowering window, and a progressive decline at later reproductive stages. Using long-term weather records (1957-2010) and modelled water stress for 104 locations, we identified three major patterns of water deficit nation wide. Environment type 1 (ET1) represents the most favourable condition, with no stress during most of the pre-flowering phase and gradual development of mild stress after flowering. Type 2 is characterised by increasing water deficit between 400 degree-days before flowering and 200 degree-days after flowering and rainfall that relieves stress late in the season. Type 3 represents the more stressful condition with increasing water deficit between 400 degree-days before flowering and maturity. Across Australia, the frequency of occurrence was 24% for ET1, 32% for ET2 and 43% for ET3, highlighting the dominance of the most stressful condition. Actual yield averaged 2.2 t/ha for ET1, 1.9 t/ha for ET2 and 1.4 t/ha for ET3, and the frequency of each pattern varied substantially among locations. Shifting from a nominal (i.e. location and season) to a quantitative (i.e. stress type) characterisation of environments could help improving breeding efficiency of field pea in Australia.