<?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.2023.1121833</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>Laboratory study on movement characteristics of a river plume using particle image velocimetry</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Xu</surname>
<given-names>Zhenshan</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/2136263"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xu</surname>
<given-names>Chiyuan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Gang</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2136336"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Jiabo</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Yongping</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1417205"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Harbor, Coastal and Offshore Engineering, Hohai University</institution>, <addr-line>Nanjing, Jiangsu</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Marine Ecological Restoration and Smart Ocean Engineering Research Center of Hebei Province, Marine Ceological Resources Survey Center of Hebei Province</institution>, <addr-line>Qinhuangdao, Hebei</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Sha Lou, Tongji University, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Guoxiang Wu, Ocean University of China, China; Fan Xu, East China Normal University, China</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Zhenshan Xu, <email xlink:href="mailto:zsxu2006@hhu.edu.cn">zsxu2006@hhu.edu.cn</email>
</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Coastal Ocean Processes, a section of the journal Frontiers in Marine Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>17</day>
<month>02</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1121833</elocation-id>
<history>
<date date-type="received">
<day>12</day>
<month>12</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>06</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Xu, Xu, Wang, Zhang and Chen</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Xu, Xu, Wang, Zhang and Chen</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>Rivers are important passages for land-based materials transported to the sea, such as fresh water, sediment, pollutants, and nutrients. River discharge with land-based materials will form a river plume with a buoyancy flow. The river plume is subject to the interaction between nearshore runoff with seawater. In this study, the movement characteristics of the river plume are investigated based on particle image velocimetry (PIV) and the dye tracing method. It is found that the change of flow rate and environmental water density will shape the river plume pattern in both the plan view and the side view. The combined effect of flow rate and environmental water density could be described by the Froude number. The river plume has a free surface in the x-z plane and the horizontal velocity of the plume can be fitted with a 1/2 Gaussian distribution curve. The change of flow rate has little effect on the type of plume thickness curve, while the increase of environmental water density will change the thickness curve type of the river plume. The stable thickness of the river plume increases with the increase of the Froude number. The maximum value of turbulent kinetic energy is located in the middle layer, and the increase of the flow rate or the density difference leads to the increase in the turbulent kinetic energy.</p>
</abstract>
<kwd-group>
<kwd>river plume</kwd>
<kwd>dye tracing</kwd>
<kwd>particle image velocimetry</kwd>
<kwd>flow field</kwd>
<kwd>turbulent kinetic energy</kwd>
</kwd-group>
<counts>
<fig-count count="9"/>
<table-count count="1"/>
<equation-count count="7"/>
<ref-count count="17"/>
<page-count count="11"/>
<word-count count="4614"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Rivers are important passages for land-based materials transported to the sea. Lots of fresh water, sediment, pollutants, and nutrients are discharged into the adjacent sea area through estuaries. Land-based materials will directly affect the near-shore water environment and the aquatic ecological conditions. The impact area and the impact extent of land-based materials carried by rivers are primarily affected by the hydrodynamics of the river discharged into the sea. From the perspective of coastal environmental hydraulics, the discharge of rivers with land-based materials will form a river plume with a buoyancy flow. The river plume is the interaction between the nearshore runoff with the seawater (<xref ref-type="bibr" rid="B1">Cameron and Pritchard, 1963</xref>).</p>
<p>River plume has several characteristics that distinguishes it from pipe flow. The spatial scale of a river plume is much larger than that of a pipe flow. The river plume has a free surface and will move above the ambient water due to buoyancy. Similar to the jet or plume flow formed <italic>via</italic> a pipe, the movement and transport of river plume is driven by the momentum of the river and the dynamic conditions of the environmental water. Therefore, there are two main types of factors affecting the mixing and transport of river plume. One is the parameters related to the river plume itself, such as the flow rate, the estuary width, and the angle between river plume and coastline. The other is the parameters related to the environmental water, such as the water depth, the tidal current, and the wave. Due to the huge kinetic energy contained in the runoff, an interface of the upper freshwater layer and the lower salt water layer is formed. The interface shear is quite strong, resulting in the lower low-momentum and high-salinity water being swept and mixed by the upper flow (<xref ref-type="bibr" rid="B3">Garvine, 1975</xref>). At the same time, the various substances carried by the runoff, especially the pollutants and nutrients suspended in it, are also mixed and diluted by the mixing of freshwater with salt water.</p>
<p>Over the past decades, many researchers have conducted various studies on the river plume and analyzed the salinity change, the front structure, and dynamic process affected by it through field measurement data and the numerical simulation method. For example, <xref ref-type="bibr" rid="B7">Lentz (2004)</xref> developed a scale analysis theory to distinguish the features related to the river plume such as flow type, structure, and propagation velocity. <xref ref-type="bibr" rid="B4">Horner-Devine et&#xa0;al., (2009)</xref> provides a conceptual summary of the interaction between recirculated buoyancy flow and tides by cruise data. The tide modifies the structure of the plume in the region near the river mouth. The tidal plume flows over the top of the re-circulating plume and is typically bounded by strong fronts that may penetrate well below the re-circulating plume water and eventually spawn internal waves that mix the re-circulating plume further. In addition to the tide, the dynamic action of environmental fluid also includes wind and wave. The influence mechanism of wind on the mixing of estuarine plumes is different from that of tidal current. The wind increases the turbulence intensity on the surface of the flow, thereby enhancing the mixing effect of river plumes and ambient fluid. The wind also generates wind-driven flow on the surface of the water, affecting the movement track of river plumes (<xref ref-type="bibr" rid="B10">Nezlin and DiGiacomo, 2005</xref>; <xref ref-type="bibr" rid="B2">Feddersen et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B13">Rijnsburger et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B12">Qu and Hetland 2019</xref>). The impact of waves on river plumes is mainly that the frontal incoming waves increase the broadening of river plumes (<xref ref-type="bibr" rid="B9">Nardin et&#xa0;al., 2013</xref>).</p>
<p>However, the field observation of a river plume has certain limitations, which cannot accurately realize the motion observation of all water masses within a certain range. Therefore, laboratory experiments are required as a supplement to conduct more targeted research on the structure and movement characteristics of a river plume. <xref ref-type="bibr" rid="B14">Yuan et&#xa0;al. (2011</xref>, <xref ref-type="bibr" rid="B15">Yuan and Horner-Devine 2013</xref>) compared the river plume in the laboratory with and without lateral limits and found that lateral spreading dramatically modifies the plume structure. The spreading plume layer consists of approximately linear density and velocity profiles that extend to the surface, whereas the channelized plume profiles are uniform near the surface. <xref ref-type="bibr" rid="B16">Yuan et&#xa0;al. (2018)</xref> studied the influence of periodically changing flow on coastal buoyancy flow through the laboratory rotating platform. It was found that the bulge geometry oscillates between a circular plume structure that extends mainly in the offshore direction, and a compressed plume structure that extends mainly in the alongshore direction. The oscillations result in periodic variations in the width and depth of the bulge and the incidence angle formed where the bulge flow re-attaches with the coastal wall.</p>
<p>Despite meaningful results from previous research, laboratory investigations of river plume are still insufficient. The effect of the buoyancy and the flow rate on the movement characteristics of the river plume has been less studied. In fact, the difference in density between the upper and lower fluid layers has a great influence on the change of fluid turbulence intensity and velocity. The mixing process and movement characteristics of a river plume under the effect of different buoyancy and the plume flow rate are comprehensively studied using particle image velocimetry (PIV) technology. The main contents of this paper are organized as follows. In Section 2, the experimental setup and PIV measuring system are briefly described. The results, including the river plume qualitative description, the river plume quantitative flow field, Gaussian distribution of horizontal velocity, the river plume thickness, and the turbulent kinetic energy of the river plume, are presented in Section 3. Finally, conclusions are summarized in Section 4.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methodology</title>
<sec id="s2_1">
<label>2.1</label>
<title>Experimental setup</title>
<p>As shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>, the experiment was conducted in a rectangular tank (1.5m long, 0.6m wide, and 0.5m deep). The rectangular tank was made with transparent acrylic plates with a good light transmission. The river plume was discharged from the right side of the rectangular tank, and the left side was designed with extra length as a buffer to reduce the backflow. The source of the river plume is a constant water tank located 3.0m above the rectangular tank, which is connected to a rectangular channel through a pipe. The center of the river mouth water surface was defined as the origin of the Cartesian coordinate. x positive direction is right, representing the onshore coordinates; z positive direction is upward, representing vertical coordinates; y is alongshore coordinate. The mouth of the rectangular channel is submerged in the environmental water. The width of the rectangular channel is 5&#xa0;cm and the slope is 0.16.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Diagram of the experimental setup: qualitative plan-view experiment <bold>(A)</bold> and qualitative and quantitative side-view experiment <bold>(B)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g001.tif"/>
</fig>
<p>Two stages of experiment were designed in this study: the qualitative and the quantitative measurement. In the qualitative experiment, the plume was stained with dye, and two cameras were set above and in front of the rectangular tank, similar to the work of <xref ref-type="bibr" rid="B11">Pan et&#xa0;al. (2022)</xref>. As a result, the planar view of the plume movement at the x-y plane and the side view at the x-z plane were captured. For quantitative experiments, the velocity of the plume movement at the x-z plane was recorded by the PIV system. The size of the field of view (FOV) was 35&#xa0;cm &#xd7; 12&#xa0;cm.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>PIV system and post-processing</title>
<p>The PIV system consists of a high-speed camera with 2-Gigabyte-memory storage and a continuous laser with a wavelength of 532nm. The resolution of the high-speed camera is 2048 &#xd7; 2048 pixels. The camera framing rate was set to 200fps. The output power of the MGL-W-532 solid state laser was set to 8W. In the experiment, the incident light of the laser was parallel to the rectangular tank bottom and then was reflected by a plane mirror which had an angle of 45&#xb0; with the bottom. As a result, the laser light could pass vertically upwards through the central axis of the river plume. The thickness of the laser light was about 2mm. Hence, the river plume was illuminated by the laser light. At the same time, tracer particles with a diameter of 50&#x3bc;m were added into the river plume.</p>
<p>A multi-grid interrogation method, which can increase the spatial and temporal resolution of flow field (<xref ref-type="bibr" rid="B5">Hsieh, 2008</xref>; <xref ref-type="bibr" rid="B6">Hsieh et&#xa0;al., 2016</xref>), was used in the post-processing stage. The multi-grid interrogation method consisted of three different interrogation window sizes: 32&#xd7;32 pixels, 16&#xd7;16 pixels, and 8&#xd7;8 pixels. All passes adopted a 50% overlap between adjacent sub-windows. In traditional PIV image processing, a fixed time interval is usually used to calculate flow field. But if the velocity gradient is large (i.e., jet), adopting a fixed-time interval will inevitably result in bias error when calculating relatively small velocities. By adopting the multi-time interval method, the aforementioned bias error can be reduced. Four different time intervals (1, 3, 9, and 21) were used in this method; the most suitable time interval would be selected to reduce errors in the post-processing stage. Details of this multi-time interval algorithm can be found in <xref ref-type="bibr" rid="B6">Hsieh et&#xa0;al. (2016)</xref>.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Experimental cases</title>
<p>In this study, nine cases of river plume were carried out for both the qualitative dyeing experiments and the velocity measurements. The specific parameters of the experimental cases are shown in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>. The effects of the buoyancy change and the flow rate on the movement characteristics of the river plume were examined.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>The parameters of the experiments.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">expt</th>
<th valign="middle" align="center">Q<sub>0</sub>(Lpm)</th>
<th valign="middle" align="center">
<italic>&#x3c1;<sub>0</sub> (g/cm<sup>3</sup>)</italic>
</th>
<th valign="middle" align="center">g<sub>0</sub>&#x2019;</th>
<th valign="middle" align="center">U<sub>0</sub>(cm/s)</th>
<th valign="middle" align="center">
<italic>Fr<sub>i</sub>
</italic>
</th>
<th valign="middle" align="center">Re</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">ND-1</td>
<td valign="middle" align="center">1.50</td>
<td valign="middle" align="center">0.9997</td>
<td valign="middle" align="center">0.49</td>
<td valign="middle" align="center">3.00</td>
<td valign="bottom" align="center">1.92</td>
<td valign="bottom" align="center">7634</td>
</tr>
<tr>
<td valign="middle" align="left">ND-2</td>
<td valign="middle" align="center">2.00</td>
<td valign="middle" align="center">0.9997</td>
<td valign="middle" align="center">0.49</td>
<td valign="middle" align="center">3.80</td>
<td valign="bottom" align="center">2.43</td>
<td valign="bottom" align="center">9670</td>
</tr>
<tr>
<td valign="middle" align="left">ND-3</td>
<td valign="middle" align="center">2.50</td>
<td valign="middle" align="center">0.9997</td>
<td valign="middle" align="center">0.49</td>
<td valign="middle" align="center">4.90</td>
<td valign="bottom" align="center">3.13</td>
<td valign="bottom" align="center">12469</td>
</tr>
<tr>
<td valign="middle" align="left">HD-1</td>
<td valign="middle" align="center">1.50</td>
<td valign="middle" align="center">1.0030</td>
<td valign="middle" align="center">3.72</td>
<td valign="middle" align="center">2.90</td>
<td valign="bottom" align="center">0.67</td>
<td valign="bottom" align="center">7404</td>
</tr>
<tr>
<td valign="middle" align="left">HD-2</td>
<td valign="middle" align="center">2.00</td>
<td valign="middle" align="center">1.0030</td>
<td valign="middle" align="center">3.72</td>
<td valign="middle" align="center">3.80</td>
<td valign="bottom" align="center">0.88</td>
<td valign="bottom" align="center">9702</td>
</tr>
<tr>
<td valign="middle" align="left">HD-3</td>
<td valign="middle" align="center">2.50</td>
<td valign="middle" align="center">1.0030</td>
<td valign="middle" align="center">3.72</td>
<td valign="middle" align="center">5.10</td>
<td valign="bottom" align="center">1.18</td>
<td valign="bottom" align="center">13021</td>
</tr>
<tr>
<td valign="middle" align="left">ED-1</td>
<td valign="middle" align="center">1.50</td>
<td valign="middle" align="center">1.0060</td>
<td valign="middle" align="center">6.63</td>
<td valign="middle" align="center">3.10</td>
<td valign="bottom" align="center">0.54</td>
<td valign="bottom" align="center">7938</td>
</tr>
<tr>
<td valign="middle" align="left">ED-2</td>
<td valign="middle" align="center">2.00</td>
<td valign="middle" align="center">1.0060</td>
<td valign="middle" align="center">6.63</td>
<td valign="middle" align="center">4.00</td>
<td valign="bottom" align="center">0.69</td>
<td valign="bottom" align="center">10243</td>
</tr>
<tr>
<td valign="middle" align="left">ED-3</td>
<td valign="middle" align="center">2.50</td>
<td valign="middle" align="center">1.0060</td>
<td valign="middle" align="center">6.63</td>
<td valign="middle" align="center">5.10</td>
<td valign="bottom" align="center">0.89</td>
<td valign="bottom" align="center">12548</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The initial velocity of the plume could be calculated by using the initial flow rate. The reduced gravitational acceleration was used to describe the effect of different buoyancies. Hence, the Froude number <italic>Fr<sub>i</sub>
</italic>, which is the ratio of inertial force to the buoyancy, could be described as</p>
<disp-formula>
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:msubsup>
<mml:mi>g</mml:mi>
<mml:mn>0</mml:mn>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <italic>&#x3c1;</italic>
<sub>0</sub> is the density of ambient fluid, <italic>&#x3c1;</italic> is the density of inflow, and <italic>g</italic> is gravitational acceleration,</p>
<disp-formula>
<label>(2)</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:msubsup>
<mml:mi>g</mml:mi>
<mml:mn>0</mml:mn>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <italic>U</italic>
<sub>0,</sub> <inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:msubsup>
<mml:mi>g</mml:mi>
<mml:mn>0</mml:mn>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, and <italic>H</italic>
<sub>0</sub> represent the velocity, reduced gravity, and the maximum depth of river plume inflow. The Reynolds number of each case is also listed in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>.</p>
<disp-formula>
<label>(3)</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:mtext>Re</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mtext>&#x3bc;</mml:mtext>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where &#x3bc; represents the dynamic viscosity of the water.</p>
<p>In these cases, <italic>Fr<sub>i</sub>
</italic> ranges from 0.54 to 3.13. If <italic>Fr<sub>i</sub>
</italic> &gt; 1, the flow is supercritical and if <italic>Fr<sub>i</sub>
</italic> &lt; 1 the flow is subcritical. Re ranges from 7404 to 13021, indicating that the flow of each case is turbulent.</p>
</sec>
</sec>
<sec id="s3" sec-type="result|discussion">
<label>3</label>
<title>Result and discussion</title>
<sec id="s3_1">
<label>3.1</label>
<title>Qualitative description of river plume</title>
<p>River plume generally enters the ocean from a channel with a certain gradient. In this study, the slope is set to 0.16. The dynamic processes of a river plume movement could be divided into the near-field process and the far-field process. In the near-field area, the initial velocity of the river plume causes the velocity shear and turbulence mixing with the sea water, and the density difference between the river plume and the sea water will cause a lateral expansion of the plume. In the far-field of the plume, the thickness of the plume tends to be stable and the mixing between the two water layers is weakened. The plan view and side view of the river plume in the qualitative experiments are illustrated in <xref ref-type="fig" rid="f2">
<bold>Figures&#xa0;2</bold>
</xref>, <xref ref-type="fig" rid="f3">
<bold>3</bold>
</xref>. All sets of the views were captured when the river plume was steady.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Plan and side view of plume with an ambient water density of 0.9997g/cm<sup>3</sup>: plan-view of Case ND-1,ND-2 and ND-3 <bold>(A&#x2013;C)</bold>,side-view of Case ND-1,ND-2 and ND-3 <bold>(D&#x2013;F)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Plan and side view of plume with the flow rate of 2.5lpm: plan-view of Case ND-3, HD-3 and ED-3 <bold>(A&#x2013;C)</bold>,side-view of Case ND-3, HD-3 and ED-3 <bold>(D&#x2013;F)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g003.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref> shows the qualitative description for cases ND-1, ND-2, and ND-3. The Froude number <italic>Fr<sub>i</sub>
</italic> of these three cases are larger than those in other cases; ND-3 is the largest of all groups and the flow is supercritical. For the plan view of these three cases, the river plume appears as jet-like currents with large offshore velocities and begins to show the tendency to expand laterally after maintaining the jet-like state. The change of plume shape caused by the increase of initial velocity shows that the lateral diffusion trend of plume decreases slightly within the observation range.</p>
<p>For the side view of these three cases, it is found that the plume maintains a steady state with a relatively constant thickness after a certain downward intrusion at the estuary. The stable thickness of the river plume in case ND-1 is obviously smaller than that in case ND-2 and case ND-3, and the stable thickness of the river plume in ND-3 is slightly greater than that in case ND-2. Comparing the pattern of mixing between the river plume and the ambient fluid, the mixing in case ND-1 is the weakest blending, and the interface between the river plume and the ambient fluid was the smoothest and most stable in the three cases. The mixing between the river plume and the ambient fluid in case ND-2 and case ND-3 at the interface was significantly stronger. The Kelvin&#x2013;Helmholtz instability at the interface was more obvious in the two cases. In the observation of qualitative experiments, it is found that a more obvious vortex is formed when the front of the river plume interacts with the ambient fluid.</p>
<p>
<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref> is the plan view and side view of the river plume at three ambient fluid densities at a flow rate of 2.5 LPM, namely case ND-3, case HD-3, and case ED-3. Among the three cases, ND-3 has the largest <italic>Fr<sub>i</sub>
</italic>, followed by <italic>Fr<sub>i</sub>
</italic> of HD-3, with <italic>Fr<sub>i</sub>
</italic> of ED-3 being the smallest.</p>
<p>In the sets of plan view, unlike case ND-3(high <italic>Fr<sub>i</sub>
</italic>) which shows a strong jet-like state, flows of case HD-3 and case ED-3 only maintain the offshore direction at the central axis and within a certain width, while the edge of the plume shows a clear lateral spread trend. In other words, the alongshore movement of the two plumes are significantly stronger than that in case ND-3. The higher the ambient fluid density is, the more obvious the lateral spread trend of the river plume is, and the river plume shape changes from a &#x201c;trumpet&#x201d; shape to &#x201c;fan&#x201d; shape.</p>
<p>In the sets of side views, the plume pattern of case HD-3 and case ED-3 are quite similar and differ greatly from case ND-3, and the Kelvin&#x2013;Helmholtz instability at the interface gradually decreases with the increase of the ambient fluid density. That is to say, the increasing of the ambient fluid density limits the mixing at the interface between the river plume and the ambient fluid.</p>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Flow field of river plume</title>
<p>The flow field of the river plume in all cases is shown in <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>. It is obvious that the flow field of the river plume is distributed with the maximum velocity at the surface layer. The velocity decreases gradually along the downward water depth, and the velocity of each layer of the river plume decreases gradually along the entrance toward the offshore direction. With the increase of the flow rate, the overall velocity of the river plume increases (the average surface velocities of HD-1, HD-2, and HD-3 are 2.2, 3.4, and 4.7 cm/s respectively); the high velocity range of plume also increases, while the low velocity range changes a little.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Flow field of the river plume in all cases: Case ND-1 <bold>(A)</bold>, Case ND-2 <bold>(B)</bold>, Case ND-3 <bold>(C)</bold>,Case HD-1 <bold>(D)</bold>, Case HD-2 <bold>(E)</bold>, Case HD-3 <bold>(F)</bold>, Case ED-1 <bold>(G)</bold>, Case ED-2 <bold>(H)</bold> and Case ED-3 <bold>(I)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g004.tif"/>
</fig>
<p>Due to the slope of the rectangular channel being 0.16, in case ND-2 and case ND-3, the direction of velocity vectors at the middle and bottom layer presents a downward direction deviating from the horizontal in the near mouth, but after moving for a distance under the buoyancy effect (relatively smaller than in HD and ED groups), they gradually change to deviating from the horizontal to the upward direction.</p>
<p>For the plume in case HD and ED groups, their velocity vectors all deviate from the horizontal upward direction near the river mouth. After moving to the surface, the velocity vector turns to the horizontal direction. It indicates that the buoyancy effect on the plume is obviously greater than the effect of its inertial force in these cases, and it is consistent with the phenomena in the plan view of <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>. Under the action of buoyancy, the flow in the lower layer gradually rises to the surface, and the fluid that was originally on the surface of the plume is pushed to two flanks.</p>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Gaussian distribution of horizontal velocity</title>
<p>
<xref ref-type="bibr" rid="B17">Zhao et&#xa0;al. (2019)</xref> fitted the vertical velocity of the vertically discharged buoyant flow with a Gaussian distribution curve and the vertical velocity in different water depths conforms to the curve well. The horizontal velocity distribution along the water depth of the river plume in this study also follows the Gaussian distribution well. The river plume has a free surface in the x-z plane, and the surface velocity is maximum. Therefore, this paper selects half of the Gaussian distribution to fit the horizontal velocity of the river plume and uses the axis of the Gaussian distribution to correspond to the horizontal velocity of the free surface, obtaining a good fitting effect. <xref ref-type="bibr" rid="B16">Yuan et&#xa0;al. (2018)</xref> used the variance of Gaussian distribution to define the width of plume in plan-view (b=C&#x3c3;). According to the experimental observation and freshwater conservation, C was set to 4, which was consistent with the width of buoyancy flow. In this study, we use C=2.2 to represent the horizontal velocity distribution width of each section because half of the Gaussian distribution was fitted with the data.</p>
<p>
<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref> shows the fitting result of the horizontal velocity of each section for case ND-1, case ND-2, and case ND-3. The horizontal velocity of each section in the three cases increases significantly with the increase of flow rate and decreases with the increase of offshore distance. The distribution width of horizontal velocity also decreases along the river plume movement direction.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Gaussian distribution of horizontal velocity of ND groups: Case ND-1 <bold>(A)</bold>, Case ND-2 <bold>(B)</bold> and Case ND-3 <bold>(C)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g005.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref> shows the fitting result of the horizontal velocity of each section for case ND-3, case HD-3, and case ED-3 under the condition that the plume flow rate is 2.5 LPM. The horizontal velocities of each section in the three groups change a little with the density of the environmental fluid, while the distribution width of the horizontal velocities change greatly. The increase of environmental fluid density reduces the width of horizontal velocity distribution at the same location, but each section still conforms well to the Gaussian distribution.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Gaussian distribution of horizontal velocity of D-3 groups: Case ND-3 <bold>(A)</bold>, Case HD-3 <bold>(B)</bold> and Case ED-3 <bold>(C)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g006.tif"/>
</fig>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Thickness of the river plume</title>
<p>The water depth where the horizontal velocity drops to 1/e of the surface velocity is defined as the plume thickness in this study, which is similar to the definition of the half-width with a free jet. Starting from the maximum thickness location, the thickness scatter points of the river plume on each section were fitted with different types of curve, as shown in <xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Thickness of river plume in different cases: ND group <bold>(A)</bold> and D3 group <bold>(B)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g007.tif"/>
</fig>
<p>For the cases in the ND group, the densities of environmental water are the same. The thickness change along the river plume movement direction is similar in these cases, and the curve types do not change. All of them in these three cases are the inverse tangent function type, and the increase of flow rate only causes the change of parameter A1. The thickness in case ND-2 and case ND-3 increases first and then decreases along the river plume movement direction, and the decreasing of the flow rate makes the thickness curve tend to be flat. For the cases in the D-3 group, the change of the environmental water density causes a great change in the shape of the thickness curve. The shape of the thickness curve changes from an inverse tangent curve to a log curve with the increase of the environmental water density. The thickness curves of case HD-3 and case ED-3 are in the form of log curves. Compared with case HD-3, the environmental water density in case ED-3 is higher, resulting in the maximum thickness of plumes being smaller, the stable thickness being smaller, and the decreasing rate of the thickness from the maximum thickness to the stable thickness being significantly smaller.</p>
<p>Under the effects of the buoyancy and the inertial force, the upper and lower layer of water shear with each other and the plume thickness changes along the river plume direction. <xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref> shows the relationship between the thickness of river plume after reaching stability with the Froude number <italic>Fr<sub>i</sub>
</italic>. The Froude number of case ED and HD groups is relatively small, and the stable thickness of the river plume will be small. The scattered points lie in the lower left corner. The Froude number of case ND group is relatively large, and as a result the stable thickness of river plume is also large. An obvious linear relationship could be found between the stable thickness of the river plume with the Froude number. It indicates that the Froude number is a key parameter to describe the behaviors of the river plume.</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>The relationship between the stable thickness and the Froude number.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g008.tif"/>
</fig>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Turbulent kinetic energy</title>
<p>Turbulent kinetic energy is an important parameter reflecting the turbulent characteristics of river plume. Turbulent kinetic energy is defined as follows:</p>
<disp-formula>
<label>(4)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:mtext>k</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msubsup>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>'</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>'</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(5)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:msubsup>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>'</mml:mo>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(6)</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>U</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>'</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(7)</label>
<mml:math display="block" id="M7">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>'</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>v</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, k is turbulent kinetic energy and <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:msubsup>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mtext>&#xa0;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> is the fluctuating velocity, defined as the D-value of velocity <italic>u<sub>i</sub>
</italic> and average speed <inline-formula>
<mml:math display="inline" id="im3">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Therefore, the turbulent kinetic energy at the measuring point is <italic>k<sub>i =</sub> k<sub>Ui</sub> + k<sub>Vi</sub>
</italic>.</p>
<p>In the flow field, the turbulent kinetic energy <italic>k<sub>i</sub>
</italic> is used to represent the velocity fluctuating energy, that is, if the turbulent kinetic energy at a certain position is higher, there will be greater velocity fluctuating, and more energy transfer and dissipation.</p>
<p>We chose the section where the plume thickness is the largest as an example. The distribution of turbulent kinetic energy along the water depth in such a section was drawn in <xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9</bold>
</xref>. It is obvious that the increase of flow rate will cause the increase of turbulent kinetic energy. In each case, the maximum value of turbulent kinetic energy is located in the middle layer, and the turbulent kinetic energy gradually increases from the surface to the middle layer and then decreases to the bottom. The turbulent kinetic energy of the surface layer is obviously smaller than that of the middle layer. At the bottom boundary of the plume, the mixing of plume makes the velocity in this area reduce significantly. The increase of environmental water density will also enhance the turbulent kinetic energy. <xref ref-type="bibr" rid="B8">MacDonald and Geyer (2004)</xref> measured the velocities at three passes at the estuarine front in Fraser River, and the distribution shape of TKE production is quite similar to this study. TKE of the pass with higher salinity is significantly larger, and the peak is also closer to the upper layer. The results here verify the distribution pattern of TKE, and show the change of TKE under the influence of environmental water density and flow rate. A higher TKE production implies that the mixing between the river plume and the ocean will be increased and will enhance the exchange of substances.</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>Turbulent kinetic energy at the maximum thickness section: ND group <bold>(A)</bold>, HD group <bold>(B)</bold> and ED group <bold>(C)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1121833-g009.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="conclusions">
<label>4</label>
<title>Conclusions</title>
<p>In this study, the mixing process and movement characteristics of a river plume under the effect of different buoyancies and the plume flow rates are comprehensively studied using particle image velocimetry (PIV) technology. The main conclusions are summarized as follows.</p>
<list list-type="simple">
<list-item>
<p>(1) The change of flow rate and environmental water density will shape the river plume pattern in both the plan view and the side view. The combined effect of flow rate and environmental water density could be described by the Froude number.</p>
</list-item>
<list-item>
<p>(2) The river plume has a free surface in the x-z plane and the horizontal velocity of the plume could be fitted with a 1/2 Gaussian distribution curve.</p>
</list-item>
<list-item>
<p>(3) The thickness of the river plume increases with the increase of flow rate, and the increase in the density of ambient water changes the shape of the plume thickness curve from an arctan-type to log-type. The stable thickness of the river plume increases with the increase of the Froude number.</p>
</list-item>
<list-item>
<p>(4) The maximum value of turbulent kinetic energy is located in the middle layer, and the increase of the flow rate or the density difference leads to the increase in the turbulent kinetic energy.</p>
</list-item>
</list>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>ZX: conceptualization, methodology, formal analysis, and writing&#x2014;review and editing. CX: experiment, and writing&#x2014;original draft. GW: formal analysis and writing&#x2014;review and editing. JZ: formal analysis, visualization, and writing&#x2014;review. YC: writing&#x2014;review and editing. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>We are sincerely grateful for the support from the Key Research and Development Program of Hebei Province (21373302D), Marine Ecological Restoration and Smart Ocean Engineering Research Center of Hebei Province (HBMESO2308) and the National Natural Science Foundation of China (51979076).</p>
</sec>
<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>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Cameron</surname> <given-names>W. M.</given-names>
</name>
<name>
<surname>Pritchard</surname> <given-names>D. W.</given-names>
</name>
</person-group> (<year>1963</year>). &#x201c;<article-title>Estuaries</article-title>,&#x201d; in <source>The Sea</source>, vol. <volume>2</volume> . Ed. <person-group person-group-type="editor">
<name>
<surname>Hill</surname> <given-names>M. N.</given-names>
</name>
</person-group> (<publisher-loc>New York</publisher-loc>: <publisher-name>Wiley</publisher-name>), <fpage>306</fpage>&#x2013;<lpage>324</lpage>.</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feddersen</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Olabarrieta</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Guza</surname> <given-names>R. T.</given-names>
</name>
<name>
<surname>Winters</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Raubenheimer</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Elgar</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Observations and modeling of a tidal inlet dye tracer plume</article-title>. <source>J. Geophys. Res.: Oceans</source> <volume>121</volume> (<issue>10</issue>), <fpage>7819</fpage>&#x2013;<lpage>7844</lpage>. doi: <pub-id pub-id-type="doi">10.1002/2016JC011922</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garvine</surname> <given-names>R. W.</given-names>
</name>
</person-group> (<year>1975</year>). <article-title>The distribution of salinity and temperature in the Connecticut river estuary</article-title>. <source>J. Geophys. Res.</source> <volume>80</volume> (<issue>9</issue>), <fpage>1176</fpage>&#x2013;<lpage>1183</lpage>. doi: <pub-id pub-id-type="doi">10.1029/JC080i009p01176</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Horner-Devine</surname> <given-names>A. R.</given-names>
</name>
<name>
<surname>Jay</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Orton</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Spahn</surname> <given-names>E. Y.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>A conceptual model of the strongly tidal Columbia river plume</article-title>. <source>J. Mar. Syst.</source> <volume>78</volume> (<issue>3</issue>), <fpage>460</fpage>&#x2013;<lpage>475</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jmarsys.2008.11.025</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hsieh</surname> <given-names>S. C.</given-names>
</name>
</person-group> (<year>2008</year>). <source>Establishment of high time-resolved PIV system with application to the characteristics of a near wake flow behind a circular cylinder</source> (<publisher-loc>Taiwan</publisher-loc>: <publisher-name>National Chung Hsing University</publisher-name>).</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hsieh</surname> <given-names>S. C.</given-names>
</name>
<name>
<surname>Low</surname> <given-names>Y. M.</given-names>
</name>
<name>
<surname>Chiew</surname> <given-names>Y. M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Flow characteristics around a circular cylinder subjected to vortex-induced vibration near a plane boundary</article-title>. <source>J. Fluids Structures</source> <volume>65</volume>, <fpage>257</fpage>&#x2013;<lpage>277</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jfluidstructs.2016.06.007</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lentz</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>The response of buoyant coastal plumes to upwelling-favorable winds</article-title>. <source>J. Phys. oceanogr.</source> <volume>34</volume> (<issue>11</issue>), <fpage>2458</fpage>&#x2013;<lpage>2469</lpage>. doi: <pub-id pub-id-type="doi">10.1175/JPO2647.1</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>MacDonald</surname> <given-names>D. G.</given-names>
</name>
<name>
<surname>Geyer</surname> <given-names>W. R.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Turbulent energy production and entrainment at a highly stratified estuarine front</article-title>. <source>J. Geophys. Res.: Oceans</source> <volume>109</volume>, <fpage>C05004</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2003JC002094</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nardin</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Mariotti</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Edmonds</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Guercio</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Fagherazzi</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Growth of river mouth bars in sheltered bays in the presence of frontal waves</article-title>. <source>J. Geophys. Res.: Earth Surface</source> <volume>118</volume>, <fpage>872</fpage>&#x2013;<lpage>886</lpage>. doi: <pub-id pub-id-type="doi">10.1002/jgrf.20057</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nezlin</surname> <given-names>N. P.</given-names>
</name>
<name>
<surname>DiGiacomo</surname> <given-names>P. M.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Satellite ocean color observations of stormwater runoff plumes along the San Pedro shelf (southern California) during 1997-2003</article-title>. <source>Continental Shelf Res.</source> <volume>25</volume>, <fpage>1692</fpage>&#x2013;<lpage>1711</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.csr.2005.05.001</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yin</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y. P.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>Y. B.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>C. Y.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Z. S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>An experimental study on the evolution of a submerged berm under the effects of regular waves in low-energy conditions</article-title>. <source>Coast. Eng.</source> <volume>176</volume>, <fpage>104169</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2022.104169</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qu</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Hetland</surname> <given-names>R. D.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Temporal resolution of windforcing required for river plumesimulations</article-title>. <source>J. Geophys. Res.: Oceans</source> <volume>124</volume>, <fpage>1459</fpage>&#x2013;<lpage>1473</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2018JC014593</pub-id>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rijnsburger</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Flores</surname> <given-names>R. P.</given-names>
</name>
<name>
<surname>Pietrzak</surname> <given-names>J. D.</given-names>
</name>
<name>
<surname>Horner-Devine</surname> <given-names>A. R.</given-names>
</name>
<name>
<surname>Souza</surname> <given-names>A. J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The influence of tide and wind on the propagation of fronts in a shallow river plume</article-title>. <source>J. Geophys. Res.: Oceans</source> <volume>123</volume>, <fpage>5426</fpage>&#x2013;<lpage>5442</lpage>. doi: <pub-id pub-id-type="doi">10.1029/2017JC013422</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Avener</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Horner-Devine</surname> <given-names>A. R.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A two-color optical method for determining layer thickness in two interacting buoyant plumes</article-title>. <source>Experiments fluids</source> <volume>50</volume> (<issue>5</issue>), <fpage>1235</fpage>&#x2013;<lpage>1245</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00348-010-0969-y</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Horner-Devine</surname> <given-names>A. R.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Laboratory investigation of the impact of lateral spreading on buoyancy flux in a river plume</article-title>. <source>J. Phys. Oceanogr.</source> <volume>43</volume> (<issue>12</issue>), <fpage>2588</fpage>&#x2013;<lpage>2610</lpage>. doi: <pub-id pub-id-type="doi">10.1175/JPO-D-12-0117.1</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Horner-Devine</surname> <given-names>A. R.</given-names>
</name>
<name>
<surname>Avener</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bevan</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The role of periodically varying discharge on river plume structure and transport</article-title>. <source>Continental Shelf Res.</source> <volume>158</volume>, <fpage>15</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.csr.2018.02.009</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W.</given-names>
</name>
<name>
<surname>He</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Tan</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>The process of mass transport for a buoyant plume in linearly stratified environment</article-title>. <source>J. Shanghai Jiaotong Univ.</source> <volume>53</volume> (<issue>04</issue>), <fpage>473</fpage>&#x2013;<lpage>479</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.16183/j.cnki.jsjtu.2019.04.012</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>