A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Comprehensive Thermodynamic Theory of Stability of Irreversible Processes (Cttsip)

Total Page:16

File Type:pdf, Size:1020Kb

A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Comprehensive Thermodynamic Theory of Stability of Irreversible Processes (Cttsip) International Journal of Knowledge Engineering ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, 2012, pp.-91-94. Available online at http://www.bioinfo.in/contents.php?id=40 A THEORETICAL STUDY OF BIOLOGICAL LOTKA-VOLTERRA ECOLOGICAL MODEL USING COMPREHENSIVE THERMODYNAMIC THEORY OF STABILITY OF IRREVERSIBLE PROCESSES (CTTSIP) RAWAT S.G.1, BHALEKAR A.A.2 AND TANGDE V.M.3* 1Department of Applied Chemistry, Priyadarshini College of Engineering, Nagpur, MS, India. 2Department of Chemistry, Rashtrasant Tukadoji Maharaj Nagpur University, Nagpur, MS, India. 3Department of Applied Chemistry, Smt. Bhagwati Chaturvedi College of Engineering, Nagpur, MS, India. *Corresponding Author: Email- [email protected] Received: February 28, 2012; Accepted: March 06, 2012 Abstract- The dynamical relationship between predator and prey (Lotka-Volterra model or host-parasitoid system)[1] is one of the dominant themes in ecology. It was observed from the population data that interaction between a pair of predator-prey influences the population growth of both the species. This paper presents the study of thermodynamic stability of periodic Lotka-Volterra system against Prey popula- tion perturbation. The thermodynamic stability of representative model of Lotka-Volterra ecosystem has been investigated using proposed thermodynamic Lyapunov function in CTTSIP [2, 3, 4, 5] which follows the steps of Lyapunov's second method (also termed as direct meth- od) of stability of motion[6, 7]. The thermodynamic Lyapunov function, used herein is the excess rate of entropy production in thermodynam- ic perturbation space that conforms well with the dictates of second law of thermodynamics[8, 9]. Moreover, present study reveals the re- gions of stability, asymptotic stability and instability. Keywords- Lotka-Volterra model, perturbation, predator, prey, Lyapunov function, irreversible thermodynamics Citation: Rawat S.G., Bhalekar A.A. and Tangde V.M. (2012) A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Com- prehensive Thermodynamic Theory of Stability of Irreversible Processes (CTTSIP). International Journal of Knowledge Engineering, ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, pp.-91-94. Copyright: Copyright©2012 Rawat S.G., et al. This is an open-access article distributed under the terms of the Creative Commons Attribu- tion License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cred- ited. Introduction an appropriate thermodynamic space using the entropy source Lyapunov's second method of stability of motion [6, 7, 10, 11, 12] strength function, the latter is a positive definite quantity as per has a generality element, which is very much akin to that of ther- the second law of thermodynamics that appears in Clausius- modynamics. In this direct method there is no need to rigorously Duhem inequality[8, 9, 15], that reads as, solve the differential equation of motion. Indeed, this course is ds adopted in Lyapunov's first method [6, 10, 11, 12], wherein one rs+divJ ss = ³ 0 tackles every case individually. This then allows the investigator to dt (1) map the regions of stability and instability (if any). In essence, Lyapunov's second method involves effecting of a sufficiently r J s small perturbations and then observing the response in the corre- where, is the mass density, s is per unit mass entropy, sponding perturbation space. Previously, the stability aspects of s s some elementary chemical reactions including autocatalytic ones is the entropy flux density, is the entropy source strength [13] and industrial sulphur dioxide oxidation [14] has already been and t is the time. investigated using this generalized framework of CTTSIP[5]. Here 0 in the present case we will be using “Comprehensive thermody- Yti () namic theory of stability of irreversible processes(CTTSIP)” which Let 's be the thermodynamic coordinates of the given uses the fabric of Lyapunov's direct method of stability of motion. Yti () CTTSIP[2, 3, 4, 5] defines a thermodynamic Lyapunov function in process on the unperturbed trajectory and 's be that on International Journal of Knowledge Engineering ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, 2012 Bioinfo Publications 91 A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Comprehensive Thermodynamic Theory of Stability of Irreversible Processes (CTTSIP) i ()t i the perturbed one. The perturbation coordinates are de- punov’s theory of stability of motion20-24 ’s need to be fined as, sufficiently small the higher order terms in the expansion of 0 (t ) ( Y ( t ) Y ( t )); ( i 1,2, , n ) LtS () i i i (2) can be ignored and hence the expressions for it get and hence the equation of unperturbed trajectory is obtained as, restricted to, 0 0. 00 i 0 SS L (3) LS S S > i i iiii ()t (10) Thus i 's determine the thermodynamic perturbation space. 0 The constitutive equations of motion within the perturbation space LSi/ in general are described as , In Eq. (10) each is either positive or negative d L i f( , , , , t ) ( i 1,2, n ). The total time derivative of S then reads as, dt i 111 (4) dL L L d S S S i LS The thermodynamic Lyapunov function, is defined as, dt ti i dt (11) 0 LSSS( ( t ) ( t )) 0 0 (5) Si/0 On the other hand, if and hence 0 S S 0 Where and are the corresponding entropy LSi/0 source strengths. , and as stated above in view of the Obviously, Eq. (26) leads to the following two possibilities, either requirement of Lyapunov’s theory of stability of motion i ’s L L( , , , , t ) 0 need to be sufficiently small the third and higher order terms in the SS1 2 3 ( 6 ) Lt() or expansion of S can be ignored and hence the expressions LSS L(1 , 2 , 3 , , t ) 0 L ( 7 ) for S , get reduced to, 0 2 and on the unperturbed trajectory we have, 1 L () 0 S L L(0,0,0, , t ) 0 S S S i j SS 2 ij, ( 8 ) ij 0 L 2 S 1 L where, , is a differentiable function. In the generalized setup S ˆ 0 ij (12) 2 ij, ij of CTTSIP we expand the entropy source strength, S , and L thermodynamic Lyapunov function, S , in a Taylor expansion in L The total time derivative of S then reads as, i 00 terms of the perturbation coordinates 's as, dL1 22 d S S S i i j i 0 0 dt2 t dt 1 2 ii j j i i j L () 0 sS (13) S S S i i j ii2 i, j i j where we have used the following expression for the local time 0 3 1 LS S derivative of , namely: 6 i j k i,, j k i j k 0 (9 2 LSS1 Thus, even the bilateral perturbations can also be treated by this ij tt2 0 ij, ij (14) Si/0 formalism. However, if and hence Thermodynamic stability is then determined by finding out the sign 0 dL/ dt LSi/0 S , and as per the requirement of Lya- of on the lines of Lyapunov’s direct method of stability International Journal of Knowledge Engineering ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, 2012 Bioinfo Publications 92 Rawat S.G., Bhalekar A.A. and Tangde V.M. of motion[6,7,11,12]. Thus the unperturbed trajectory is said to be dx=- x x stable if , Perturbation in Prey population ( 0 ) The Thermodynamic stability of above model is then analyzed dLS LS .0 using CTTSIP under the perturbation effect in Prey population ( dt (15) The asymptotic stability is obtained if the result is either, ). Let us consider a case when the concentra- dL x L 0, S 0, tion of prey is perturbed by a suffciently small amount say, S dt (16) dx . We then have by the definition of perturbation coordinate, dL S dx=- x x 0 ˜ 0 LS 0, 0, (20) Or dt (17) Now, we investigate the effect of such perturbation on the thermo- dynamic stability of process below. Let us consider the population where, , is a strictly positive number that van- rate equation of prey as ishes only at the origin. The stability under constantly acting small disturbances, as per Malkin’s theorem [6], is obtained if in addition dx =-ax bxy dt LSi/ (21) to Eq. (16) or (17) each is finite. Following the track of CTTSIP let us directly move for deriving a suitable thermodynamic Lyapunov function for the present case, System under investigation as In this study, we investigate the thermodynamic stability of Lotka- Volterra ecosystem model (Prey-Predator interaction). A L()=-a by0 d x In ecosystem, one species do not occur in isolation of another. s T They compete for resources, assists, excludes and kills one an- (22) other. In ecology, predation describes a biological interaction and on solving the steps through-out we finally reach to the total where a predator feeds on its prey. The predator-prey relationship L is important in maintaining the balance between the different spe- time derivative of s as cies. Without prey, there would be no predators and vice-versa. e. d (L ) A g. Cheetahs and gazelles, foxes and rabbits, lion and foxes, foxes s =-2 (a by02 ) d x and goat, birds and insects,goat and grass, etc. The Lotka- dt T Volterra Prey-Predator model involves coupled equations, one (23) which describes how the prey population changes and second Conclusion which describes how the predator population changes, namely: Analytical conclusions regarding the thermodynamic stability of above system against the perturbation in prey population have put forward in following different situations, namely: dx =-ax bxy dt 1.
Recommended publications
  • An Assessment of Marine Ecosystem Damage from the Penglai 19-3 Oil Spill Accident
    Journal of Marine Science and Engineering Article An Assessment of Marine Ecosystem Damage from the Penglai 19-3 Oil Spill Accident Haiwen Han 1, Shengmao Huang 1, Shuang Liu 2,3,*, Jingjing Sha 2,3 and Xianqing Lv 1,* 1 Key Laboratory of Physical Oceanography, Ministry of Education, Ocean University of China, Qingdao 266100, China; [email protected] (H.H.); [email protected] (S.H.) 2 North China Sea Environment Monitoring Center, State Oceanic Administration (SOA), Qingdao 266033, China; [email protected] 3 Department of Environment and Ecology, Shandong Province Key Laboratory of Marine Ecology and Environment & Disaster Prevention and Mitigation, Qingdao 266100, China * Correspondence: [email protected] (S.L.); [email protected] (X.L.) Abstract: Oil spills have immediate adverse effects on marine ecological functions. Accurate as- sessment of the damage caused by the oil spill is of great significance for the protection of marine ecosystems. In this study the observation data of Chaetoceros and shellfish before and after the Penglai 19-3 oil spill in the Bohai Sea were analyzed by the least-squares fitting method and radial basis function (RBF) interpolation. Besides, an oil transport model is provided which considers both the hydrodynamic mechanism and monitoring data to accurately simulate the spatial and temporal distribution of total petroleum hydrocarbons (TPH) in the Bohai Sea. It was found that the abundance of Chaetoceros and shellfish exposed to the oil spill decreased rapidly. The biomass loss of Chaetoceros and shellfish are 7.25 × 1014 ∼ 7.28 × 1014 ind and 2.30 × 1012 ∼ 2.51 × 1012 ind in the area with TPH over 50 mg/m3 during the observation period, respectively.
    [Show full text]
  • Vegetation Demographics in Earth System Models: a Review of Progress and Priorities
    Lawrence Berkeley National Laboratory Recent Work Title Vegetation demographics in Earth System Models: A review of progress and priorities. Permalink https://escholarship.org/uc/item/3912p4m3 Journal Global change biology, 24(1) ISSN 1354-1013 Authors Fisher, Rosie A Koven, Charles D Anderegg, William RL et al. Publication Date 2018 DOI 10.1111/gcb.13910 Peer reviewed eScholarship.org Powered by the California Digital Library University of California Received: 11 April 2017 | Revised: 12 August 2017 | Accepted: 17 August 2017 DOI: 10.1111/gcb.13910 RESEARCH REVIEW Vegetation demographics in Earth System Models: A review of progress and priorities Rosie A. Fisher1 | Charles D. Koven2 | William R. L. Anderegg3 | Bradley O. Christoffersen4 | Michael C. Dietze5 | Caroline E. Farrior6 | Jennifer A. Holm2 | George C. Hurtt7 | Ryan G. Knox2 | Peter J. Lawrence1 | Jeremy W. Lichstein8 | Marcos Longo9 | Ashley M. Matheny10 | David Medvigy11 | Helene C. Muller-Landau12 | Thomas L. Powell2 | Shawn P. Serbin13 | Hisashi Sato14 | Jacquelyn K. Shuman1 | Benjamin Smith15 | Anna T. Trugman16 | Toni Viskari12 | Hans Verbeeck17 | Ensheng Weng18 | Chonggang Xu4 | Xiangtao Xu19 | Tao Zhang8 | Paul R. Moorcroft20 1National Center for Atmospheric Research, Boulder, CO, USA 2Lawrence Berkeley National Laboratory, Berkeley, CA, USA 3Department of Biology, University of Utah, Salt Lake City, UT, USA 4Los Alamos National Laboratory, Los Alamos, NM, USA 5Department of Earth and Environment, Boston University, Boston, MA, USA 6Department of Integrative Biology,
    [Show full text]
  • Emergent Biogeography of Microbial Communities in a Model Ocean
    REPORTS germ insects not only uncovers those features es- mRNA localization indeed appears to be an sential to this developmental mode but also sheds important component of long-germ embryogene- light on how the bcd-dependent anterior patterning sis, perhaps even playing a role in the transition program might have evolved. Through analysis of from the ancestral short-germ to the derived long- the regulation of the trunk gap gene Kr in Dro- germ fate. sophila and Nasonia,wehavebeenabletodem- onstrate that anterior repression of Kr is essential References and Notes for head and thorax formation and is a common 1. G. K. Davis, N. H. Patel, Annu. Rev. Entomol. 47, 669 (2002). feature of long-germ patterning. Both insects 2. T. Berleth et al., EMBO J. 7, 1749 (1988). accomplish this task through maternal, anteriorly 3. W. Driever, C. Nusslein-Volhard, Cell 54, 83 (1988). localized factors that either indirectly (Drosophila) 4. J. Lynch, C. Desplan, Curr. Biol. 13, R557 (2003). or directly (Nasonia) repress Kr and, hence, trunk 5. J. A. Lynch, A. E. Brent, D. S. Leaf, M. A. Pultz, C. Desplan, Nature 439, 728 (2006). fates. In Drosophila, the terminal system and bcd 6. J. Savard et al., Genome Res. 16, 1334 (2006). regulate expression of gap genes, including Dm-gt, 7. G. Struhl, P. Johnston, P. A. Lawrence, Cell 69, 237 (1992). that repress Dm-Kr. Nasonia’s bcd-independent 8. A. Preiss, U. B. Rosenberg, A. Kienlin, E. Seifert, long-germ embryos must solve the same problem, H. Jackle, Nature 313, 27 (1985). Fig. 4.
    [Show full text]
  • Meta-Ecosystems: a Theoretical Framework for a Spatial Ecosystem Ecology
    Ecology Letters, (2003) 6: 673–679 doi: 10.1046/j.1461-0248.2003.00483.x IDEAS AND PERSPECTIVES Meta-ecosystems: a theoretical framework for a spatial ecosystem ecology Abstract Michel Loreau1*, Nicolas This contribution proposes the meta-ecosystem concept as a natural extension of the Mouquet2,4 and Robert D. Holt3 metapopulation and metacommunity concepts. A meta-ecosystem is defined as a set of 1Laboratoire d’Ecologie, UMR ecosystems connected by spatial flows of energy, materials and organisms across 7625, Ecole Normale Supe´rieure, ecosystem boundaries. This concept provides a powerful theoretical tool to understand 46 rue d’Ulm, F–75230 Paris the emergent properties that arise from spatial coupling of local ecosystems, such as Cedex 05, France global source–sink constraints, diversity–productivity patterns, stabilization of ecosystem 2Department of Biological processes and indirect interactions at landscape or regional scales. The meta-ecosystem Science and School of perspective thereby has the potential to integrate the perspectives of community and Computational Science and Information Technology, Florida landscape ecology, to provide novel fundamental insights into the dynamics and State University, Tallahassee, FL functioning of ecosystems from local to global scales, and to increase our ability to 32306-1100, USA predict the consequences of land-use changes on biodiversity and the provision of 3Department of Zoology, ecosystem services to human societies. University of Florida, 111 Bartram Hall, Gainesville, FL Keywords 32611-8525,
    [Show full text]
  • Chapter 4 – Steady State Models
    CHAPTER 4 Steady State Models X.-Z. Kong, F.-L. Xu1, W. He, W.-X. Liu, B. Yang Peking University, Beijing, China 1Corresponding author: E-mail: xufl@urban.pku.edu.cn OUTLINE 4.1 Steady State Model: Ecopath as an 4.2.4.2 Changes in Ecosystem Example 66 Functioning 79 4.1.1 Steady State Model 66 4.2.4.3 Collapse in the Food 4.1.2 Ecopath Model 66 Web: Differences in 4.1.3 Future Perspectives 68 Structure 81 4.2.4.4 Toward an Immature 4.2 Ecopath Model for a Large Chinese but Stable Ecosystem 83 Lake: A Case Study 68 4.2.4.5 Potential Driving 4.2.1 Introduction 68 Factors and 4.2.2 Study Site 70 Underlying 4.2.3 Model Development 73 Mechanisms 84 4.2.3.1 Model Construction 4.2.4.6 Hints for Future Lake and Parameterization 73 Fishery and 4.2.3.2 Evaluation of Ecosystem Restoration 85 Functioning 74 4.2.3.3 Determination of 4.3 Conclusions 86 Trophic Level 76 References 86 4.2.4 Results and Discussion 76 4.2.4.1 Basic Model Performance 76 Ecological Model Types, Volume 28 Ó 2016 Elsevier B.V. http://dx.doi.org/10.1016/B978-0-444-63623-2.00004-9 65 All rights reserved. 66 4. STEADY STATE MODELS 4.1 STEADY STATE MODEL: ECOPATH AS AN EXAMPLE 4.1.1 Steady State Model Steady state ecological models are established to describe conditions in which the modeled components (mass or energy) are stable, i.e., do not change over time (Jørgensen and Fath, 2011).
    [Show full text]
  • Models for an Ecosystem Approach to Fisheries
    ISSN 0429-9345 FAO FISHERIES 477 TECHNICAL PAPER 477 Models for an ecosystem approach to fisheries Models for an ecosystem approach to fisheries This report reviews the methods available for assessing the impacts of interactions between species and fisheries and their implications for marine fisheries management. A brief description of the various modelling approaches currently in existence is provided, highlighting in particular features of these models that have general relevance to the field of ecosystem approach to fisheries (EAF). The report concentrates on the currently available models representative of general types such as bionergetic models, predator-prey models and minimally realistic models. Short descriptions are given of model parameters, assumptions and data requirements. Some of the advantages, disadvantages and limitations of each of the approaches in addressing questions pertaining to EAF are discussed. The report concludes with some recommendations for moving forward in the development of multispecies and ecosystem models and for the prudent use of the currently available models as tools for provision of scientific information on fisheries in an ecosystem context. FAO Cover: Illustration by Elda Longo FAO FISHERIES Models for an ecosystem TECHNICAL PAPER approach to fisheries 477 by Éva E. Plagányi University of Cape Town South Africa FOOD AND AGRICULTURE AND ORGANIZATION OF THE UNITED NATIONS Rome, 2007 The designations employed and the presentation of material in this information product do not imply the expression of any opinion whatsoever on the part of the Food and Agriculture Organization of the United Nations concerning the legal or development status of any country, territory, city or area or of its authorities, or concerning the delimitation of its frontiers or boundaries.
    [Show full text]
  • Can We Detect Ecosystem Critical Transitions and Signals of Changing Resilience from Paleo-Ecological Records? Zofia E
    Can we detect ecosystem critical transitions and signals of changing resilience from paleo-ecological records? Zofia E. Taranu, Stephen R. Carpenter, Victor Frossard, Jean-Philippe Jenny, Zoe Thomas, Jesse C. Vermaire, Marie-Elodie Perga To cite this version: Zofia E. Taranu, Stephen R. Carpenter, Victor Frossard, Jean-Philippe Jenny, Zoe Thomas, etal.. Can we detect ecosystem critical transitions and signals of changing resilience from paleo-ecological records?. Ecosphere, Ecological Society of America, 2018, 9 (10), 10.1002/ecs2.2438. hal-01959637 HAL Id: hal-01959637 https://hal.archives-ouvertes.fr/hal-01959637 Submitted on 18 Dec 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution| 4.0 International License Can we detect ecosystem critical transitions and signals of changing resilience from paleo-ecological records? 1, 2 3 3,4 € 5 ZOFIA E. TARANU, STEPHEN R. CARPENTER, VICTOR FROSSARD, JEAN-PHILIPPE JENNY, ZOE THOMAS, 6 7 JESSE C. VERMAIRE, AND MARIE-ELODIE PERGA 1Department of Biology, University
    [Show full text]
  • Models of the World's Large Marine Ecosystems: GEF/LME Global
    Intergovernmental Oceanographic Commission technical series 80 Models of the World’s Large Marine Ecosystems UNESCO Intergovernmental Oceanographic Commission technical series 80 Models of the World’s Large Marine Ecosystems* GEF/LME global project Promoting Ecosystem-based Approaches to Fisheries Conservation and Large Marine Ecosystems UNESCO 2008 * As submitted to IOC Technical Series, UNESCO, 22 October 2008 IOC Technical Series No. 80 Paris, 20 October 2008 English only The designations employed and the presentation of the material in this publication do not imply the expression of any opinion whatsoever on the part of the Secretariats of UNESCO and IOC concerning the legal status of any country or territory, or its authorities, or concerning the delimitation of the frontiers of any country or territory. For bibliographic purposes, this document should be cited as follows: Models of the World’s Large Marine Ecosystems GEF/LME global project Promoting Ecosystem-based Approaches to Fisheries Conservation and Large Marine Ecosystems IOC Technical Series No. 80. UNESCO, 2008 (English) Editors: Villy Christensen1, Carl J. Walters1, Robert Ahrens1, Jackie Alder2, Joe Buszowski1, Line Bang Christensen1, William W.L. Cheung1, John Dunne3, Rainer Froese4, Vasiliki Karpouzi1, Kristin Kastner5, Kelly Kearney6, Sherman Lai1, Vicki Lam1, Maria L.D. Palomares1,7, Aja Peters-Mason8, Chiara Piroddi1, Jorge L. Sarmiento6, Jeroen Steenbeek1, Rashid Sumaila1, Reg Watson1, Dirk Zeller1, and Daniel Pauly1. Technical Editor: Jair Torres 1 Fisheries Centre,
    [Show full text]
  • Mapping Fisheries Onto Marine Ecosystems for Regional, Oceanic and Global Integrations
    Large Marine Ecosystems of the World G. Hempel and K. Sherman (Editors) 0 2003 Elsevier B. V. All rights reserved 375 Mapping Fisheries onto Marine Ecosystems for Regional, Oceanic and Global Integrations Reg Watson, Daniel Pauly, Villy Christensen,Rainer Froese,Alan Longhurst, Trevor Platt, Shubha Sathyendranath,Kenneth Sherman,John O'Reilly, and Peter Celone ABSTRACT Research on ecosystem-based fisheries management, marine biodiversity conservation,and other marine fields requires appropriate maps of the major natural regions of the oceans,and their ecosystems.A global ocean classification system proposed by T. Platt and S. Sathyendranathand implemented by A.R. Longhurst, defined largely by physical parameters that subdivide the oceansinto four 'biomes' and 57 'biogeochemicalprovinces' (BGCPs),is merged with the systemof 64 Large Marine Ecosystems (LMEs) identified by K. Sherman and colleagues as transboundarygeographic coastaland watershedunits. This arrangementenhances each of the systems,and renders them mutually compatible. LMEs are ecologically defined to serve as a framework for the assessmentand managementof coastal fisheries and environments including watersheds,while the BGCPshave physical definitions, including borders defined by natural features, and extend over open oceanregions. The combined mapping will, for example, allow the computation of GI5-derived properties such as temperature,primary production, and their analysis in relation to fishery abundancedata for any study area in the combinedsystem. A further useful aspectof the integration is that it allows for the quantification, even within the EEZs of various countries, of the distribution of marine features (e.g. primary production, coral reef areas)so far not straightforwardly associatedwith different coastalstates. Applications to shelf,coral reef and oceanicfisheries, and to the mapping of marine biodiversity are briefly discussed.
    [Show full text]
  • Parameter Optimisations of a Marine Ecosystem Model in the Southern Ocean
    PARAMETER OPTIMISATIONS OF A MARINE ECOSYSTEM MODEL IN THE SOUTHERN OCEAN by Mehera Kidston A thesis submitted in fulfilment of the requirements for the degree of Doctor of Philosophy within the University of New South Wales, Australia 2010 i Abstract The Southern Ocean plays a central role in global biological production (Sarmiento et al., 2004). Quantifying the mechanisms that regulate magnitude and variability in plankton biomass in the Southern Ocean is a key challenge in understanding controls on global climate (Busalacchi, 2004). Numerical models are an integral means of understanding large scale views of seasonal plankton cycles (Fasham et al., 1990). Marine biogeochemical models are characterised by non-linear dynamics and the annual model trajectories are highly sensitive to the parameters used to run them. Experiments were performed to select a stochastic inverse method to objectively estimate the parameters of a simple four component nitrogen based mixed-layer marine ecosystem model for the Southern Ocean. Twin experiments using the Metropolis-Hastings algorithm and simulated annealing show that simulated annealing holds promise as a standard means of assigning the parameters of marine biogeochemical models in a way that improves the model agreement to available observations. SeaWiFS surface chlorophyll estimates in the Sub-Antarctic Zone show low concentrations south west of Tasmania and high concentrations south east of Tasmania. Simulated annealing was used to estimate the model parameters at two locations in the Sub-Antartic Zone (station P1 at 140 E , 46.5 S and station P3 at 152 E , 45.5 S ) through assimilation of SeaWiFS chlorophyll observations. Model parameter estimates were compared to in situ parameter estimates from the SAZ-Sense (Sub-Antarctic Zone Sensitivity to environmental change) project stations P1 and P3 in the austral summer of 2007.
    [Show full text]
  • Cascading Ecological Effects of Eliminating Fishery Discards
    ARTICLE Received 19 Nov 2013 | Accepted 15 Apr 2014 | Published 13 May 2014 DOI: 10.1038/ncomms4893 OPEN Cascading ecological effects of eliminating fishery discards Michael R. Heath1, Robin M. Cook1, Angus I. Cameron1, David J. Morris1 & Douglas C. Speirs1 Discarding by fisheries is perceived as contrary to responsible harvesting. Legislation seeking to end the practice is being introduced in many jurisdictions. However, discarded fish are food for a range of scavenging species; so, ending discarding may have ecological consequences. Here we investigate the sensitivity of ecological effects to discarding policies using an ecosystem model of the North Sea—a region where 30–40% of trawled fish catch is currently discarded. We show that landing the entire catch while fishing as usual has conservation penalties for seabirds, marine mammals and seabed fauna, and no benefit to fish stocks. However, combining landing obligations with changes in fishing practices to limit the capture of unwanted fish results in trophic cascades that can benefit birds, mammals and most fish stocks. Our results highlight the importance of considering the broader ecosystem consequences of fishery management policy, since species interactions may dissipate or negate intended benefits. 1 Department of Mathematics and Statistics, University of Strathclyde, Livingstone Tower, 26 Richmond Street, Glasgow G1 1XH, UK. Correspondence and requests for materials should be addressed to M.R.H. (email: [email protected]). NATURE COMMUNICATIONS | 5:3893 | DOI: 10.1038/ncomms4893 | www.nature.com/naturecommunications 1 & 2014 Macmillan Publishers Limited. All rights reserved. ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4893 ood subsidies to wildlife as a result of human activity are nitrogen mass in groups of functionally similar taxa and recognized as having an important effect on terrestrial and materials31 rather than as individual species32, but spans the Faquatic ecosystems1.
    [Show full text]
  • Numerical Modeling of Bioluminescence Distributions in the Coastal Ocean*
    1060 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 20 Numerical Modeling of Bioluminescence Distributions in the Coastal Ocean* IGOR SHULMAN Naval Research Laboratory, Stennis Space Center, Mississippi STEVEN H. D. HADDOCK Monterey Bay Aquarium Research Institute, Moss Landing, California DENNIS J. MCGILLICUDDY JR. Woods Hole Oceanographic Institution, Woods Hole, Massachusetts JEFFREY D. PADUAN Naval Postgraduate School, Monterey, California W. P AUL BISSETT Florida Environmental Research Institute, Tampa, Florida (Manuscript received 30 April 2002, in ®nal form 31 January 2003) ABSTRACT Bioluminescence (BL) predictability experiments (predictions of the intensity, depth, and distance offshore of the BL maximum) were conducted using an advective±diffusive tracer model with velocities and diffusivities from a ®ne-resolution model of the Monterey Bay, California, area. For tracer initialization, observations were assimilated into the tracer model while velocities and diffusivities were taken from the hydrodynamic model and kept unchanged during the initialization process. This dynamic initialization procedure provides an equilibrium tracer distribution that is balanced with the velocity and diffusivity ®elds from the hydrodynamic model. This equilibrium BL distribution was used as the initial BL ®eld for 3 days of prognostic calculations. Two cross-shore surveys of bioluminescence data conducted at two locations (north of the bay and inside the bay) were used in four numerical experiments designed to estimate the limits of bioluminescence predictions by tracers. The cross-shore sections extended to around 25 km offshore, they were around 30 m deep, and on average they were approximately 35 km apart from each other. Bioluminescence predictability experiments demonstrated a strong utility of the tracer model (combined with limited bioluminescence observations and with the output from a circulation model) in predicting (over a 72-h period and over 25±35-km distances) the location and intensity of the BL maximum.
    [Show full text]