<?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.1323271</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>Deconstructing co-occurring marine heatwave and phytoplankton bloom events in the Arkona Sea in 2018</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Cahill</surname>
<given-names>Bronwyn</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/2256101"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chrysagi</surname>
<given-names>Evridiki</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/771348"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Vortmeyer-Kley</surname>
<given-names>Rahel</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/645484"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gr&#xe4;we</surname>
<given-names>Ulf</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/572585"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Physical Oceanography and Instrumentation, Leibniz Institute for Baltic Sea Research</institution>, <addr-line>Warnem&#xfc;nde</addr-line>, <country>Germany</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Institute for Chemistry and Biology of the Marine Environment (ICBM), Carl von Ossietzky Universit&#xe4;t Oldenburg</institution>, <addr-line>Oldenburg</addr-line>, <country>Germany</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Gemma Kulk, Plymouth Marine Laboratory, United Kingdom</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Yizhen Li, Woods Hole Oceanographic Institution, United States</p>
<p>Ulrike L&#xf6;ptien, University of Kiel, Germany</p>
<p>Inia Soto Ramos, Morgan State University, United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Bronwyn Cahill, <email xlink:href="mailto:bronwyn.cahill@io-warnemuende.de">bronwyn.cahill@io-warnemuende.de</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>03</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>11</volume>
<elocation-id>1323271</elocation-id>
<history>
<date date-type="received">
<day>17</day>
<month>10</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>02</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Cahill, Chrysagi, Vortmeyer-Kley and Gr&#xe4;we</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Cahill, Chrysagi, Vortmeyer-Kley and Gr&#xe4;we</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>Between May and August 2018, two separate marine heatwaves (MHWs) occurred in the Arkona Sea in the western Baltic Sea. These heatwaves bookended an extended period of phytoplankton growth in the region. Data from the Ocean and Land Colour Instrument (OLCI) on board the European Sentinel-3 satellite revealed an eddy-like structure containing high chlorophyll a (Chl-a) concentrations (ca. 25 mg.m<sup>-3</sup>) persisting for several days at the end of May in the Arkona Sea. Combining ocean colour observations, a coupled bio-optical ocean model and a particle tracking model, we examined the three dimensional relationship between these co-occurring MHW and phytoplankton bloom events. We find that the onset of the MHW in May provided the optimal conditions for phytoplankton growth, i.e. sufficient light and nutrients. Wind-driven surface eddy circulation, geostrophic eddy stirring and transient submesoscale dynamics along the edges of the eddy provided a transport path for nutrient fluxes and carbon export, and helped to sustain the phytoplankton bloom. The bloom may have indirectly had an enhancing effect on the MHW, through the impact of water constituent-induced heating rates on air-sea energy fluxes. The subsurface signature of the MHW plays a critical role in de-coupling surface and subsurface dynamics and terminating the phytoplankton bloom. Subsurface temperature anomalies of up to 8&#xb0;C between 15 and 20 m depth are found to persist up to 15 days after the surface signature of the MHW has disappeared. The study reveals how surface and subsurface dynamics of MHWs and phytoplankton blooms are connected under different environmental conditions. It extends our knowledge on surface layer processes obtained from satellite data.</p>
</abstract>
<kwd-group>
<kwd>marine heatwaves</kwd>
<kwd>phytoplankton blooms</kwd>
<kwd>light and nutrient availability</kwd>
<kwd>biooptical modelling</kwd>
<kwd>ocean colour</kwd>
<kwd>mesoscale eddy stirring</kwd>
<kwd>submesoscale dynamics</kwd>
</kwd-group>
<counts>
<fig-count count="13"/>
<table-count count="2"/>
<equation-count count="1"/>
<ref-count count="80"/>
<page-count count="20"/>
<word-count count="9039"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Ocean Observation</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Satellite and ocean data reveal a marked increase in the Earth&#x2019;s heating rate (<xref ref-type="bibr" rid="B54">Loeb et&#xa0;al., 2021</xref>), with the Earth trapping nearly twice as much heat as it did in 2005. This trend is likely to continue in the near term due to global warming and the increase in frequency and magnitude of heat waves (<xref ref-type="bibr" rid="B41">IPCC, 2023</xref>). There is also evidence of a growing trend in the frequency and duration of marine heatwaves (MHWs) in the global oceans (<xref ref-type="bibr" rid="B64">Oliver et&#xa0;al., 2018</xref>). MHWs, defined as periods where the surface temperature of the ocean exceeds the 90th percentile of the 30 year local mean for longer than 5 days, can have severe and destructive consequences on marine species, ecosystems and biogeochemical processes (<xref ref-type="bibr" rid="B72">Smale et&#xa0;al., 2019</xref>). Marginal seas have warmed faster than the global ocean, with the Baltic Sea warming at a rate up to four times the global mean warming rate (<xref ref-type="bibr" rid="B3">Belkin, 2009</xref>). As of 2020, the summer of 2018 was the warmest on instrumental record in Europe, and the warmest summer in the past 30 years in the southern Baltic Sea with surface-water temperatures 4-5&#xb0;C above the 1990-2018 long-term mean (<xref ref-type="bibr" rid="B62">Naumann et&#xa0;al., 2019</xref>) and bottom water temperatures of 20.5&#xb0;C recorded at 32 m at the Tv&#xe4;rminne Zoological Station (TZS) in southern Finland (<xref ref-type="bibr" rid="B39">Humborg et&#xa0;al., 2019</xref>).</p>
<p>A recent statistical analysis by <xref ref-type="bibr" rid="B55">Lorenz (2019)</xref> of historical sea surface temperature (SST) and surface chlorophyll (Chl-a) satellite data products showed that sea surface temperature (SST) and MHW are meaningful parameters for the initiation and development of phytoplankton spring blooms in the Baltic Sea and North Sea, and that there is a relationship between co-occurring MHWs and bloom events. During the last 20 years, a trend towards an earlier spring bloom start has developed which is significantly stronger in the Baltic Sea (<xref ref-type="bibr" rid="B77">Wasmund et&#xa0;al., 2019a</xref>). In addition, there is some evidence of a positive trend in bloom sum and peak, and therefore, magnitude in the Baltic Sea which may be connected to climate variability and eutrophication (<xref ref-type="bibr" rid="B42">Jaanus et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B44">Kahru et&#xa0;al, 2016</xref>; <xref ref-type="bibr" rid="B55">Lorenz, 2019</xref>). Lorenz found significant positive relationships between MHW and bloom indices of magnitude indicating one event may have an enhancing effect on the other by creating a positive feedback.</p>
<p>Spatio-temporal studies of MHW typically examine the response of the sea surface using blended satellite data products (e.g. NOAA OI SST V2, <xref ref-type="bibr" rid="B38">Huang et&#xa0;al., 2021</xref>). While these provide the spatial and temporal surface coverage needed to investigate the frequency and duration of MHWs, they do not provide any diagnostic information on the drivers of MHWs or any prognostic information on surface and subsurface processes which may be impacted by MHWs, i.e. stratification, vertical mixing, light and nutrient availability, phytoplankton growth and optically significant water constituent concentrations. A number of recent studies discuss the subsurface response of marginal and shallow shelf seas to MHWs. An observation-based study by <xref ref-type="bibr" rid="B20">Elzahaby and Schaeffer (2019)</xref> highlight how MHWs occurring in shallow seas (&lt; 150m), occur predominantly during the stratified season (summer/autumn) in various mesoscale structures (cyclonic, anticyclonic or no eddies). They are characterized by stratified and fresher surface waters and the depth to which they extend is correlated with the SST anomaly. The surface origin is likely a response to air-sea flux forcing whereby anomalous solar radiation and decreased wind stress act on the latent heat-flux to reduce evaporation (<xref ref-type="bibr" rid="B17">Chen et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B7">Bond et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B16">Chen et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B4">Benthuysen et&#xa0;al., 2018</xref>). <xref ref-type="bibr" rid="B71">Schaeffer and Roughan (2017)</xref> show that MHWs regularly extend over the full depth of the water column in coastal waters off southeastern Australia, and their maximum intensity occurs below the surface and can persist long after the surface signature of the MHW has disappeared. <xref ref-type="bibr" rid="B32">Hayashida et&#xa0;al. (2020)</xref> use daily output from a near-global ocean physical-biogeochemical model to explore how background nutrient concentrations determine the response of co-occurring phytoplankton blooms and MHWs in regional seas.</p>    <p>In this paper, we examine two different late spring and mid-summer MHW events which bookmark the evolution and decay of a phytoplankton bloom in the Arkona Sea in the Western Baltic Sea. We use a combination of satellite data, a coupled bio-optical ocean model and a particle tracking model, in order to better understand the full three-dimensional impact MHWs have on bloom dynamics. Our objective is to explore the relationship between co-occurring MHW and phytoplankton bloom events with the following specific questions in mind:</p>
<list list-type="simple">
<list-item>
<p>1. How, and under which circumstances, do MHWs contribute to the initiation of phytoplankton blooms?</p>
</list-item>
<list-item>
<p>2. Which dynamics play a role in sustaining the bloom?</p>
</list-item>
<list-item>
<p>3. How deep and for how long is the impact of the MHW felt?</p>
</list-item>
<list-item>
<p>4. What role, if any, do MHWs play in terminating a phytoplankton bloom?</p>
</list-item>
<list-item>
<p>5. Do phytoplankton blooms have an enhancing effect on MHWs by creating a positive feedback from water constituent&#x2013;induced surface heating?</p>
</list-item>
</list>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Study site</title>    <p>Our study site is the Arkona Sea located in the Western Baltic Sea (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). This site was selected because two significant MHW events took place here in May and July 2018 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>). During the peak surface heating period of the May event, the sea surface temperature (SST) anomaly on 29 May 2018 shows a distinct eddy-like structure of warmer water (up to 6&#xb0;C) in the Arkona Sea (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3A</bold>
</xref>). Coincident satellite data from the Ocean and Land Colour Instrument (OLCI) on board the European Sentinel-3 series satellites also shows a similar eddy-like structure containing high Chl-a concentrations (c. 25 mg m<sup>-3</sup>) in the Arkona Sea (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4A</bold>
</xref>). During the July MHW, surface Chl-a concentrations were much lower (c. 3 mg m<sup>-3</sup>), with no distinct structure visible in the satellite data (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4H</bold>
</xref>). Prevailing winds show weak northeasterly and southeasterly winds converging along 55&#xb0; N during the May event, with wind speeds on the order of 5 m s<sup>-1</sup>, while in July, prevailing winds were ca. 2 m s<sup>-1</sup> with northerly winds west of 13&#xb0;30&#x2019; E, and southerly east of 13&#xb0;30&#x2019; E (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5A, B</bold>
</xref>). The Arkona Sea is characterized by a basin which covers an area of approximately 18,700 km<sup>2</sup>, and has a maximum depth of 47 m. A semi-permanent halocline separates the fresher surface water (6 &#x2013; 8 PSU) from the more saline deep water (12 &#x2013; 14 PSU) between 20 and 40 m depth. Large variability in baroclinic circulation caused by the imbalance between filling and emptying the dense bottom water pool results in a range of reported Rossby radii (2.5 to 6 km) and baroclinic phase velocities (0.27 to 0.75 ms<sup>-1</sup>) for the first baroclinic mode (<xref ref-type="bibr" rid="B23">Fennel et&#xa0;al., 1991</xref>; <xref ref-type="bibr" rid="B49">Lass and Mohrholz, 2003</xref>). Seasonal hypoxia can also occur. The mean stratification can be disturbed by geostrophic eddies with a characteristic radius equal to the first baroclinic Rossby radius and up- and downwelling occurring along the rim of the eddy. <xref ref-type="bibr" rid="B75">Vortmeyer-Kley et&#xa0;al. (2019)</xref> found in a modelling case study of the surface velocity field in the Western Baltic Sea from May 1 to October 31, 2010 about 28 000 eddies, while <xref ref-type="bibr" rid="B68">Rei&#xdf;mann (2005)</xref> detected in a CTD measurement campaign 5 to 18 eddies as three dimensional isolated anomalies in pressure in the Arkona Basin in October 1999. These structures are found at a mean depth of about 15 - 22 m in <xref ref-type="bibr" rid="B68">Rei&#xdf;mann (2005)</xref>. In their satellite image study for 2009 - 2011, <xref ref-type="bibr" rid="B45">Karimova and Gade (2016)</xref> found the Arkona Basin as region of eddies that might be caused by sharp thermal gradients. Nonlinear eddies are important for biological production because they trap fluid, phytoplankton and nutrients within them (<xref ref-type="bibr" rid="B15">Chelton et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B14">Chelton et&#xa0;al., 2011a</xref>, <xref ref-type="bibr" rid="B13">Chelton et&#xa0;al., 2011b</xref>; <xref ref-type="bibr" rid="B58">McGillicuddy, 2016</xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>
<bold>(A)</bold> Global distribution of chlorophyll-a, as seen by MODIS, May 2018 (source: <uri xlink:href="https://earthobservatory.nasa.gov/global-maps/MY1DMM_CHLORA">https://earthobservatory.nasa.gov/global-maps/MY1DMM_CHLORA</uri>). Black rectangle indicates location of our study region. <bold>(B)</bold> Western Baltic Sea model domain bathymetry (m). Black rectangle shows the location of the Arkona Sea, the red dot indicates the location of Marnet Arkona Buoy long term mooring site (13&#xb0;52&#x2019; E; 54&#xb0;53&#x2019; N) and the dashed lines show the location of the transects used in the analysis.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g001.tif"/>
</fig>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>    <p>
<bold>(A)</bold> modelled 2018 temperature anomaly at Marnet Arkona Buoy location (see <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>), 2018 temperature minus 90<sup>th</sup> percentile 40 year climatology (1979 &#x2013; 2019); <bold>(B)</bold> surface temperature at Marnet Arkona Buoy location: modelled (2018 ROMS), observed (2018 Marnet Buoy) and 30 year mean climatology (CLIM) and 90<sup>th</sup> percentile (90<sup>th</sup>) [using the NOAA OI SST V2 High Resolution Data Set (<xref ref-type="bibr" rid="B38">Huang et&#xa0;al., 2021</xref>)] (ROMS vs Marnet Buoy statistics: r<sup>2</sup>: 0.99, RMSE: 0.016, BIAS: -0.0010). Pink shaded areas indicate timing of MHW-1 and MHW-2, respectively; <bold>(C)</bold> modelled 2018 surface Chl-a (blue) and optically significant water constituent-induced surface heating rate (orange) at the Marnet Arkona Buoy location (dashed horizontal line indicates threshold value for onset of phytoplankton bloom (see text); <bold>(D)</bold> modelled 2018 water column Chl-a at the Marnet Arkona Buoy location. (grey dashed vertical lines indicate onset and end of MHW events, solid grey vertical line indicates day of detailed analysis during MHW event, green vertical dashed lines indicate start and end of phytoplankton blooms).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>SST anomaly on 29 May 2018 <bold>(A)</bold> and 26 July 2018 <bold>(B)</bold> in the Western Baltic Sea; <bold>(C)</bold> cumulative contribution of the air-sea heat flux (blue line), horizontal advection (dashed green line) to the temperature anomaly and the mixed layer depth temperature anomaly (orange line) in the Arkona Sea in 2018. The SST Reanalysis (<uri xlink:href="https://doi.org/10.48670/moi-00156">https://doi.org/10.48670/moi-00156</uri>) described in <xref ref-type="bibr" rid="B37">H&#xf8;yer and She (2007)</xref> and <xref ref-type="bibr" rid="B36">H&#xf8;yer and Karagali (2016)</xref> was used to calculate the SST anomalies shown in <bold>(A, B)</bold>, while the Baltic Sea Physics Reanalysis (<uri xlink:href="https://doi.org/10.48670/moi-00013">https://doi.org/10.48670/moi-00013</uri>) and the ERA-5 Global Reanalysis (<uri xlink:href="https://doi.org/10.1002/qj.3803">https://doi.org/10.1002/qj.3803</uri>) described in <xref ref-type="bibr" rid="B33">Hersbach et&#xa0;al. (2020)</xref> were used to calculate the temperature anomalies shown in <bold>(C)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g003.tif"/>
</fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Sentinel-3 Ocean Land Cover Instrument (OLCI) L3 300 m resolution Chl-a on selected dates during the study period <bold>(A-H)</bold> 29, 30 May, 2, 6, 7, 9, 30 June, 26 July 2018, respectively). Black dot indicates the location of the Marnet Arkona Buoy long term mooring site.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g004.tif"/>
</fig>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Rossby number at 1m and 15m on May 29 and July 26, 2018 <bold>(A&#x2013;D)</bold>. DWD-ICON 3-hourly surface wind vectors are overlaid on <bold>(A, B)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g005.tif"/>
</fig>
<p>The Arkona Sea is an optically complex and biologically productive region, influenced by significant inputs of terrestrial organic matter from neighbouring rivers, especially colour dissolved organic matter (CDOM) during months of intensive mixing and high riverine discharge, March, April and November (<xref ref-type="bibr" rid="B46">Kowalczuk et&#xa0;al., 2006</xref>). In recent years, the growing season of phytoplankton in the Western Baltic Sea has extended from 159 days in the period 1988 to 1992 to 284 days in the period 2014 to 2017 (<xref ref-type="bibr" rid="B77">Wasmund et&#xa0;al., 2019a</xref>) in response to climate change. In 2018, <xref ref-type="bibr" rid="B76">Wasmund et&#xa0;al. (2019b)</xref> observed a shift in the spring bloom peak in the region to May with a prolonged period of moderate phytoplankton growth in the Arkona Sea primarily dominated by diatoms and dinoflagellates in May and June 2018.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Modelling the bio-optical ocean state in the Arkona Sea in 2018</title>
<p>We use the Regional Ocean Modelling System, ROMS, which drives the physics and the advection and diffusion of tracers, coupled with the Ecosim/Bio-Optic module (herein referred to as ROMS-BioOptic) which drives the ecosystem and a spectrally-resolved underwater light field. This setup is used to simulate the bio-optical ocean state in the western Baltic Sea for the year 2018 and is described in detail in <xref ref-type="bibr" rid="B9">Cahill et&#xa0;al. (2023)</xref> and references therein. Here we summarize important features. ROMS, is widely used for shelf circulation (e.g. <xref ref-type="bibr" rid="B31">Haidvogel et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B78">Wilkin et&#xa0;al., 2011</xref>) and coupled physical-biological applications (e.g. <xref ref-type="bibr" rid="B25">Fennel et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B10">Cahill et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B26">Fennel et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B24">Fennel and Wilkin, 2009</xref>; <xref ref-type="bibr" rid="B11">Cahill et&#xa0;al., 2016</xref>, <xref ref-type="bibr" rid="B9">Cahill et&#xa0;al., 2023</xref>). Ecosim is a carbon-based, ecological/optical modelling system (<xref ref-type="bibr" rid="B5">Bissett et&#xa0;al., 1999a</xref>, <xref ref-type="bibr" rid="B6">Bissett et&#xa0;al., 1999b</xref>) which was developed for simulations of carbon cycling and biological productivity. Ecosim simulates up to four phytoplankton functional groups each with a characteristic pigment suite which varies with the group carbon-to-chlorophyll-a ratio, C:Chl-a. Each groups&#x2019; C:Chl-a ratio varies between some maximum and minimum value, as a function of light or nutrient limitation. The properties of each functional group evolve over time as a function of light and nutrient conditions (i.e. NO<sub>3</sub>, NH<sub>4</sub>, PO<sub>4</sub>, SiO and FeO). The maximum phytoplankton growth is modulated by temperature (<xref ref-type="bibr" rid="B22">Eppley, 1972</xref>). Loss processes are represented by grazing and excretion. Grazing accounts for the majority of the biomass sink in the model and is considered the closure term of the phytoplankton equations (<xref ref-type="bibr" rid="B73">Steele and Henderson, 1992</xref>). It is modelled as a Michaelis-Menten function based on the functional groups&#x2019; biomass (<xref ref-type="bibr" rid="B5">Bissett et&#xa0;al., 1999a</xref>). Marine and riverine sources of dissolved organic carbon (DOC and CDOC) are accounted for and explicitly resolved into labile (e.g. available for biological and photo-degradation) and relict (e.g. available for photo-degradation) forms. Dissolved inorganic carbon (DIC) is also accounted for. Riverine sources of carbon and nutrients are introduced via point sources. The underwater light field is spectrally-resolved at 5nm intervals between 400 and 700 nm. This allows for differential growth of different phytoplankton groups that have unique pigment complements.</p>
<p>Ecosim&#x2019;s daylight module explicitly calculates the in-water spectrally-resolved absorption coefficients for phytoplankton, detritus and CDOM, the scattering and backscattering coefficients for phytoplankton and detritus, the average cosine, downwelling irradiance attenuation coefficient, K<sub>d</sub>, in addition to the scalar, E<sub>0</sub>, and downward, E<sub>d</sub>, irradiance fields following <xref ref-type="bibr" rid="B61">Morel (1988)</xref>. <xref ref-type="bibr" rid="B9">Cahill et&#xa0;al. (2023)</xref> recently updated the K<sub>d</sub> formulation following <xref ref-type="bibr" rid="B50">Lee et&#xa0;al., 2005</xref> which accounts for some of the optical complexity found in coastal waters. The spectrally-resolved underwater light field drives the evolution of all the water constituents in the ecosystem model (phytoplankton, detritus and CDOM), while the water constituents in turn determine the evolution of the light field in each layer by absorption and scattering of the light. This means that their contribution to the divergence of the heat flux (<xref ref-type="bibr" rid="B61">Morel, 1988</xref>) can be accounted for within the full hydrodynamic solution. Furthermore, water constituent-induced heating rates can be assessed (<xref ref-type="bibr" rid="B9">Cahill et&#xa0;al., 2023</xref>) and their impact on the ocean sea surface temperature can be communicated to the bulk flux formulation of the atmosphere in the modelling system.</p>    <p>The ROMS-BioOptic model was configured as described in <xref ref-type="bibr" rid="B9">Cahill et&#xa0;al., 2023</xref> except with a higher resolution 600m horizontal grid in the Western Baltic Sea (see <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). A bulk flux atmosphere was forced with DWD-ICON output (<xref ref-type="bibr" rid="B79">Z&#xe4;ngl et&#xa0;al., 2015</xref>) and river forcing including runoff and biogeochemistry (NO<sub>3</sub>, NH<sub>4</sub>, PO<sub>4</sub>, SiO, DOC, CDOM and DIC) from 9 rivers which influence the region was derived from HELCOM PLC (Pollution Load Compilation) data (Neumann, pers. comm). Open boundaries to the north and east were forced with output from GETM physics using a combination of Chapman/Flather conditions for u and v velocities and transports, and Radiation + Nudging for temperature and salinity. This 3D setup is based on an existing GETM physics setup which has been previously evaluated and published (<xref ref-type="bibr" rid="B29">Gr&#xe4;we et&#xa0;al., 2015a</xref>, <xref ref-type="bibr" rid="B30">Gr&#xe4;we et&#xa0;al., 2015b</xref>). The Ecosim/Bio-Optic module was configured with four phytoplankton functional groups representative of small and large diatoms, large dinoflagellates and cyanobacteria. Initial conditions for the ecosystem model were obtained from previously evaluated and published output from the Ecological Regional Ocean Model (ERGOM) (<xref ref-type="bibr" rid="B63">Neumann et&#xa0;al., 2022</xref>). Our simulation period was 1 January to 31 December 2018. Daily averages and snapshots were output for the entire year. Hourly averages and snapshots were output for selected periods during both MHW events.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Configuration of ROMS-BioOptic Western Baltic Sea Application.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" colspan="2" align="left">ROMS-BioOptic</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">
<bold>Application Name</bold>
</td>
<td valign="middle" align="left">3D Western Baltic Sea</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Model Grid</bold>
</td>
<td valign="middle" align="left">853 x 505 (600 m), 30 sigma levels</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Simulation Period</bold>
</td>
<td valign="middle" align="left">2018</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Boundary Conditions</bold>
</td>
<td valign="middle" align="left">Chapman for zeta, Flather for ubar and vbar; Radiation + Nudging for temperature and salinity</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Bulk Flux Atmosphere</bold>
</td>
<td valign="middle" align="left">DWD-ICON 3-hourly</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>River Forcing</bold>
</td>
<td valign="middle" align="left">HELCOM PLC (Pollution Load Compilation), Neumann (pers. comm.)</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Initial Conditions</bold>
</td>
<td valign="middle" align="left">GETM/ERGOM</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Time Step</bold>
</td>
<td valign="middle" align="left">DT = 30s; NDTFAST = 20s</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Ecosim</bold>
</td>
<td valign="middle" align="left">4 phytoplankton groups (small and large diatoms, large dinoflagellates &amp; cyanobacteria)</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Spectral Resolution</bold>
</td>
<td valign="middle" align="left">5 nm intervals between 400 and 700 nm</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Detecting MHW events</title>
<p>Following <xref ref-type="bibr" rid="B34">Hobday et&#xa0;al. (2018)</xref>, surface temperature data from the BSH (Bundesamt f&#xfc;r Seeschifffahrt und Hydrographie) MARNET Arkona Buoy, surface temperature output from our 600m ROMS-BioOptic simulation of the Western Baltic Sea and sea surface temperature data from the NOAA OI SST V2 High Resolution Dataset (<xref ref-type="bibr" rid="B38">Huang et&#xa0;al., 2021</xref>) were used to diagnose the May and July 2018 MHW events (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>; <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). The heat budget diagnostic approach summarized in <xref ref-type="bibr" rid="B21">Elzahaby et&#xa0;al. (2021)</xref> and based on <xref ref-type="bibr" rid="B16">Chen et&#xa0;al. (2015)</xref> and <xref ref-type="bibr" rid="B8">Bowen et&#xa0;al. (2017)</xref> was used to diagnose whether the MHWs were atmospheric-driven or oceanic-process driven (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). Herein, the respective contributions of advection and air-sea heat flux anomalies to the mixed layer temperature tendency anomaly during each event are used to classify the drivers of MHW drivers. For this purpose, we used the Baltic Sea Physics Reanalysis (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.48670/moi-00013">https://doi.org/10.48670/moi-00013</ext-link>) and the ERA-5 Global Reanalysis (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1002/qj.3803">https://doi.org/10.1002/qj.3803</ext-link>) described in <xref ref-type="bibr" rid="B33">Hersbach et&#xa0;al. (2020)</xref> to calculate the temperature anomalies.</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Summary of MHW indices in the Arkona Sea in May and July 2018 (after <xref ref-type="bibr" rid="B34">Hobday et&#xa0;al., 2018</xref>).</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">MHW #</th>
<th valign="middle" align="center">Start</th>
<th valign="middle" align="center">Day</th>
<th valign="middle" align="center">End</th>
<th valign="middle" align="center">Day</th>
<th valign="middle" align="center"># Days</th>
<th valign="middle" align="center">I<sub>mean</sub>
</th>
<th valign="middle" align="center">I<sub>max</sub>
</th>
<th valign="middle" align="center">Category</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">
<bold>MHW-1</bold>
</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">13/05/2018</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">133</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">19/06/2018</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">170</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">38</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">4.0&#xb0;C</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">5.3&#xb0;C</td>
<td valign="middle" align="center" style="background-color:#ffffff">23 days moderate</td>
</tr>
<tr>
<td valign="middle" align="center" style="background-color:#ffffff">15 days strong</td>
</tr>
<tr>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">
<bold>MHW-2</bold>
</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">25/07/2018</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">206</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">09/08/2018</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">221</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">16</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">4.4&#xb0;C</td>
<td valign="middle" rowspan="2" align="center" style="background-color:#ffffff">5.2&#xb0;C</td>
<td valign="middle" align="center" style="background-color:#ffffff">15 days moderate</td>
</tr>
<tr>
<td valign="middle" align="center" style="background-color:#ffffff">1 day strong</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Detecting phytoplankton blooms events</title>
<p>We define a bloom event following a threshold approach as used by <xref ref-type="bibr" rid="B74">Thomalla et&#xa0;al. (2011)</xref> and applied by <xref ref-type="bibr" rid="B55">Lorenz (2019)</xref> whereby the initiation of the bloom is understood to be the period of the year which registers a relative increase in chlorophyll concentration, irrelevant of the actual value. For our purposes, we define the threshold as the first day that surface Chl-a rises 15% above the annual median as follows:</p>
<p>CHL<sub>S</sub> = CHL<sub>MEDIAN</sub> + 0.15 * CHL<sub>MEDIAN</sub>
</p>
<p>Thus, the bloom start condition is:</p>
<p>CHL<sub>t</sub> &lt; CHL<sub>S</sub> and CHL<sub>t</sub>+1 &gt; CHL<sub>S</sub>
</p>
<p>and, the bloom end condition is:</p>
<p>CHL<sub>t</sub> &gt; CHL<sub>S</sub> and CHL<sub>t+1</sub> &lt; CHL<sub>S</sub>
</p>
<p>In the run up to a bloom event, Chl-a tends to pulsate, often exceeding the threshold CHL<sub>S</sub> for a short period of time (<xref ref-type="bibr" rid="B67">Racault et&#xa0;al., 2015</xref>). Therefore, events which are shorter than 7 days are not considered blooms and blooms which are less than 6 days apart are considered as one event.</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Detecting and tracking eddy-like finite-time coherent structures</title>
<p>Building on the ideas applied by <xref ref-type="bibr" rid="B75">Vortmeyer-Kley et&#xa0;al. (2019)</xref>, we use our modelled velocity output to search for three dimensional, eddy-like coherent structures in the Arkona Sea during the time of the MHW events. To do this, we apply a particle tracking method combined with the concept of finite-time coherent sets by <xref ref-type="bibr" rid="B27">Froyland and Junge (2018)</xref>. In general, eddies can be described as separated waterbodies that minimally mix with or leak into their neighbourhood. This links the idea of eddies to the concept of finite-time coherent sets (<xref ref-type="bibr" rid="B27">Froyland and Junge, 2018</xref>; <xref ref-type="bibr" rid="B28">Froyland et&#xa0;al., 2019</xref>).</p>
<p>
<xref ref-type="bibr" rid="B27">Froyland and Junge (2018)</xref> and <xref ref-type="bibr" rid="B28">Froyland et&#xa0;al. (2019)</xref> formulate the definition of a finite-time coherent set as the solution of the weak eigenproblem of the dynamic Laplacian <xref ref-type="disp-formula" rid="eq1">Equation 1</xref>:</p>
<disp-formula id="eq1">
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
</mml:munder>
<mml:msup>
<mml:mi>D</mml:mi>
<mml:mi>t</mml:mi>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>with u as set of eigenvectors, whose information correspond to the finite-time coherent set, &#x3bb; as corresponding eigenvalues, D as the stiffness matrix and M as the mass matrix. The entries in the mass matrix M can be interpreted as the volume of all tetrahedrons built from a three dimensional triangulation from the tracer position in space at time t if the tracers are advected by a flow field for the time period &#x3c4;. The entries in the stiffness matrix D correspond to the change of the shape of all tetrahedrons (details cf. <xref ref-type="bibr" rid="B27">Froyland and Junge, 2018</xref>). The information of the coherent sets in the eigenvectors is often not cleanly represented in the eigenvectors, so we use the SEBA algorithm (<xref ref-type="bibr" rid="B28">Froyland et&#xa0;al., 2019</xref>) to disentangle the information and save them into SEBA vectors whose entries represent the coherence level of the finite-time coherent sets.</p>
<p>The backbone of the above mentioned algorithm is a particle tracking that provides the trajectory calculations needed to describe the tracer positions in space in the time period &#x3c4; properly. Here we use ROMSpath (<xref ref-type="bibr" rid="B40">Hunter et&#xa0;al., 2022</xref>; code at <ext-link ext-link-type="uri" xlink:href="https://github.com/imcslatte/ROMSPath/tree/V1.0.0">https://github.com/imcslatte/ROMSPath/tree/V1.0.0</ext-link>) and integrate massless tracer trajectories on a longitude-latitude-depth grid for 24h starting tracers from the same grid every hour in the period May 27, 2018 00:30am to June 2, 2018 01:30am and July 22, 2018 05:30am to July 29, 2018 11:30pm using the velocity fields from ROMS-BioOptic simulations.</p>
<p>We apply the approach by <xref ref-type="bibr" rid="B27">Froyland and Junge (2018)</xref> and <xref ref-type="bibr" rid="B28">Froyland et&#xa0;al. (2019)</xref> using their matlab scripts (<ext-link ext-link-type="uri" xlink:href="https://github.com/gaioguy/FEMDL">https://github.com/gaioguy/FEMDL</ext-link> and <ext-link ext-link-type="uri" xlink:href="https://github.com/gfroyland/SEBA">https://github.com/gfroyland/SEBA</ext-link>) to calculate eddy-like finite-time coherent sets for each hour from the calculated tracer trajectories. We requested 40 eigenvectors and chose the first 25 SEBA vectors to find the largest and most coherent sets. From this assemblage of hourly detected eddy-like coherent sets, we build consecutive tracks of coherent sets by searching in the assemblage for sets of the same type (positive or negative 24h-mean relative vorticity in the centre of the set) that are spatially close to each other in successive time steps. (Spatially close to each other means that their inner product is larger than 0.75.) We define the 0.7 coherence level as the outer shapes of the eddy-like coherent sets in space. We take into account tracks of structures with a lifetime larger than 11 h, to consider only structures that live long enough to have an ecological impact.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Characteristics of the May and July MHW events</title>
<p>Following the approach described in section 2.3, two MHW events were identified using observations from the Marnet Arkona Buoy in 2018 and the blended NOAA OI SST V2 product from <xref ref-type="bibr" rid="B38">Huang et&#xa0;al. (2021)</xref> (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>). Coincident modelled SST was compared to the Marnet Arkona Buoy data, to ensure the fitness for purpose of the physical model for the analysis of MHW surface and subsurface dynamics. The spatial extent of the MHWs on 29 May and 26 July 2018 was calculated as the SST anomaly using the SST Reanalysis (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.48670/moi-00156">https://doi.org/10.48670/moi-00156</ext-link>) described in <xref ref-type="bibr" rid="B37">H&#xf8;yer and She (2007)</xref> and <xref ref-type="bibr" rid="B36">H&#xf8;yer and Karagali (2016)</xref> (<xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3A, B</bold>
</xref>). The MHW events are described in terms of their duration, mean intensity, I<sub>mean</sub>, maximum intensity, I<sub>max</sub>, and category (moderate, strong, severe or extreme) (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). The May MHW (herein referred to as MHW-1) lasted 38 days, in which 23 are classified as moderate and 15 as strong. MHW-1 I<sub>mean</sub> is 4.0&#xb0;C and I<sub>max</sub> is 5.3&#xb0;C. The July MHW (herein referred to as MHW-2) lasted 17 days, in which 15 are classified as moderate and 2 as strong. MHW-2 I<sub>mean</sub> is 4.4&#xb0;C and I<sub>max</sub> is 5.2&#xb0;C. The heat budget diagnostic (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>) reveals that MHW-1 is an atmospheric-driven event, while MHW-2 is a mixed atmospheric- and horizontal advection-driven event.</p>
<p>There is very good agreement between our modelled ROMS-BioOptic sea surface temperature at Arkona Sea in 2018 and the observed sea surface temperature in 2018 (r<sup>2 =</sup> 0.99, RMSE = 0.016, Bias = -0.001) (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>). The timing and surface characteristics of both MHW events are also captured very well in the model. This gives us confidence to examine the sub-surface properties of the 2018 modelled temperature anomaly using a 40 year modelled climatology (1979 &#x2013; 2019) derived from General Estuarine Transport Model (GETM) simulations (<xref ref-type="bibr" rid="B29">Gr&#xe4;we et&#xa0;al., 2015a</xref>) (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>).</p>
<p>The impact of MHW-1 extends to ca. 15 m during the first half of the event, after which it deepens over the full extent of the water column (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>). The maximum temperature anomaly (ca. 6&#xb0;C) occurs a few days after the surface signature of MHW-1 has disappeared at depths between 15 and 25m. This maximum subsurface anomaly persists for about 3 days, after which it relaxes to ca. 4&#xb0;C, but remains positive until the onset of MHW-2.</p>
<p>The subsurface impact of MHW-2 is much more pronounced compared to MHW-1. MHW-2 extends from the surface to ca. 15 m for the duration of the event with a temperature anomaly of ca. 4&#xb0;C in the surface layer. Between 15 m and 30 m, a negative temperature anomaly (ca. 2&#xb0;C) lies over a positive bottom temperature anomaly (ca. 3&#xb0;C). The maximum temperature anomaly (ca. 8&#xb0;C) is concentrated between 15 m and 20m for 15 days after the surface signature of MHW-2 has disappeared. Its extent then deepens from 15 m to about 30 m for further 10 days, following which it decays over a period of about 10 days, finally relaxing at the end of September.</p>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Dynamics of phytoplankton blooms and MHWs</title>
<p>Following the threshold approach described in section 2.4, three distinct phytoplankton bloom periods were identified in the modelled Chl-a in 2018 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2C</bold>
</xref>). The first bloom developed on 29<sup>th</sup> March and persisted until 1<sup>st</sup> May. Peak Chl-a concentration of 4.95 mg m<sup>-3</sup> occurred on 14<sup>th</sup> April. A second bloom developed on 20<sup>th</sup> May, 7 days after the onset of MHW-1. This persisted until 24<sup>th</sup> July, one day before the onset of MHW-2. Peak Chl-a concentration of 5.4 mg m<sup>-3</sup> occurred on 29<sup>th</sup> May and coincided with maximum water constituent-induced surface heating rates (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2C</bold>
</xref>). Modelled phytoplankton, CDOM and detrital absorption at 442 nm show that phytoplankton absorption dominates the diffuse attenuation coefficient at 492 nm during the second bloom event (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6C&#x2013;F</bold>
</xref>) and thus contributes most to the surface heating rates. A third bloom event developed on 12<sup>th</sup> August, two days after the end of MHW-2 event and persisted until 14<sup>th</sup> September. Peak Chl-a concentration of 3.7 mg m<sup>-3</sup> occurred on 4<sup>th</sup> September. We focus our attention on the evolution and decay of the 2<sup>nd</sup> phytoplankton bloom event which is bookended by the MHW-1 and MHW-2.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>
<bold>(A)</bold> Temperature, <bold>(B)</bold> salinity, <bold>(C-E)</bold> phytoplankton, detrital and CDOM absorption at 442 nm, respectively, and <bold>(F)</bold> the diffuse attenuation coefficient at 492 nm [K<sub>d</sub>(492)] in Arkona Sea 2018. (Grey dashed vertical lines indicate onset and end of two MHW events, solid grey vertical line indicates day of detailed analysis during MHW event, green vertical dashed lines indicate start and end of the prolonged phytoplankton bloom event discussed in the text.).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g006.tif"/>
</fig>
<p>Our analysis is centered on two separate days within the MHW-1 and MHW-2 events, 29 May and 26 July 2018. These days were selected for a number of reasons: they exhibit different bio-optical and biogeochemical responses to the MHW events; 29 May 2018 coincides with the maximum surface Chl-a concentrations and peak water constituent-induced heating rates during MHW-1 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2C</bold>
</xref>); 26 July 2018 coincides with the end of the 2<sup>nd</sup> phytoplankton bloom event, a steady decline in the water constituent-induced surface heating rates and the onset of MHW-2. A sequence of selected cloud-free Chl-a OLCI satellite data (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>) starting from 29 May 2018 show a surface Chl-a signature persisting within the eddy-like structure for at least 10 days. As will be shown below, a similar eddy pattern was found in the model results in the same location at the end of May, although we do not expect the positions of the observed and simulated eddies to agree completely, because trajectories of eddies typically contain a stochastic element.</p>
<p>We evaluated the simulated surface Chl-a with coincident OLCI data. Ocean colour instruments receive most of their in-water signal from the surface down to one optical depth. When we refer to the simulated surface chlorophyll, we actually refer to the mean of the simulated chlorophyll over the first optical depth, rendering such satellite-derived and simulated chlorophyll concentrations comparable. In our area of interest, the OLCI data characterize the diffuse attenuation coefficient at 490 nm, K<sub>d</sub>(490), to be about 0.5 m<sup>-1</sup> on 29 May 2018 and 0.25 m<sup>-1</sup> on 26 July 2018, resulting in a remotely sensed layer at 490 nm of about 2 m and 4 m, respectively. The evaluation shows that the model captures the surface dynamics of the phytoplankton bloom but highlights the difficulty in capturing the magnitude of the surface bloom event that is seen in the satellite data in May in the model (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figures&#xa0;1A-D</bold>
</xref>). However, we would not expect the model to necessarily reproduce the magnitude of the bloom event observed by OLCI in May as these type of events are very difficult to reproduce in a model without data assimilation. Good agreement is seen between the model and OLCI data in July. Previous evaluations of modelled Chl-a and other water constituents using OLCI data also found good agreement between the model and OLCI data background values in the region (<xref ref-type="bibr" rid="B9">Cahill et&#xa0;al., 2023</xref>). The structure of the simulated surface Chl-a concentrations captures the extent of the eddy feature observed by OLCI in May (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;1A</bold>
</xref>). This supports our application of the model to explore the relationship between co-occurring MHWs and phytoplankton bloom events.</p>
<p>In order to explore the physical differences between the two MHWs further, we examine horizontal cross sections across the area of interest at 1 m and 15 m, and vertical transects across 13&#xb0;58&#x2019; E and 55&#xb0;8&#x2019; N. We use the Rossby number, Ro (defined as balance between the vertical component of the relative vorticity and planetary vorticity) to situate the flow regime (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5</bold>
</xref>, <xref ref-type="fig" rid="f7">
<bold>7</bold>
</xref>). This will be &lt;&lt; 1 in mesoscale regimes, where planetary rotation constrains the flow and vertical stratification dominates, and O(1) in submesoscale flow regimes, where relative vorticity becomes important. Vertical sections of the horizontal (u and v) and vertical (w) velocity components (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>) provide more information on the structure of the flow.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Rossby number along 13&#xb0;58&#x2019; E (TE) and 55&#xb0;8&#x2019; N (TN) on May 29 and July 26, 2018 <bold>(A&#x2013;D)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g007.tif"/>
</fig>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>U, V, W velocities along 13&#xb0;58&#x2019; E (TE) and 55&#xb0;8&#x2019; N (TN) on May 29 and July 26, 2018.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g008.tif"/>
</fig>
<sec id="s3_2_1">
<label>3.2.1</label>
<title>29<sup>th</sup> May 2018 (MHW-1)</title>
<p>Horizontal cross sections of Ro at 1 m and 15 m (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5A, C</bold>
</xref>) show a cyclonic eddy structure centered around 13&#xb0;58&#x2019; E and 55&#xb0;8&#x2019; N present on 29 May 2018. This structure extends from the surface to approximately 27 m depth (<xref ref-type="fig" rid="f7">
<bold>Figures&#xa0;7A, C</bold>
</xref>). The cyclonic nature of the structure is clear in the horizontal velocity components (<xref ref-type="fig" rid="f8">
<bold>Figures&#xa0;8A, C</bold>
</xref>) and relatively strong vertical velocities on the order of 10<sup>-4</sup> ms<sup>-1</sup> are seen on the northern and western flanks of the structure (<xref ref-type="fig" rid="f8">
<bold>Figures&#xa0;8E, G</bold>
</xref>) where Ro is +/- 0.8 (<xref ref-type="fig" rid="f7">
<bold>Figures&#xa0;7A, C</bold>
</xref>). The strong vertical velocities coincide with sharp lateral density gradients, a doming of isopycnals and upwelling of cooler subsurface water (<xref ref-type="fig" rid="f9">
<bold>Figures&#xa0;9A, C</bold>
</xref>). Converging north-southeasterly winds prevailed during this time (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5A</bold>
</xref>).</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>
<bold>(A-D)</bold> Temperature (&#xb0;C), <bold>(E-H)</bold> K<sub>d</sub>(492) and <bold>(I-L)</bold> Chl-a along 13&#xb0;58&#x2019; E (TE) and 55&#xb0;8&#x2019; N (TN) on May 29 and July 26, 2018. Density contours are plotted as grey lines.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g009.tif"/>
</fig>
<p>Horizontal cross sections of the diffuse attenuation coefficient at 492 nm, K<sub>d</sub>(492) and Chl-a concentration at 1 m and 15 m on 29 May 2018, show higher subsurface values for both quantities compared to the surface values (<xref ref-type="fig" rid="f10">
<bold>Figures&#xa0;10A, C, E, G</bold>
</xref>). A subsurface Chl-a maximum accumulates on the southern flank of the eddy (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9I</bold>
</xref>). Sharp lateral gradients in both quantities also coincide with the sharp lateral density gradients (<xref ref-type="fig" rid="f9">
<bold>Figures&#xa0;9E, G, I, K</bold>
</xref>) and there are strong upward and downward vertical fluxes of Chl-a (ca. 60 mg m<sup>-2</sup> d<sup>-1</sup>) along the northern, southern and western flanks of the eddy (<xref ref-type="fig" rid="f11">
<bold>Figures&#xa0;11A, C</bold>
</xref>).</p>
<fig id="f10" position="float">
<label>Figure&#xa0;10</label>
<caption>
<p>
<bold>(A-D)</bold> K<sub>d</sub>(492) and <bold>(E-H)</bold> Chl-a at 1m and 15m on May 29 and July 26, 2018.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g010.tif"/>
</fig>
<fig id="f11" position="float">
<label>Figure&#xa0;11</label>
<caption>
<p>Vertical flux of Chl-a (w*Chl-a, mg m<sup>-2</sup> d<sup>-1</sup>) along 13&#xb0;58&#x2019; E (TE) and 55&#xb0;8&#x2019; N (TN) on May 29 and July 26, 2018 <bold>(A-D)</bold>. Density contours are plotted as grey lines; <bold>(E)</bold> simplified schematic of coupled surface-deep layer dynamics driven by cyclonic eddy at the peak of MHW-1.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g011.tif"/>
</fig>
</sec>
<sec id="s3_2_2">
<label>3.2.2</label>
<title>26<sup>th</sup> July 2018 (MHW-2)</title>
<p>The characteristics of Ro at 1 m and 15 m (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5B, D</bold>
</xref>, <xref ref-type="fig" rid="f7">
<bold>7B, D</bold>
</xref>) on 26<sup>th</sup> July 2018 are quite different from those seen in May. The values are much smaller (ca. +/- 0.2) and no clear eddy-like structures are evident in the domain. Weaker northerly winds prevail east of 13&#xb0;58&#x2019; E while southerly winds prevail west of 13&#xb0;58&#x2019; E (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5B</bold>
</xref>). Along the southern Swedish coast, east of 13&#xb0;30&#x2019; E, surface velocities are directed onshore (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8D</bold>
</xref>) and an overturning circulation is evident characterized by downwelling near the coast and upwelling offshore at 55&#xb0;8&#x2019; N (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8F</bold>
</xref>). Some doming of isopycnals and upwelling of cooler subsurface water (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9B</bold>
</xref>) can also be seen at 55&#xb0;12&#x2019; N. This coincides with a lateral gradient in the diffuse attenuation coefficient at 492nm, K<sub>d</sub>(492) (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9F</bold>
</xref>) and Chl-a concentrations (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9J</bold>
</xref>) at this location. Strong downward and upward vertical fluxes of Chl-a (ca. 50 mg m<sup>-2</sup> d<sup>-1</sup>) are also seen near the coast and offshore at 55&#xb0;12&#x2019; N, respectively (<xref ref-type="fig" rid="f11">
<bold>Figure&#xa0;11B</bold>
</xref>). Horizontal cross sections of Kd(492) and Chl-a concentration at 1 m and 15 m on 26 July 2018, also reveal higher subsurface values for both quantities compared to surface values (<xref ref-type="fig" rid="f10">
<bold>Figures 10B, D, F, H</bold>
</xref>). However, lateral gradients in temperature, Kd(492) and Chl-a values are absent along the 55&#xb0; 8' N transect (<xref ref-type="fig" rid="f11">
<bold>Figures 9D, H, L</bold>
</xref>) and vertical fluxes of Chl-a are weak (<xref ref-type="fig" rid="f11">
<bold>Figure 11D</bold>
</xref>).</p>
</sec>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Coherent structures found during the MHW events</title>
<p>During MHW-1, a long-living, eddy-like coherent structure according to the methodology described in Section 2.5 is detected below the surface eddy structure visible in the satellite and modelled Chl-a and K<sub>d</sub>(490/492) fields (shown in <xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4A</bold>
</xref>, <xref ref-type="fig" rid="f10">
<bold>10C, G</bold>
</xref>). The structure emerges on May 28, 2018 at 04:30am at 13&#xb0;55&#x2019; E and 55&#xb0;6&#x2019; N and dies out on June 02, 2018 at 01:30am at 13&#xb0;47&#x2019; E and 55&#xb0;0.6&#x2019; N, travelling about 13 km (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>). The structure extends during its life time from about 24 m to about 44 m depth and has a mean volume of 0.4 km<sup>3</sup> (an equivalent diameter of about 4 km if the structure is assumed to be cylindrical). During its lifetime, the tracers inside the structure show tendencies of upward motion (<xref ref-type="fig" rid="f12">
<bold>Figures&#xa0;12A, B</bold>
</xref>) which coincides with a doming of isopycnals above the structure (<xref ref-type="fig" rid="f9">
<bold>Figures&#xa0;9A, C</bold>
</xref>). Tracers seeded along a rectangular slice at 13&#xb0;58&#x2019; E and 55&#xb0;8&#x2019; N in May 29, 2018 at 11:30am show a strong semi-circular shaped south-eastward displacement after 24h of integration (<xref ref-type="fig" rid="f12">
<bold>Figures&#xa0;12A, B</bold>
</xref>) in the surface layer which correlates with an increase in Chl-a at the eddy boundaries (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>). During MHW-2, the tracer displacement starting from the same slice on July 26, 2018 at 11:30am shows a less dynamic displacement into different directions with a small trend towards the north in the surface layers (<xref ref-type="fig" rid="f12">
<bold>Figures&#xa0;12C, D</bold>
</xref>), consistent with the onshore displacement of surface waters and overturning circulation seen earlier in the vertical velocity fields (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8F</bold>
</xref>), and the vertical fluxes of Chl-a (<xref ref-type="fig" rid="f12">
<bold>Figure&#xa0;12B</bold>
</xref>).</p>
<fig id="f12" position="float">
<label>Figure&#xa0;12</label>
<caption>
<p>
<bold>(A)</bold> 24h-mean chlorophyll along a 13&#xb0;58&#x2019; E slice <bold>(B)</bold> 55&#xb0;8&#x2019; N slice) and the tracer position (black dots) after 24h of particle tracking starting on a grid in the black rectangle at May 29, 2018 11:30am. The gray structure corresponds to the position of the coherent structure at that time. <bold>(C)</bold> 24h-mean chlorophyll along a 13&#xb0;58&#x2019; E slice <bold>(D)</bold> 55&#xb0;8&#x2019; N slice) and the tracer position (black dots) after 24h of particle tracking starting on a grid in the black rectangle at July 26, 2018 11:30am.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g012.tif"/>
</fig>
<p>Comparing the spatial dynamics of the May coherent structure with the Chl-a content in the water column for May 28 to May 31 (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>), we found that the structure correlates with a low Chl-a patch during the whole period of time and weak vertical fluxes of Chl-a. The long lifetime of the deeper coherent structure in May 2018 appears to be linked to the persistence of the eddy in the top 30m and the dynamics of the deeper layers are more coupled to the surface layers during MHW-1. In July, a stronger thermocline de-couples the dynamics between the surface and deeper layers. Combined with onshore winds in July, the tendency is rather to support onshore displacement of surface waters toward the Swedish coast, and a coastal overturning circulation cell.</p>
<p>The high Chl-a patch that bends around the southwestern edge of the eddy in May (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10G</bold>
</xref>) can be interpreted as the impact of a &#x201c;sticking&#x201d; (unstable) manifold (<xref ref-type="bibr" rid="B51">Lehahn et&#xa0;al., 2007</xref>) on the distribution of particles in the flow. Stable and unstable manifolds act as organizing structures of the flow separating regions of different dynamical behaviour (<xref ref-type="bibr" rid="B65">Prants, 2013</xref>). Due to the attracting properties of the unstable manifold, nutrients and plankton are collected along it. This &#x201c;sticking&#x201d; manifold is visible as a dark blue singular line in <xref ref-type="fig" rid="f13">
<bold>Figure&#xa0;13A</bold>
</xref> which also correspond to the high K<sub>d</sub> and Chl-a patch seen in <xref ref-type="fig" rid="f10">
<bold>Figures&#xa0;10C, G</bold>
</xref>. The manifold extends down into the upper water column, as seen in <xref ref-type="fig" rid="f13">
<bold>Figure&#xa0;13B</bold>
</xref> and provides a three-dimensional transport barrier around the southern, southeastern and eastern edge of the eddy. This mechanism could also explain the surface Chl-a patterns seen in May in the OLCI data (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<fig id="f13" position="float">
<label>Figure&#xa0;13</label>
<caption>
<p>Sticking manifold calculated as the trajectory length of backward integrated trajectories of the velocity field (<xref ref-type="bibr" rid="B66">Prants et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B65">Prants, 2013</xref>; <xref ref-type="bibr" rid="B43">Jimenez Madrid and Mancho, 2009</xref>; <xref ref-type="bibr" rid="B60">Mendoza and Mancho, 2010</xref>). The values of the trajectory length are assigned to the starting points of the trajectories. All these values make up a three-dimensional map. <bold>(A)</bold> Map of trajectory length of for 24h backward integrated trajectories starting at May 29, 2018 11:30am. The sticking manifold is displayed as the dark blue singular line. <bold>(B)</bold> Three-dimensional view of the sticking manifold.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1323271-g013.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>    <p>We examined two different MHW events in 2018 which bookmark the evolution and decay of a phytoplankton bloom in the Arkona Sea. Our objective was to explore the relationship between co-occurring MHW and phytoplankton bloom events with the following specific questions in mind:</p>
<list list-type="simple">
<list-item>
<p>1. How, and under which circumstances, do MHWs contribute to the initiation of phytoplankton blooms?</p>
</list-item>
<list-item>
<p>2. Which dynamics play a role in sustaining the bloom?</p>
</list-item>
<list-item>
<p>3. How deep and for how long is the impact of the MHW felt?</p>
</list-item>
<list-item>
<p>4. What role, if any, do MHWs play in terminating a phytoplankton bloom?</p>
</list-item>
<list-item>
<p>5. Do phytoplankton blooms have an enhancing effect on MHWs by creating a positive feedback from water constituent&#x2013;induced surface heating?</p>
</list-item>
</list>
<p>Marine heatwaves are typically associated with shallower mixed layer depths (<xref ref-type="bibr" rid="B32">Hayashida et&#xa0;al., 2020</xref>). <xref ref-type="bibr" rid="B19">Cook et&#xa0;al. (2022)</xref> show that atmospheric pressure systems, wind speed and latent heat fluxes are important contributing factors to the generation and decline of MHWs. Surface flux-driven MHWs are shallower and occur predominantly in summer (<xref ref-type="bibr" rid="B21">Elzahaby et&#xa0;al., 2021</xref>). MHW-1 develops mid-May in parallel with a strengthening and shoaling of the seasonal thermocline (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6A</bold>
</xref>). The heat budget diagnostic confirms that MHW-1 is an atmospheric-driven event. Surface salinity also decreases during the onset of MHW-1 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6B</bold>
</xref>), consistent with <xref ref-type="bibr" rid="B20">Elzahaby and Schaeffer (2019)</xref> whereby increased anomalous solar radiation and decreased wind stress act on the latent heat flux to reduce evaporation. Given sufficient light and supply of nutrients, phytoplankton growth will occur. Increased phytoplankton biomass in the surface will increase the surface temperature due to absorption of light by phytoplankton in the surface layer. A thermal structure is established in the water column which will impact the growth, transport and fate of phytoplankton biomass. The availability of light below the productive layer will be strongly reduced. In the absence of other physical transport processes on timescales which are relevant for phytoplankton growth (<xref ref-type="bibr" rid="B47">Kuhn et&#xa0;al., 2019</xref>), nutrients will become depleted in the surface layer, the supply of nutrients from deeper waters will be inhibited by the stronger thermocline mid-summer, and phytoplankton growth will be expected to decrease.</p>
<p>A model study by <xref ref-type="bibr" rid="B32">Hayashida et&#xa0;al. (2020)</xref> shows that background nutrient conditions will determine the response of phytoplankton blooms co-occurring with MHWs. Generally, they find that in nutrient poor waters, blooms are weaker during MHWs, whereas in nutrient-rich waters, blooms are stronger during MHWs. However, transport dynamics play a critical role in the supply of nutrients to the surface mixed layer. In our scenario, wind driven surface circulation underpins the development of a series of cyclonic eddies in May and June over the Arkona Basin. The eddy we focus on late May has Ro O(1), an indication that relative vorticity is the same order of magnitude as the Coriolis force, that circulation has departed from geostrophy, and that ephemeral submesoscale dynamics may be at play. Indeed, large vertical velocities (ca. 35 m d<sup>-1</sup>) along the edges of the eddy are seen, giving rise to upwelling and downwelling along the boundaries of the eddy (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8E, G</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>). Vertical fluxes of Chl-a are between 40 and 60 mg m<sup>-2</sup> d<sup>-1</sup>, along the north- and south-western boundaries of the eddy (<xref ref-type="fig" rid="f11">
<bold>Figures&#xa0;11A, C</bold>
</xref>). Upwelling vertical fluxes may also transport nutrients to the surface contributing to sustained phytoplankton growth. This potential coupling of surface and subsurface dynamics is illustrated schematically in <xref ref-type="fig" rid="f11">
<bold>Figure&#xa0;11E</bold>
</xref>.</p>
<p>A number of mechanisms relate the role of submesocales and mesoscale eddies and their impact on the horizontal and vertical distribution of Chl-a. Submesoscale phytoplankton patchiness is often visible in satellite images (<xref ref-type="bibr" rid="B48">Lapeyre and Klein, 2006</xref>; <xref ref-type="bibr" rid="B52">Levy et&#xa0;al., 2012</xref>), whereas in the Baltic Sea intense submesoscale activity has been observed through its imprint on cyanobacteria blooms (<xref ref-type="bibr" rid="B59">McWilliams, 2016</xref>). Several dynamical mechanisms that occur in the submesoscale regime, e.g., frontal subduction and submesoscale restratification (<xref ref-type="bibr" rid="B18">Chrysagi et&#xa0;al., 2021</xref>), have been proposed to explain not only the surface but also the subsurface biogeochemistry signals (<xref ref-type="bibr" rid="B35">Hosegood et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B53">Levy et&#xa0;al., 2018</xref>). Eddy stirring which describes the direction of rotational flow of the eddy field, is known to be a source of phytoplankton patchiness (<xref ref-type="bibr" rid="B1">Abraham, 1998</xref>; <xref ref-type="bibr" rid="B57">Martin, 2003</xref>) and will determine the position of the chlorophyll anomaly and the direction of propagation of the eddy relative to the ambient chlorophyll field (<xref ref-type="bibr" rid="B58">McGillicuddy, 2016</xref>). A cyclonic eddy in the northern hemisphere will result in a positive anomaly in the southwest quadrant, a negative anomaly in the northeast quadrant and westward propagation. We see a Chl-a maximum accumulate on the southern flank of the eddy at about 15 m depth (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9I</bold>
</xref>).</p>
<p>According to <xref ref-type="bibr" rid="B56">Mahadevan (2016)</xref>, the vertical transport of nutrients into the euphotic zone can be achieved either by the vertical movement of the nutrient-rich isopycnal layers or by an advective flux. The former occurs mainly in the presence of internal waves or within mesoscale eddies but at different timescales. Eddies typically mix properties along isopycnals, and diapycnal mixing is considered to be weak. On the other hand, the advective flux of nutrients dominates in submesoscale features and tends to occur along the vertically tilted isopycnals. <xref ref-type="bibr" rid="B69">Rei&#xdf;mann et&#xa0;al., 2009</xref>, show that in the Baltic Sea, mesoscale eddies, known as Beddies may contribute to vertical mixing through different mechanisms. In particular, Beddies can contribute to the diapycnal mixing, inside the permanent halocline region, either through the vertical displacement of water and isopycnals, or through their decay. Another potential yet indirect mechanism through which Beddies may impact the halocline mixing, is through their interaction with internal waves, although this remains to be verified. Nevertheless, the exact dynamical mechanisms by which eddies might affect vertical mixing, are out of the scope of this study, since here we focus mainly on the co-occurence of MHWs and phytoplankton blooms.</p>
<p>We see indications of an unstable &#x201c;sticking&#x201d; manifold (<xref ref-type="fig" rid="f13">
<bold>Figure&#xa0;13A</bold>
</xref>), arising from geostrophic eddy stirring, acting as both a horizontal transport barrier for Chl-a and a facilitator of phytoplankton growth mediated by nutrient upwelling, as observed by <xref ref-type="bibr" rid="B51">Lehahn et&#xa0;al. (2007)</xref>. The manifold extends down into the upper water column (<xref ref-type="fig" rid="f13">
<bold>Figure&#xa0;13B</bold>
</xref>) and provides a three-dimensional transport barrier around the southern, southeastern and eastern edge of the eddy. The emergence of the long-living, eddy-like coherent structure below the surface eddy structure illustrates the role of the surface eddy in isolating a cold, dense, sub-surface, low Chl-a water body below the surface eddy pointing at a coupling of surface and subsurface dynamics (<xref ref-type="fig" rid="f12">
<bold>Figures&#xa0;12A, B</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure&#xa0;1</bold>
</xref>). Outcropping of isopycnals and upwelling of colder water along the southwestern edge of the eddy (as seen in <xref ref-type="fig" rid="f9">
<bold>Figures&#xa0;9A, C</bold>
</xref>) produce cold sea surface temperature anomalies which would tend to draw heat into the ocean from the atmosphere, further increasing stratification in those features relative to ambient waters. The complexity of eddy-induced transport mechanisms and biophysical interactions is reviewed in depth in <xref ref-type="bibr" rid="B58">McGillicuddy (2016)</xref> and references therein).</p>
<p>According to <xref ref-type="bibr" rid="B53">Levy et&#xa0;al. (2018)</xref>, the combination of the seasonal distribution of light and the vertical supply of nutrients through ephemeral fronts are essential ingredients for sustained phytoplankton growth. The timing of the onset of MHW-1 in May provides optimal light conditions which, given sufficient supply of nutrients, may contribute to the initiation of the phytoplankton bloom. The series of eddies which ensue, support both geostrophic eddy stirring and transient submesoscale dynamics along the edges of the eddies which may in turn provide both an upward and downward transport path for nutrient fluxes and carbon export (<xref ref-type="bibr" rid="B12">Callbeck et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B70">Ruiz et&#xa0;al., 2019</xref>).</p>
<p>It is important to note that the surface signature of both MHWs, does not represent the deeper signature and impact of the MHWs. <xref ref-type="bibr" rid="B71">Schaeffer and Roughan (2017)</xref> suggest that vertical mixing, in combination with downwelling favourable winds, can weaken stratification and enable MHWs to extend deeper than the surface mixed layer, thus homogenizing the water column. We see a deepening of the thermocline take place in MHW-1 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6A</bold>
</xref>) shortly after the peak in Chl-a concentration and associated water constituent-induced heating rate at the end of May (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<p>Water constituent-induced surface heating increases with the onset the bloom, and peaks in the middle of MHW-1 at 0.7 K d<sup>-1</sup> (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). This increase in water constituent-induced surface heating is clearly a response to increases in Chl-a concentration and is directly related to the absorption of light by phytoplankton (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6C</bold>
</xref>) and to a lesser extent, detrital and CDOM absorption (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6D, E</bold>
</xref>). MHW-1 contributes to the initiation of the bloom, which in turn contributes to an increase in water constituent-induced heating rates. A recent study by <xref ref-type="bibr" rid="B9">Cahill et&#xa0;al. (2023)</xref> found that, in 2018, water constituent-induced surface heating rates in the Western Baltic Sea could reach up to 0.4 to 0.8&#x2009;K&#x2009;d<sup>&#x2212;1</sup> during the period April to September. Moreover, this water constituent-induced surface warming resulted in a mean loss of heat (ca. 5&#x2009;W&#x2009;m<sup>&#x2212;2</sup>) from the sea to the atmosphere, primarily in the form of latent and sensible heat fluxes. Thus, the phytoplankton bloom may have an enhancing effect on the MHW as a consequence of its contribution to surface warming but air-sea energy flux exchange plays a role in regulating this.</p>
<p>MHW-2 develops late-July, and coincides with a re-stratification of the water column, a stronger, shallower thermocline and the end of the phytoplankton bloom. The Rossby number is small, indicating that there is no significant submesoscale activity. The end of the phytoplankton bloom is preceded by a period where the temperature anomaly extended the full depth of the water column. <xref ref-type="bibr" rid="B80">Zhan et&#xa0;al (2023)</xref> examine the roles of atmospheric forcing-driven and oceanic processes-driven MHWs in driving changes in Chl-a concentrations and phytoplankton biomass. They show that atmospheric forcing-driven MHWs (like MHW-1) tend to increase Chl-a concentrations and phytoplankton biomass, while oceanic processes-driven MHWs tend to decrease Chl-a concentrations and phytoplankton biomass. The heat budget diagnostic shows MHW-2 to be a mixed atmospheric-and horizontal advection-driven event. Indeed, advection of higher salinity, warmer bottom water precedes MHW-2 (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6A, B</bold>
</xref>). Changes in the seasonality of saltwater inflows from the North Sea to the Baltic Sea has been shown to cause exceptional warming trends in the Western Baltic Sea (<xref ref-type="bibr" rid="B2">Barghorn et&#xa0;al., 2023</xref>). The stronger thermocline in MHW-2 de-couples the dynamics between the surface and deeper layers and any supply of nutrients to the surface from the deeper waters is cut off. Only along the Swedish coast do we see evidence of some downwelling and upwelling (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8F</bold>
</xref>) and strong vertical fluxes of Chl-a (ca. 40 mg m<sup>-2</sup> d<sup>-1</sup>) (<xref ref-type="fig" rid="f11">
<bold>Figure&#xa0;11B</bold>
</xref>).</p>
<p>The persistence (up to 35 days between 15 m and 20 m) of the subsurface maximum temperature anomaly after the surface signature of MHW-2 has disappeared is remarkable. This underscores the importance of considering subsurface hydrography in order to fully understand the impact of MHWs on biological production (<xref ref-type="bibr" rid="B71">Schaeffer and Roughan, 2017</xref>). Moreover, the vertical extent of the subsurface temperature anomaly will play a role in the distribution of horizontal gradients in density and intensification of fronts, and thus determine transport pathways for nutrients and carbon.</p>
<p>In summary, we find that in the shallow Arkona Sea, the timing of atmospheric driven MHWs can contribute to the initiation of a phytoplankton bloom by providing the optimal conditions for phytoplankton growth. The vertical component of an eddy-like structure&#x2019;s vorticity balance determines the strength of coupling between the surface and subsurface dynamics, and seems to provide a vertical transport pathway for nutrients. These coupled dynamics might in turn contribute to sustaining phytoplankton growth. Depth-integrated phytoplankton will be restricted within a shallower mixed layer which will in turn increase surface heating. Thus, phytoplankton blooms may have an enhancing effect on MHWs as a consequence of their contribution to surface warming but air-sea energy flux exchange will play a role in regulating this exchange. The subsurface signature of MHWs is often stronger and persists for longer than its surface signature. Moreover, the subsurface signature of the MHW plays a critical role in de-coupling surface and subsurface dynamics and terminating a phytoplankton bloom.</p>
<p>Remotely sensed ocean colour data provides a window into the spatial complexity of optically active constituents in surface waters and how these constituents are transported by surface circulation. Knowledge of coincident surface winds, can provide some clues as to what may be occurring subsurface but used in tandem with 3D biogeochemical ocean models, it is possible to see how surface and subsurface dynamics are really connected and extend our knowledge on surface layer processes obtained from satellite data. In the last decade, ocean colour observations have been recognized as essential climate variables (ECVs) and become an integral part of ocean observing systems. While the volume of this data set increases, it remains underexploited in operational biogeochemical modelling and forecasting. Moreover, the availability of new satellites, automated measurement systems such as Bio-Argo floats, drones, and computational resources presents both a challenge and an opportunity to advance integrated observing systems which combine optical observations, including remotely sensed ocean colour, with biogeochemical models to monitor and predict the impact of extreme events on biogeochemical cycles and ecosystem functioning.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>Publicly available datasets were analyzed in this study. This data can be found here: NOAA OI SST V2 High Resolution Dataset: <uri xlink:href="https://psl.noaa.gov/data/gridded/data.noaa.oisst.v2.highres.html">https://psl.noaa.gov/data/gridded/data.noaa.oisst.v2.highres.html</uri> OLCI Level 3 300m Baltic Sea Ocean Colour Plankton, Transparency and Optics NRT daily observations, CMEMS: <uri xlink:href="https://doi.org/10.48670/moi-00294">https://doi.org/10.48670/moi-00294</uri>. The ROMS-BioOptic model code can be accessed at <uri xlink:href="https://www.myroms.org">https://www.myroms.org</uri>. The code of ROMSpath can be found at <uri xlink:href="https://github.com/imcslatte/ROMSPath/tree/V1.0.0">https://github.com/imcslatte/ROMSPath/tree/V1.0.0</uri>. For the calculation of the coherent sets the 3d version of the matlab scripts <uri xlink:href="https://github.com/gaioguy/FEMDL">https://github.com/gaioguy/FEMDL</uri> and <uri xlink:href="https://github.com/gfroyland/SEBA">https://github.com/gfroyland/SEBA</uri> were used. Data of the tracer trajectories and the coherent set tracks upon request by RVK.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>BC: Writing &#x2013; review &amp; editing, Writing &#x2013; original draft, Visualization, Validation, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. EC: Writing &#x2013; review &amp; editing, Visualization, Investigation, Conceptualization. RV-K: Writing &#x2013; review &amp; editing, Visualization, Investigation, Formal analysis, Conceptualization. UG: Writing &#x2013; review &amp; editing, Resources.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. BC was supported by funding from the German Research Foundation (Grant No. CA 1347/2-1, 2018 to 2021, Temporary Position for Principal Investigator). EC was supported by funding from the CONMAR (CONcepts for conventional MArine Munition Remediation in the German North and Baltic Sea) project. RV-K would like to thank DFG Project VO 2508/1-1 for funding. UG was supported by the Federal Ministry of Education and Research (BMBF, funding code 03F0860H), within the project: ECAS-Baltic project: Strategies of ecosystem-friendly coastal protection and ecosystem-supporting coastal adaptation for the German Baltic Sea Coast. The authors gratefully acknowledge the computing time granted by the Resource Allocation Board and provided on the supercomputer Lise and Emmy at NHR@ZIB and NHR@G&#xf6;ttingen as part of the NHR (North German Supercomputing Alliance) infrastructure. The calculations for this research were conducted with computing resources under the project ID bek00027.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>BC and RV-K would like to thank John Wilkin, Rutgers University for his advice on setting up the ROMSpath simulations.</p>
</ack>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The 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="s9" 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="s10" 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.1323271/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmars.2024.1323271/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.pdf" id="SM1" mimetype="application/pdf"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abraham</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>The generation of plankton patchiness by turbulent stirring</article-title>. <source>Nature</source> <volume>391</volume>, <fpage>577</fpage>&#x2013;<lpage>580</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1038/35361</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barghorn</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Meier</surname> <given-names>H. E. M.</given-names>
</name>
<name>
<surname>Radtke</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>). Changes in seasonality of saltwater inflows caused by exceptional warming trends in the Western Baltic Sea</article-title>. <source>Geophysical Res. Lett.</source> <volume>50</volume>, <elocation-id>e2023GL103853</elocation-id>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2023GL103853</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Belkin</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Rapid warming of large marine ecosystems</article-title>. <source>Prog. Oceanography</source> <volume>81</volume>, <fpage>207</fpage>&#x2013;<lpage>213</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/J.POCEAN.2009.04.011</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benthuysen</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Oliver</surname> <given-names>E. C. J.</given-names>
</name>
<name>
<surname>Feng</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Marshall</surname> <given-names>A. G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Extreme marine warming across tropical Australia during austral summer 2015&#x2013;2016</article-title>. <source>J. Geophys. Res. Oceans</source> <volume>123</volume>, <fpage>1301</fpage>&#x2013;<lpage>1326</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2017JC013326</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bissett</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Walsh</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Dieterle</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Carder</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>1999</year>a). <article-title>Carbon cycling in the upper waters of the Sargasso Sea: I. Numerical simulation of differential carbon and nitrogen fluxes</article-title>. <source>Deep Sea Res.</source> <volume>46</volume>, <fpage>205</fpage>&#x2013;<lpage>269</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0967-0637(98)00062-4</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bissett</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Walsh</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Dieterle</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Carder</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>1999</year>b). <article-title>Carbon cycling in the upper waters of the Sargasso Sea: II Numerical simulation of apparent and inherent optical properties</article-title>. <source>Deep Sea Res.</source> <volume>46</volume>, <fpage>271</fpage>&#x2013;<lpage>317</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0967-0637(98)00063-6</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bond</surname> <given-names>N. A.</given-names>
</name>
<name>
<surname>Cronin</surname> <given-names>M. F.</given-names>
</name>
<name>
<surname>Freeland</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Mantua</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Causes and impacts of the 2014 warm anomaly in the NE Pacific. Geophys</article-title>. <source>Res. Lett.</source> <volume>42</volume>, <fpage>414</fpage>&#x2013;<lpage>420</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2015GL063306</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bowen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Markham</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Sutton</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Shears</surname> <given-names>N. T.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Interannual variability of sea surface temperature in the southwest Pacific and the role of ocean dynamics</article-title>. <source>J. Climate</source>. <volume>30</volume>, <page-range>7481&#x2013;7492</page-range> doi:&#xa0;<pub-id pub-id-type="doi">10.1175/JCLI-D-16-0852.1</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cahill</surname> <given-names>B. E.</given-names>
</name>
<name>
<surname>Kowalczuk</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Kritten</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Gr&#xe4;we</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fischer</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Estimating the seasonal impact of optically significant water constituents on surface heating rates in the western Baltic Sea</article-title>. <source>Biogeosciences</source> <volume>20</volume>, <fpage>2743</fpage>&#x2013;<lpage>2768</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.5194/bg-20-2743-2023</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cahill</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Schofield</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Chant</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Hunter</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Glenn</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>). <article-title>Dynamics of turbid buoyant plumes and the feedbacks on near-shore biogeochemistry and physics</article-title>. <source>Geo. Res. Lett.</source> <volume>35</volume>, <fpage>1</fpage>&#x2013;<lpage>6</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.29/2008GL033595</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cahill</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fennel</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Vandemark</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Friedrichs</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Interannual and seasonal variabilities in air-sea CO2 fluxes along the U.S eastern continental shelf and their sensitivity to increasing air temperatures and variable winds</article-title>. <source>J. Geophys. Res.</source> <volume>121</volume>, <fpage>295</fpage>&#x2013;<lpage>311</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2015JG002939</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Callbeck</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Lavik</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Stramma</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Kuypers</surname> <given-names>M. M. M.</given-names>
</name>
<name>
<surname>Bristow</surname> <given-names>L. A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Enhanced nitrogen loss by eddy-induced vertical transport in the offshore Peruvian oxygen minimum zone</article-title>. <source>PLoSONE</source> <volume>12</volume>, <elocation-id>e0170059</elocation-id>doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0170059</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chelton</surname> <given-names>D. B.</given-names>
</name>
<name>
<surname>Gaube</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Schlax</surname> <given-names>M. G.</given-names>
</name>
<name>
<surname>Early</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Samelson</surname> <given-names>R. M.</given-names>
</name>
</person-group> (<year>2011</year>b). <article-title>The Influence of Nonlinear mesoscale eddies on near-surface oceanic chlorophyll</article-title>. <source>Science</source> <volume>334</volume>, <fpage>328</fpage>&#x2013;<lpage>332</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1208897</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chelton</surname> <given-names>D. B.</given-names>
</name>
<name>
<surname>Schlax</surname> <given-names>M. G.</given-names>
</name>
<name>
<surname>Samelson</surname> <given-names>R. M.</given-names>
</name>
</person-group> (<year>2011</year>a). <article-title>Global observations of nonlinear mesoscale eddies</article-title>. <source>Prog. Oceanography</source> <volume>91</volume>, <fpage>167</fpage>&#x2013;<lpage>216</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.pocean.2011.01.002</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chelton</surname> <given-names>D. B.</given-names>
</name>
<name>
<surname>Schlax</surname> <given-names>M. G.</given-names>
</name>
<name>
<surname>Samelson</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>de Szoeke</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Global observations of large oceanic eddies</article-title>. <source>Ge. Res. Lett.</source> <volume>34</volume>, <fpage>L15606</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2007GL030812</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Gawarkiewicz</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Kwon</surname> <given-names>Y. O.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W. G.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The role of atmospheric forcing versus ocean advection during the extreme warming of the Northeast U.S. continental shelf in 2012</article-title>. <source>J. Geophys. Res. Oceans</source> <volume>120</volume>, <fpage>4324</fpage>&#x2013;<lpage>4339</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2014JC010547</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Gawarkiewicz</surname> <given-names>G. G.</given-names>
</name>
<name>
<surname>Lentz</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Bane</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Diagnosing the warming of the Northeastern U.S. Coastal Ocean in 2012: a linkage between the atmospheric jet stream variability and ocean response</article-title>. <source>J. Geophys. Res. Oceans</source> <volume>119</volume>, <fpage>218</fpage>&#x2013;<lpage>227</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2013JC009393</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chrysagi</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Umlauf</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Holtermann</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Klingbeil</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Burchard</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>High-resolution simulations of submesoscale processes in the Baltic Sea: The role of storm events</article-title>. <source>J. Geophysical Research: Oceans</source>, <fpage>126</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2020JC016411</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cook</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>R. O.</given-names>
</name>
<name>
<surname>Roughan</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Cullen</surname> <given-names>N. J.</given-names>
</name>
<name>
<surname>Shears</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Bowen</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Marine heatwaves in shallow coastal ecosystems are coupled with the atmosphere: Insights from half a century of daily in <italic>situ</italic> temperature records</article-title>. <source>Front. Climate</source> <volume>4</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fclim.2022.1012022</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elzahaby</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Schaeffer</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Observational insight into the subsurface anomalies of marine heatwaves</article-title>. <source>Front. Mar. Sci.</source> <volume>6</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmars.2019.00745</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elzahaby</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Schaeffer</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Roughan</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Delaux</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Oceanic circulation drives the deepest and longest marine heatwaves in the East Australian Current System</article-title>. <source>Geophysical Res. Lett.</source> <volume>48</volume>, <elocation-id>e2021GL094785</elocation-id>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2021GL0947785</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eppley</surname> <given-names>R. W.</given-names>
</name>
</person-group> (<year>1972</year>). <article-title>Temperature and phytoplankton growth in the sea, Fish</article-title>. <source>Bull.</source> <volume>70</volume>, <fpage>1063</fpage>&#x2013;<lpage>1085</lpage>.</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fennel</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Seifert</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kayser</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Rossby radii and phase speeds in the Baltic Sea</article-title>. <source>Continental Shelf Res. 11</source> <volume>1</volume>, <fpage>23</fpage>&#x2013;<lpage>236</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0278-4343(91)90032-2</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fennel</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Quantifying biological carbon export for the northwest North Atlantic continental shelves</article-title>. <source>Geo. Res. Lett.</source> <volume>36</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2009GL039818</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fennel</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Levin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Moisan</surname> <given-names>J.</given-names>
</name>
<name>
<surname>O&#x2019;Reilly</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Haidvogel</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Nitrogen cycling in the Middle Atlantic Bight: Results from a three dimensional model and implications for the North Atlantic nitrogen budget</article-title>. <source>Glob. Biogeo. Cycles</source> <volume>20</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2005GB002456</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fennel</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Previdi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Najjar</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Denitrification effects on air-sea CO<sub>2</sub> flux in the coastal ocean: Simulations for the northwest North Atlantic</article-title>. <source>Geo. Res. Lett.</source> <volume>35</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2008GL036147</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Froyland</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Junge</surname> <given-names>O.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Robust FEM-based extraction of finite-time coherent sets using scattered, sparse, and incomplete trajectories</article-title>. <source>SIAM J. Appl. Dynamical Syst. 17</source> <volume>2</volume>, <fpage>1891</fpage>&#x2013;<lpage>1924</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1137/17M1129738</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Froyland</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Rock</surname> <given-names>C. P.</given-names>
</name>
<name>
<surname>Sakellariou</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Sparse eigenbasis approximation: Multiple feature extraction across spatiotemporal scales with application to coherent set identification</article-title>. <source>Commun. Nonlinear Sci. Numerical Simulation</source> <volume>77</volume>, <fpage>81</fpage>&#x2013;<lpage>107</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cnsns.2019.04.012</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gr&#xe4;we</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Holtermann</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Klingbeil</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Burchard</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>a). <article-title>Advantages of vertically adaptive coordinates in numerical models of stratified shelf seas</article-title>. <source>Ocean Model.</source> <volume>92</volume>, <fpage>56</fpage>&#x2013;<lpage>68</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ocemod.2015.05.008</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gr&#xe4;we</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Naumann</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Mohrholz</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Burchard</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>b). <article-title>Anatomizing one of the largest saltwater inflows into the Baltic Sea in December 2014</article-title>. <source>J. Geophysical Research: Oceans</source> <volume>120</volume>, <fpage>7676</fpage>&#x2013;<lpage>7697</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2015JC011269</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Haidvogel</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Arango</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Budgell</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Cornuelle</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Curchitser</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Lorenzo</surname> <given-names>E.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>). <article-title>Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System</article-title>. <source>J. Comp. Phys.</source> <volume>227</volume>, <fpage>3595 3624</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jcp.2007.06.016</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hayashida</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Matear</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Strutton</surname> <given-names>P. G.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Background nutrient concentration determines phytoplankton bloom response to marine heatwaves</article-title>. <source>Global Change Biol.</source> <volume>26</volume>, <fpage>4800</fpage>&#x2013;<lpage>4811</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1111/gcb.15255</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hersbach</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Bell</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Berrisford</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Hirahara</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Horanyi</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Munoz-Sabater</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>The ERA5 global reanalysis</article-title>. <source>Q. J. R. Meteorological Soc.</source> <volume>146</volume>, <fpage>1999</fpage>&#x2013;<lpage>2049</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/qj.3803</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hobday</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Oliver</surname> <given-names>E. C. J.</given-names>
</name>
<name>
<surname>Gupta</surname> <given-names>A. S.</given-names>
</name>
<name>
<surname>Benthuysen</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Burrows</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Donat</surname> <given-names>M. G.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Categorizing and naming marine heatwaves</article-title>. <source>Oceanography Special Issue Ocean Warming</source> <volume>31</volume>, <fpage>162</fpage>&#x2013;<lpage>173</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5670/oceanog</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hosegood</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Gregg</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Alford</surname> <given-names>M. H.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Wind-driven submesoscale subduction at the north Pacific subtropical front, J. Geophys</article-title>. <source>Res. Oceans</source> <volume>118</volume>, <fpage>5333</fpage>&#x2013;<lpage>5352</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/jgrc.20385</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>H&#xf8;yer</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Karagali</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Sea surface temperature climate data record for the North Sea and Baltic Sea</article-title>. <source>J. Climate</source> <volume>29</volume>, <fpage>2529</fpage>&#x2013;<lpage>2541</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1175/JCLI-D-15-0663.1</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>H&#xf8;yer</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>She</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Optimal interpolation of sea surface temperature for the North Sea and Baltic Sea, J. Mar. Sys., Vol 65, 1-4, pp H&#xf8;yer, J. L. and She, J., Optimal interpolation of sea surface temperature for the North Sea and Baltic Sea</article-title>. <source>J. Mar. Sys. Vol</source> <volume>65</volume>, <fpage>1</fpage>&#x2013;<lpage>4</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jmarsys.2005.03.008</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Banzon</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Freeman</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Graham</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Hankins</surname> <given-names>B.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Improvments of the daily optimum interpoaltion sea surface temperature (DOISST) version 2.1</article-title>. <source>J. Climate</source> <volume>348</volume>, <page-range>2923&#x2013;2939</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1175/JCLI-D-20-0166.1</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Humborg</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Geibel</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>X.</given-names>
</name>
<name>
<surname>McCrackin</surname> <given-names>M.</given-names>
</name>
<name>
<surname>M&#xf6;rth</surname> <given-names>C.-M.</given-names>
</name>
<name>
<surname>Stranne</surname> <given-names>C.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>High emissions of carbon dioxide and methane from the coastal Baltic Sea at the end of a summer heat wave</article-title>. <source>Front. Mar. Sci.</source> <volume>6</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmars.2019.00493</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hunter</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Fuchs</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Wilkin</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Gerbi</surname> <given-names>G. P.</given-names>
</name>
<name>
<surname>Chant</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Garwood</surname> <given-names>J. C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>ROMSPath v1.0: offline particle tracking for the Regional Ocean Modeling System (ROMS)</article-title>. <source>Geoscientific Model. Dev. 15</source> <volume>11</volume>, <fpage>4297</fpage>&#x2013;<lpage>4311</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.5194/gmd-15-4297-2022</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>IPCC</collab>
</person-group> (<year>2023</year>). &#x201c;<article-title>Summary for policymakers. In: climate change 2023: synthesis report. A report of the intergovernmental panel on climate change</article-title>,&#x201d; in <source>Contribution of working groups I, II and III to the sixth assessment report of the intergovernmental panel on climate change [CoreWriting team</source>, vol. <volume>36</volume> . Eds. <person-group person-group-type="editor">
<name>
<surname>Lee</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Romero</surname> <given-names>J.</given-names>
</name>
</person-group> (<publisher-name>Switzerland</publisher-name>, <publisher-loc>IPCC, Geneva</publisher-loc>).</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jaanus</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Andersson</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Olenina</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Toming</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kaljurand</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Changes in phytoplankton communities along a north-south gradient in the Baltic Sea between 1990 and 2008</article-title>. <source>Boreal Environ. Res.</source> <volume>16</volume>, <fpage>191</fpage>&#x2013;<lpage>208</lpage>.</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jim&#xe9;nez Madrid</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Mancho</surname> <given-names>A. M.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Distinguished trajectories in time dependent vector fields</article-title>. <source>Chaos 1;</source> <volume>19</volume>, <fpage>013111</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1063/1.3056050</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kahru</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Elmgren</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Sachuk</surname> <given-names>O.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Changing seasonality of the baltic sea</article-title>. <source>Biogeosciences.</source> <volume>113</volume>, <fpage>1009</fpage>&#x2013;<lpage>1018</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.5194/bg-13-1009-2016</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Karimova</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Gade</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Improved statistics of sub-mesoscale eddies in the Baltic Sea retrieved from SAR imagery</article-title>. <source>Int. J. Remote Sens.</source> <volume>37</volume>, <fpage>2394</fpage>&#x2013;<lpage>2414</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/01431161.2016.1145367</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kowalczuk</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Stedmon</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Markager</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Modelling absorption by CDOM in the Baltic Sea from season, salinity and chlorophyll</article-title>. <source>Mar. Chem.</source> <volume>101</volume>, <fpage>1</fpage>&#x2013;<lpage>11</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.marchem.2005.12.005</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kuhn</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Dutkiewicz</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Jahn</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Clayton</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Rynearson</surname> <given-names>T. S.</given-names>
</name>
<name>
<surname>Mazloff</surname> <given-names>M. R.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Temperoal and spatial scales of correlation in marine phytoplankton communities</article-title>. <source>J. Geophysical Research: Oceans</source> <volume>124</volume>, <fpage>9417</fpage>&#x2013;<lpage>9438</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2019JC015331</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lapeyre</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Klein</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Impact of the small-scale elongated filaments on the ocean vertical pump</article-title>. <source>J. Mar. Res.</source> <volume>64</volume>, <fpage>835</fpage>&#x2013;<lpage>851</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1357/002224006779698369</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lass</surname> <given-names>H. U.</given-names>
</name>
<name>
<surname>Mohrholz</surname> <given-names>V.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>On the dynamics and mixing of inflowing saltwater in the Arkona Sea</article-title>. <source>J. Geophys. Res.</source> <volume>108</volume>, <fpage>3042</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2002JC001465</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Du</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Arnone</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>A model for the diffuse attenuation coefficient of downwelling irradiance</article-title>. <source>J. Geophysical Research Oceans</source> <volume>110</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2004JC002275</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lehahn</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>d&#x2019;Ovidio</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Levy</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Heifetz</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Stirring of the northeast Atlantic spring bloom: A Lagrangian analysis based on multisatellite data</article-title>. <source>J. Geophys. Res.</source> <volume>112</volume>, <fpage>C08005</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2006JC003927</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Levy</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ferrari</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Franks</surname> <given-names>P. J. S.</given-names>
</name>
<name>
<surname>Martin</surname> <given-names>A. P.</given-names>
</name>
<name>
<surname>Rivi&#xe8;re</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Bringing physics to life at the submesoscale</article-title>. <source>Geophys. Res. Lett.</source> <volume>39</volume>, <fpage>L14602</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2012GL052756</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Levy</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Franks</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The role of submesoscale currents in structuring marine ecosystems</article-title>. <source>Nat. Commununications</source> <volume>9</volume>, <fpage>4758</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-018-07059-3</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Loeb</surname> <given-names>N. G.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>G. C.</given-names>
</name>
<name>
<surname>Thorsen</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Lyman</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Rose</surname> <given-names>F. G.</given-names>
</name>
<name>
<surname>Kato</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Satellite and ocean data reveal marked increase in Earth&#x2019;s heating rate</article-title>. <source>Geophysical Res. Lett.</source> <volume>48</volume>, <elocation-id>e2021GL093047</elocation-id>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2021GL093047</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Lorenz</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Marine heatwaves and the phytoplankton spring bloom: a comparative survey of marine heatwaves, spring bloom events and their relationship in the North Sea and Baltic Sea</source> (<publisher-name>University of Potsdam</publisher-name>), <fpage>86</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.5281/zenodo.10194802</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mahadevan</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The impact of submesoscale physics on primary productivity of plankton</article-title>. <source>Annu. Rev. Mar. Sci.</source> <volume>8</volume>, <fpage>161</fpage>&#x2013;<lpage>184</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev-marine-010814-015912</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Martin</surname> <given-names>A. P.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Phytoplankton patchiness: the role of lateral stirring and mixing</article-title>. <source>Prog. Oceanogr.</source> <volume>57</volume>, <fpage>125</fpage>&#x2013;<lpage>174</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0079-6611(03)00085-5</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McGillicuddy</surname> <given-names>D. J.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Mechanisms of physical-biological-biogeochemical interaction at the oceanic mesoscale</article-title>. <source>Annu. Rev. Mar. Sci.</source> <volume>8</volume>, <fpage>125</fpage>&#x2013;<lpage>159</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev-marine-010814-015606</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McWilliams</surname> <given-names>J. C.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Submesoscale currents in the ocean</article-title>. <source>Proc. R. Soc. A: Mathematical Phys. Eng. Sci.</source> <volume>472</volume>, <fpage>20160117</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1098/rspa.2016.0117</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mendoza</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Mancho</surname> <given-names>A. M.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Hidden geometry of ocean flows</article-title>. <source>Phys. Rev. Lett.</source> <volume>105</volume>, <elocation-id>38501</elocation-id>doi:&#xa0;<pub-id pub-id-type="doi">10.1103/PhysRevLett.105.038501</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Morel</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Optical modelling of the upper ocean in relation to its biogenous matter content (Case I waters)</article-title>. <source>J. Geophys. Res.</source> <volume>93</volume>, <fpage>749</fpage>&#x2013;<lpage>768</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1029/JC093iC09p10749</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Naumann</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Gr&#xe4;we</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Mohrholz</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Kuss</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Waniek</surname> <given-names>J. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Hydrographic-hydrochemical assessment of the baltic sea 2018</article-title>. <source>Meereswiss. Ber. Warnem&#xfc;nde</source> <volume>110</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.12754/msr-2019-0110</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Neumann</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Radtke</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Cahill</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Schmidt and G. Rehder</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Non-Redfieldian carbon model for the Baltic Sea (ERGOM version 1.2) &#x2013; implementation and budget estimates</article-title>. <source>Geoscientific Model. Dev.</source> <volume>15</volume>, <fpage>8473</fpage>&#x2013;<lpage>8540</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.5194/gmd-15-8473-2022</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oliver</surname> <given-names>E. C. J.</given-names>
</name>
<name>
<surname>Donat</surname> <given-names>M. G.</given-names>
</name>
<name>
<surname>Burrows</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Moore</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Smale</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Alexander</surname> <given-names>L. V.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Longer and more frequent marine heatwaves over the past century</article-title>. <source>Nat. Commun.</source> <volume>9</volume>, <fpage>1324</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-018-03732-9</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prants</surname> <given-names>S. V.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Dynamical systems theory methods to study mixing and transport in the ocean</article-title>. <source>Phys. Scr.</source> <volume>87</volume>, <elocation-id>38115</elocation-id>doi:&#xa0;<pub-id pub-id-type="doi">10.1088/0031-8949/87/03/038115</pub-id>
</citation>
</ref>
<ref id="B66">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prants</surname> <given-names>S. V.</given-names>
</name>
<name>
<surname>Uleysky</surname> <given-names>M. Y.</given-names>
</name>
<name>
<surname>Budyansky</surname> <given-names>M. V.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Numerical simulation of propagation of radioactive pollution in the ocean from the Fukushima Dai-ichi nuclear power plant</article-title>. <source>Dokl. Earth Sc.</source> <volume>439</volume>, <fpage>1179</fpage>&#x2013;<lpage>1182</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1134/S1028334X11080277</pub-id>
</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Racault</surname> <given-names>M.-F.</given-names>
</name>
<name>
<surname>Raitsos</surname> <given-names>D. E.</given-names>
</name>
<name>
<surname>Berumen</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Brewin</surname> <given-names>R. J. W.</given-names>
</name>
<name>
<surname>Platt</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Sathyendranath</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Phytoplankton phenology indices in coral reef ecosystems: Application to ocean-color observations in the Red Sea</article-title>. <source>Remote Sens. Environ.</source> <volume>160</volume>, <fpage>222</fpage>&#x2013;<lpage>234</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rse.2015.01.019</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rei&#xdf;mann</surname> <given-names>J. H.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>An algorithm to detect isolated anomalies in three-dimensional stratified data fields with an application to density fields from four deep basins of the Baltic Sea</article-title>. <source>J. Geophysical Research: Oceans</source> <volume>110</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2005JC002885</pub-id>
</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rei&#xdf;mann</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Burchard</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Feistel</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Hagen</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Lass</surname> <given-names>H. U.</given-names>
</name>
<name>
<surname>Mohrholz</surname> <given-names>V.</given-names>
</name>
<etal/>
</person-group>. (<year>2009</year>). <article-title>Vertical mixing in the Baltic Sea and consequences for eutrophication &#x2013; A review</article-title>. <source>Prog. Oceanography</source> <volume>821</volume>, <fpage>47</fpage>&#x2013;<lpage>80</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.pocean.2007.10.004</pub-id>
</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ruiz</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Claret</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Pascual</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Olita</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Troupin</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Capet</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Effects of oceanic mesoscale and submesoscale frontal processes on the vertical transport of phytoplankton</article-title>. <source>J. Geophysical Research Oceans</source> <volume>124</volume>, <fpage>5999</fpage>&#x2013;<lpage>6014</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2019JC015034</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schaeffer</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Roughan</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Subsurface intensification of marine heatwaves off southeastern Australia: the tole of stratification and local winds</article-title>. <source>Geophysical Res. Lett.</source> <volume>44</volume>, <fpage>5025</fpage>&#x2013;<lpage>5033</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2017GL073714</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Smale</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Wernberg</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Oliver</surname> <given-names>E. C. J.</given-names>
</name>
<name>
<surname>Thomsen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Harvey</surname> <given-names>B. P.</given-names>
</name>
<name>
<surname>Straub</surname> <given-names>S. C.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Marine heatwaves threaten global biodiversity and the provision of ecosystem services</article-title>. <source>Nat. Climate Change.</source> <volume>9</volume>, <fpage>206</fpage>&#x2013;<lpage>312</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41558-019-0412-1</pub-id>
</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Steele</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Henderson</surname> <given-names>E. W.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>The role of predation in plankton models</article-title>. <source>J. Plankton Res.</source> <volume>14</volume>, <fpage>157</fpage>&#x2013;<lpage>172</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1093/plankt/14.1.157</pub-id>
</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thomalla</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Fauchereau</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Swart</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Monteiro</surname> <given-names>P. M. S.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Regional scale characteristics of the seasonal cycle of chlorophyll in the Southern Ocean</article-title>. <source>Biogeosciences</source> <volume>8</volume>, <fpage>2849</fpage>&#x2013;<lpage>2866</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.5194/bg-8-2849-2011</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vortmeyer-Kley</surname> <given-names>R.</given-names>
</name>
<name>
<surname>L&#xfc;nsmann</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Berthold</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Gr&#xe4;we</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Feudel</surname> <given-names>U.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Eddies: fluid dynamical niches or transporters?&#x2013;A case study in the western baltic sea</article-title>. <source>Front. Mar. Sci.</source> <volume>6</volume>doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmars.2019.00118</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wasmund</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Dutz</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Kremp</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Zettler</surname> <given-names>M. L.</given-names>
</name>
</person-group> (<year>2019</year>b). <source>Biological assessment of the baltic sea 2018</source> (<publisher-loc>Warnemuende</publisher-loc>: <publisher-name>Meereswiss. Ber.</publisher-name>), <fpage>112</fpage>doi:&#xa0;<pub-id pub-id-type="doi">10.12754/msr-2019-0112</pub-id>
</citation>
</ref>
<ref id="B77">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wasmund</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Nausch</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Gerth</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Busch</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Burmeister</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Hansen</surname> <given-names>R.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>a). <article-title>Extension of the growing season of phytoplankton in the western Baltic Sea in response to climate change</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>622</volume>, <fpage>1</fpage>&#x2013;<lpage>16</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps12994</pub-id>
</citation>
</ref>
<ref id="B78">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wilkin</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W. G.</given-names>
</name>
<name>
<surname>Cahill</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Chant</surname> <given-names>R. C.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>Integrating coastal models and observations for studies of ocean dynamics, observing systems and forecasting</article-title>,&#x201d; in <source>Operational oceanography in the 21st century</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Schiller</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Brassington</surname>
</name>
</person-group> (<publisher-name>Springer</publisher-name>)doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-94-007-0332-2_19</pub-id>
</citation>
</ref>
<ref id="B79">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Z&#xe4;ngl</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Reinert</surname> <given-names>D.</given-names>
</name>
<name>
<surname>R&#xed;podas</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Baldauf</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The ICON (ICOsahedral Non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core</article-title>. <source>Quar. J. R. Met. Soc.</source> <volume>141</volume>, <fpage>563</fpage>&#x2013;<lpage>579</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1002/qj.2378</pub-id>
</citation>
</ref>
<ref id="B80">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhan</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>He</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Zhan</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Shifting responses of phytoplankton to atmospheric and oceanic forcing in a prolonged marine heatwave</article-title>. <source>Limnology Oceanography</source>. <volume>68</volume>, <page-range>1821&#x2013;1834</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/lno.12388</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>