<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2024.1374869</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Hindcasting estuary ecological states using sediment cores, modelled historic nutrient loads, and a Bayesian network</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Hale</surname>
<given-names>Rachel</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/886848"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zeldis</surname>
<given-names>John</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1528372"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Dudley</surname>
<given-names>Bruce D.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Haddadchi</surname>
<given-names>Arman</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Plew</surname>
<given-names>David</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/548683"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shankar</surname>
<given-names>Ude</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Swales</surname>
<given-names>Andrew</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2672610"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Roberts</surname>
<given-names>Keryn</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>O&#x2019;Connell-Milne</surname>
<given-names>Sorrel</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Verburg</surname>
<given-names>Piet</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>National Institute of Water &amp; Atmospheric Research (NIWA) Nelson</institution>, <addr-line>Nelson</addr-line>, <country>New Zealand</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>National Institute of Water &amp; Atmospheric Research (NIWA) Christchurch</institution>, <addr-line>Christchurch</addr-line>, <country>New Zealand</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>National Institute of Water &amp; Atmospheric Research (NIWA) Hamilton</institution>, <addr-line>Hamilton</addr-line>, <country>New Zealand</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Salt Ecology</institution>, <addr-line>Nelson</addr-line>, <country>New Zealand</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>National Institute of Water &amp; Atmospheric Research (NIWA) Wellington</institution>, <addr-line>Wellington</addr-line>, <country>New Zealand</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Angel P&#xe9;rez-Ruzafa, University of Murcia, Spain</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Autumn Oczkowski, United States Environmental Protection Agency (EPA), United States</p>
<p>Savanna Barry, University of Florida, United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Rachel Hale, <email xlink:href="mailto:rachel.hale@niwa.co.nz">rachel.hale@niwa.co.nz</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>21</day>
<month>03</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>11</volume>
<elocation-id>1374869</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>03</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Hale, Zeldis, Dudley, Haddadchi, Plew, Shankar, Swales, Roberts, O&#x2019;Connell-Milne and Verburg</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Hale, Zeldis, Dudley, Haddadchi, Plew, Shankar, Swales, Roberts, O&#x2019;Connell-Milne and Verburg</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Estuaries are impacted by catchment land use changes, driving degradation associated with eutrophication and alterations in sediment dynamics. Estuarine ecological monitoring has typically covered only recent periods, so magnitudes and trajectories of degradation are often poorly described. Here, we develop a multi-method approach to hindcast historic estuary condition. We determined geochemical properties in dated sediment cores from two sites in the eutrophic New River Estuary (NRE), Aotearoa New Zealand. Mud, organic matter, carbon, and nitrogen (N) contents increased from ~1847 to 2019 (~20 to 100%, ~1 to ~7%, ~0.2 to ~3%, and ~0.03 to ~0.37% respectively). Historic land use modelling shows increasing catchment loads and water column total N rising from 165 to 663 mg m<sup>-3</sup>, with more rapid change in the late 20<sup>th</sup> and early 21<sup>st</sup> centuries. Strong correlations (R<sup>2 =</sup> 0.93-0.95) were found between sediment N content and water column N concentrations. Trophic indicators hindcasted by a Bayesian Belief Network worsened from the mid-20<sup>th</sup> to 21<sup>st</sup> centuries. This degradation is associated with increased sheep, beef and dairy farming in the catchment. Our approach identifies nutrient and sediment loadings associated with &#x2018;tipping points&#x2019; for indicators and can inform policy by identification of hazardous levels of contaminant loading for other estuaries.</p>
</abstract>
<kwd-group>
<kwd>estuaries</kwd>
<kwd>sedimentation</kwd>
<kwd>eutrophication</kwd>
<kwd>carbon cycling</kwd>
<kwd>nitrogen cycling</kwd>
<kwd>sediment nutrients</kwd>
<kwd>sediment stable isotopes</kwd>
<kwd>nutrient load modelling</kwd>
</kwd-group>
<contract-sponsor id="cn001">Ministry for Business Innovation and Employment<named-content content-type="fundref-id">10.13039/501100004629</named-content>
</contract-sponsor>
<counts>
<fig-count count="8"/>
<table-count count="5"/>
<equation-count count="0"/>
<ref-count count="105"/>
<page-count count="19"/>
<word-count count="10860"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Marine Ecosystem Ecology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Estuaries are among earth&#x2019;s most anthropogenically impacted ecosystems and threatened habitats (<xref ref-type="bibr" rid="B43">Kennish, 2002</xref>; <xref ref-type="bibr" rid="B4">Barbier et&#xa0;al., 2011</xref>). Globally, estuarine ecological degradation has accelerated over the past 150 to 300 years (<xref ref-type="bibr" rid="B55">Lotze et&#xa0;al., 2006</xref>) as a result of anthropogenic overexploitation, habitat transformation and pollution (<xref ref-type="bibr" rid="B41">Jackson et&#xa0;al., 2001</xref>; <xref ref-type="bibr" rid="B56">Lotze and Milewski, 2004</xref>; <xref ref-type="bibr" rid="B57">Lotze et&#xa0;al., 2005</xref>), particularly from increased diffuse-source catchment inputs of sediment and nutrients (<xref ref-type="bibr" rid="B18">Conley, 1999</xref>; <xref ref-type="bibr" rid="B53">Le Moal et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B94">Vieillard et&#xa0;al., 2020</xref>). The effects of these contaminant load increases include eutrophication (phytoplankton and macroalgal blooms; <xref ref-type="bibr" rid="B7">Barr et&#xa0;al., 2020</xref>), and increased sedimentation and organic content in sediments (<xref ref-type="bibr" rid="B28">Engelsen et&#xa0;al., 2008</xref>), with attendant ecological disruption including increased light attenuation, hypoxia and sulfide toxicity (<xref ref-type="bibr" rid="B39">Howarth and Marino, 2006</xref>; <xref ref-type="bibr" rid="B30">Fowler et&#xa0;al., 2013</xref>). These stressors can result in loss of important macrofaunal communities and habitats including seagrass beds and their ecosystem services (<xref ref-type="bibr" rid="B32">Green et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B70">Pratt et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B89">Sutula et&#xa0;al., 2014</xref>).</p>
<p>Estuaries in Aotearoa New Zealand have been increasingly impacted by catchment land use changes causing excessive nutrient and sediment inputs in freshwater flows (<xref ref-type="bibr" rid="B54">Lohrer et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B66">Morrison et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B35">Handley et&#xa0;al., 2020</xref>), especially since European colonization (<xref ref-type="bibr" rid="B102">Zeldis et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B84">Snelder et&#xa0;al., 2018</xref>). Increased sedimentation due to sediment losses from land has been attributed to deforestation for agriculture (<xref ref-type="bibr" rid="B40">Hume et&#xa0;al., 1992</xref>), and agricultural land uses and urbanization have increased nutrient inputs to freshwater (<xref ref-type="bibr" rid="B47">Larned et&#xa0;al., 2016</xref>) and estuaries (<xref ref-type="bibr" rid="B24">Dudley et&#xa0;al., 2020</xref>). Major changes in landcover in New Zealand resulted from large-scale catchment deforestation and conversion to pastoral agriculture during the late 19<sup>th</sup> and early 20<sup>th</sup> centuries. Technological advances, such as fertilizer top-dressing, increased agricultural production in the mid-20<sup>th</sup> century, and increases in irrigation, stock densities and fertilizer application to pasture after the mid-20<sup>th</sup> century (<xref ref-type="bibr" rid="B62">Ministry for the Environment, 2007</xref>; <xref ref-type="bibr" rid="B38">Howard-Williams et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B65">Ministry for the Environment and Stats NZ, 2019</xref>) allowing the intensification of agricultural land use and an increase in the area of urban land cover and human populations have resulted in increases of estuarine contaminant loads.</p>
<p>To manage these impacts, there are current New Zealand National Planning Framework (NPF) requirements that local governments develop policy statements (e.g. the National Policy Statement for Freshwater Management; <xref ref-type="bibr" rid="B64">Ministry for the Environment, 2023</xref>) that consider receiving waters (including estuaries) when setting freshwater and marine environmental limits and targets for restoration. Toward this, there is demand for management tools that quantify relationships between nutrient and sediment loads and estuary ecological condition (<xref ref-type="bibr" rid="B91">Townsend and Lohrer, 2015</xref>; <xref ref-type="bibr" rid="B68">Plew et&#xa0;al., 2020</xref>). Recent work in New Zealand has delivered a series of tools to link contaminant loads to estuarine trophic state (<xref ref-type="bibr" rid="B69">Plew et&#xa0;al., 2018</xref>, <xref ref-type="bibr" rid="B68">2020</xref>; <xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>), although these have only been applied in the context of modern estuarine conditions. Furthermore, quantitative monitoring data describing the progression and impact of nutrient/sediment stressors in New Zealand&#x2019;s estuaries was usually begun only within the last two decades, often after degradation of estuary condition had commenced or was well advanced (<xref ref-type="bibr" rid="B36">Hewitt et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B24">Dudley et&#xa0;al., 2020</xref>). This has made it difficult to ascertain historical estuary condition and trajectories of change.</p>
<p>Sediments deposited within estuaries contain a wealth of information about the history and time course of changes of estuary biogeochemistry and ecology (<xref ref-type="bibr" rid="B83">Savage et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B98">Wigand et&#xa0;al., 2014</xref>). Here, we provide an approach to determine historic estuarine condition by combining estimates of historical nitrogen (N) loading to an estuary with sediment core biogeochemical properties and an ecological model of responses to sedimentation and nutrient loading. We hypothesize that these tools will provide an estuary-specific time series of ecological response to contaminant load changes and allow examination of estuarine condition tipping points (e.g., sediment and/or nutrient inputs at which seagrass or macrobenthic habitats and health were lost). Such knowledge could be used to determine realistic restoration objectives and their associated contaminant input levels, thereby providing management targets for restoration and catchment limit settings.</p>
<p>To test this hypothesis, we relate hindcasts of estuarine nutrient loading derived from catchment land use information to sediment accumulation rates (SAR) and geochemical and biological attributes in dated strata from sediment cores in a eutrophic estuary: New River Estuary (NRE) in Southland, New Zealand (<xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>). This estuary has a history of increasing nutrient and sediment loading corresponding with land clearance for sheep and beef farming from the early 20<sup>th</sup> century, with intensive dairy production in the late 20<sup>th</sup> century. It has also received urban-sourced wastewater nutrient loading from the adjacent Invercargill city, all of which have led to ecological change (<xref ref-type="bibr" rid="B25">Dudley et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B72">Roberts et&#xa0;al., 2022a</xref>). We use contents of mud, photosynthetic pigments, carbon (C) and N isotopes, and bulk C and N from dated sediment core layers to describe estuary geochemical and biological conditions through time, back to 1847. We then estimate historic and present-day nutrient loading to the NRE by applying archival GIS-based land use layers within a nutrient runoff model (<xref ref-type="bibr" rid="B27">Elliott et&#xa0;al., 2016</xref>), and estuary nutrient concentrations arising from those loads using an estuary mixing model (<xref ref-type="bibr" rid="B69">Plew et&#xa0;al., 2018</xref>). Finally, these inputs are combined with other known drivers of estuary condition within a Bayesian network model (<xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>) to hindcast values of a range of estuary ecological health indicators and an overall estuary health score, from 1847 to the present day.</p>
<p>To date, Bayesian modelling has only been applied to make present-day descriptions of ecological condition for New Zealand estuaries. However, by combining sediment core results with historic nutrient load and concentration estimates, the modelling can show &#x2018;from where we have come&#x2019; in terms of estuary ecological condition over time compared to an earlier, pre-impacted state. Examining those historic states (including the pre-European colonization less impacted state), perspectives may be gained on when important changes in ecological condition occurred, along with co-occurring land use, and nutrient and sediment load forcings. By providing views of past ecological states and their drivers, we can inform efforts to set management targets and accompanying load limits for estuary remediation.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Sediment core collection</title>
<p>Replicate sediment cores (n = 3) were collected from two depositional sites (Upper: 46&#xb0;25&#x2019;51.7&#x201d;S 168&#xb0;19&#x2019;35.5&#x201d;E, Central: 46&#xb0;26&#x2019;28.4&#x201d;S 168&#xb0;19&#x2019;39.8&#x201d;E, WGS84) in the Waihopai Arm of the NRE, Southland, Aotearoa New Zealand (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>), a region of the estuary where the most significant sedimentation has occurred with observed eutrophic changes (increased macroalgal cover and reduced dissolved oxygen). Sediment cores were collected using an Aquatic Research Large Bore Universal Sediment Corer fitted with 1.2 m x 102 mm OD (95 mm ID) polycarbonate core barrels with sharpened ends. Sediment cores of 1.0 - 1.2 m were collected by driving the cores into the sediment using a 3 kg percussion hammer attached to the top of a polyethylene and stainless-steel core head unit. The head unit is fitted with a one-way check valve that seals automatically, creating a vacuum when retrieving the core. The sediment around the core was dug out and removed using spades to a depth of 600 - 700 mm to facilitate easier retrieval. Polyethylene end caps and EVA inserts were fitted after retrieval and taped in place with duct tape, prior to shipping. Each core was allocated to be subsampled for layer dating, physical and biogeochemical analysis (particle size, phytopigment analysis, organic matter content, carbon content, nitrogen content, &#x3b4;<sup>13</sup>C, &#x3b4;<sup>15</sup>N, and sediment metal content), or stored (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Location of the Upper and Central sediment coring sites within <bold>(A)</bold> Aotearoa/New Zealand showing the location of the Waihopai Estuary and its upstream catchment within Southland, <bold>(B)</bold> the catchment river network, with nutrient load monitoring sites and point sources, and <bold>(C)</bold> the lower New River Estuary and its main tributaries, the Waihopai River and the &#x14c;reti River, and location of the Invercargill Wastewater Treatment Plant (WWTP) discharge at Clifton.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g001.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Details of the sediment cores collected from the two locations (Upper, Central) in the Waihopai Estuary, core length retrieved (cm), estimated compression during collection, and the analyses undertaken.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Location</th>
<th valign="top" align="left">Core</th>
<th valign="top" align="left">Length (cm)</th>
<th valign="top" align="left">Compression (cm)</th>
<th valign="top" align="left">Analyses</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Upper</td>
<td valign="top" align="left">UW1A</td>
<td valign="top" align="left">95.5</td>
<td valign="top" align="left">17</td>
<td valign="top" align="left">Dating refinement<break/>Physical and biogeochemical</td>
</tr>
<tr>
<td valign="top" align="left">Upper</td>
<td valign="top" align="left">UW1B</td>
<td valign="top" align="left">83</td>
<td valign="top" align="left">21</td>
<td valign="top" align="left">Stored</td>
</tr>
<tr>
<td valign="top" align="left">Upper</td>
<td valign="top" align="left">UW1C</td>
<td valign="top" align="left">102.5</td>
<td valign="top" align="left">14</td>
<td valign="top" align="left">Dating</td>
</tr>
<tr>
<td valign="top" align="left">Central</td>
<td valign="top" align="left">CW3A</td>
<td valign="top" align="left">no data</td>
<td valign="top" align="left">15</td>
<td valign="top" align="left">Stored</td>
</tr>
<tr>
<td valign="top" align="left">Central</td>
<td valign="top" align="left">CW3B</td>
<td valign="top" align="left">91.5</td>
<td valign="top" align="left">13</td>
<td valign="top" align="left">Dating refinement<break/>Physical and biogeochemical</td>
</tr>
<tr>
<td valign="top" align="left">Central</td>
<td valign="top" align="left">CW3C</td>
<td valign="top" align="left">102</td>
<td valign="top" align="left">12</td>
<td valign="top" align="left">Dating</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Geochronology</title>
<p>To determine the age of sediment preserved in the deposited layers, the cores were sampled in 1-cm thick increments (40 &#x2013; 60 g dry sediment sample) every 5 cm downcore from one replicate core per site (UW1C, n =13, CW3C, n = 13, <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>) down core and dated using gamma spectroscopy. The activity of excess <sup>210</sup>Pb and <sup>137</sup>Cs (number of disintegrations per second, Bq kg<sup>-1</sup>) was counted at the ESR National Radiation Laboratory for 23 hours using a Canberra Model BE5030 hyper-pure germanium detector. Excess <sup>210</sup>Pb was determined from the <sup>226</sup>Ra (t<sub>&#xbd;</sub> 1622 y) activity and uncertainty determined following <xref ref-type="bibr" rid="B90">Swales et&#xa0;al. (2021)</xref>. Excess <sup>210</sup>Pb was used to determine time-averaged SAR using a regression analysis of natural log-transformed data (<xref ref-type="bibr" rid="B90">Swales et&#xa0;al., 2021</xref>). The maximum depth of <sup>137</sup>Cs was used to verify the excess <sup>210</sup>Pb dating estimates. <sup>137</sup>Cs deposition was first detected in 1953 in New Zealand and therefore provides an approximate baseline for appearance in the sediment layer, with peak atmospheric deposition in 1963/64 (<xref ref-type="bibr" rid="B59">Matthews, 1989</xref>; <xref ref-type="bibr" rid="B90">Swales et&#xa0;al., 2021</xref>). Following initial dating, further refinement was provided by additional analyses of layers subsampled from a second replicate (UW1A n = 3, CW3B n = 2) to confirm the <sup>210</sup>Pb SAR inflection points in the regression analysis. However, the SAR estimates carry a certain amount of uncertainty. Fitting the unsupported <sup>210</sup>Pb data underpinning the SAR calculations involves some subjectivity when selecting the most likely break points in their trajectory. When estimated ages of adjacent core slices are widely spread, this can introduce error in the estimation of dates of the break points, which may have contributed to the differences found in the estimated break point dates between the Upper (1983) and Central (2006) cores. Also, naturally occurring frequent deposition and erosion of the surface sediment layers (such as by high runoff or storm events) mean that long-term sediment accumulation captures only a small fraction of the sediment deposited at any one time (<xref ref-type="bibr" rid="B10">Bentley et&#xa0;al., 2014</xref>). Additionally, bioturbation of sediments, resulting in mixing of the sediment surface with deeper layers, results in disruption of clear strata formation until sediment is buried below the maximum mixing depth (estimated to be approximately 15 cm at another similar New Zealand estuarine site (<xref ref-type="bibr" rid="B10">Bentley et&#xa0;al., 2014</xref>). Below this horizon, we assumed that C and N were no longer remineralised and lost across the sediment water interface, preserving the C and N contents and their ratios within the deeper sediment strata (<xref ref-type="bibr" rid="B14">Canfield, 1994</xref>). We also assumed that these deeper sediments were mixed during the period when they were within the bioturbation horizon, which may disrupt sediment core strata. Finally, core compression during collection adds uncertainty to age estimations of the sediment layers. We assumed the compression occurred predominantly in the upper layers of the core that had lower bulk density, meaning recent age estimations were more affected, while the age determination of the older layers was more accurate.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Physical and biogeochemical analyses</title>
<p>Sediment cores were sliced into 1 cm samples based on the calculated age model of the layers to obtain a time series of core sediment properties. For particle size analysis, a subsample of each slice (10 g) was taken, lyophilised and ultrasounded for 4 minutes to break up flocculants. The mean particle size (&#xb5;m) and sediment mud content (%) were determined via volume particle diameter using an EyeTech time-of-transit (obscuration time) instrument (AmbiValue B.V., Dussen, Netherlands). A subsample (10 g) was taken for phytopigment analysis for chlorophyll-<italic>a</italic> and phaeophytin-<italic>a</italic> content (&#xb5;g g<sup>-1</sup>) following <xref ref-type="bibr" rid="B81">Sartory (1985)</xref>. Pigments were extracted in 95% ethanol and determined spectrophotometrically. A preservation index (the ratio of chlorophyll-<italic>a</italic>:phaeophytin-<italic>a</italic>) was determined as a proxy for organic matter preservation (<xref ref-type="bibr" rid="B52">Le Guitton et&#xa0;al., 2015</xref>). Chlorophyll-<italic>a</italic> degrades more slowly in low oxygen conditions (<xref ref-type="bibr" rid="B87">Sun et&#xa0;al., 1991</xref>, <xref ref-type="bibr" rid="B88">1993</xref>), therefore a greater preservation index can indicate either high deposition rates of chlorophyll-<italic>a</italic> or hypoxic/anoxic conditions of the overlying water (<xref ref-type="bibr" rid="B71">Rabalais et&#xa0;al., 2004</xref>). Organic matter content (%) was determined on a subsample (25 g) using loss on ignition. Contents of carbon (C, %), nitrogen (N, %) and stable isotopes &#x3b4;<sup>13</sup>C and &#x3b4;<sup>15</sup>N (&#x2030;) were determined on subsamples (5 g) using a DELTA V Plus continuous flow isotope ratio mass spectrometer linked to a Flash 2000 elemental analyser using a MAS 200 R autosampler (Thermo-Fisher Scientific, Bremen, Germany) at the NIWA Environmental and Ecological Stable Isotope Facility in Wellington, New Zealand. For sediment metal concentrations, subsamples (50 g) were air dried at 35<sup>&#xb0;</sup>C, acid digested (nitric/hydrochloric acid: US EPA 200.2) and assayed using ICP-MS mass spectrometer.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Modelling historic nutrient load to the New River Estuary</title>
<p>Point sources (e.g., discharges from waste-water treatment plants either directly into the estuary or into tributaries of the estuary) and diffuse sources (run-off from land) both contribute substantial nutrient loads to the NRE.</p>
<sec id="s2_4_1">
<label>2.4.1</label>
<title>Point sources of nutrients</title>
<p>Loads of the four major point sources within the catchment (i.e., the Invercargill City Council Wastewater Treatment Plant (ICC WWTP), and three meat processing plants: Prime Range Meats, Alliance Lorneville, and Alliance Makarewa) were determined for all years with available data by summing daily loads derived from multiplication of total nitrogen (TN) concentration data and flow data. This method is reliant on accurate reporting of flow on all days when discharges occur; however, unreported flows have been observed on some dates between 2012 and 2013 (<xref ref-type="bibr" rid="B78">Robertson and Stevens, 2013a</xref>). Because it is not possible to accurately quantify unreported losses through time, we calculate loads using reported discharges and note that the loads from point sources may be higher than those stated here.</p>
<p>Wastewater loads from the ICC WWTP to NRE in years for which no data were available were calculated by assuming a linear relationship between wastewater TN and TP discharges and human population within the Invercargill urban area. This relationship was calculated based on population and ICC WWTP discharge records between 2004 and 2014. This period spanned the earliest date that daily consent data were available (2004) until the diversion of the Prime Range Meats meat processing plant flow to the ICC WWTP in 2015. Population data were sourced from Statistics New Zealand, and populations for non-census years were calculated by linear interpolations between census years.</p>
<p>Discharges from meat processing within the catchment in years for which no data were available were calculated from relationships between total meat processing plant TN and TP discharges and stock units (number of livestock) within the Southland Region. This relationship was calculated from discharge and stock data between 1999 and 2011. Annual Southland stock data were taken from Statistics New Zealand records.</p>
</sec>
<sec id="s2_4_2">
<label>2.4.2</label>
<title>Diffuse (non-point) sources of nutrients</title>
<p>Diffuse source nitrogen loads to the NRE between 1840 and 2020 were estimated using the Catchment Land Use for Environmental Sustainability (CLUES) model (<xref ref-type="bibr" rid="B27">Elliott et&#xa0;al., 2016</xref>). The CLUES model calculates annual riverine loads of nitrogen (as TN), based on land use within sub-catchments covered by the area of one or more of 19 land use classes representing key rural enterprises, native and exotic forest, and urban areas.</p>
<p>The land use together with animal stock numbers, physical catchment characteristics, rainfall and mean flows were used in a regression model to predict the total loads. Land use input data for the pre-European colonisation and post-M&#x101;ori arrival period for the year 1840 were calculated from the vegetation layers of <xref ref-type="bibr" rid="B49">Leathwick et&#xa0;al. (2012)</xref>, adapted from the work of <xref ref-type="bibr" rid="B48">Leathwick (2001)</xref> and <xref ref-type="bibr" rid="B50">Leathwick et&#xa0;al. (2003)</xref> which predicted potential vegetation cover composition from regressions relating the distributions of major canopy trees to environmental variables. This potential vegetation cover was adjusted to account for clearance of lowland indigenous forest prior to European arrival, using the 1840 grassland extent map of <xref ref-type="bibr" rid="B3">Ausseil et&#xa0;al. (2011)</xref>. Land use for 1935 was derived by digitising the map of grasslands surveyed by <xref ref-type="bibr" rid="B37">Hilgendorf (1935)</xref>. The vegetation cover from the first edition of New Zealand Land Resource Inventory was used as a 1978 Land use map (<xref ref-type="bibr" rid="B46">Landcare Research Ltd., 2018</xref>). Land use maps for the years between 2000 and 2020 were from Land Cover Database maps (LCDB), ranging from LCDB1 to LCDB5 (<xref ref-type="bibr" rid="B45">Landcare Research Ltd., 2013</xref>), intersected with Agribase land use boundaries (<xref ref-type="bibr" rid="B2">AsureQuality, 2008</xref>). Due to the lack of Agribase data for the year 2012, the 2012 CLUES model run was based on the 2008 land use/agribase layer but with updated coefficients using 2012 stock intensity data. Slope data for the catchment was acquired from the Land Environments of New Zealand classification (LENZ, <xref ref-type="bibr" rid="B51">Leathwick et&#xa0;al., 2002</xref>). Generated land use maps from 1840 to 2020 are shown in <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>. Agribase and LENZ were used to split the LCDB classification of grassland into pastoral land uses for different stock types that are characterised by different contaminant yields (e.g., dairy, deer, other animals and lowland intensive, hill county and high-country sheep and beef farming). Other land use classes such as cropland, orchards and perennial crops were further subdivided into specific classes, such as apples, kiwifruit and grapes, or for crops such as onions, potatoes, or maize. The rest of the LCDB classes were mapped directly without subdividing (e.g., native and exotic forests, scrub, tussock, urban areas, rivers and other water bodies).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Derived land use maps of the New River Estuary catchment at seven timepoints from the pre-European (1840) period to 2020.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g002.tif"/>
</fig>
<p>The CLUES model uses two metrics to adjust nutrient loads generated on each land use class: a proportional adjustment for stock density, and an adjustment for land use intensity, both relative to the 2008 baseline used in model development (<xref ref-type="bibr" rid="B27">Elliott et&#xa0;al., 2016</xref>). Stock numbers (used to calculate stock density) were obtained from Agribase (<xref ref-type="bibr" rid="B2">AsureQuality, 2008</xref>) and summed accordingly for combined classes (e.g., sheep and beef). The maps were clipped to the New River catchment boundary and added to the CLUES GIS application directly and run as different scenarios to predict the loads. In the absence of stock density data prior to the year 2000, stock density for &#x2018;sheep and beef&#x2019;, &#x2018;dairy&#x2019;, and &#x2018;deer&#x2019; land use classes were approximated using changes in stock numbers (or stock unit equivalent in the case of sheep and beef) in the Southland region (excluding Fiordland and the islands which have negligible stock), relative to the 2008 baseline. Adjustments for land use intensity for sheep and beef, and deer land classes were made using annual excretal N loads for these stocks taken from the New Zealand Greenhouse Gas Inventory 1990-2019 (<xref ref-type="bibr" rid="B63">Ministry for the Environment, 2021</xref>). Excretal N loads for the years prior to 1990 use the 1990 values. Adjustments for the dairy land class were made using milk solids per cow statistics based on relationships between milk solids output per cow and nutrient inputs per cow to dairy farmland (<xref ref-type="bibr" rid="B58">Luo and Ledgard, 2021</xref>). Area of land use by type, livestock density, and estimated nitrogen load are shown in <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>New River Estuary catchment <bold>(A)</bold> land use, <bold>(B)</bold> livestock density, and <bold>(C)</bold> estimated nitrogen load input.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g003.tif"/>
</fig>
</sec>
<sec id="s2_4_3">
<label>2.4.3</label>
<title>Calibration of CLUES model results</title>
<p>Recent observed nutrient data from the six monitoring sites in rivers draining into the NRE were used to calibrate CLUES model results. These include daily flow and monthly concentration data from Makarewa River at Wallacetown (2007-2021), &#x14c;reti River at Lumsden Bridge (2006-2020), &#x14c;reti River at Three Kings (2007-2020), &#x14c;reti River at Wallacetown (2006-2020), Otepuni Creek at Nith Street (2007-2020), and Waihopai River upstream of Queens Drive (2007-2020). Daily and annual TN concentrations at these locations were calculated by interpolating between monthly values for each site. Daily loads were then calculated by multiplying daily flow, nutrient concentration values and time steps. A linear regression was made between observed annual loads (t y<sup>-1</sup>) and CLUES model results for years with land use maps and observed data. A separate regression relationship was made for sites on the &#x14c;reti and Makarewa Rivers and two other sites at Waihopai River and Otepuni Creek. There was a reasonable agreement between observed and modelled TN with R<sup>2</sup> values of 0.88 and 0.97 for combined TN load from &#x14c;reti and Makarewa Rivers sites and Waihopai and Otepuni sites, respectively. The uncertainty associated with the regression relationships between observed annual TN load and model estimates was assessed using the standard error of the regression. This measure provides insight into the variability or margin of error in the predictions made by the model. The standard error between observed annual TN load and CLUES model estimates were 304 t y<sup>-1</sup> for the &#x14c;reti and Makarewa River reaches, and 31 t yr<sup>-1</sup> for the Waihopai and Otepuni. These high regression values and relatively low error estimates, compared with the average of observed TN values underscore the effectiveness of the model in accurately predicting TN levels in the study area. Coefficients from these regression relationships were then used to calibrate CLUES model results for years without observed data. To calculate annual NRE TN load, calibrated CLUES results for all terminal river reaches which drain into the NRE were summed.</p>
<p>A separate linear regression relationship was developed to estimate TN loads for years between 1840 and 2020 without available land use maps, using dairy stock excreta N load for the Southland region (for the years 2000-2020) and total stock excreta N (for the years 1840-1999) as the predictors. Stock excreta data from Ministry for the Environment, and stock numbers from Statistics New Zealand were used for this calculation (<xref ref-type="bibr" rid="B63">Ministry for the Environment, 2021</xref>).</p>
<p>To investigate the sensitivity of the CLUES model estimates to changes in stock density values, the data for the year 1935 was run in the CLUES model using 0.5 and 2 times the original stock density values, but with the same land use intensities. By reducing the stock density values to 0.5 of the original values, the calibrated TN load was reduced by 8.2% (123 t y<sup>-1</sup>). When doubling the stock density values, the TN load increased by 9.6% (143 t y<sup>-1</sup>). We note that the difference between modelled TN loads caused by these changes in stock density was less than the uncertainty in the CLUES model estimates.</p>
</sec>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Dilution modelling</title>
<p>Total nitrogen concentration in NRE was calculated as &#x2018;potential TN&#x2019;, which is defined as the concentration in the estuary the absence of uptake by algae or changes due to non-conservative processes such as denitrification (<xref ref-type="bibr" rid="B69">Plew et&#xa0;al., 2018</xref>). As such, potential TN represents the nitrogen available for primary production. Past and present estuary potential TN concentrations were estimated using estuary hydrological dilution models developed within ETI Tool 1 (<xref ref-type="bibr" rid="B69">Plew et&#xa0;al., 2018</xref>, <xref ref-type="bibr" rid="B68">2020</xref>), which enables users to assess susceptibility of estuaries to eutrophication based on their nutrient loads and their flushing/dilution characteristics. Potential TN concentrations were calculated based on the relative proportions of riverine and marine water within the estuary, using salinity to tune a modified tidal prism model (<xref ref-type="bibr" rid="B69">Plew et&#xa0;al., 2018</xref>). For NRE, the tidal prism model was tuned using volume-averaged salinity extracted from a calibrated three-dimensional hydrodynamic model (<xref ref-type="bibr" rid="B61">Measures, 2016</xref>). The potential TN concentrations represent a time- (annual) and volume-averaged concentration for the entire estuary.</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>Ecological state modelling</title>
<p>Past and present estuary ecological states were estimated using Tool 3 of the online Estuary Trophic Index (ETI) toolset (<xref ref-type="bibr" rid="B105">Zeldis et&#xa0;al., 2017</xref>). ETI Tool 3 is a Bayesian Belief Network (BBN, <xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>) that uses knowledge of the ecological connections between drivers of estuary trophic condition (e.g., estuary type, nutrient and sediment loads, flushing rate, etc.) and responses of estuary health indicators to calculate scores for the individual indicators, and a combined overall estuary health (ETI) score. Indicators are distinguished as &#x2018;primary&#x2019;, being macroalgae and phytoplankton that have demonstrably increasing responses to nutrient enrichment, and &#x2018;secondary&#x2019;, being ecological attributes that show symptoms of estuary health impairment (e.g., high sediment carbon, reduced sediment oxygen, impacted seagrass and macrobenthic communities, in response to nutrient loads and sedimentation. The ETI BBN combines the primary and secondary indicator scores into the final ETI score (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>), where the score of the most impaired primary indicator is combined with the average of the secondary indicator scores and is discretised into 10 bands (<xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>). &#x2018;Intermediate calculation&#x2019; nodes do not directly input to the scores but are used to calculate some of the indicators (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Schematic of ETI Tool 3 BBN. Information on &#x2018;drivers&#x2019; (blue nodes) are input by the user and the BBN calculates values of primary and secondary trophic indicators (yellow and pink nodes). Primary and secondary indicator values are used to produce ETI primary and secondary indicator scores, which are then combined to give the final ETI &#x2018;performance&#x2019; score and band (red nodes). Intermediate calculation nodes (grey nodes) are used to estimate linked indicator values. Sediment accumulation rate, % mud, and % TOC were set using the sediment core results (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). Because these are the downstream &#x2018;children&#x2019; of the sediment load driver node (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>), that node was not set. All other driver nodes were set (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>), including hindcasted potential TN concentrations obtained from ETI Tool 1 (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). Adapted from <xref ref-type="bibr" rid="B104">Zeldis and Plew (2022)</xref> under <uri xlink:href="https://creativecommons.org/licenses/by/4.0">CC-BY 4.0</uri>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g004.tif"/>
</fig>
<p>In NZ water quality policy statements (e.g., the National Policy Statement for Freshwater Management; NPS-FM; <xref ref-type="bibr" rid="B64">Ministry for the Environment, 2023</xref>), freshwater quality health variables are scored using a four-level banding (A-D). Therefore, the ETI score can be recast as four bands (A-D) if needed, to make comparison with the 4-state banding used within the NPS-FM. The estuary health conditions associated with the 4-state bandings of the ETI final score are described in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>.</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Estuary health conditions associated with the bandings (A &#x2013; D) and ETI scores 0 to 1 of the ETI Tool 3 BBN (<xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>).</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="center">Band A</th>
<th valign="top" align="center">Band B</th>
<th valign="top" align="center">Band C</th>
<th valign="top" align="center">Band D</th>
</tr>
<tr>
<th valign="top" align="center">0 &#x2013; 0.25</th>
<th valign="top" align="center">&gt; 0.25 &#x2013; 0.5</th>
<th valign="top" align="center">&gt; 0.5 &#x2013; 0.75</th>
<th valign="top" align="center">&gt; 0.75 &#x2013; 1</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Ecological communities are healthy and resilient.<break/>Primary producers: dominated by seagrasses and microalgae or in coastal lakes, macrophytes.<break/>Water column: high clarity, well-oxygenated.<break/>Sediment: well oxygenated, low organic matter, low sulphides and ammonium, diverse macrobenthic community with low abundance of enrichment tolerant species.</td>
<td valign="top" align="left">Ecological communities are slightly impacted by additional algal growth arising from nutrient levels that are elevated.<break/>Primary producers: seagrass/microalgae/macrophytes still present but increasing biomass of opportunistic macroalgae.<break/>Water column: moderate clarity, can have moderate - poor DO especially at depth.<break/>Sediment: moderate oxygenation, organic matter, and sulphides, diverse macrobenthic community with increasing abundance of enrichment tolerant species.</td>
<td valign="top" align="left">Ecological communities are highly impacted by macroalgal and/or phytoplankton biomass elevated well above natural conditions. Conditions likely to affect habitat available for seagrass/macrophytes.<break/>Primary producers: opportunistic macroalgal biomass high, seagrass cover low. Increasing phytoplankton where residence time long.<break/>Water column: reduced water clarity may affect deep seagrass beds. Low DO, especially at depth.<break/>Sediment: poor oxygenation, high organic matter and sulphides, macrobenthos dominated by high abundance of enrichment tolerant species.</td>
<td valign="top" align="left">Ecological communities extremely impacted. Excessive macroalgal and/or phytoplankton growth making ecological communities at high risk of undergoing a regime shift to a persistent, degraded state without macrophyte/seagrass cover.<break/>Excessive macroalgal growth making ecological communities at high risk of nuisance algal bloom situation. Phytoplankton community may be dominated by nuisance types (e.g., cyanobacteria, picoplankton).<break/>Water column: low clarity, deoxygenated at depth.<break/>Sediment: anoxic, very high organic matter and sulphides, subsurface macrobenthos very limited or absent.</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In default applications of the ETI Tool 3 BBN, all inputs are provided solely by values of the driver nodes (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). However, the BBN can also be informed by inputting values of indicator or intermediate calculation nodes, if these are known <italic>a priori</italic>, to improve its predictions (<xref ref-type="bibr" rid="B92">Uusitalo, 2007</xref>; <xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>). Accordingly, in this application, intermediate calculation node values for sediment accumulation rate and % mud, and secondary node values for sediment % TOC (all derived from the sediment core results: <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>), were input <italic>a priori</italic> to the BBN model runs. Because these are the downstream &#x2018;children&#x2019; of the sediment load driver node (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>), we did not set the value of the sediment load node in the BBN model runs. All other driver nodes were input, including hindcasted potential TN concentrations obtained from ETI Tool 1 dilution modelling (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>).</p>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Values used to set potential TN (mg m<sup>-3</sup>), % mud, sediment accumulation rate (SAR) (mm y<sup>-1</sup>) and % TOC nodes in the BBN, for years which had values for % mud and/or % TOC derived from the Central and Upper cores.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" colspan="5" align="left">Central core</th>
<th valign="middle" colspan="5" align="left">Upper core</th>
</tr>
<tr>
<th valign="middle" align="center">Year</th>
<th valign="middle" align="center">Potential TN</th>
<th valign="middle" align="center">% mud</th>
<th valign="middle" align="center">SAR</th>
<th valign="middle" align="center">% TOC</th>
<th valign="middle" align="left">Year</th>
<th valign="middle" align="center">Potential TN</th>
<th valign="middle" align="center">% mud</th>
<th valign="middle" align="center">SAR</th>
<th valign="middle" align="center">% TOC</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">2019</td>
<td valign="middle" align="center">663</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">3.1</td>
<td valign="middle" align="left">2019</td>
<td valign="middle" align="center">663</td>
<td valign="middle" align="center">88</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">3.0</td>
</tr>
<tr>
<td valign="middle" align="left">2017</td>
<td valign="middle" align="center">669</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.5</td>
<td valign="middle" align="left">2016</td>
<td valign="middle" align="center">656</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">2.8</td>
</tr>
<tr>
<td valign="middle" align="left">2016</td>
<td valign="middle" align="center">656</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.6</td>
<td valign="middle" align="left">2013</td>
<td valign="middle" align="center">596</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">2.8</td>
</tr>
<tr>
<td valign="middle" align="left">2013</td>
<td valign="middle" align="center">596</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.5</td>
<td valign="middle" align="left">2009</td>
<td valign="middle" align="center">551</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">*</td>
</tr>
<tr>
<td valign="middle" align="left">2011</td>
<td valign="middle" align="center">595</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.5</td>
<td valign="middle" align="left">2007</td>
<td valign="middle" align="center">476</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">2.4</td>
</tr>
<tr>
<td valign="middle" align="left">2010</td>
<td valign="middle" align="center">572</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.5</td>
<td valign="middle" align="left">2004</td>
<td valign="middle" align="center">441</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">2.3</td>
</tr>
<tr>
<td valign="middle" align="left">2008</td>
<td valign="middle" align="center">662</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="left">2001</td>
<td valign="middle" align="center">414</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">2.1</td>
</tr>
<tr>
<td valign="middle" align="left">2006</td>
<td valign="middle" align="center">457</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">30.8</td>
<td valign="middle" align="center">2.1</td>
<td valign="middle" align="left">1998</td>
<td valign="middle" align="center">414</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">2.1</td>
</tr>
<tr>
<td valign="middle" align="left">1998</td>
<td valign="middle" align="center">414</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">2.0</td>
<td valign="middle" align="left">1995</td>
<td valign="middle" align="center">408</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">1.9</td>
</tr>
<tr>
<td valign="middle" align="left">1993</td>
<td valign="middle" align="center">404</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">1.9</td>
<td valign="middle" align="left">1989</td>
<td valign="middle" align="center">395</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">1.5</td>
</tr>
<tr>
<td valign="middle" align="left">1985</td>
<td valign="middle" align="center">385</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">1.3</td>
<td valign="middle" align="left">1983</td>
<td valign="middle" align="center">381</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">16.7</td>
<td valign="middle" align="center">1.5</td>
</tr>
<tr>
<td valign="middle" align="left">1980</td>
<td valign="middle" align="center">374</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">1.8</td>
<td valign="middle" align="left">1966</td>
<td valign="middle" align="center">301</td>
<td valign="middle" align="center">43</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="center">*</td>
</tr>
<tr>
<td valign="middle" align="left">1966</td>
<td valign="middle" align="center">301</td>
<td valign="middle" align="center">61</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">1.6</td>
<td valign="middle" align="left">1934</td>
<td valign="middle" align="center">220</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="center">0.8</td>
</tr>
<tr>
<td valign="middle" align="left">1953</td>
<td valign="middle" align="center">245</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">1.3</td>
<td valign="middle" align="left">1925</td>
<td valign="middle" align="center">210</td>
<td valign="middle" align="center">13</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="center">*</td>
</tr>
<tr>
<td valign="middle" align="left">1943</td>
<td valign="middle" align="center">234</td>
<td valign="middle" align="center">30</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">0.7</td>
<td valign="middle" align="left">1892</td>
<td valign="middle" align="center">188</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="center">0.8</td>
</tr>
<tr>
<td valign="middle" align="left">1929</td>
<td valign="middle" align="center">215</td>
<td valign="middle" align="center">100</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="left">1880</td>
<td valign="middle" align="center">182</td>
<td valign="middle" align="center">8</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="center">*</td>
</tr>
<tr>
<td valign="middle" align="left">1927</td>
<td valign="middle" align="center">212</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">0.7</td>
<td valign="middle" align="left">1847</td>
<td valign="middle" align="center">166</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">2.2</td>
<td valign="middle" align="center">0.3</td>
</tr>
<tr>
<td valign="middle" align="left">1916</td>
<td valign="middle" align="center">201</td>
<td valign="middle" align="center">33</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">0.6</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
</tr>
<tr>
<td valign="middle" align="left">1900</td>
<td valign="middle" align="center">193</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">0.4</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
</tr>
<tr>
<td valign="middle" align="left">1874</td>
<td valign="middle" align="center">180</td>
<td valign="middle" align="center">*</td>
<td valign="middle" align="center">3.8</td>
<td valign="middle" align="center">0.2</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Potential TN values were hindcasted using methods of Section 2.5. Asterisks indicate missing values.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T4" position="float">
<label>Table&#xa0;4</label>
<caption>
<p>Values used to set driver nodes in the BBN, for all years and both cores.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="bottom" align="center">Estuary Type</th>
<th valign="bottom" align="center">Salinity</th>
<th valign="bottom" align="center">Intertidal area (%)</th>
<th valign="bottom" align="center">Flushing time (d)</th>
<th valign="bottom" align="center">Seasonality</th>
<th valign="bottom" align="center">Stratification</th>
<th valign="bottom" align="center">Closure</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="bottom" align="center">Tidal lagoon</td>
<td valign="bottom" align="center">26.1</td>
<td valign="bottom" align="center">67%</td>
<td valign="bottom" align="center">3.9</td>
<td valign="bottom" align="center">0.695</td>
<td valign="bottom" align="center">No</td>
<td valign="bottom" align="center">Open</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>These values were available from ETI Tool 1 output.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>We provide resulting scores for four eutrophication and sedimentation associated indicators: Macroalgal ecological quality rating (EQR) is provided with thresholds of 0.8, 0.6 and 0.4 between A&#x2013;B (minimal-moderate), B&#x2013;C (moderate-high) and C&#x2013;D (high&#x2013;very high) bands of macroalgal eutrophication respectively (<xref ref-type="bibr" rid="B76">Robertson et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B68">Plew et&#xa0;al., 2020</xref>); Apparent Redox Potential Discontinuity (aRPD) depth, the depth of the boundary between oxic near-surface sediment and the underlying suboxic or anoxic sediment, where shallower scores are related to reduced benthic substrate oxygen penetration and habitat quality; Seagrass health driven by SAR and macroalgal growth (EQR) where larger scores indicate greater deleterious impacts; and a macrobenthos NZ AZTI Marine Benthic Index (AMBI, <xref ref-type="bibr" rid="B12">Borja et&#xa0;al., 2000</xref>) ranging from &#x2018;Normal&#x2019; (band B: 1.2 &#x2013; 3.3) to &#x2018;Transitional to pollution&#x2019; (band C: 3.3 &#x2013; 4.3), to &#x2018;Polluted&#x2019; (band D: &gt; 4.3).</p>
<p>The relationships between nodes of the ETI Tool 3 BBN are detailed in <xref ref-type="bibr" rid="B104">Zeldis and Plew (2022)</xref> and its supplementary information file. These relationships and modifications to the BBN since publication of <xref ref-type="bibr" rid="B104">Zeldis and Plew (2022)</xref> are provided in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Datasheet 1</bold>
</xref>, <xref ref-type="supplementary-material" rid="SM1">
<bold>Section 1</bold>
</xref> with modifications also in the &#x2018;Change Log&#x2019; file available in the online ETI Tool 3.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Sediment accumulation rate</title>
<p>Sediments at the Upper site were dated back to 1847 and at the Central site to 1874. At the Upper site, the SAR from 1847 (91 cm below surface) to 1983 (60 cm below surface) was 2.4 mm y<sup>-1</sup> and from 1983 to the present day it was 16.7 mm y<sup>-1</sup> (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Datasheet 1</bold>
</xref>, <xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S1</bold>
</xref>). At the Central site, SAR from 1874 (91 cm below surface) to 2006 (40 cm below surface) was 3.8 mm y<sup>-1</sup> and from 2006 to the present-day SAR was 33.6 mm y<sup>-1</sup> (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Datasheet 1</bold>
</xref>, <xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S1</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Log-linear regression graphical interpretations of unsupported <sup>210</sup>Pb profiles with sediment core depth for <bold>(A)</bold> the Upper (orange) and <bold>(B)</bold> the Central (blue) sites, and <bold>(C)</bold> calculated sediment accumulation from the modelled log-liner regression fitted ages at the Upper (from 1847) and the Central (from 1874) locations in the New River Estuary.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g005.tif"/>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Sediment physical and biogeochemical properties</title>
<p>Sediment organic matter and C content increased from the oldest dated layers to the present. Sediment OM ranged from 1.1% to 6.7% at the Upper site and from 1.4% to 8.2% at the Central site (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6A</bold>
</xref>) and % C ranged from 0.32% to 2.95% at the Upper site and from 0.16% to 3.12% at the Central site (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6B</bold>
</xref>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Sediment core physical and biogeochemical properties at the Upper (orange) and Central (blue) sites in the New River Estuary: <bold>(A)</bold> organic matter content (%), <bold>(B)</bold> carbon (C) content (%), <bold>(C)</bold> &#x3b4;<sup>13</sup>C (&#x2030;), <bold>(D)</bold> mud content (%), <bold>(E)</bold> nitrogen (N) content, <bold>(F)</bold> &#x3b4;<sup>15</sup>N (&#x2030;), <bold>(G)</bold> C:N molar ratio, <bold>(H)</bold> chlorophyll-<italic>a</italic> content (&#xb5;g g-1), <bold>(I)</bold> organic matter preservation index (chlorophyll-<italic>a</italic>:phaeophytin-<italic>a</italic>).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g006.tif"/>
</fig>
<p>&#x3b4;<sup>13</sup>C values, corrected for the Suess effect (the decrease in atmospheric &#x3b4;<sup>13</sup>C caused by CO<sub>2</sub> emissions of fossil fuel depleted in <sup>13</sup>C; <xref ref-type="bibr" rid="B93">Verburg, 2007</xref>), decreased from -24.47&#x2030; to -26.64&#x2030; in the Upper core between 1847-1983 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6C</bold>
</xref>). In the Central core there was little variation (range from -25.27&#x2030; to -26.16&#x2030;) and no clear trend between 1874 and 1985 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6C</bold>
</xref>). However, sediment laid down in recent decades (i.e., post-1985) showed a steady increase in &#x3b4;<sup>13</sup>C values at both locations up to -23.67&#x2030; at the Upper site and -23.38&#x2030; at the Central site. The conversion of the &#x3b4;<sup>13</sup>C data to correct for the Suess effect was crucial: without it &#x3b4;<sup>13</sup>C in the first two decades of this century was lower than in the early 20th century, suggesting a decrease in dissolved carbon demand and lower primary productivity this century compared with a century ago. This would not be consistent with the multi-pronged evidence for eutrophication from other sediment characteristics.</p>
<p>Mud content had steep increases in both the Upper and Central cores from the older sediments (1880 and 1916 at the Upper and Central sites respectively) to the present (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6D</bold>
</xref>, <xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). While the increase in mud content at the Central site was somewhat variable prior to the 1980&#x2019;s, its subsequent persistent 100% values (at both core sites) indicated formation of a very muddy environment in recent decades.</p>
<p>Sediment core N content increased from 0.03% to 0.34% at the Upper site and from 0.01% to 0.36% at the Central site (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6E</bold>
</xref>) from oldest dated layers sampled to the present day. From ~ 1850 to 1990, &#x3b4;<sup>15</sup>N values were within the range of marine-derived N sources in southern New Zealand of ~ 5-6&#x2030; (<xref ref-type="bibr" rid="B19">Cornelisen et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B26">Durante et&#xa0;al., 2021</xref>) and slightly below the 6.6 &#x2013; 8.8&#x2030; range of &#x3b4;<sup>15</sup>N values given for macroalgal (<italic>Ulva</italic> spp) tissue collected from unimpacted, exposed coastal sites around New Zealand <xref ref-type="bibr" rid="B5">Barr et&#xa0;al., 2013</xref>; <xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6F</bold>
</xref>). Increased &#x3b4;<sup>15</sup>N values, above the marine derived 5-6&#x2030; range, are apparent from ~ 1990 and both cores show particularly sharp increases in &#x3b4;<sup>15</sup>N values in sediments deposited since the early 2000&#x2019;s. C:N ratios (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6G</bold>
</xref>) largely followed inverse trends to those seen in &#x3b4;<sup>13</sup>C and &#x3b4;<sup>15</sup>N values, with steady decreases from ~11 to 9 from the early 2000&#x2019;s in both cores.</p>
<p>Sediment chlorophyll-<italic>a</italic> concentrations were relatively stable in both cores over the 20th century (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6H</bold>
</xref>), but from ~2000 onwards, sharp increases occurred at both core sites. The preservation index (chlorophyll-<italic>a</italic>:phaeophytin-<italic>a</italic> ratio) was relatively invariant until the early 2000&#x2019;s, when it underwent sharp increases at both sites (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6I</bold>
</xref>).</p>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Total N load and potential TN time series</title>
<p>CLUES-modelled TN loads to the estuary from freshwater and point sources increased from 947 t y<sup>-1</sup> in 1840, up to a maximum of 5715 t y<sup>-1</sup> in 2007 (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>). Hindcasted estuary potential TN concentration (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7B</bold>
</xref>) derived from dilution modelling increased ~4-fold from 1840 to the present, consistent with the increasing total load to the estuary. During the more recent period of rapid increase from 2000-2005 to 2015-2019 (5-year averages to reduce interannual variability), potential TN increased from 449 to 713 mg m<sup>-3</sup>. There were very high correlations between the potential TN and N in the sediments at corresponding times in their respective time series (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7C</bold>
</xref>: Upper core R<sup>2</sup> = 0.949, TN = 110 + 1751 (% N); Central core R<sup>2 </sup>= 0.926, TN = 112 + 1916 (% N)).</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>
<bold>(A)</bold> CLUES modelled (white circles) and interpolated values for TN load (t y<sup>-1</sup>), <bold>(B)</bold> potential N (mg m<sup>-3</sup>) in the New River Estuary from 1847 to 2019 at the Upper (orange) and Central (blue sites), and the <bold>(C)</bold> relationship between sediment core nitrogen content (%) in the dated sediment core layers and the calculated potential N (mg m<sup>-3</sup>). Upper core R<sup>2</sup> = 0.949 (y = 1751.40x - 110.20), Central core R<sup>2</sup> = 0.926 (y = 1916.36x - 112.00).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g007.tif"/>
</fig>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Ecological state modelling</title>
<p>NRE ecological state over the time series was described by BBN indicator and ETI scores (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>) that resulted when the BBN was informed by the driver settings and observed sediment core values in <xref ref-type="table" rid="T3">
<bold>Tables&#xa0;3</bold>
</xref>, <xref ref-type="table" rid="T4">
<bold>4</bold>
</xref>. Hindcasted macroalgal EQR was stable at band A (healthy: greater than or near 0.8) for much of the time series prior to the 1950&#x2019;s (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8A</bold>
</xref>). It worsened through band B (0.6 to &lt; 0.8) to band C (0.4 to &lt; 0.6) by the late 1980&#x2019;s, coincident with gradually increasing potential TN (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>), but then declined more steeply starting in the late 1990&#x2019;s in response to sharply increasing potential TN, reaching band D (&lt; 0.4) by ~2010.</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Responses of trophic indicator scores <bold>(A)</bold> macroalgae ecological quality rating (EQR), <bold>(B)</bold> sediment aRPD, <bold>(C)</bold> seagrass impact, <bold>(D)</bold> macrobenthic NZ AMBI, and <bold>(E)</bold> ETI score, with associated bandings (for <bold>(A&#x2013;D)</bold> <italic>c.f.</italic> <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1374869-g008.tif"/>
</fig>
<p>BBN-predicted aRPD (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8B</bold>
</xref>) was relatively stable prior to the 1930&#x2019;s at band A (&gt; 4 cm depth), after which it started shallowing initially in response to increasing % TOC (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6B</bold>
</xref>; <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>) and then declining macroalgal EQR (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8A</bold>
</xref>) after the 1950&#x2019;s. It reached band D (&lt; 1 cm depth) in the early 2000&#x2019;s, driven both by increasing % TOC and steeply declining macroalgae EQR. Overall, the shallowing aRPD indicated that the depth of the boundary between oxic near-surface sediment and the underlying suboxic/anoxic sediment was likely to have shallowed to very near the surface of the estuary bed.</p>
<p>Seagrass decline (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8C</bold>
</xref>) was predicted to be minor to moderate (bands A-B: scores &lt; 0.5) until ~1983 for the Upper core area and ~ 2006 for the Central core area, after which they underwent severe decline to band D (&gt; 0.8) at both sites. These dates corresponded with inflection points in SAR at the respective sites (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5C</bold>
</xref>), although the transition was less well defined for the Upper core because of the lack of SAR results between 1966 and 1983 (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). For both cores, D-band status corresponded with steep declines in macroalgal EQR (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8A</bold>
</xref>).</p>
<p>Macrobenthic NZ AMBI index (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8D</bold>
</xref>) indicated &#x2018;normal&#x2019; (&#x2264; 3.3) status until 1966 and 1953 for the Upper and Central cores respectively, after which they assumed &#x2018;transitional to pollution&#x2019; (&#x2264; 4.3) status, eventually reaching &#x2018;polluted&#x2019; status (&gt; 4.3) in 1998 at both sites. Their worsening scores occurred in response to increasing % TOC and % mud (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6B, D</bold>
</xref>) and finally, in response to steeply worsening macroalgae EQR (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8A</bold>
</xref>).</p>
<p>Overall NRE estuary health, indexed by the ETI score (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8E</bold>
</xref>), was in band A (&#x2264; 0.25: <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>) until the early 20<sup>th</sup> century. It degraded slowly through band B (&gt; 0.25 to &#x2264; 0.50) from the 1930&#x2019;s to the 1980&#x2019;s, but then worsened rapidly through band C (&gt; 0.50 to &#x2264; 0.75), reaching band D (&gt; 0.75) by 2007. The degradation arose as an interaction of steadily increasing SAR, % mud, % TOC and potential TN, but its acceleration in the late 1990&#x2019;s was associated with rapidly worsening macroalgae EQR state, driven by the sudden increase of potential TN.</p>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>In this study, we reconstructed the historic trophic condition of the NRE. We used historical land use information to calculate nutrient loading to the estuary and determine estuary-wide time series of water column TN concentrations. We compared TN concentration time series to biogeochemical indices in dated sediment core layers at two locations in the estuary. Both TN concentration model results and sediment core chemistry indicated catchment and estuary N loads increased greatly from the year 2000 to the present day, and further modelling using a BBN (<xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>) provided a time series of ecological state values. Our results indicate that the currently eutrophic NRE (<xref ref-type="bibr" rid="B80">Robertson et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B29">Forrest et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B72">Roberts et&#xa0;al., 2022a</xref>), whilst historically muddy, was not necessarily eutrophic in the past with a much lower SAR. While overall, the NRE showed increasing degradation, it showed differing responses to differing land use profiles through time (i.e., sheep and beef, then dairy), firstly with increased muddiness, and then macroalgal eutrophication.</p>
<sec id="s4_1">
<label>4.1</label>
<title>Temporal trajectories of sediment and nutrient loading</title>
<p>The sediment core results from the upper NRE showed clear sign of increased sedimentation with increased muddiness and % TOC commencing ca. early-20<sup>th</sup> century. SAR increased steadily at relatively low rates until inflection points appeared in the late 20<sup>th</sup> and early 21<sup>st</sup> centuries, with sharply increased rates thereafter. These rates concur with sedimentation estimates from sediment cores collected in 2007 by <xref ref-type="bibr" rid="B77">Robertson and Stevens (2007)</xref> and in 2017 by <xref ref-type="bibr" rid="B13">Brown (2019)</xref> in depositional zones in the NRE. Sedimentation rate in the Upper Waihopai Arm was estimated to have increased from 3 mm y<sup>-1</sup> from 1906 to 1967 to 28 mm y<sup>-1</sup> from 2001 to 2007 by <xref ref-type="bibr" rid="B77">Robertson and Stevens (2007)</xref>, from 13.3 mm y<sup>-1</sup> from 1915 to 1935 to 22.4 mm y<sup>-1</sup> from 2009 to 2017 by <xref ref-type="bibr" rid="B13">Brown (2019)</xref>, and from 5.9 mm y<sup>-1</sup> from 1956 to 2017 to 17.5 mm y<sup>-1</sup> from 2007 to 2017 in the Lower Waihopai Arm by <xref ref-type="bibr" rid="B13">Brown (2019)</xref>. Sediment plate monitoring conducted in in Lower Waihopai Arm estimated sedimentation rates of 30-40 mm y<sup>-1</sup> (including both fine sediment an macroalgae entrained within the sediments) from 2007 to 2020 (<xref ref-type="bibr" rid="B74">Roberts and Ward, 2021</xref>) and 18.1-23.2 mm y<sup>-1</sup> from 2007 to 2022 (<xref ref-type="bibr" rid="B29">Forrest et&#xa0;al., 2022</xref>).</p>
<p>Increases in several core biogeochemical indicators (% OM, % C and % N) since the 1800&#x2019;s, and at a more rapid rate from mid-20<sup>th</sup> century, indicated gradually accelerating sedimentation and eutrophication in the upper NRE as catchment stock numbers and human population increased (<xref ref-type="bibr" rid="B2">AsureQuality, 2008</xref>; <xref ref-type="bibr" rid="B31">Grant, 2023</xref>). Additionally, and in contrast to % OM, % C, and % N, between about 1980 and the early 21<sup>st</sup> century, core sediments showed a sudden switch from little to no change in sediment &#x3b4;<sup>13</sup>C, &#x3b4;<sup>15</sup>N, chlorophyll-<italic>a</italic>, and C:N ratios since the 1800&#x2019;s to steeply increasing sediment &#x3b4;<sup>13</sup>C and &#x3b4;<sup>15</sup>N values, large increases in chlorophyll-<italic>a</italic> content, and decreasing C:N ratios. C:N ratio decreases in organic material in mixohaline systems are indicative of increased contributions of marine-derived algal material relative to contributions of vascular plant material, and/or an increase in N availability to primary producers (<xref ref-type="bibr" rid="B17">Cloern et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B6">Barr and Rees, 2003</xref>; <xref ref-type="bibr" rid="B83">Savage et&#xa0;al., 2010</xref>). These changes taken together indicate increases in autochthonous algal primary production rates and deposition, as N loading accelerated during recent agricultural intensification in the catchment (<xref ref-type="bibr" rid="B15">Cifuentes et&#xa0;al., 1988</xref>; <xref ref-type="bibr" rid="B83">Savage et&#xa0;al., 2010</xref>). Sediment indicators of recent eutrophication mirrored the evidence seen in regular ecological monitoring of the intertidal zone of NRE since 2001, showing increases in algal biomass and changes in faunal communities (<xref ref-type="bibr" rid="B79">Robertson and Stevens, 2013b</xref>; <xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>).</p>
<p>In addition to increases in N availability from the estuary catchment, increases in sediment &#x3b4;<sup>15</sup>N in the mid-late 20<sup>th</sup> century likely indicate contributions of N from treated sewage. Hydrodynamic modelling of the NRE suggests that N from the Invercargill WWTP main municipal wastewater discharge (in its modern location at Clifton) likely forms a sizable proportion of total water column N at the locations of both cores (<xref ref-type="bibr" rid="B61">Measures, 2016</xref>). Given this, it is possible that some of the increases in sediment &#x3b4;<sup>15</sup>N resulted from the commencement of effluent discharge from the WWTP (in 1969), construction of a secondary effluent treatment processing facilities (in 1992), and addition of maturation ponds and wetlands (in 2004) that likely increased N removal from effluent. Typically, the nitrification-denitrification processes that cause gaseous losses of N during secondary and tertiary effluent processing result in &#x3b4;<sup>15</sup>N-enrichment of the remaining effluent-N (<xref ref-type="bibr" rid="B82">Savage, 2005</xref>; <xref ref-type="bibr" rid="B1">Alldred et&#xa0;al., 2023</xref>). Also, a range of N-cycling processes within catchments exhibit concentration-dependent expression of discrimination against &#x3b4;<sup>15</sup>N, leading to progressively greater &#x3b4;<sup>15</sup>N-enrichment of inorganic N pools passing to estuaries and coastal waters, as catchment land use (and N application) intensifies (<xref ref-type="bibr" rid="B8">Bedard-Haughn et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B67">Needoba et&#xa0;al., 2004</xref>; <xref ref-type="bibr" rid="B23">Dudley et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B44">Kritee et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B97">Wells and Eyre, 2019</xref>). As a result, sediments and primary producers in eutrophic, N-enriched estuaries and coastal environments influenced by rivers tend to have higher &#x3b4;<sup>15</sup>N values (<xref ref-type="bibr" rid="B60">McClelland et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B95">Voss et&#xa0;al., 2000</xref>).</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Impact of land use on nutrient and sediment loads</title>
<p>Increased nutrient and sediment loading to the NRE since the early-to-mid-20<sup>th</sup> century resulted largely from changing rural land use and land use practices. Human populations within the Invercargill urban area have increased only marginally between the early 1970&#x2019;s (50,681 in 1971) and recent times (54,204 at the 2018 census). We estimate that N releases via urban wastewater discharge to the estuary also increased little over this time. In contrast, our results indicate that inputs of N from agricultural land within the NRE catchment approximately doubled over the same period. This increase in N loss from the catchment mostly resulted from agricultural intensification, particularly via increases in stock numbers per unit area of farmland, and conversion of land used for less intensive agriculture such as sheep and beef farming, to more intensive use, particularly dairy farming. Note that while the area occupied by dairy farming in the 2010&#x2019;s was similar to or less than that in 1935 and 1978, there were far more dairy cows after 2010. These patterns of agricultural intensification seen in Southland since the mid-20<sup>th</sup> century mirror those seen in other areas of New Zealand (<xref ref-type="bibr" rid="B99">Wynyard, 2016</xref>; <xref ref-type="bibr" rid="B64">Ministry for the Environment, 2023</xref>).</p>
<p>Land use intensification can also lead to an increase in sediment generation from the catchment, a trend that has been observed in various catchments across New Zealand (<xref ref-type="bibr" rid="B33">Haddadchi and Hicks, 2020</xref>). Changes in hydrology as a result of land reclamation in the Waihopai Arm of the NRE in the early 1900s have been suggested to be a cause of increased sedimentation in this region, however we find that sedimentation following the reclamations remains low compared to recent SAR and that changes in land use are the likely cause of increased sedimentation observed since the 1980s. SAR increases in more recent times likely indicate a positive feedback often seen between fine sediment and nutrient loading increases, with muddy sediments and high nutrient loading favouring nutrient retention and macroalgal outgrowth (<xref ref-type="bibr" rid="B75">Robertson and Savage, 2018</xref>) and with macroalgae, in turn, trapping more mud (<xref ref-type="bibr" rid="B103">Zeldis et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>).</p>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Using sediment cores to determine historic N loads and potential N concentration</title>
<p>We found very high correlations between the time series of estuary-wide potential TN in the water column (determined from N load and dilution modelling) and sediment N content in the dated sediment core strata. This suggested that historic potential TN water column concentrations are predictable from sediment core data. However, the potential TN concentrations represent an average value over the whole estuary, while the cores were taken from the impacted depositional zone in the Lower Waihopai Arm of the estuary. Spatially-distributed nutrient sampling (<xref ref-type="bibr" rid="B22">Dudley et&#xa0;al., 2022</xref>) in NRE and hydrodynamic modelling using tracers as a proxy for nutrients (<xref ref-type="bibr" rid="B61">Measures, 2016</xref>) show that the inner parts of the estuary, such as the Waihopai Arm, have higher nutrient concentrations than &#x2018;estuary-averaged&#x2019; values. Additionally, the best/high-temporal resolution sedimentary records come from high-SAR environments. Sediment mixing by infauna is less effective as SAR increases so that physical processes dominate. This suggests that a historic time series of &#x2018;local&#x2019; potential TN predicted for the inner estuary and other locations with different SAR rates would have a different relation with sediment N content, likely with a higher slope, than we show here for &#x2018;estuary averaged&#x2019; values. While recognising this caveat, the fact that estuary-wide historic potential TN appears predictable from profiles of sediment core N content in NRE indicates that similar coring work on other suitable (depositional) sites in NRE and other estuaries would be beneficial, to determine if the strong relationship found here between modelled TN loads and sediment core N content is a general phenomenon.</p>
</sec>
<sec id="s4_4">
<label>4.4</label>
<title>BBN hindcasts of NRE ecological health</title>
<p>The BBN resolved a complex interaction of effects, operating at different times and time scales, within the overall pattern of degradation. The ETI BBN showed that NRE ecological health has deteriorated across all the ecological indicators calculated, particularly from the mid-20<sup>th</sup> century to the present day. Macroalgal EQR deteriorated mostly from the late 1980&#x2019;s (through the C band to D band) at both core sites, associated with the strong increase in potential TN occurring at that time (exceeding 400 mg m<sup>-3</sup> TN). The EQR results from both sites showed identical trajectories of change, due to the estuary-wide value of potential TN driving the EQR node. Monitoring of the area of intertidal flats of the estuary showed that estuary area covered by &gt; 50% macroalgae increased from 43 ha in 2001 to 364 ha in 2016 (<xref ref-type="bibr" rid="B80">Robertson et&#xa0;al., 2017</xref>), to 479.7 ha in 2022 (<xref ref-type="bibr" rid="B72">Roberts et&#xa0;al., 2022a</xref>), and a decrease in calculated EQR score for the Waihopai Arm from 0.418 to 0.308 from 2012 to 2022 (<xref ref-type="bibr" rid="B72">Roberts et&#xa0;al., 2022a</xref>). This was consistent with, although slightly lower than, our calculated decrease in macroalgal EQR from 0.30 to 0.21 over the same period, following the rapid decrease in EQR after 2008. Monitoring showed that the EQR decline was driven primarily by expansion in substrate coverage of the red alga species <italic>Agarophyton</italic> spp. (formerly <italic>Gracilaria</italic> spp.), and to a lesser extent by <italic>Ulva</italic> spp., which were limited in extent in the NRE before 2000 (<xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>). These macroalgae are typically associated with increases in sediment mud content, loss of indigenous seagrass, decreased sediment oxygenation, and impacted macrobenthos (<xref ref-type="bibr" rid="B77">Robertson and Stevens, 2007</xref>; <xref ref-type="bibr" rid="B79">Robertson and Stevens, 2013b</xref>; <xref ref-type="bibr" rid="B80">Robertson et&#xa0;al., 2017</xref>). Since 2007, <italic>Agarophyton</italic> spp. has become dominant in the depositional areas of the upper NRE (<xref ref-type="bibr" rid="B29">Forrest et&#xa0;al., 2022</xref>), underpinning the declines in ecological conditions seen in the sediment core and BBN indicators retrieved from that location.</p>
<p>The modelled BBN aRPD showed a steady decrease from a &#x2018;low impact&#x2019; value of 4.4 cm (A band) to 0.5 cm (D band), below the &#x2018;exhaustion threshold&#x2019; levels of <xref ref-type="bibr" rid="B89">Sutula et&#xa0;al. (2014)</xref>. Measurements of redox potential in the western Waihopai Arm from 2013 show an aRPD of &lt; 1 cm (<xref ref-type="bibr" rid="B80">Robertson et&#xa0;al., 2017</xref>), consistent with the present day values yielded by the BBN. Further measurements from 2016 showed the sediment was anoxic at 3 cm depth (redox potential &lt; 150mV, the depth at which the primary indicator measurement for assessment is taken) in the Waihopai arm, with worsening of the general condition of the sediments observed, including surface anoxia rich in sulphide (<xref ref-type="bibr" rid="B80">Robertson et&#xa0;al., 2017</xref>). Thus, both the BBN predictions and the present observations of aRPD implied greatly reduced benthic habitat quality and volume for macrobenthic communities (<xref ref-type="bibr" rid="B32">Green et&#xa0;al., 2014</xref>).</p>
<p>The macrobenthic taxon abundance and richness index NZ AMBI steadily increased (worsened) over the time series after ~1960, reflecting its dependence on bulk properties (% TOC and % mud: <xref ref-type="bibr" rid="B76">Robertson et&#xa0;al., 2016</xref>). Its sudden increase after ~1998 resulted from steeply worsening macroalgal conditions, interacting with % TOC and % mud (<xref ref-type="bibr" rid="B32">Green et&#xa0;al., 2014</xref>). Macrobenthic state in NRE was not routinely measured prior to 2001, but after 2001 it was shown to have &#x2018;polluted&#x2019; status at many sites (NZ AMBI &gt; 4.3) (<xref ref-type="bibr" rid="B76">Robertson et&#xa0;al., 2016</xref>, <xref ref-type="bibr" rid="B80">2017</xref>; <xref ref-type="bibr" rid="B29">Forrest et&#xa0;al., 2022</xref>) at sites with high mud content and macroalgal cover. The worsening in NZ AMBI was consistent with findings that New Zealand estuaries with less mud tend to be more favourable for healthy macrobenthic communities (<xref ref-type="bibr" rid="B16">Clark et&#xa0;al., 2020</xref>).</p>
<p>Sporadic aerial imagery shows intertidal seagrass beds present on the central flats of the NRE (south of the Waihopai) in the 1940&#x2019;s (first images in 1946) through to the 60&#x2019;s however these beds were not recorded on imagery in the 80&#x2019;s (there is no whole estuary imagery for the 70&#x2019;s and imagery from 1996 was taken at high tide so inconclusive, K. Roberts, unpubl. data March 2023). In 2001 seagrass was mapped as part of the National Estuary Monitoring Protocol development and was localised to the Waihopai Arm, sheltered areas on the eastern margin and at the entrance (<xref ref-type="bibr" rid="B73">Roberts et&#xa0;al., 2022b</xref>). Seagrass was then subsequently lost over the period 2001-2022, the largest losses occurred between 2001 and 2012, with a 85.2% decrease (-83.1 ha) in seagrass areas with &gt;50% cover in the NRE between 2001 and 2022 with much of this decline occurring in the Waihopai Arm (<xref ref-type="bibr" rid="B72">Roberts et&#xa0;al., 2022a</xref>, <xref ref-type="bibr" rid="B73">b</xref>). It is likely that these historic seagrass beds were affected by reductions in water clarity associated with increased SAR, when SAR transitioned from 2.4 and 3.8 mm y<sup>-1</sup> to 17 and 34 mm y<sup>-1</sup> at the Upper and Central sites, respectively. The latter values greatly exceeded the 5 mm y<sup>-1</sup> SAR threshold described by <xref ref-type="bibr" rid="B100">Zabarte-Maeztu et&#xa0;al. (2021)</xref> and <xref ref-type="bibr" rid="B9">Benham et&#xa0;al. (2019)</xref> for <italic>Z. muelleri</italic> decline. These declines were also coincident with heavy macroalgal growth, particularly <italic>Agarophyton</italic> spp., and eutrophication (<xref ref-type="bibr" rid="B85">Stevens and Forrest, 2020</xref>; <xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>) driven by the steeply increasing potential TN.</p>
<p>The BBN ETI scores slightly overestimate the worsening condition of the NRE after 2012 compared to ETI scores calculated by <xref ref-type="bibr" rid="B72">Roberts et&#xa0;al. (2022a)</xref> in 2022 obtained via direct estuary monitoring in the Waihopai Arm (<xref ref-type="table" rid="T5">
<bold>Table&#xa0;5</bold>
</xref>), however, we can still see the trend in estuarine trophic quality degradation in both the direct estuarine monitoring-based estimates and the BBN estimates.</p>
<table-wrap id="T5" position="float">
<label>Table&#xa0;5</label>
<caption>
<p>Estuarine Trophic Index (ETI) scores calculated at the two sediment coring sites (Upper and Central) and Waihopai Arm scores calculated by <xref ref-type="bibr" rid="B72">Roberts et&#xa0;al. (2022a)</xref> for 2012, 2016, and 2022.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Year</th>
<th valign="top" align="left">Waihopai Arm Score</th>
<th valign="top" align="left">Upper</th>
<th valign="top" align="left">Central</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">2012</td>
<td valign="top" align="left">0.774</td>
<td valign="top" align="left">0.82 (2013)</td>
<td valign="top" align="left">0.82 (2011)</td>
</tr>
<tr>
<td valign="top" align="left">2016</td>
<td valign="top" align="left">0.782</td>
<td valign="top" align="left">0.82</td>
<td valign="top" align="left">0.82</td>
</tr>
<tr>
<td valign="top" align="left">2022</td>
<td valign="top" align="left">0.797</td>
<td valign="top" align="left">0.86 (2019)</td>
<td valign="top" align="left">0.86 (2019)</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Where a particular yearly calculated score is unavailable, the available closest year ETI score is provided with the year in brackets.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Another potential reason for the discrepancy between the monitoring-based ETI estimates and the BBN estimates is that the BBN ETI scores may approach, but not reach, the theoretical limits of the scores (0 and 1). As described by <xref ref-type="bibr" rid="B104">Zeldis and Plew (2022)</xref> the distribution of ETI scores predicted by the BBN become compressed at the upper and lower limits, with achieved maximum and minimum values of ~0.9 and ~0.2. This is an artifact of the probability distributions built into BBNs generally, which can spread predictions for indicators, and their resulting scores, over several bands (<xref ref-type="bibr" rid="B92">Uusitalo, 2007</xref>). Mean values for each indicator and the ETI score were calculated from those probability distributions (using the centre of each band weighted by the probability distribution across each band), tending to bring the ETI scores away from the limits of 0 or 1. Confidence intervals were calculated in a similar manner. The consequence of this is that confidence intervals were wide (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>) due, in part, to the discretised structure of BBNs (<xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>).</p>
<p>The EQR values calculated within the BBN can be assumed to represent the whole of the estuary (<xref ref-type="bibr" rid="B104">Zeldis and Plew, 2022</xref>), because of the use of estuary-wide potential TN estimates and spatially representative EQR (<xref ref-type="bibr" rid="B96">Water Framework Directive - United Kingdom Technical Advisory Group, 2014</xref>). However, scores for other BBN indicators were sensitive to the inputs acquired from the localised area where the cores were collected (Waihopai Arm) including % TOC and % mud. Waihopai Arm is recognised to be in relatively poor ecological condition in the present-day estuary (<xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>), so this needs to be considered when interpreting results. The indicators predicted from those inputs (i.e., aRPD, seagrass, macrobenthos) were likely &#x2018;more impacted&#x2019; than would have been estimated if &#x2018;estuary average&#x2019; core-derived values had been available. It would therefore be beneficial to expand this study using cores from other suitable areas in the NRE recognised to represent a wider range of potential impact to provide a better estuary-wide comparison.</p>
</sec>
<sec id="s4_5">
<label>4.5</label>
<title>Historic ecological health and prospects for remediation</title>
<p>The coring, hindcasted potential TN and BBN results show &#x2018;from where we have come&#x2019; in terms of NRE ecological conditions dating from mid-19<sup>th</sup> to early 21<sup>st</sup> centuries. Severe ecological impact has arisen from both sediment and TN loading, with strong impacts from sediments appearing earlier (mid-20<sup>th</sup> century) than those from TN (late 20<sup>th</sup>/early 21<sup>st</sup> century). These differing impacts are associated with differing land use profiles over time (land clearance, grazing from sheep and beef farming eliciting the earlier sedimentation and moderate TN increases, accelerated N runoff from intensified dairy farming accounting for recent rapid TN increases). The correlations of sediment properties with runoff and associations with environmental degradation resembled responses in Jamaica Bay Estuary (<xref ref-type="bibr" rid="B98">Wigand et&#xa0;al., 2014</xref>), although in that case the runoff variability was driven largely by human population size changes, rather than agricultural intensification.</p>
<p>Because the ecological responses in NRE show temporal associations with agricultural land use profiles, we believe that the results may give guidance in designating restoration aspirations and targets for limit-setting, by identifying times of pivotal changes in estuary state and associated catchment sediment and nutrient inputs. For example, identifying the sediment loading conditions that triggered turbidity-driven subtidal seagrass and macrobenthic declines prior to ~1990 could indicate the sediment runoff conditions required to re-establish those communities, or those required to prevent deterioration of similar, but currently less impacted, estuaries. The large decrease in seagrass cover (in the early 21<sup>st</sup> century) and continued decline of macrobenthic health could indicate the potential TN levels required to prevent macroalgal outbreaks and their overgrowth of intertidal seagrass beds and impacts on macrobenthos in this and other estuaries.</p>
<p>Although our study indicates previous estuary ecological states and their prevailing drivers, the compounding effects of sedimentation and nutrient pressures also indicate that reversion to previous states may not follow a linear recovery trajectory, and that interactions between the pressures could be difficult to predict. There could be hysteresis in the recovery trajectory (<xref ref-type="bibr" rid="B21">Duarte et&#xa0;al., 2009</xref>, <xref ref-type="bibr" rid="B20">2015</xref>; <xref ref-type="bibr" rid="B34">Hale et&#xa0;al., 2023</xref>), potentially following the &#x2018;hysteresis&#x2019; or &#x2018;hysteresis with threshold&#x2019; patterns of <xref ref-type="bibr" rid="B42">Kemp et&#xa0;al. (2009)</xref>, wherein degradation and recovery do not follow similar but opposite trajectories (<xref ref-type="bibr" rid="B20">Duarte et&#xa0;al., 2015</xref>). The compounding stressors would likely have different time scales of remediation. If catchment sediment loading and SAR were significantly reduced, there could still be a long time-lag before the deep deposits of mud (with its high C and N contents) were washed naturally from NRE. Also, macroalgae trap mud in NRE (<xref ref-type="bibr" rid="B86">Stevens et&#xa0;al., 2022</xref>), so unless their excessive biomass was reduced, sediment winnowing could be slow. Reducing TN loading to reduce macroalgal eutrophication would likely be beneficial, but the benefit could be slowed by persistence of sediment-bound N in the muddy deposits (<xref ref-type="bibr" rid="B28">Engelsen et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B11">Borja et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B103">Zeldis et&#xa0;al., 2019</xref>) and a response lag in <italic>Agarophyton</italic> spp. reduction due to utilization of internal nutrient storage pools (<xref ref-type="bibr" rid="B22">Dudley et&#xa0;al., 2022</xref>). High porewater DIN concentrations (~160 &#xb5;mol) occur in NRE sediments (<xref ref-type="bibr" rid="B75">Robertson and Savage, 2018</xref>), and these may be important (along with water column DIN) in supporting its eutrophic biomasses of macroalgae. The situation contrasts with the significant, rapid recovery of ecological condition in Avon-Heathcote Estuary (Christchurch, NZ) upon reduction of its N load (<xref ref-type="bibr" rid="B7">Barr et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B101">Zeldis et&#xa0;al., 2020</xref>), which was facilitated by the naturally coarse (sandy), well-flushed sediments of that estuary.</p>
<p>This discussion is not intended to discourage remediation efforts, but to indicate possible outcomes. There would undoubtedly be benefits accruing from better runoff and wastewater management in NRE, including efforts to simultaneously reduce sediment and nutrient loads, in terms of realising improved ecological states. These results are also informative by illustrating historic sediment and nutrient runoff conditions that could lead to NRE-like conditions in other similar, but currently less impacted, estuaries.</p>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Summary and conclusion</title>
<p>In this work we have provided reconstructed NRE ecological health records based on historical biogeochemical properties derived from a novel combination of aged sediment core data, models of hindcasted nutrient loads and estuary potential N concentrations, and BBN predictions of current and historical ecological conditions. Our results suggest that catchment N load (including diffuse and point source loads) is predictable from sediment core N and, more broadly, that the combination of methods used is effective in hindcasting estuary health over centuries. However, work is needed to further validate the ecological state hindcasts, and to establish how transferable the results are, both to other parts of NRE and to other estuaries.</p>
<p>Using sediment cores to derive historic estuarine condition, in the absence of historic time series measurements or catchment land use information, could inform nutrient load limits for other estuaries, to protect them from impacts of excessive nutrient and sediment loads. Runoff conditions corresponding to when significant shifts in estuarine health occurred, and the corresponding states of land use, could be used to propose targets for nutrient and sediment input reduction and associated limit-setting in the catchment, or determination of upper limits to prevent further estuarine decline. The combination of sediment core sampling with modelling tools demonstrated here shows &#x2018;from where we have come&#x2019; in terms of history of change.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Datasheet 2</bold>
</xref>. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>RH: Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. JZ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. BD: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. AH: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. DP: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. US: Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing &#x2013; review &amp; editing. AS: Formal analysis, Funding acquisition, Investigation, Validation, Visualization, Writing &#x2013; review &amp; editing. KR: Investigation, Resources, Writing &#x2013; review &amp; editing. SO: Investigation, Resources, Writing &#x2013; review &amp; editing. PV: Formal analysis, Investigation, Methodology, Writing &#x2013; review &amp; editing.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was funded under the Coasts and Estuaries platform by the National Institute for Water &amp; Atmospheric Research Ltd.&#x2019;s Strategic Science Investment Fund (FWCE2106, FWCE2206, CECE2305, CECE2402) from New Zealand&#x2019;s Ministry for Business, Innovation and Employment core funding.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>Thanks to Environment Southland for in-kind support with fieldwork and provision of historic maps for digitization. Thanks to Invercargill City Council and Adrian Cocker for provision of water quality monitoring data. Stable Isotope Analysis was carried out by the NIWA Environmental and Ecological Stable Isotope Facility in Wellington, New Zealand: Sample acidification and sample weighing was carried out by Josette Delgado and Anna Kilimnik. Standards weighing was carried out by Anna Kilimnik and Rahul Peethambaran. Stable isotope analysis was completed by Josette Delgado and Rahul Peethambaran. Daily stable isotope data were compiled by Josette Delgado. Data processing and calculations were carried out by Julie Brown. Sample processing, data analyses and calculations were overseen, checked, and amended where necessary by Sarah Bury. We thank Amandine Sabadel for assistance with interpretation of stable isotopes. Radio-dating was carried out by ESR labs, and we thank Michael Lechermann and Levi Bourke for assistance. We thank Ron Ovendon for their assistance with sediment core slicing and NIWA Hamilton Water Quality Laboratory for sediment sample analysis.</p>
</ack>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Authors KR and SO were employed by the company Salt Ecology.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmars.2024.1374869/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmars.2024.1374869/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="DataSheet_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="DataSheet_2.xlsx" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alldred</surname> <given-names>F. C.</given-names>
</name>
<name>
<surname>Gr&#xf6;cke</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Leung</surname> <given-names>C. Y.</given-names>
</name>
<name>
<surname>Wright</surname> <given-names>L. P.</given-names>
</name>
<name>
<surname>Banfield</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Diffuse and concentrated nitrogen sewage pollution in island environments with differing treatment systems</article-title>. <source>Sci. Rep.</source> <volume>13</volume>, <fpage>4838</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-023-32105-6</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>AsureQuality</collab>
</person-group> (<year>2008</year>) <source>Agribase land use boundaries</source>. Available online at: <uri xlink:href="https://www.asurequality.com/services/agribase/">https://www.asurequality.com/services/agribase/</uri> (Accessed <access-date>2022</access-date>).</citation>
</ref>
<ref id="B3">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Ausseil</surname> <given-names>A.-G. E.</given-names>
</name>
<name>
<surname>Dymond</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Weeks</surname> <given-names>E. S.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>Provision of Natural Habitat for Biodiversity: Quantifying Recent Trends in New Zealand</article-title>,&#x201d; in <source>Biodiversity Loss in a Changing Planet</source>. eds. <person-group person-group-type="editor">
<name>
<surname>Grillo</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Venora</surname> <given-names>G.</given-names>
</name>
</person-group> (<publisher-loc>Rijeka</publisher-loc>: <publisher-name>IntechOpen</publisher-name>), Ch. <fpage>9</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5772/24969</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barbier</surname> <given-names>E. B.</given-names>
</name>
<name>
<surname>Hacker</surname> <given-names>S. D.</given-names>
</name>
<name>
<surname>Kennedy</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Koch</surname> <given-names>E. W.</given-names>
</name>
<name>
<surname>Stier</surname> <given-names>A. C.</given-names>
</name>
<name>
<surname>Silliman</surname> <given-names>B. R.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>The value of estuarine and coastal ecosystem services</article-title>. <source>Ecol. Monogr.</source> <volume>81</volume>, <fpage>169</fpage>&#x2013;<lpage>193</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/10-1510.1</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barr</surname> <given-names>N. G.</given-names>
</name>
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Rogers</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Cornelisen</surname> <given-names>C. D.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Broad-scale patterns of tissue-&#x3b4;15N and tissue-N indices in frondose Ulva spp.; developing a national baseline indicator of nitrogen-loading for coastal New Zealand</article-title>. <source>Mar. pollut. Bull.</source> <volume>67</volume>, <fpage>203</fpage>&#x2013;<lpage>216</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.marpolbul.2012.11.033</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barr</surname> <given-names>N. G.</given-names>
</name>
<name>
<surname>Rees</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Nitrogen status and metabolism in the green seaweed Enteromorpha intestinalis: an examination of three natural populations</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>249</volume>, <fpage>133</fpage>&#x2013;<lpage>144</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps249133</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barr</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Scheuer</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Schiel</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Macroalgal bioindicators of recovery from eutrophication in a tidal lagoon following wastewater diversion and earthquake disturbance</article-title>. <source>Estuaries Coasts</source> <volume>43</volume>, <fpage>240</fpage>&#x2013;<lpage>255</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-019-00654-7</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bedard-Haughn</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Van Groenigen</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Van Kessel</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Tracing 15 N through landscapes: potential uses and precautions</article-title>. <source>J. Hydrology</source> <volume>272</volume>, <fpage>175</fpage>&#x2013;<lpage>190</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0022-1694(02)00263-9</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benham</surname> <given-names>C. F.</given-names>
</name>
<name>
<surname>Beavis</surname> <given-names>S. G.</given-names>
</name>
<name>
<surname>Jackson</surname> <given-names>E. L.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Tolerance of tropical seagrasses <italic>Zostera muelleri</italic> and <italic>Halophila ovalis</italic> to burial: Toward an understanding of threshold effects</article-title>. <source>Estuarine Coast. Shelf Sci.</source> <volume>218</volume>, <fpage>131</fpage>&#x2013;<lpage>138</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecss.2018.11.005</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bentley</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Swales</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Pyenson</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Dawe</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Sedimentation, bioturbation, and sedimentary fabric evolution on a modern mesotidal mudflat: A multi-tracer study of processes, rates, and scales</article-title>. <source>Estuarine Coast. Shelf Sci.</source> <volume>141</volume>, <fpage>58</fpage>&#x2013;<lpage>68</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecss.2014.02.004</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Borja</surname> <given-names>&#xc1;.</given-names>
</name>
<name>
<surname>Dauer</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Elliott</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Simenstad</surname> <given-names>C. A.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Medium-and long-term recovery of estuarine and coastal ecosystems: patterns, rates and restoration effectiveness</article-title>. <source>Estuaries Coasts</source> <volume>33</volume>, <fpage>1249</fpage>&#x2013;<lpage>1260</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-010-9347-5</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Borja</surname> <given-names>&#xc1;.</given-names>
</name>
<name>
<surname>Franco</surname> <given-names>J.</given-names>
</name>
<name>
<surname>P&#xe9;rez</surname> <given-names>V.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>A marine biotic index to establish the ecological quality of soft-bottom benthos within European estuarine and coastal environments</article-title>. <source>Mar. pollut. Bull.</source> <volume>40</volume>, <fpage>1100</fpage>&#x2013;<lpage>1114</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0025-326X(00)00061-8</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Brown</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Geochemistry and isotopic composition of sediment cores to understand the lithological and anthropogenic controls on eutrophication in the New River Estuary, Southland, New Zealand</source> (<publisher-loc>Ontario, Canada</publisher-loc>: <publisher-name>Masters, Laurentian University Sudbury</publisher-name>).</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Canfield</surname> <given-names>D. E.</given-names>
</name>
</person-group> (<year>1994</year>). <article-title>Factors influencing organic carbon preservation in marine sediments</article-title>. <source>Chem. Geology</source> <volume>114</volume>, <fpage>315</fpage>&#x2013;<lpage>329</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0009-2541(94)90061-2</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cifuentes</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Sharp</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fogel</surname> <given-names>M. L.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Stable carbon and nitrogen isotope biogeochemistry in the Delaware estuary</article-title>. <source>Limnology oceanography</source> <volume>33</volume>, <fpage>1102</fpage>&#x2013;<lpage>1115</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.1988.33.5.1102</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Clark</surname> <given-names>D. E.</given-names>
</name>
<name>
<surname>Hewitt</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Pilditch</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Ellis</surname> <given-names>J. I.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The development of a national approach to monitoring estuarine health based on multivariate analysis</article-title>. <source>Mar. pollut. Bull.</source> <volume>150</volume>, <fpage>110602</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.marpolbul.2019.110602</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cloern</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Canuel</surname> <given-names>E. A.</given-names>
</name>
<name>
<surname>Harris</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Stable carbon and nitrogen isotope composition of aquatic and terrestrial plants of the San Francisco Bay estuarine system</article-title>. <source>Limnology Oceanography</source> <volume>47</volume>, <fpage>713</fpage>&#x2013;<lpage>729</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.2002.47.3.0713</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Conley</surname> <given-names>D. J.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Biogeochemical nutrient cycles and nutrient management strategies</article-title>. In: <person-group person-group-type="editor">
<name>
<surname>Garnier</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Mouchel</surname> <given-names>J. M.</given-names>
</name>
</person-group> (eds) <source>Man and River Systems. Developments in Hydrobiology</source>. (<publisher-loc>Dordrecht</publisher-loc>: <publisher-name>Springer</publisher-name>) vol <volume>146</volume>. doi: <pub-id pub-id-type="doi">10.1007/978-94-017-2163-9_10</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cornelisen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Wing</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Clark</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Bowman</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Frew</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Hurd</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Patterns in the d13C and d15N signature of Ulva pertusa: Interaction between physical gradients and nutrient source pools</article-title>. <source>Limnology Oceanography</source> <volume>52</volume>, <fpage>820</fpage>&#x2013;<lpage>832</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.2007.52.2.0820</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Duarte</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Borja</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Carstensen</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Elliott</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Krause-Jensen</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Marb&#xe0;</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Paradigms in the recovery of estuarine and coastal ecosystems</article-title>. <source>Estuaries Coasts</source> <volume>38</volume>, <fpage>1202</fpage>&#x2013;<lpage>1212</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-013-9750-9</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Duarte</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Conley</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Carstensen</surname> <given-names>J.</given-names>
</name>
<name>
<surname>S&#xe1;nchez-Camacho</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Return to Neverland: shifting baselines affect eutrophication restoration targets</article-title>. <source>Estuaries Coasts</source> <volume>32</volume>, <fpage>29</fpage>&#x2013;<lpage>36</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-008-9111-2</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Barr</surname> <given-names>N. G.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Scheuer</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Experiments to parametrise a growth and nutrient storage model for Agarophyton spp</article-title>. <source>Estuarine Coast. Shelf Sci.</source> <volume>264</volume>, <fpage>107660</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecss.2021.107660</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Barr</surname> <given-names>N. G.</given-names>
</name>
<name>
<surname>Shima</surname> <given-names>J. S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Influence of light intensity and nutrient source on &#x3b4;13C and &#x3b4;15N signatures in Ulva pertusa</article-title>. <source>Aquat. Biol.</source> <volume>9</volume>, <fpage>85</fpage>&#x2013;<lpage>93</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/ab00241</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Burge</surname> <given-names>O. R.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Effects of agricultural and urban land cover on New Zealand&#x2019;s estuarine water quality</article-title>. <source>New Z. J. Mar. Freshw. Res.</source> <volume>54</volume>, <fpage>372</fpage>&#x2013;<lpage>392</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00288330.2020.1729819</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Dudley</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Burge</surname> <given-names>O.</given-names>
</name>
</person-group> (<year>2017</year>). <source>New Zealand Coastal Water Quality Assessment. NIWA Report No 2016093CH. Prepared for Ministry for the Environment</source>. (<publisher-loc>Christchurch, NZ</publisher-loc>: <publisher-name>NIWA</publisher-name>). <fpage>84</fpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Durante</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>R. O.</given-names>
</name>
<name>
<surname>Kolodzey</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Mcmullin</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Salmond</surname> <given-names>N. H.</given-names>
</name>
<name>
<surname>Schlieman</surname> <given-names>C. D.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Oceanographic transport along frontal zones forms carbon, nitrogen, and oxygen isoscapes on the east coast of New Zealand: Implications for ecological studies</article-title>. <source>Continental Shelf Res.</source> <volume>216</volume>, <fpage>104368</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.csr.2021.104368</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elliott</surname> <given-names>A. H.</given-names>
</name>
<name>
<surname>Semadeni-Davies</surname> <given-names>A. F.</given-names>
</name>
<name>
<surname>Shankar</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Wheeler</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>A national-scale GIS-based system for modelling impacts of land use on water quality</article-title>. <source>Environ. Model. Software</source> <volume>86</volume>, <fpage>131</fpage>&#x2013;<lpage>144</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.envsoft.2016.09.011</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Engelsen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Hulth</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Pihl</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Sundb&#xe4;ck</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Benthic trophic status and nutrient fluxes in shallow-water sediments</article-title>. <source>Estuarine Coast. Shelf Sci.</source> <volume>78</volume>, <fpage>783</fpage>&#x2013;<lpage>795</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecss.2008.02.018</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Forrest</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Roberts</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Rabel</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). &#x201c;<article-title>Fine scale intertidal monitoring of new river estuary: synthesis of survey data 2001-2022</article-title>,&#x201d; in <source>Salt Ecology Report 100, prepared for Environment Southland</source>. (<publisher-loc>Nelson, NZ</publisher-loc>: <publisher-name>Salt Ecology</publisher-name>) <fpage>47</fpage> p.</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fowler</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Coyle</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Skiba</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Sutton</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Cape</surname> <given-names>J. N.</given-names>
</name>
<name>
<surname>Reis</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>The global nitrogen cycle in the twenty-first century</article-title>. <source>Philos. Trans. R. Soc. B: Biol. Sci.</source> <volume>368</volume>, <fpage>20130164</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1098/rstb.2013.0164</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Grant</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2023</year>) <source>Southland places - Invercargill', Te Ara - the Encyclopedia of New Zealand</source>. Available online at: <uri xlink:href="http://www.TeAra.govt.nz/en/southland-places/page-1">http://www.TeAra.govt.nz/en/southland-places/page-1</uri> (Accessed <access-date>28 August 2023</access-date>).</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Green</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Sutula</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Fong</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>How much is too much? Identifying benchmarks of adverse effects of macroalgae on the macrofauna in intertidal flats</article-title>. <source>Ecol. Appl.</source> <volume>24</volume>, <fpage>300</fpage>&#x2013;<lpage>314</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/13-0524.1</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Haddadchi</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Hicks</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Understanding the effect of catchment characteristics on suspended sediment dynamics during flood events</article-title>. <source>Hydrological Processes</source> <volume>34</volume>, <fpage>1558</fpage>&#x2013;<lpage>1574</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/hyp.13682</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hale</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Depree</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Broekhuizen</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Simulating fish farm enrichment and fallowing impacts reveals unequal biogeochemical recovery of benthic variables</article-title>. <source>Aquaculture Environ. Interact.</source> <volume>15</volume>, <fpage>115</fpage>&#x2013;<lpage>131</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/aei00456</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Handley</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Horrocks</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Willis</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Bradley</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Bury</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Brown</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Historical analyses of coastal marine sediments reveal land-based impacts on ecosystem productivity</article-title>. <source>New Z. J. Ecol.</source> <volume>44</volume>, <fpage>3415</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.20417/nzjecol.44.18</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hewitt</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Bell</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Costello</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Cummings</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Currie</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Ellis</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2014</year>). &#x201c;<article-title>Development of a national marine environment monitoring programme (MEMP) for New Zealand</article-title>,&#x201d; in <source>NNew Zealand Aquatic Environment and Biodiversity Report No. 141</source>. (<publisher-loc>Wellington, NZ</publisher-loc>: <publisher-name>Ministry for Primary Industries</publisher-name>).</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hilgendorf</surname> <given-names>H. W.</given-names>
</name>
</person-group> (<year>1935</year>). <article-title>The grasslands of the South Island of New Zealand: An ecological survey</article-title>. <source>DSIR Bull.</source> <volume>47</volume>, <fpage>24</fpage>.</citation>
</ref>
<ref id="B38">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Howard-Williams</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Davies-Colley</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Rutherford</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Wilcock</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>Diffuse pollution and freshwater degradation: New Zealand perspectives</article-title>&#x201d;, in <person-group person-group-type="editor">
<name>
<surname>Bochove van</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Vanrolleghem</surname> <given-names>P. A.</given-names>
</name>
<name>
<surname>Chambers</surname> <given-names>P. A.</given-names>
</name>
<name>
<surname>Th&#xe9;riault</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Novotn&#xe1;</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Burkart</surname> <given-names>M. R.</given-names>
</name>
<etal/>
</person-group> (eds) <conf-name>Issues and Solutions to Diffuse Pollution, selected papers from the 14th international conference of the IWA Diffuse Pollution Specialist Group</conf-name>, <publisher-name>DIPCON 2010</publisher-name>.</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Howarth</surname> <given-names>R. W.</given-names>
</name>
<name>
<surname>Marino</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Nitrogen as the limiting nutrient for eutrophication in coastal marine ecosystems: Evolving views over three decades</article-title>. <source>Limnology Oceanography</source> <volume>51</volume>, <fpage>364</fpage>&#x2013;<lpage>376</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.2006.51.1_part_2.0364</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hume</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>Bell</surname> <given-names>R. G.</given-names>
</name>
<name>
<surname>De Lange</surname> <given-names>W. P.</given-names>
</name>
<name>
<surname>Healy</surname> <given-names>T. R.</given-names>
</name>
<name>
<surname>Hicks</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Kirk</surname> <given-names>R. M.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>Coastal oceanography and sedimentology in New Zealand 1967&#x2013;91</article-title>. <source>New Z. J. Mar. Freshw. Res.</source> <volume>26</volume>, <fpage>1</fpage>&#x2013;<lpage>36</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00288330.1992.9516497</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jackson</surname> <given-names>J. B. C.</given-names>
</name>
<name>
<surname>Kirby</surname> <given-names>M. X.</given-names>
</name>
<name>
<surname>Berger</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Bjorndal</surname> <given-names>K. A.</given-names>
</name>
<name>
<surname>Botsford</surname> <given-names>L. W.</given-names>
</name>
<name>
<surname>Bourque</surname> <given-names>B. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2001</year>). <article-title>Historical overfishing and the recent collapse of coastal ecosystems</article-title>. <source>Science</source> <volume>293</volume>, <fpage>629</fpage>&#x2013;<lpage>637</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1059199</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kemp</surname> <given-names>W. M.</given-names>
</name>
<name>
<surname>Testa</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Conley</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Gilbert</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Hagy</surname> <given-names>J. D.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Temporal responses of coastal hypoxia to nutrient loading and physical controls</article-title>. <source>Biogeosciences</source> <volume>6</volume>, <fpage>2985</fpage>&#x2013;<lpage>3008</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5194/bg-6-2985-2009</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kennish</surname> <given-names>M. J.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Environmental threats and environmental future of estuaries</article-title>. <source>Environ. Conserv.</source> <volume>29</volume>, <fpage>78</fpage>&#x2013;<lpage>107</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1017/S0376892902000061</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kritee</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Sigman</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Granger</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ward</surname> <given-names>B. B.</given-names>
</name>
<name>
<surname>Jayakumar</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Deutsch</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Reduced isotope fractionation by denitrification under conditions relevant to the ocean</article-title>. <source>Geochimica Cosmochimica Acta</source> <volume>92</volume>, <fpage>243</fpage>&#x2013;<lpage>259</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.gca.2012.05.020</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Landcare Research Ltd</collab></person-group>. (<year>2013</year>). <source>Land Cover Database</source>. Available at: <uri xlink:href="https://lris.scinfo.org.nz/layer/104400-lcdb-v50-land-cover-database-version-50-mainland-new-zealand/">https://lris.scinfo.org.nz/layer/104400-lcdb-v50-land-cover-database-version-50-mainland-new-zealand/</uri>.</citation>
</ref>
<ref id="B46">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Landcare Research Ltd</collab>
</person-group>. (<year>2018</year>). <source>New Zealand Land Resource Inventory</source>. Available at: <uri xlink:href="https://lris.scinfo.org.nz/data/?q=NZLRI">https://lris.scinfo.org.nz/data/?q=NZLRI</uri>.</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Larned</surname> <given-names>S. T.</given-names>
</name>
<name>
<surname>Snelder</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Unwin</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Mcbride</surname> <given-names>G. B.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Water quality in New Zealand rivers: current state and trends</article-title>. <source>New Z. J. Mar. Freshw. Res.</source> <volume>50</volume>, <fpage>389</fpage>&#x2013;<lpage>417</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00288330.2016.1150309</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leathwick</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>New Zealand's potential forest pattern as predicted from current species-environment relationships</article-title>. <source>New Z. J. Bot.</source> <volume>39</volume>, <fpage>447</fpage>&#x2013;<lpage>464</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/0028825X.2001.9512748</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Leathwick</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Mcglone</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Walker</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2012</year>)<article-title>New Zealand's potential vegetation pattern</article-title>. In: <source>Landcare Research</source>. Available online at: <uri xlink:href="https://lris.scinfo.org.nz/layer/48289-potential-vegetation-of-new-zealand/">https://lris.scinfo.org.nz/layer/48289-potential-vegetation-of-new-zealand/</uri> (Accessed <access-date>25/05/2018</access-date>).</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leathwick</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Overton</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Mcleod</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>An environmental domain classification of New Zealand and its use as a tool for biodiversity management</article-title>. <source>Conserv. Biol.</source> <volume>17</volume>, <fpage>1612</fpage>&#x2013;<lpage>1623</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1523-1739.2003.00469.x</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Leathwick</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wilson</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Rutledge</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Wardle</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Morgan</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Johnston</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2002</year>). <source>Land environments of New Zealand</source> (<publisher-loc>Auckland, New Zealand</publisher-loc>: <publisher-name>Bateman</publisher-name>).</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Le Guitton</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Soetaert</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Damst&#xe9;</surname> <given-names>J. S. S.</given-names>
</name>
<name>
<surname>Middelburg</surname> <given-names>J. J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Biogeochemical consequences of vertical and lateral transport of particulate organic matter in the southern North Sea: A multiproxy approach</article-title>. <source>Estuarine Coast. Shelf Sci.</source> <volume>165</volume>, <fpage>117</fpage>&#x2013;<lpage>127</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecss.2015.09.010</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Le Moal</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Gascuel-Odoux</surname> <given-names>C.</given-names>
</name>
<name>
<surname>M&#xe9;nesguen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Souchon</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>&#xc9;trillard</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Levain</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Eutrophication: A new wine in an old bottle</article-title>? <source>Sci. Total Environ.</source> <volume>651</volume>, <fpage>1</fpage>&#x2013;<lpage>11</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.scitotenv.2018.09.139</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lohrer</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Hewitt</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Thrush</surname> <given-names>S. F.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Assessing far-field effects of terrigenous sediment loading in the coastal marine environment</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>315</volume>, <fpage>13</fpage>&#x2013;<lpage>18</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps315013</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lotze</surname> <given-names>H. K.</given-names>
</name>
<name>
<surname>Lenihan</surname> <given-names>H. S.</given-names>
</name>
<name>
<surname>Bourque</surname> <given-names>B. J.</given-names>
</name>
<name>
<surname>Bradbury</surname> <given-names>R. H.</given-names>
</name>
<name>
<surname>Cooke</surname> <given-names>R. G.</given-names>
</name>
<name>
<surname>Kay</surname> <given-names>M. C.</given-names>
</name>
<etal/>
</person-group>. (<year>2006</year>). <article-title>Depletion, degradation, and recovery potential of estuaries and coastal seas</article-title>. <source>Science</source> <volume>312</volume>, <fpage>1806</fpage>&#x2013;<lpage>1809</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1128035</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lotze</surname> <given-names>H. K.</given-names>
</name>
<name>
<surname>Milewski</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Two centuries of multiple human impacts and successive changes in a North Atlantic food web</article-title>. <source>Ecol. Appl.</source> <volume>14</volume>, <fpage>1428</fpage>&#x2013;<lpage>1447</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/03-5027</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lotze</surname> <given-names>H. K.</given-names>
</name>
<name>
<surname>Reise</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Worm</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Van Beusekom</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Busch</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ehlers</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2005</year>). <article-title>Human transformations of the Wadden Sea ecosystem through time: a synthesis</article-title>. <source>Helgoland Mar. Res.</source> <volume>59</volume>, <fpage>84</fpage>&#x2013;<lpage>95</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10152-004-0209-z</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ledgard</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>New Zealand dairy farm systems and key environmental effects</article-title>. <source>Front. Agr. Sci. Eng.</source> <volume>8</volume>, <fpage>148</fpage>&#x2013;<lpage>158</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.15302/J-FASE-2020372</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Matthews</surname> <given-names>K. M.</given-names>
</name>
</person-group> (<year>1989</year>). <source>Radioactive fallout in the South Pacific: a history Part 1: deposition in New Zealand</source> (<publisher-loc>Christchurch New Zealand</publisher-loc>: <publisher-name>National Radiation Lab</publisher-name>).</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McClelland</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Valiela</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Michener</surname> <given-names>R. H.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Nitrogen-stable isotope signatures in estuarine food webs: A record of increasing urbanization in coastal watersheds</article-title>. <source>Limnology Oceanography</source> <volume>42</volume>, <fpage>930</fpage>&#x2013;<lpage>937</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.1997.42.5.0930</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Measures</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2016</year>). &#x201c;<article-title>New River Estuary Hydrodynamic Modelling. Model build, calibration and tracer simulations</article-title>,&#x201d; in <source>NIWA Client Report No: 2016077CH</source>. (<publisher-loc>Christchurch, NZ</publisher-loc>: <publisher-name>NIWA</publisher-name>).</citation>
</ref>
<ref id="B62">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Ministry for the Environment</collab>
</person-group> (<year>2007</year>). <source>Environment New Zealand 2007</source> (<publisher-loc>Wellington, New Zealand</publisher-loc>: <publisher-name>Ministry for the Environment</publisher-name>).</citation>
</ref>
<ref id="B63">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Ministry for the Environment</collab>
</person-group> (<year>2021</year>). <source>New Zealand&#x2019;s Greenhouse Gas Inventory 1990&#x2013;2019</source> (<publisher-loc>Wellington, New Zealand</publisher-loc>: <publisher-name>Ministry for the Environment</publisher-name>).</citation>
</ref>
<ref id="B64">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Ministry for the Environment</collab>
</person-group> (<year>2023</year>). <source>National Policy Statement for Freshwater Management 2020</source> (<publisher-loc>Wellington, New Zealand</publisher-loc>: <publisher-name>Ministry for the Environment</publisher-name>).</citation>
</ref>
<ref id="B65">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Ministry for the Environment, and Stats Nz</collab>
</person-group>. (<year>2019</year>). <source>New Zealand&#x2019;s Environmental Reporting Series: Environment Aotearoa 2019</source>. <publisher-loc>Wellington, NZ</publisher-loc>: <publisher-name>Ministry for the Environment and Stats NZ</publisher-name>.</citation>
</ref>
<ref id="B66">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Morrison</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lowe</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Parsons</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Usmar</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Mcleod</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2009</year>). &#x201c;<article-title>A review of land-based effects on coastal fisheries and supporting biodiversity in New Zealand</article-title>,&#x201d; in <source>New Zealand Aquatic Environment and Biodiversity Report 37</source> (<publisher-name>New Zealand Ministry of Fisheries</publisher-name>, <publisher-loc>Wellington</publisher-loc>).</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Needoba</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Sigman</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Harrison</surname> <given-names>P. J.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>The mechanism of isotope fractionation during algal nitrate assimilation as illuminated by the 15n/14n of intracellular nitrate</article-title>. <source>J. Phycology</source> <volume>40</volume>, <fpage>517</fpage>&#x2013;<lpage>522</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1529-8817.2004.03172.x</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Whitehead</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Robertson</surname> <given-names>B. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Assessing the eutrophic susceptibility of New Zealand estuaries</article-title>. <source>Estuaries Coasts</source> <volume>43</volume>, <fpage>2015</fpage>&#x2013;<lpage>2033</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-020-00729-w</pub-id>
</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Shankar</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Elliott</surname> <given-names>A. H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Using simple dilution models to predict New Zealand estuarine water quality</article-title>. <source>Estuaries Coasts</source> <volume>41</volume>, <fpage>1643</fpage>&#x2013;<lpage>1659</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-018-0387-6</pub-id>
</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pratt</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Lohrer</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Pilditch</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Thrush</surname> <given-names>S. F.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Changes in ecosystem function across sedimentary gradients in estuaries</article-title>. <source>Ecosystems</source> <volume>17</volume>, <fpage>182</fpage>&#x2013;<lpage>194</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10021-013-9716-6</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rabalais</surname> <given-names>N. N.</given-names>
</name>
<name>
<surname>Atilla</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Normandeau</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Eugene Turner</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Ecosystem history of Mississippi River-influenced continental shelf revealed through preserved phytoplankton pigments</article-title>. <source>Mar. pollut. Bull.</source> <volume>49</volume>, <fpage>537</fpage>&#x2013;<lpage>547</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.marpolbul.2004.03.017</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Roberts</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Scott-Simmonds</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Forrest</surname> <given-names>B. M.</given-names>
</name>
</person-group> (<year>2022</year>a). <source>Broad Scale Intertidal Habitat Mapping of New River Estuary</source>. Salt Ecology Report 097. (<publisher-name>Environment Southland</publisher-name>: <publisher-name>Salt Ecology, Nelson, Nz</publisher-name>).</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roberts</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Forrest</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Shankar</surname> <given-names>U.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>b). <article-title>Use of a multi-metric macroalgal index to track changes in response to nutrient loads, New River Estuary</article-title>. In <source>Salt Ecology Report 085/ NIWA Client Report 2022153CH</source>. (<publisher-loc>Salt Ecology, Nelson, NZ</publisher-loc>: <publisher-name>Environment Southland by Salt Ecology and NIWA under an Envirolink medium advice grant</publisher-name>). <fpage>50</fpage>p.</citation>
</ref>
<ref id="B74">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Roberts</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Ward</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2021</year>). <source>Proposed reference conditions in Southland estuaries: review of historical data and literature</source> (<publisher-loc>Invercargill</publisher-loc>: <publisher-name>Environment southland</publisher-name>).</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Savage</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Mud-entrained macroalgae utilise porewater and overlying water column nutrients to grow in a eutrophic intertidal estuary</article-title>. <source>Biogeochemistry</source> <volume>139</volume>, <fpage>53</fpage>&#x2013;<lpage>68</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10533-018-0454-x</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>B. P.</given-names>
</name>
<name>
<surname>Savage</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Gardner</surname> <given-names>J. P. A.</given-names>
</name>
<name>
<surname>Robertson</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Optimising a widely-used coastal health index through quantitative ecological group classifications and associated thresholds</article-title>. <source>Ecol. Indic.</source> <volume>69</volume>, <fpage>595</fpage>&#x2013;<lpage>605</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecolind.2016.04.003</pub-id>
</citation>
</ref>
<ref id="B77">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2007</year>). <source>Broad Scale Habitat Mapping and Sedimentation Rate</source> (<publisher-loc>Nelson, NZ</publisher-loc>: <publisher-name>Wriggle Ltd</publisher-name>.).</citation>
</ref>
<ref id="B78">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2013</year>a). <source>New River Estuary Preliminary Nutrient and Sediment Load Estimates 2012/13</source> (<publisher-loc>Nelson, NZ</publisher-loc>: <publisher-name>Wriggle Ltd</publisher-name>.).</citation>
</ref>
<ref id="B79">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2013</year>b). <source>New River Estuary: fine scale monitoring of highly eutrophic arms 2012/2013</source> (<publisher-loc>Nelson, New Zealand</publisher-loc>: <publisher-name>Report prepared by Wriggle Coastal Management for Environment Southland</publisher-name>).</citation>
</ref>
<ref id="B80">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Ward</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Robertson</surname> <given-names>B. P.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Condition of Southland&#x2019;s shallow, intertidal dominated estuaries in relation to eutrophication and sedimentation: output 1: data analysis and technical assessment - habitat mapping, vulnerability assessment and monitoring recommendations related to issues of eutrophication and sedimentation</article-title>,&#x201d; in <source>Report prepared for Environment Southland</source> (<publisher-name>Wriggle Coastal Management</publisher-name>, <publisher-loc>Nelson</publisher-loc>).</citation>
</ref>
<ref id="B81">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sartory</surname> <given-names>D. P.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>The determination of algal chlorophyllous pigments by high performance liquid chromatography and spectrophotometry</article-title>. <source>Water Res.</source> <volume>19</volume>, <fpage>605</fpage>&#x2013;<lpage>610</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0043-1354(85)90066-1</pub-id>
</citation>
</ref>
<ref id="B82">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Savage</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Tracing the influence of sewage nitrogen in a coastal ecosystem using stable nitrogen isotopes</article-title>. <source>Ambio</source> <volume>34</volume>, <fpage>145</fpage>&#x2013;<lpage>150</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1579/0044-7447-34.2.145</pub-id>
</citation>
</ref>
<ref id="B83">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Savage</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Leavitt</surname> <given-names>P. R.</given-names>
</name>
<name>
<surname>Elmgren</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Effects of land use, urbanization, and climate variability on coastal eutrophication in the Baltic Sea</article-title>. <source>Limnology Oceanography</source> <volume>55</volume>, <fpage>1033</fpage>&#x2013;<lpage>1046</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.2010.55.3.1033</pub-id>
</citation>
</ref>
<ref id="B84">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Snelder</surname> <given-names>T. H.</given-names>
</name>
<name>
<surname>Larned</surname> <given-names>S. T.</given-names>
</name>
<name>
<surname>Mcdowell</surname> <given-names>R. W.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Anthropogenic increases of catchment nitrogen and phosphorus loads in New Zealand</article-title>. <source>New Z. J. Mar. Freshw. Res.</source> <volume>52</volume>, <fpage>336</fpage>&#x2013;<lpage>361</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00288330.2017.1393758</pub-id>
</citation>
</ref>
<ref id="B85">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Forrest</surname> <given-names>B. M.</given-names>
</name>
</person-group> (<year>2020</year>). <source>Macroalgal Monitoring of New River Estuary</source>. Salt Ecology Report 054. (<publisher-loc>Salt Ecology, Nelson, NZ</publisher-loc>: <publisher-name>Environment Southland</publisher-name>).</citation>
</ref>
<ref id="B86">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stevens</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Forrest</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Dudley</surname> <given-names>B. D.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Zeldis</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Shankar</surname> <given-names>U.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Use of a multi-metric macroalgal index to document severe eutrophication in a New Zealand estuary</article-title>. <source>N. Z. J. Mar. Freshwater Res.</source> <volume>56</volume> (<issue>3</issue>), <fpage>410</fpage>&#x2013;<lpage>429</lpage>. doi: <pub-id pub-id-type="doi">10.1080/00288330.2022.2093226</pub-id>
</citation>
</ref>
<ref id="B87">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Aller</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Early diagenesis of chlorophyll-a in Long Island Sound sediments: A measure of carbon flux and particle reworking</article-title>. <source>J. Mar. Res.</source> <volume>49</volume>, <fpage>379</fpage>&#x2013;<lpage>401</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1357/002224091784995927</pub-id>
</citation>
</ref>
<ref id="B88">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname> <given-names>M.-Y.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Aller</surname> <given-names>R. C.</given-names>
</name>
</person-group> (<year>1993</year>). <article-title>Anoxic and oxic degradation of 14C-labeled chloropigments and a 14C-labeled diatom in Long Island Sound sediments</article-title>. <source>Limnology Oceanography</source> <volume>38</volume>, <fpage>1438</fpage>&#x2013;<lpage>1451</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4319/lo.1993.38.7.1438</pub-id>
</citation>
</ref>
<ref id="B89">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sutula</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Green</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Cicchetti</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Detenbeck</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Fong</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Thresholds of adverse effects of Macroalgal abundance and sediment organic matter on benthic habitat quality in estuarine intertidal flats</article-title>. <source>Estuaries Coasts</source> <volume>37</volume>, <fpage>1532</fpage>&#x2013;<lpage>1548</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-014-9796-3</pub-id>
</citation>
</ref>
<ref id="B90">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Swales</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Gibbs</surname> <given-names>M. M.</given-names>
</name>
<name>
<surname>Handley</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Olsen</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Ovenden</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Wadhwa</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). &#x201c;<article-title>Sources of fine sediment and contribution to sedimentation in the inner Pelorus Sound/Te Hoiere</article-title>&#x201d;. (<source>NIWA Client report prepared for Marlborough District Council, 2021291HN, MDC17201</source>.). (<publisher-loc>Christchurch, NZ</publisher-loc>: <publisher-name>NIWA</publisher-name>).</citation>
</ref>
<ref id="B91">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Townsend</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lohrer</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2015</year>). <source>ANZECC Guidance for Estuary Sedimentation</source>. Prepared for Ministry for the Environment, September 2015. NIWA Client Report No HAM2015-096. (Published by <publisher-name>Hamilton, NZ</publisher-name>: <publisher-name>NIWA</publisher-name>). <volume>45</volume>.</citation>
</ref>
<ref id="B92">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Uusitalo</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Advantages and challenges of Bayesian networks in environmental modelling</article-title>. <source>Ecol. Model.</source> <volume>203</volume>, <fpage>312</fpage>&#x2013;<lpage>318</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecolmodel.2006.11.033</pub-id>
</citation>
</ref>
<ref id="B93">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Verburg</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>The need to correct for the Suess effect in the application of d13C as a productivity proxy in sediment of autotrophic Lake Tanganyika in the anthropocene</article-title>. <source>J. Paleolimnology</source> <volume>37</volume>, <fpage>591</fpage>&#x2013;<lpage>602</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10933-006-9056-z</pub-id>
</citation>
</ref>
<ref id="B94">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vieillard</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Newell</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Thrush</surname> <given-names>S. F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Recovering from bias: A call for further study of underrepresented tropical and low-nutrient estuaries</article-title>. <source>J. Geophysical Research: Biogeosciences</source> <volume>125</volume>, <elocation-id>e2020JG005766</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2020JG005766</pub-id>
</citation>
</ref>
<ref id="B95">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Voss</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Larsen</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Leivuori</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Vallius</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Stable isotope signals of eutrophication in Baltic Sea sediments</article-title>. <source>J. Mar. Syst.</source> <volume>25</volume>, <fpage>287</fpage>&#x2013;<lpage>298</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0924-7963(00)00022-1</pub-id>
</citation>
</ref>
<ref id="B96">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Water Framework Directive - United Kingdom Technical Advisory Group</collab>
</person-group> (<year>2014</year>). <source>UKTAG River Assessment Method Macrophytes and Phytobenthos. Macrophytes (River LEAFPACS2)</source>. (<publisher-loc>Stirling, Scotland</publisher-loc>: <publisher-name>United Kingdom Advisory Group</publisher-name>) (WFD-UKTAG).</citation>
</ref>
<ref id="B97">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wells</surname> <given-names>N. S.</given-names>
</name>
<name>
<surname>Eyre</surname> <given-names>B. D.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>&#x3b4;15N patterns in three subtropical estuaries show switch from nitrogen &#x201c;reactors&#x201d; to &#x201c;pipes&#x201d; with increasing degradation</article-title>. <source>Limnology Oceanography</source> <volume>64</volume>, <fpage>860</fpage>&#x2013;<lpage>876</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/lno.11080</pub-id>
</citation>
</ref>
<ref id="B98">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wigand</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Roman</surname> <given-names>C. T.</given-names>
</name>
<name>
<surname>Davey</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Stolt</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Hanson</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2014</year>). <article-title>Below the disappearing marshes of an urban estuary: historic nitrogen trends and soil structure</article-title>. <source>Ecol. Appl.</source> <volume>24</volume>, <fpage>633</fpage>&#x2013;<lpage>649</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/13-0594.1</pub-id>
</citation>
</ref>
<ref id="B99">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wynyard</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2016</year>). <source>The Price of Milk: Primitive accumulation and the New Zealand Dairy Industry 1814-2014</source>. PhD thesis. (<publisher-name>University of Auckland</publisher-name>).</citation>
</ref>
<ref id="B100">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zabarte-Maeztu</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Matheson</surname> <given-names>F. E.</given-names>
</name>
<name>
<surname>Manley-Harris</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Davies-Colley</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Hawes</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Fine sediment effects on seagrasses: A global review, quantitative synthesis and multi-stressor model</article-title>. <source>Mar. Environ. Res.</source> <volume>171</volume>, <fpage>105480</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.marenvres.2021.105480</pub-id>
</citation>
</ref>
<ref id="B101">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeldis</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Depree</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Gongol</surname> <given-names>C.</given-names>
</name>
<name>
<surname>South</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Marriner</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Schiel</surname> <given-names>D. R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Trophic indicators of ecological resilience in a tidal lagoon estuary following wastewater diversion and earthquake disturbance</article-title>. <source>Estuaries Coasts</source> <volume>43</volume>, <fpage>223</fpage>&#x2013;<lpage>239</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12237-019-00637-8</pub-id>
</citation>
</ref>
<ref id="B102">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Zeldis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Hicks</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Trustrum</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Orpin</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Nodder</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Probert</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2010</year>). &#x201c;<article-title>New zealand continental margin fluxes</article-title>,&#x201d; in <source>Carbon and nutrient fluxes in continental margins</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Liu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Atkinson</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Qui&#xf1;ones</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Talaue-Mcmanus</surname> <given-names>L.</given-names>
</name>
</person-group> (<publisher-name>Springer-Verlag</publisher-name>, <publisher-loc>Berlin Heidelberg</publisher-loc>), <fpage>273</fpage>&#x2013;<lpage>287</lpage>.</citation>
</ref>
<ref id="B103">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Zeldis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Measures</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Stevens</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Matheson</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Dudley</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Remediation Options for Southland Estuaries</source> (<publisher-loc>NIWA, Christchurch, NZ</publisher-loc>: <publisher-name>Environment Southland</publisher-name>).</citation>
</ref>
<ref id="B104">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeldis</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D. R.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Predicting and scoring estuary ecological health using a Bayesian belief network</article-title>. <source>Front. Mar. Sci.</source> <volume>9</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmars.2022.898992</pub-id>
</citation>
</ref>
<ref id="B105">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Zeldis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Storey</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Plew</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Whitehead</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Madarasz-Smith</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Oliver</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>) <source>The New Zealand Estuary Trophic Index (ETI) Tools: Tool 3 - Assessing Estuary Trophic State using a Bayesian Belief Network</source>. Available online at: <uri xlink:href="https://shiny.niwa.co.nz/Estuaries-Screening-Tool-3/">https://shiny.niwa.co.nz/Estuaries-Screening-Tool-3/</uri>.</citation>
</ref>
</ref-list>
</back>
</article>