Effects of trawling on sessile megabenthos in the Great Barrier Reef and evaluation of the efficacy of management strategies

A series of related research studies over 15 years assessed the effects of prawn trawling on sessile megabenthos in the Great Barrier Reef, to support management for sustainable use in the World Heritage Area. These large-scale studies estimated impacts on benthos (particularly removal rates per trawl pass), monitored subsequent recovery rates, measured natural dynamics of tagged megabenthos, mapped the regional distribution of seabed habitats and benthic species, and integrated these results in a dynamic modelling framework together with spatio-temporal fishery effort data and simulated management. Typical impact rates were between 5 and 25% per trawl, recovery times ranged from several years to several decades, and most sessile megabenthos were naturally distributed in areas where little or no trawling occurred and so had low exposure to trawling. The model simulated trawl impact and recovery on the mapped species distributions, and estimated the regional scale cumulative changes due to trawling as a time series of status for megabenthos species. The regional status of these taxa at time of greatest depletion ranged from ∼77% relative to pre-trawl abundance for the worst case species, having slow recovery with moderate exposure to trawling, to ∼97% for the least affected taxon. The model also evaluated the expected outcomes for sessile megabenthos in response to major management interventions implemented between 1999 and 2006, including closures, effort reductions, and protected areas. As a result of these interventions, all taxa were predicted to recover (by 2-14% at 2025); the most affected species having relatively greater recovery. Effort reductions made the biggest positive contributions to benthos status for all taxa, with closures making smaller contributions for some taxa. The results demonstrated that management actions have arrested and reversed previous unsustainable trends for all taxa assessed, and have led to a prawn trawl fishery with improved environmental sustainability. © 2015 International Council for the Exploration of the Sea 2015. All rights reserved. For Permissions, please email: journals.permissions@oup.com.

FAM49 - A novel regulator of the protrusive behaviour and motility of cells

Most eukaryotic cell motility relies on plasma membrane protrusions, which depend on the actin cytoskeleton and its tight regulation. The SCAR/WAVE complex, a pentameric assembly comprising SCAR/WAVE, Nap1, CYFIP/Pir121, Abi and HSPC300, is a key driver of actin-based protrusions such as pseudopods. SCAR/WAVE is thought to activate the Arp2/3 complex, a crucial actin nucleator, after being itself activated by upstream signals such as active Rac1. Despite recent progress on the study of the SCAR/WAVE complex, its regulation is still incompletely understood, with Nap1’s role being particularly enigmatic. Upon screening for potential Nap1 binding partners in the social amoeba Dictyostelium discoideum – a well established model organism in the study of the actin cytoskeleton and cell motility – we found FAM49, a ~36 kDa protein of unknown function which is highly conserved in Metazoa (animals) and evolutionarily closer species such as D. discoideum. Interestingly, D. discoideum’s FAM49 and its homologs contain a DUF1394 domain, which is also predicted in CYFIP/Pir121 proteins and most likely involved in their direct binding to active Rac1, which in turn contributes to SCAR/WAVE’s activation. FAM49’s unknown role, apparent high degree of conservation and potential connections to SCAR/WAVE and Rac1 persuaded us to start investigating its function and biological relevance in D. discoideum, leading to the work presented in this thesis. Several pieces of our data collectively support a function for FAM49 in modulating the protrusive behaviour, and ultimately motility, of D. discoideum cells, as well as a regulatory link between FAM49 and Rac1. FAM49’s involvement in protrusion regulation was first hinted at by our observation that GFP-tagged FAM49 is enriched in pseudopods. The possibility of a link with Rac1 was then strengthened by two additional observations: first, pseudopodial GFP-FAM49 is substantially co-enriched with active Rac, both showing fairly comparable spatio-temporal accumulation dynamics; second, when dominant-active (G12V) Rac1 is expressed in cells, it triggers the recruitment and persistent accumulation of GFP-FAM49 at the plasma membrane, where both become highly co-enriched. We subsequently determined that fam49 KO cells differ from wild-type cells in the way they protrude and move, as assessed in under-agarose chemotaxis assays. In particular, our data indicate that fam49 KO cells tend to display a lower degree of global protrusive activity, their protrusions extend more slowly and are less discrete, and the cells end up moving at lower speeds and with higher directional persistence. This phenotype was substantially rescued by FAM49 re-expression. While re-expressing FAM49 in fam49 KO cells we generated putative FAM49 overexpressor cells; compared to wild-type cells, they displayed atypically thin pseudopods and what seemed to be an excessively dynamic, and perhaps less coordinated, protrusive behaviour. Additional data in our study suggest that pseudopods made by fam49 KO cells are still driven by SCAR/WAVE, which is clearly not being replaced by WASP (as is now known to be the case in D. discoideum cells lacking a functional SCAR/WAVE complex). Nonetheless, the peculiar dynamics of those pseudopods imply that SCAR/WAVE’s activity is regulated differently when FAM49 is lost, though it remains to be determined how. This thesis is the first report of a dedicated study on FAM49 and lays the foundation for future research on it.

What can exploratory modelling tell us about the ecobiology of European sea bass ( Dicentrarchus labrax ): a comprehensive overview

European sea bass, Dicentrarchus labrax, is a highly valuable species in Europe, both for aquaculture in the Mediterranean Sea and for commercial and recreational fisheries in the North East Atlantic Ocean. Subjected to increasing fishing pressure, the wild population has recently experienced significant recruitment fluctuation as well as a northward extension of its distribution area in the North Sea. While the nature of the ecological and/or physiological processes involved remains unresolved, ontogenetic habitat shifts and adult site fidelity could increase the species’ vulnerability to climate change and overfishing. As managers look for expert information to propose management scenarios leading to sustainable exploitation, exploratory modelling appears to be a cost-efficient approach to enhance the understanding of recruitment dynamics and the spatio-temporal scales over which fish populations function. A conceptual modelling framework and its specific data requirements are discussed to tackle some sound ecological questions regarding this species. We consequently provide an updated review of current knowledge on bass population structure, biology and ecology. This paper will hence be particularly valuable to develop spatially-explicit models of European sea bass dynamics under environmental and anthropogenic forcing. Knowledge gaps requiring further research efforts are also reported.

Instabilities in threshold-diffusion equations with delay

The introduction of delays into ordinary or partial differential equation models is well known to facilitate the production of rich dynamics ranging from periodic solutions through to spatio-temporal chaos. In this paper we consider a class of scalar partial differential equations with a delayed threshold nonlinearity which admits exact solutions for equilibria, periodic orbits and travelling waves. Importantly we show how the spectra of periodic and travelling wave solutions can be determined in terms of the zeros of a complex analytic function. Using this as a computational tool to determine stability we show that delays can have very different effects on threshold systems with negative as opposed to positive feedback. Direct numerical simulations are used to confirm our bifurcation analysis, and to probe some of the rich behaviour possible for mixed feedback.

Rare sex or out of reach equilibrium? The dynamics of F-IS in partially clonal organisms

Background: Partially clonal organisms are very common in nature, yet the influence of partial asexuality on the temporal dynamics of genetic diversity remains poorly understood. Mathematical models accounting for clonality predict deviations only for extremely rare sex and only towards mean inbreeding coefficient (F-IS) over bar < 0. Yet in partially clonal species, both F-IS < 0 and F-IS > 0 are frequently observed also in populations where there is evidence for a significant amount of sexual reproduction. Here, we studied the joint effects of partial clonality, mutation and genetic drift with a state-and-time discrete Markov chain model to describe the dynamics of F-IS over time under increasing rates of clonality. Results: Results of the mathematical model and simulations show that partial clonality slows down the asymptotic convergence to F-IS = 0. Thus, although clonality alone does not lead to departures from Hardy-Weinberg expectations once reached the final equilibrium state, both negative and positive F-IS values can arise transiently even at intermediate rates of clonality. More importantly, such "transient" departures from Hardy Weinberg proportions may last long as clonality tunes up the temporal variation of F-IS and reduces its rate of change over time, leading to a hyperbolic increase of the maximal time needed to reach the final mean (F-IS,F-infinity) over bar value expected at equilibrium. Conclusion: Our results argue for a dynamical interpretation of F-IS in clonal populations. Negative values cannot be interpreted as unequivocal evidence for extremely scarce sex but also as intermediate rates of clonality in finite populations. Complementary observations (e.g. frequency distribution of multiloci genotypes, population history) or time series data may help to discriminate between different possible conclusions on the extent of clonality when mean (F-IS) over bar values deviating from zero and/or a large variation of F-IS over loci are observed.

Mode locking in a spatially extended neuron model: active soma and compartmental tree

Understanding the mode-locked response of excitable systems to periodic forcing has important applications in neuroscience. For example it is known that spatially extended place cells in the hippocampus are driven by the theta rhythm to generate a code conveying information about spatial location. Thus it is important to explore the role of neuronal dendrites in generating the response to periodic current injection. In this paper we pursue this using a compartmental model, with linear dynamics for each compartment, coupled to an active soma model that generates action potentials. By working with the piece-wise linear McKean model for the soma we show how the response of the whole neuron model (soma and dendrites) can be written in closed form. We exploit this to construct a stroboscopic map describing the response of the spatially extended model to periodic forcing. A linear stability analysis of this map, together with a careful treatment of the non-differentiability of the soma model, allows us to construct the Arnol'd tongue structure for 1:q states (one action potential for q cycles of forcing). Importantly we show how the presence of quasi-active membrane in the dendrites can influence the shape of tongues. Direct numerical simulations confirm our theory and further indicate that resonant dendritic membrane can enlarge the windows in parameter space for chaotic behavior. These simulations also show that the spatially extended neuron model responds differently to global as opposed to point forcing. In the former case spatio-temporal patterns of activity within an Arnol'd tongue are standing waves, whilst in the latter they are traveling waves.

Abundance of sedentary consumers and sessile organisms along the wave exposure gradient of subtropical rocky shores of the south-west Atlantic

Sedentary consumers play an important role on populations of prey and, hence, their patterns of abundance, distribution and coexistence on shores are important to evaluate their potential influence on ecosystem dynamics. Here, we aimed to describe their spatio-temporal distribution and abundance in relation to wave exposure in the intertidal rocky shores of the south-west Atlantic to provide a basis for further understanding of ecological processes in this system. The abundance and composition of the functional groups of sessile organisms and sedentary consumers were taken by sampling the intertidal of sheltered and moderately exposed shores during a period of one year. The sublittoral fringe of sheltered areas was dominated by macroalgae, while the low midlittoral was dominated by bare rock and barnacles. In contrast, filter-feeding animals prevailed at exposed shores, probably explaining the higher abundance of the predator Stramonita haemastoma at these locations. Limpets were more abundant at the midlittoral zone of all shores while sea urchins were exclusively found at the sublittoral fringe of moderately exposed shores, therefore, adding grazing pressure on these areas. The results showed patterns of coexistence, distribution and abundance of those organisms in this subtropical area, presumably as a result of wave action, competition and prey availability. It also brought insights on the influence of top-down and bottom-up processes in this area.

Rogue waves generation via nonlinear soliton collision in multiple-soliton state of a mode-locked fiber laser

We report for the first time, rogue waves generation in a mode-locked fiber laser that worked in multiple-soliton state in which hundreds of solitons occupied the whole laser cavity. Using real-time spatio-temporal intensity dynamics measurements, it is unveiled that nonlinear soliton collision accounts for the formation of rogue waves in this laser state. The nature of interactions between solitons are also discussed. Our observation may suggest similar formation mechanisms of rogue waves in other systems.

Modelling and predictive control techniques for building heating systems

Model predictive control (MPC) has often been referred to in literature as a potential method for more efficient control of building heating systems. Though a significant performance improvement can be achieved with an MPC strategy, the complexity introduced to the commissioning of the system is often prohibitive. Models are required which can capture the thermodynamic properties of the building with sufficient accuracy for meaningful predictions to be made. Furthermore, a large number of tuning weights may need to be determined to achieve a desired performance. For MPC to become a practicable alternative, these issues must be addressed. Acknowledging the impact of the external environment as well as the interaction of occupants on the thermal behaviour of the building, in this work, techniques have been developed for deriving building models from data in which large, unmeasured disturbances are present. A spatio-temporal filtering process was introduced to determine estimates of the disturbances from measured data, which were then incorporated with metaheuristic search techniques to derive high-order simulation models, capable of replicating the thermal dynamics of a building. While a high-order simulation model allowed for control strategies to be analysed and compared, low-order models were required for use within the MPC strategy itself. The disturbance estimation techniques were adapted for use with system-identification methods to derive such models. MPC formulations were then derived to enable a more straightforward commissioning process and implemented in a validated simulation platform. A prioritised-objective strategy was developed which allowed for the tuning parameters typically associated with an MPC cost function to be omitted from the formulation by separation of the conflicting requirements of comfort satisfaction and energy reduction within a lexicographic framework. The improved ability of the formulation to be set-up and reconfigured in faulted conditions was shown.

Airborne particulate matter at a bus station : concentration levels and governing parameters

Traffic emissions are an important contributor to ambient air pollution, especially in large cities featuring extensive and high density traffic networks. Bus fleets represent a significant part of inner city traffic causing an increase in exposure to general public, passengers and drivers along bus routes and at bus stations. Limited information is available on quantification of the levels, and governing parameters affecting the air pollution exposure at bus stations. The presented study investigated the bus emissions-dominated ambient air in a large, inner city bus station, with a specific focus on submicrometer particles. The study’s objectives were (i) quantification of the concentration levels; (ii) characterisation of the spatio-temporal variation; (iii) identification of the parameters governing the emissions levels at the bus station and (iv) assessment of the relationship between particle concentrations measured at the street level (background) and within the bus station. The results show that up to 90% of the emissions at the station are ultrafine particles (smaller than 100 nm), with the concentration levels up to 10 times the value of urban ambient air background (annual) and up to 4 times the local ambient air background. The governing parameters affecting particle concentration at the station were bus flow rate and meteorological conditions (wind velocity). Particle concentration followed a diurnal trend, with an increase in the morning and evening, associated with traffic rush hours. Passengers’ exposure could be significant compared to the average outdoor and indoor exposure levels.

Hierarchical models for 2D presence/absence data having ambiguous zeroes: With a biogeographical case study on dingo behaviour

This dissertation is primarily an applied statistical modelling investigation, motivated by a case study comprising real data and real questions. Theoretical questions on modelling and computation of normalization constants arose from pursuit of these data analytic questions. The essence of the thesis can be described as follows. Consider binary data observed on a two-dimensional lattice. A common problem with such data is the ambiguity of zeroes recorded. These may represent zero response given some threshold (presence) or that the threshold has not been triggered (absence). Suppose that the researcher wishes to estimate the effects of covariates on the binary responses, whilst taking into account underlying spatial variation, which is itself of some interest. This situation arises in many contexts and the dingo, cypress and toad case studies described in the motivation chapter are examples of this. Two main approaches to modelling and inference are investigated in this thesis. The first is frequentist and based on generalized linear models, with spatial variation modelled by using a block structure or by smoothing the residuals spatially. The EM algorithm can be used to obtain point estimates, coupled with bootstrapping or asymptotic MLE estimates for standard errors. The second approach is Bayesian and based on a three- or four-tier hierarchical model, comprising a logistic regression with covariates for the data layer, a binary Markov Random field (MRF) for the underlying spatial process, and suitable priors for parameters in these main models. The three-parameter autologistic model is a particular MRF of interest. Markov chain Monte Carlo (MCMC) methods comprising hybrid Metropolis/Gibbs samplers is suitable for computation in this situation. Model performance can be gauged by MCMC diagnostics. Model choice can be assessed by incorporating another tier in the modelling hierarchy. This requires evaluation of a normalization constant, a notoriously difficult problem. Difficulty with estimating the normalization constant for the MRF can be overcome by using a path integral approach, although this is a highly computationally intensive method. Different methods of estimating ratios of normalization constants (N Cs) are investigated, including importance sampling Monte Carlo (ISMC), dependent Monte Carlo based on MCMC simulations (MCMC), and reverse logistic regression (RLR). I develop an idea present though not fully developed in the literature, and propose the Integrated mean canonical statistic (IMCS) method for estimating log NC ratios for binary MRFs. The IMCS method falls within the framework of the newly identified path sampling methods of Gelman & Meng (1998) and outperforms ISMC, MCMC and RLR. It also does not rely on simplifying assumptions, such as ignoring spatio-temporal dependence in the process. A thorough investigation is made of the application of IMCS to the three-parameter Autologistic model. This work introduces background computations required for the full implementation of the four-tier model in Chapter 7. Two different extensions of the three-tier model to a four-tier version are investigated. The first extension incorporates temporal dependence in the underlying spatio-temporal process. The second extensions allows the successes and failures in the data layer to depend on time. The MCMC computational method is extended to incorporate the extra layer. A major contribution of the thesis is the development of a fully Bayesian approach to inference for these hierarchical models for the first time. Note: The author of this thesis has agreed to make it open access but invites people downloading the thesis to send her an email via the 'Contact Author' function.

Spatial analysis of notified dengue fever infections

This study aimed to investigate the spatial clustering and dynamic dispersion of dengue incidence in Queensland, Australia. We used Moran’s I statistic to assess the spatial autocorrelation of reported dengue cases. Spatial empirical Bayes smoothing estimates were used to display the spatial distribution of dengue in postal areas throughout Queensland. Local indicators of spatial association (LISA) maps and logistic regression models were used to identify spatial clusters and examine the spatio-temporal patterns of the spread of dengue. The results indicate that the spatial distribution of dengue was clustered during each of the three periods of 1993–1996, 1997–2000 and 2001–2004. The high-incidence clusters of dengue were primarily concentrated in the north of Queensland and low-incidence clusters occurred in the south-east of Queensland. The study concludes that the geographical range of notified dengue cases has significantly expanded in Queensland over recent years.

Hierarchical Bayesian models for estimating the extent of plant pest invasions

Plant biosecurity requires statistical tools to interpret field surveillance data in order to manage pest incursions that threaten crop production and trade. Ultimately, management decisions need to be based on the probability that an area is infested or free of a pest. Current informal approaches to delimiting pest extent rely upon expert ecological interpretation of presence / absence data over space and time. Hierarchical Bayesian models provide a cohesive statistical framework that can formally integrate the available information on both pest ecology and data. The overarching method involves constructing an observation model for the surveillance data, conditional on the hidden extent of the pest and uncertain detection sensitivity. The extent of the pest is then modelled as a dynamic invasion process that includes uncertainty in ecological parameters. Modelling approaches to assimilate this information are explored through case studies on spiralling whitefly, Aleurodicus dispersus and red banded mango caterpillar, Deanolis sublimbalis. Markov chain Monte Carlo simulation is used to estimate the probable extent of pests, given the observation and process model conditioned by surveillance data. Statistical methods, based on time-to-event models, are developed to apply hierarchical Bayesian models to early detection programs and to demonstrate area freedom from pests. The value of early detection surveillance programs is demonstrated through an application to interpret surveillance data for exotic plant pests with uncertain spread rates. The model suggests that typical early detection programs provide a moderate reduction in the probability of an area being infested but a dramatic reduction in the expected area of incursions at a given time. Estimates of spiralling whitefly extent are examined at local, district and state-wide scales. The local model estimates the rate of natural spread and the influence of host architecture, host suitability and inspector efficiency. These parameter estimates can support the development of robust surveillance programs. Hierarchical Bayesian models for the human-mediated spread of spiralling whitefly are developed for the colonisation of discrete cells connected by a modified gravity model. By estimating dispersal parameters, the model can be used to predict the extent of the pest over time. An extended model predicts the climate restricted distribution of the pest in Queensland. These novel human-mediated movement models are well suited to demonstrating area freedom at coarse spatio-temporal scales. At finer scales, and in the presence of ecological complexity, exploratory models are developed to investigate the capacity for surveillance information to estimate the extent of red banded mango caterpillar. It is apparent that excessive uncertainty about observation and ecological parameters can impose limits on inference at the scales required for effective management of response programs. The thesis contributes novel statistical approaches to estimating the extent of pests and develops applications to assist decision-making across a range of plant biosecurity surveillance activities. Hierarchical Bayesian modelling is demonstrated as both a useful analytical tool for estimating pest extent and a natural investigative paradigm for developing and focussing biosecurity programs.

Hierarchical Bayesian modelling of early detection surveillance for plant pest invasions

Early detection surveillance programs aim to find invasions of exotic plant pests and diseases before they are too widespread to eradicate. However, the value of these programs can be difficult to justify when no positive detections are made. To demonstrate the value of pest absence information provided by these programs, we use a hierarchical Bayesian framework to model estimates of incursion extent with and without surveillance. A model for the latent invasion process provides the baseline against which surveillance data are assessed. Ecological knowledge and pest management criteria are introduced into the model using informative priors for invasion parameters. Observation models assimilate information from spatio-temporal presence/absence data to accommodate imperfect detection and generate posterior estimates of pest extent. When applied to an early detection program operating in Queensland, Australia, the framework demonstrates that this typical surveillance regime provides a modest reduction in the estimate that a surveyed district is infested. More importantly, the model suggests that early detection surveillance programs can provide a dramatic reduction in the putative area of incursion and therefore offer a substantial benefit to incursion management. By mapping spatial estimates of the point probability of infestation, the model identifies where future surveillance resources can be most effectively deployed.