<?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.1105112</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>Numerical study of sediment suspension affected by rigid cylinders under unidirectional and combined wave&#x2013;current flows</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Lou</surname>
<given-names>Sha</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1934736"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Xiaolan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhou</surname>
<given-names>Shengyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ma</surname>
<given-names>Gangfeng</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/268216"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Shuguang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/268216"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Radnaeva</surname>
<given-names>Larisa Dorzhievna</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2169164"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Nikitina</surname>
<given-names>Elena</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2147933"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fedorova</surname>
<given-names>Irina Viktorovna</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Hydraulic Engineering, Tongji University</institution>, <addr-line>Shanghai</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Key Laboratory of Yangtze River Water Environment, Ministry of Education, Tongji University</institution>, <addr-line>Shanghai</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Department of Civil and Environmental Engineering, Old Dominion University</institution>, <addr-line>Norfolk, VA</addr-line>, <country>United States</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Laboratory of Chemistry of Natural Systems, Baikal Institute of Nature Management of Siberian Branch of the Russian Academy of Sciences</institution>, <addr-line>Ulan-Ude</addr-line>, <country>Russia</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Institute of Earth Sciences, Saint Petersburg State University</institution>, <addr-line>St Petersburg</addr-line>, <country>Russia</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Wei-Bo Chen, National Science and Technology Center for Disaster Reduction(NCDR), Taiwan</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Peter M.J. Herman, Delft University of Technology, Netherlands; Pankaj Raushan, Indian Institute of Engineering Science and Technology, Shibpur, India</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Sha Lou, <email xlink:href="mailto:lousha@tongji.edu.cn">lousha@tongji.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>06</day>
<month>02</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1105112</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>11</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>24</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Lou, Chen, Zhou, Ma, Liu, Radnaeva, Nikitina and Fedorova</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Lou, Chen, Zhou, Ma, Liu, Radnaeva, Nikitina and Fedorova</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>Sediment transport modeling for flows with cylinders is very challenging owing to the complicated flow&#x2013;cylinder&#x2013;sediment interactions, especially under the combined wave-current flows. In this paper, an improved formulation for incipient sediment suspension considering the effect of cylinder density (i.e., solid volume fraction) is employed to simulate the bottom sediment flux in the flow with cylinders. The proposed model is calibrated and validated using laboratory measurements under unidirectional and combined wave-current flows in previous studies. It is proved that the effects of cylinders on sediment suspension can be accounted for through a modified critical Shields number, and the proposed model is capable of simulating sediment suspension under both unidirectional and combined wave&#x2013;current flows reasonably well with the average the coefficients of determination and model skills greater than 0.8 and 0.64.</p>
</abstract>
<kwd-group>
<kwd>sediment suspension</kwd>
<kwd>modified critical Shields number</kwd>
<kwd>cylinder</kwd>
<kwd>unidirectional flows</kwd>
<kwd>combined wave&#x2013;current flows</kwd>
</kwd-group>
<contract-num rid="cn002">20230742500, 22ZR1464200</contract-num>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Shanghai Technology Innovation Project<named-content content-type="fundref-id">10.13039/501100019445</named-content>
</contract-sponsor>
<contract-sponsor id="cn003">Fundamental Research Funds for the Central Universities<named-content content-type="fundref-id">10.13039/501100012226</named-content>
</contract-sponsor>
<counts>
<fig-count count="9"/>
<table-count count="4"/>
<equation-count count="18"/>
<ref-count count="50"/>
<page-count count="12"/>
<word-count count="6683"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>In estuarine and coastal areas, various types of flows involving cylinders exist. These cylinders, such as aquatic vegetation in wetlands (mangrove, reed, and seaweed/kelp), masts in aquafarm for laver or oysters, and piles foundations of offshore windfarms, can significantly affect flow structure (<xref ref-type="bibr" rid="B25">Nepf and Vivoni, 2000</xref>; <xref ref-type="bibr" rid="B24">Nepf, 2012</xref>; <xref ref-type="bibr" rid="B15">Lou et&#xa0;al., 2018</xref>), wave propagation (<xref ref-type="bibr" rid="B22">Moller et&#xa0;al., 1999</xref>; <xref ref-type="bibr" rid="B21">Mendez and Losada, 2004</xref>; <xref ref-type="bibr" rid="B1">Augustin et&#xa0;al., 2009</xref>), and turbulence development (<xref ref-type="bibr" rid="B23">Nepf, 1999</xref>; <xref ref-type="bibr" rid="B5">Ghisalberti and Nepf, 2002</xref>; <xref ref-type="bibr" rid="B26">Neumeier, 2007</xref>). The cylinders reduce mean flow and convert mean kinetic energy to turbulent kinetic energy (<italic>TKE</italic>), which have great influence on sediment suspension and transport (<xref ref-type="bibr" rid="B2">Balke et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B29">Ros et&#xa0;al., 2014</xref>). Furthermore, waves and currents usually coexist in coastal zones. Wave-current interactions play an important role in coastal ocean dynamics and nearshore sediment transport (<xref ref-type="bibr" rid="B41">Van Hoften and Karaki, 1976</xref>). With the presence of cylinders, sediment suspension under the combined wave-current flow become more complicated, which makes it more difficult to simulate sediment transport in flow with cylinders.</p>
<p>Field observations indicated that the presence of cylinders reduced resuspension and promoted particle retention, for example, studies in Lake Hiidenvesi (<xref ref-type="bibr" rid="B9">Horppila and Nurminen, 2003</xref>) and Lake Taihu (<xref ref-type="bibr" rid="B50">Zhu et&#xa0;al., 2015</xref>). On the other hand, laboratory studies provided contrary results and showed that sediment resuspension was enhanced in the flows with cylinders comparing to bare bed condition (<xref ref-type="bibr" rid="B38">Tinoco and Coco, 2016</xref>; <xref ref-type="bibr" rid="B39">Tinoco and Coco, 2018</xref>; <xref ref-type="bibr" rid="B14">Lou et&#xa0;al., 2022</xref>). These contrasting findings reflected the influence of cylinders on the flow at two scales: array- and stem- scales (<xref ref-type="bibr" rid="B49">Zhang and Nepf, 2019</xref>). In the field, the expanse of cylinders (e.g. the length of region covered by vegetation in wetlands in the order of hundred meters) can lead to the damping of flow and sediment resuspension in the array-scale. In the laboratory, the cylinders enhance sediment resuspension due to stem-generated turbulence. It is found out that cylinder density significantly affects the incipient motion of the sediment (<xref ref-type="bibr" rid="B38">Tinoco and Coco, 2016</xref>; <xref ref-type="bibr" rid="B39">Tinoco and Coco, 2018</xref>; <xref ref-type="bibr" rid="B35">Tang et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B48">Yang and Nepf, 2019</xref>). Although the flow velocity is significantly damped with dense cylinders, the amount of sediment lifted into suspension increases owing to the increasing turbulent kinetic energy under unidirectional currents (<xref ref-type="bibr" rid="B38">Tinoco and Coco, 2016</xref>; <xref ref-type="bibr" rid="B48">Yang and Nepf, 2019</xref>) and oscillatory flows (<xref ref-type="bibr" rid="B29">Ros et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B35">Tang et&#xa0;al., 2019</xref>). The incipient velocity of sediment motion decreases within cylinder canopy (<xref ref-type="bibr" rid="B36">Tang et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B47">Yang et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B39">Tinoco and Coco, 2018</xref>).</p>
<p>Although field investigations and experimental studies provide a certain basis for the establishment of numerical models, sediment transport modeling for flows with cylinders is more challenging owing to complicated flow&#x2013;cylinder&#x2013;sediment interactions. Sediment suspension determined by the bottom sediment flux was related to the turbulent flow structure altered by cylinders. <xref ref-type="bibr" rid="B13">Lopez and Garcia (1998)</xref> established a two-equation turbulence model to determine the mean flow and turbulence structure of open channels through cylinders and used this model to estimate sediment transport processes, assuming that the average sediment deposition near the bed and the entrainment rate were in equilibrium. A three-dimensional turbulence model was proposed by <xref ref-type="bibr" rid="B11">Liu and Shen (2008)</xref> by introducing cylinder density and drag force into the control equations of flow in the presence of cylinders, and an empirical relation proposed by <xref ref-type="bibr" rid="B42">Van Rijn (1984)</xref> was employed to solve the sediment transport equations in flow affected by the cylinders. With an enhanced Manning&#x2019;s roughness coefficient and cylinder-induced drag force representing the effect of the cylinders on flow, sediment flux was computed with the formula of <xref ref-type="bibr" rid="B43">Van Rijn (1993)</xref>, disregarding the effect of the cylinders on sediment suspension (<xref ref-type="bibr" rid="B10">Hu et&#xa0;al., 2018</xref>). The equilibrium sediment transport rates by Wu&#x2019;s formula (<xref ref-type="bibr" rid="B46">Wu and Wang, 2007</xref>) were used for bed load, and the suspended sediment flux was modeled by bed shear stress estimated from mean flow velocity with the cylinder effects implemented by additional terms of the cylinder drag and water volume reduction caused by the cylinders (<xref ref-type="bibr" rid="B8">He et&#xa0;al., 2017</xref>). A similar method was used in a sediment transport model established by <xref ref-type="bibr" rid="B19">Ma et&#xa0;al. (2013)</xref>, in which the sediment flux was computed by bed shear stress estimated from mean flow velocity with cylinder effects implemented by additional terms of the cylinder drag. However, in flow with cylinders, the generation of turbulence in the wakes of cylinders exceeds that associated with bed shear, such that the previously used bare bed sediment transport models based on bed shear stress need to be improved to better describe sediment suspension affected by cylinders.</p>
<p>In our previous work (<xref ref-type="bibr" rid="B16">Lou et&#xa0;al., 2021</xref>), a vegetated flow and sediment transport model was developed for investigating flow&#x2013;cylinder(vegetation)&#x2013;sediment interactions. In the model, the vegetation-induced turbulence was simulated by an improved <italic>k-&#x3f5;</italic> turbulence closure, in which shear and stem-scale wake turbulent kinetic energies were solved separately due to the different length scales of turbulence in vegetated flow. It simulated both unidirectional and oscillatory flows as well as turbulence structures reasonably well. However, one of the disadvantages of the model was that it included more coefficients for calibration compared to traditional <italic>k-&#x3f5;</italic> turbulence closure, which introduced uncertainties to model results and restrained its applications.</p>
<p>In this paper, we propose a modified critical Shields number to account for the effects of cylinders on bottom sediment suspension, which will be implemented into the Non-Hydrostatic WAVE (NHWAVE) model (<xref ref-type="bibr" rid="B20">Ma et&#xa0;al., 2012</xref>). Purpose of the proposed model is to simplify the numerical simulation of sediment transport in flow with cylinders. With fewer parameters, the sediment initiation can be identified and simulated, which aims to simplify the numerical simulation of sediment transport in flow with cylinders, especially for the combined wave-current flow. In the proposed model, the modified critical Shields number affected by cylinders is determined by the cylinder features (solid volume fraction of cylinders <italic>&#x3d5;</italic> and the cylinder induced-drag coefficient) and the bottom drag coefficient, and the traditional <italic>k-&#x3f5;</italic> turbulence closure with a cylinder-induced turbulence production term including both shear and stem-scale turbulence is adopted. The model will be validated by using the experimental measurements collected in a wave-current flume at Tongji University in our previous study (<xref ref-type="bibr" rid="B3">Chen et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B14">Lou et&#xa0;al., 2022</xref>) and a dataset presented by <xref ref-type="bibr" rid="B17">Lu (2008)</xref>. The performance of the model will be evaluated by computing the model skill (<italic>M<sub>S</sub>
</italic>), the determination coefficient (<italic>R<sup>2</sup>
</italic>) and the Root Mean Square Error (<italic>RMSE</italic>). Effects of cylinder densities on sediment suspension will be analyzed and the interactions between flow&#x2013;cylinder&#x2013;sediment will be examined. The main purpose of this paper is to study sediment suspension mechanisms affected by cylinders under both unidirectional and the combined wave-current flows, and to estimate the accuracy of mathematical descriptions of the bottom sediment flux considering the effects of cylinders on sediment suspension threshold condition.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Model formulations</title>
<p>The hydrodynamic and sediment transport model in this study is based on the Non-Hydrostatic WAVE (NHWAVE) model, which was originally developed by <xref ref-type="bibr" rid="B20">Ma et&#xa0;al. (2012)</xref>. It has been well validated and widely applied to study coastal wave processes, including wave/flow&#x2013;cylinder&#x2013;sediment interactions (<xref ref-type="bibr" rid="B19">Ma et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B18">Ma et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B45">Wu et&#xa0;al., 2016</xref>).</p>
<sec id="s2_1">
<label>2.1</label>
<title>NHWAVE model</title>    <p>The vegetation module in NHWAVE is based on the Volume-averaged and Reynolds-averaged Navier-Stokes equations. The effects of the cylinders on the mean flow were modeled by drag forces in the momentum equations. This treatment has been proven successful in simulating turbulent structures in vegetated flows by other modeling efforts (<xref ref-type="bibr" rid="B30">Sheng et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B45">Wu et&#xa0;al., 2016</xref>). The governing equations, numerical schemes, and boundary conditions were introduced in <xref ref-type="bibr" rid="B19">Ma et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B20">Ma et&#xa0;al. (2012</xref>). The drag force is given by</p>
<disp-formula>
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mi>a</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mi>u</mml:mi> <mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>C<sub>D</sub>
</italic> is the cylinder drag coefficient, <italic>a</italic> is the frontal area of cylinders per unit volume, <italic>u</italic> is the velocity vector.</p>
<p>A nonlinear <italic>k</italic>&#x2013;<italic>&#x3f5;</italic> turbulence model was implemented in the NHWAVE model to simulate turbulent flows. According to <xref ref-type="bibr" rid="B23">Nepf (1999)</xref>, the large-scale mean kinetic energy is converted to small-scale turbulent kinetic energy in flows with cylinders, and turbulent diffusivity is reduced owing to a downward shift in turbulent-length scale. In the NHWAVE model, these processes were modeled by introducing a cylinder-induced turbulence production term, <italic>P<sub>v</sub>
</italic>, in the <italic>k</italic>&#x2013;<italic>&#x3f5;</italic> equations (<xref ref-type="bibr" rid="B20">Ma et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B19">Ma et&#xa0;al., 2013</xref>). The modified <italic>k</italic>&#x2013;<italic>&#x3f5;</italic> equations are written as</p>
<disp-formula>
<label>(2)</label>
<mml:math display="block" id="M2">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mo>&#x2207;</mml:mo>
<mml:mo>&#xb7;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>k</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>=</mml:mo>
<mml:mo>&#x2207;</mml:mo>
<mml:mo>&#xb7;</mml:mo>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>&#x2207;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo stretchy="false">]</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>+</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mi>&#x3f5;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<disp-formula>
<label>(3)</label>
<mml:math display="block" id="M3">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mo>&#x2207;</mml:mo>
<mml:mo>&#xb7;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>=</mml:mo>
<mml:mo>&#x2207;</mml:mo>
<mml:mo>&#xb7;</mml:mo>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>&#x2207;</mml:mo>
<mml:mi>&#x3f5;</mml:mi>
<mml:mo stretchy="false">]</mml:mo>
<mml:mo>+</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mi>&#x3f5;</mml:mi>
<mml:mi>k</mml:mi>
</mml:mfrac>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<disp-formula>
<label>(4)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mi>u</mml:mi> <mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(5)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>&#x3bc;</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mi>&#x3f5;</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>k</italic> is the turbulence kinetic energy; <italic>&#x3f5;</italic> is the energy dissipation; <italic>D</italic> is the total water depth; <italic>P<sub>s</sub>
</italic> and <italic>P<sub>b</sub>
</italic> are shear and buoyancy productions, respectively; <italic>&#x3c5;<sub>t</sub>
</italic> is the turbulent eddy viscosity and <italic>&#x3c5;</italic> is the kinematic water viscosity; <italic>&#x3c3;<sub>k</sub>
</italic>(=1.0), <italic>&#x3c3;<sub>&#x3f5;</sub>
</italic>(=1.3), <italic>C<sub>1&#x3f5;</sub>
</italic>(=1.44), <italic>C<sub>2&#x3f5;</sub>
</italic>(=1.92), and <italic>C<sub>&#x3bc;</sub>
</italic>(=0.09) are empirical coefficients (<xref ref-type="bibr" rid="B28">Rodi, 1987</xref>); the coefficient <italic>C<sub>3&#x3f5;</sub>
</italic> is selected as zero for a stably stratified water column (<xref ref-type="bibr" rid="B31">Snyder and Hsu, 2011</xref>); <italic>C<sub>fk</sub>
</italic>(=1.0) and <italic>C<sub>f&#x3f5;</sub>
</italic>(=1.33) are drag-related coefficients based on <xref ref-type="bibr" rid="B13">Lopez and Garcia (1998)</xref>.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Sediment transport model</title>
<p>To study sediment dynamics in flows with cylinders, suspended sediment concentration was computed from the advection&#x2013;diffusion equation formulated in <italic>&#x3c3;</italic> coordinates in the NHWAVE model (<xref ref-type="bibr" rid="B18">Ma et&#xa0;al., 2015</xref>).</p>
<disp-formula>
<label>(6)</label>
<mml:math display="block" id="M6">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">]</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">]</mml:mo>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>D</mml:mi>
</mml:mfrac>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">[</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3c5;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">]</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>where <italic>C</italic> is the suspended sediment concentration; <italic>&#x3c9;<sub>s</sub>
</italic> is the sediment settling velocity; <italic>&#x3c3;<sub>h</sub>
</italic> and <italic>&#x3c3;<sub>v</sub>
</italic> are the horizontal and vertical Schmidt numbers of the sediment, respectively. To solve the sediment transport equation, a no-flux boundary condition was employed at the free surface.</p>
<disp-formula>
<label>(7)</label>
<mml:math display="block" id="M7">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>D</mml:mi>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mi>C</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</disp-formula>
<p>At the bottom, the mass exchange of the suspended sediment was considered by the erosion and deposition fluxes:</p>
<disp-formula>
<label>(8)</label>
<mml:math display="block" id="M8">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3c5;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c5;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>D</mml:mi>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mi>C</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>E</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>E</italic> and <italic>D</italic> are sediment fluxes for erosion and deposition, respectively.</p>
<p>For non-cohesive sediments, a pickup function (<xref ref-type="bibr" rid="B42">Van Rijn, 1984</xref>) is typically used to simulate the fluxes, expressed as</p>
<disp-formula>
<label>(9)</label>
<mml:math display="block" id="M9">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi>E</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0.00033</mml:mn>
<mml:mo>&#xb7;</mml:mo>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mn>50</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>0.8</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#xb7;</mml:mo>
<mml:msup>
<mml:mo stretchy="false">(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
</mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
</mml:mfrac>
</mml:msup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mn>0.6</mml:mn>
<mml:mo>&#xb7;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msup>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mn>06</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#xb7;</mml:mo>
<mml:msup>
<mml:mi>&#x3c5;</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#xb7;</mml:mo>
<mml:msup>
<mml:mo stretchy="false">(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:msup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mrow>
<mml:mn>1.5</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>where <italic>D</italic>
<sub>50</sub> is the median grain size, <italic>&#x3c1;<sub>s</sub>
</italic> is the sediment density, <italic>&#x3c1;</italic> is the water density, <italic>&#x3c8;</italic> is the Shields number given by</p>
<disp-formula>
<label>(10)</label>
<mml:math display="block" id="M10">
<mml:mrow>
<mml:mi>&#x3c8;</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mn>50</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the bottom shear stress, <italic>&#x3c8;<sub>c</sub>
</italic> is the critical Shields number for incipient sediment motion.</p>
<p>For the combined wave&#x2013;current flow, the bottom shear stress <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is expressed as <inline-formula>
<mml:math display="inline" id="im3">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, calculated using Eq. (11) considering the nonlinear interactions between waves and currents (<xref ref-type="bibr" rid="B33">Soulsby, 1995</xref>; <xref ref-type="bibr" rid="B34">Soulsby, 1997</xref>).</p>
<disp-formula>
<label>(11)</label>
<mml:math display="block" id="M11">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>1.2</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>w</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>w</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>3.2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>f<sub>w</sub>
</italic> is the wave friction coefficient, <italic>f<sub>w</sub>
</italic>=1.39(<italic>A/z</italic>
<sub>0</sub>)<sup>-0.52</sup>, in which <inline-formula>
<mml:math display="inline" id="im4">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im5">
<mml:mrow>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2.5</mml:mn>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mn>50</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>30</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, <italic>T</italic> is the wave period, and <italic>U<sub>w</sub>
</italic> is the velocity at the wave boundary layer. <inline-formula>
<mml:math display="inline" id="im6">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the bottom shear stress under the unidirectional flow, calculated by</p>
<disp-formula>
<label>(12)</label>
<mml:math display="block" id="M12">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mtext>&#xa0;</mml:mtext>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow> <mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>c<sub>d</sub>
</italic> is the bottom friction coefficient, and <italic>U<sub>c</sub>
</italic> is the mean velocity near the bottom. <inline-formula>
<mml:math display="inline" id="im7">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the bottom shear stress under a pure wave condition, computed as</p>
<disp-formula>
<label>(13)</label>
<mml:math display="block" id="M13">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>w</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>w</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>In a flow without cylinders, the critical Shields number, <italic>&#x3c8;<sub>c</sub>
</italic>, is considered to be a constant which is related to the sediment characteristics, which can be calculated by</p>
<disp-formula>
<label>(14)</label>
<mml:math display="block" id="M14">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>0.15</mml:mn>
<mml:msup>
<mml:mo stretchy="false">(</mml:mo>
<mml:msup>
<mml:mi>S</mml:mi>
<mml:mo>*</mml:mo>
</mml:msup>
</mml:msup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.3</mml:mn>
</mml:mrow>
<mml:mo>+</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mn>0.045</mml:mn>
<mml:mi>exp</mml:mi>
<mml:mo stretchy="false">[</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>35</mml:mn>
<mml:msup>
<mml:mo stretchy="false">(</mml:mo>
<mml:msup>
<mml:mi>S</mml:mi>
<mml:mo>*</mml:mo>
</mml:msup>
</mml:msup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.59</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>with the dimensionless viscosity <inline-formula>
<mml:math display="inline" id="im8">
<mml:mrow>
<mml:msup>
<mml:mi>S</mml:mi>
<mml:mo>*</mml:mo>
</mml:msup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">/</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo stretchy="false">]</mml:mo>
<mml:mi>g</mml:mi>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mn>50</mml:mn>
<mml:mn>3</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>However, it has been reported that the incipient velocity decreases under both unidirectional flows and wave conditions affected by cylinders (<xref ref-type="bibr" rid="B36">Tang et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B29">Ros et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B35">Tang et&#xa0;al., 2019</xref>). The turbulence, instead of mean flow, determines sediment resuspension through arrays of rigid vegetation (<xref ref-type="bibr" rid="B38">Tinoco and Coco, 2016</xref>; <xref ref-type="bibr" rid="B47">Yang et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B39">Tinoco and Coco, 2018</xref>). Hence, the sediment transport models based on bed shear stress developed for bare bed do not work well for flows with cylinders.</p>
<p>Based on experimental measurements, <xref ref-type="bibr" rid="B38">Tinoco and Coco (2016)</xref> proposed a function to estimate the critical Shields number for unidirectional flow with cylinders (Eq. 15), which was modified for wave conditions (Eq. 16) in <xref ref-type="bibr" rid="B39">Tinoco and Coco (2018)</xref> based on the study in <xref ref-type="bibr" rid="B47">Yang et&#xa0;al. (2016)</xref>. These functions were obtained under the conditions of vegetation solid volume fraction <italic>&#x3a6;</italic>=0.008-0.079, and stem Reynolds number <italic>R<sub>ed</sub>
</italic>=200-20000 for unidirectional flow, and wave period <italic>T</italic>=2.5s, wave height <italic>H</italic>=0.01&#x2013;0.10m under wave conditions.</p>
<disp-formula>
<label>(15)</label>
<mml:math display="block" id="M15">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>37</mml:mn>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(16)</label>
<mml:math display="block" id="M16">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>0.9</mml:mn>
<mml:mo>&#xb7;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>&#x3d5;</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#xb7;</mml:mo>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">/</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>&#x3c8;<sub>c,</sub>
</italic>
<sub>0</sub> is the critical Shields number in the flows without cylinders, <italic>&#x3c8;<sub>c,&#x3d5;</sub>
</italic> is the modified critical Shields number affected by cylinders, and <italic>&#x3d5;</italic> is the solid volume fraction of cylinders.</p>
<p>Ratios of obstructed to unobstructed critical Shields values in the flows with cylinders using Eq. 15 and Eq. 16 are displayed in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>. It is shown that the values of <italic>&#x3c8;</italic>
<sub>
<italic>c</italic>,<italic>&#x3d5;</italic>
</sub>/<italic>&#x3c8;</italic>
<sub>
<italic>c</italic>,0</sub> are approximately the same when the solid volume fraction <italic>&#x3d5;</italic> is smaller than 0.03.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Ratios of obstructed to unobstructed critical Shields values in the flow with cylinders.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g001.tif"/>
</fig>
<p>To demonstrate the effects of cylinders on sediment suspension, numerical simulations were also carried out using the traditional critical Shield number without cylinders (Eq. 14). The predicted suspended sediment concentrations from the improved model (Eq. 15 and 16) and the traditional model (Eq. 14) under both unidirectional and combined wave&#x2013;current flows were compared.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Assessment of model performance</title>
<p>The developed model was validated using the experimental measurements, and the model performance was assessed by calculating the model skill (<italic>M<sub>S</sub>
</italic>), the coefficient of determination (<italic>R<sup>2</sup>
</italic>) and the Root Mean Square Error (<italic>RMSE</italic>). The model skill (<xref ref-type="bibr" rid="B44">Wilmott, 1981</xref>) describes the degree to which the measured deviations about the measured mean correspond to the simulated deviations about the measured mean, as shown in Eq. 17. Perfect agreement between the model results and measurements yields <italic>M<sub>S</sub>
</italic> of 1.0 whereas complete disagreement yields <italic>M<sub>S</sub>
</italic> of 0.</p>
<disp-formula>
<label>(17)</label>
<mml:math display="block" id="M17">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">M</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mi>&#x0020;</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>&#x0020;</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x0020;</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<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>I</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow> <mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<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>I</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow> <mml:mo>|</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
</mml:mrow>
</mml:mstyle>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow> <mml:mo>|</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>X<sub>S</sub>
</italic> and <italic>X<sub>M</sub>
</italic> are simulated and measured variables being compared, <inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> is its sample mean, <italic>I</italic> is the number of the data.</p>
<p>The coefficient of determination (<italic>R<sup>2</sup>
</italic>) is a number that indicates how well the data fit a statistical model (<xref ref-type="bibr" rid="B4">Gao et&#xa0;al., 2018</xref>). The coefficient of determination for a linear regression model with one independent variable is given by</p>
<disp-formula>
<label>(18)</label>
<mml:math display="block" id="M18">
<mml:mrow>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mtext>&#xa0;=&#xa0;[</mml:mtext>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>N</mml:mi>
</mml:mfrac>
<mml:mo>&#xb7;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<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:munderover>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow> <mml:mo>|</mml:mo>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
<mml:mrow>
<mml:mo>|</mml:mo> <mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow> <mml:mo>|</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#xb7;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mo stretchy="false">]</mml:mo>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>&#x3c3;<sub>xs</sub>
</italic> is the standard deviation of <italic>X<sub>S</sub>
</italic>, and <italic>&#x3c3;</italic>
<sub>
<italic>X</italic>
<sub>
<italic>M</italic>
</sub>
</sub> is the standard deviation of <italic>X<sub>M</sub>
</italic>. A <italic>R<sup>2</sup>
</italic> of 1 indicates that the regression line perfectly fits the data, while a <italic>R<sup>2</sup>
</italic> of 0 indicates that the line does not fit the data at all.</p>
<p>The root mean square error (<italic>RMSE</italic>), which is a typical metric for measuring the differences between samples and estimators, was also calculated as an additional assessment for the vertically varying density case.</p>
</sec>
</sec>
<sec id="s3">
<label>3</label>
<title>Model calibration and validation</title>
<p>To evaluate the performance of NHWAVE with the improved sediment transport model, a number of experimental datasets were chosen for model calibration and validation. Experiments in our previous studies (<xref ref-type="bibr" rid="B3">Chen et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B14">Lou et&#xa0;al., 2022</xref>) were adopted to calibrate the model, including conditions of unidirectional flow (Test 1) and combined wave-current flows (Test 2-4). Another dataset collected in the experiments of sediment suspension under directional flow affected by rigid submerged cylinders (<xref ref-type="bibr" rid="B17">Lu, 2008</xref>) was used to further validate the model, which has also been applied to validate the model in <xref ref-type="bibr" rid="B12">Li et&#xa0;al. (2020)</xref>. The hydraulic conditions and characteristics of cylinders in these experiments were listed in <xref ref-type="table" rid="T1">
<bold>Tables&#xa0;1</bold>
</xref>&#x2013;<xref ref-type="table" rid="T3">
<bold>3</bold>
</xref>. The cylinders were arrayed regularly in the experiments.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Hydrodynamic conditions in experiments conducted by <xref ref-type="bibr" rid="B3">Chen et&#xa0;al. (2020)</xref> and <xref ref-type="bibr" rid="B14">Lou et&#xa0;al. (2022)</xref> .</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Run number</th>
<th valign="middle" align="center">Hydrodynamic state</th>
<th valign="middle" align="center">Hydrodynamic parameters</th>
<th valign="middle" align="center">Water depths</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">Test 1</td>
<td valign="middle" align="center">Unidirectional currents</td>
<td valign="middle" align="center">
<italic>Q</italic> = 140 L/s</td>
<td valign="middle" align="center">0.8 m</td>
</tr>
<tr>
<td valign="middle" align="center">Test 2</td>
<td valign="middle" rowspan="3" align="center">Combined wave&#x2013;current flows</td>
<td valign="middle" align="center">Combination of Test 1 and wave<break/>(<italic>H</italic> = 0.06&#xa0;m, <italic>T</italic> = 1.8 s)</td>
<td valign="middle" align="center">0.8 m</td>
</tr>
<tr>
<td valign="middle" align="center">Test 3</td>
<td valign="middle" align="center">Combination of Test 1 and wave<break/>(<italic>H</italic> = 0.08&#xa0;m, <italic>T</italic> = 1.8 s)</td>
<td valign="middle" align="center">0.8 m</td>
</tr>
<tr>
<td valign="middle" align="center">Test 4</td>
<td valign="middle" align="center">Combination of Test 1 and wave<break/>(<italic>H</italic> = 0.10&#xa0;m, <italic>T</italic> = 1.8 s)</td>
<td valign="middle" align="center">0.8 m</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>H: wave height; T: wave period; Q: discharge.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Cylinder configurations in experiments conducted by <xref ref-type="bibr" rid="B3">Chen et&#xa0;al. (2020)</xref> and <xref ref-type="bibr" rid="B14">Lou et&#xa0;al. (2022)</xref>.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Cases</th>
<th valign="middle" align="center">Diameter <italic>d</italic> (mm)</th>
<th valign="middle" align="center">Height <italic>h<sub>v</sub>
</italic> (m)</th>
<th valign="middle" align="center">Stem number</th>
<th valign="middle" align="center">Density <italic>n</italic> (stems/m<sup>2</sup>)</th>
<th valign="middle" align="center">Solid volume fraction <italic>&#x3d5;</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">Case-a</td>
<td valign="middle" rowspan="6" align="center">8</td>
<td valign="middle" rowspan="2" align="center">0.4 (uniform heights and vertically uniform density)</td>
<td valign="middle" align="center">806</td>
<td valign="middle" align="center">168</td>
<td valign="middle" align="center">0.008</td>
</tr>
<tr>
<td valign="middle" align="center">Case-b</td>
<td valign="middle" align="center">1599</td>
<td valign="middle" align="center">333</td>
<td valign="middle" align="center">0.017</td>
</tr>
<tr>
<td valign="middle" rowspan="4" align="center">Case-c*</td>
<td valign="middle" rowspan="4" align="center">0.2; 0.4; 0.6; 0.8 (four cylinder heights causing the vertically varying density)</td>
<td valign="middle" align="center">369</td>
<td valign="middle" rowspan="4" align="center">333 (0&#x2264;<italic>z &#x2264;</italic> 0.2)<break/>256 (0.2&lt;<italic>z &#x2264;</italic> 0.4)<break/>166 (0.6&lt;<italic>z &#x2264;</italic> 0.8)<break/>77 (0.4&lt;<italic>z &#x2264;</italic> 0.6)</td>
<td valign="middle" rowspan="4" align="center">0.017 (0&#x2264;<italic>z &#x2264;</italic> 0.2)<break/>0.013 (0.2&lt;<italic>z &#x2264;</italic> 0.4)<break/>0.008 (0.4&lt;<italic>z &#x2264;</italic> 0.6)<break/>0.004 (0.6&lt;<italic>z &#x2264;</italic> 0.8)</td>
</tr>
<tr>
<td valign="middle" align="center">431</td>
</tr>
<tr>
<td valign="middle" align="center">430</td>
</tr>
<tr>
<td valign="middle" align="center">369</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>*Total stem number in Case-c is 1599, which is the same with that in Case-b.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Hydrodynamic conditions and cylinder configurations in experiments conducted by <xref ref-type="bibr" rid="B17">Lu (2008)</xref>.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Run<break/>number</th>
<th valign="middle" align="center">Discharge <italic>Q</italic> (L/s)</th>
<th valign="middle" align="center">Water depths (m)</th>
<th valign="middle" align="center">Cylinder Diameter <italic>d</italic> (mm)</th>
<th valign="middle" align="center">Cylinder Height <italic>h<sub>v</sub>
</italic> (m)</th>
<th valign="middle" align="center">Density <italic>n</italic> (stems/m<sup>2</sup>)</th>
<th valign="middle" align="center">Solid volume fraction <italic>&#x3d5;</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">C15-1</td>
<td valign="middle" align="center">21.00</td>
<td valign="middle" rowspan="2" align="center">0.15</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">0.06</td>
<td valign="middle" align="center">2000</td>
<td valign="middle" align="center">0.023</td>
</tr>
<tr>
<td valign="middle" align="center">C15-2</td>
<td valign="middle" align="center">19.07</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">1000</td>
<td valign="middle" align="center">0.011</td>
</tr>
<tr>
<td valign="middle" align="center">C18-1</td>
<td valign="middle" align="center">26.50</td>
<td valign="middle" rowspan="2" align="center">0.18</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">2000</td>
<td valign="middle" align="center">0.019</td>
</tr>
<tr>
<td valign="middle" align="center">C18-2</td>
<td valign="middle" align="center">23.83</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">1000</td>
<td valign="middle" align="center">0.009</td>
</tr>
<tr>
<td valign="middle" align="center">C24-1</td>
<td valign="middle" align="center">30.67</td>
<td valign="middle" rowspan="2" align="center">0.24</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">2000</td>
<td valign="middle" align="center">0.014</td>
</tr>
<tr>
<td valign="middle" align="center">C24-2</td>
<td valign="middle" align="center">29.55</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">1000</td>
<td valign="middle" align="center">0.007</td>
</tr>
<tr>
<td valign="middle" align="center">C30-1</td>
<td valign="middle" align="center">38.80</td>
<td valign="middle" rowspan="2" align="center">0.30</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">2000</td>
<td valign="middle" align="center">0.011</td>
</tr>
<tr>
<td valign="middle" align="center">C30-2</td>
<td valign="middle" align="center">36.42</td>
<td valign="middle" align="center">6</td>
<td valign="top" align="center">0.06</td>
<td valign="middle" align="center">1000</td>
<td valign="middle" align="center">0.006</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>To reproduce the vegetated flow in the experiments, the computational domain was discretized by 64 cells in the streamwise direction with intervals of 0.1&#xa0;m and 50 vertical layers. Periodic boundary conditions were imposed at the streamwise boundaries (<xref ref-type="bibr" rid="B19">Ma et&#xa0;al., 2013</xref>). The flow was driven by an external pressure gradient determined by the flow condition in the laboratory. The time step was adjusted during the simulation based on the stability restriction. Cylinder-induced bulk drag coefficient <italic>C<sub>D</sub>
</italic> is a key parameter for the simulation of flows with cylinders, which is also the main parameter of model calibration. According to laboratory studies, the value of <italic>C<sub>D</sub>
</italic> increases with cylinder density and decreases with the increasing stem Reynolds number (<xref ref-type="bibr" rid="B40">Tinoco and Cowen, 2013</xref>; <xref ref-type="bibr" rid="B27">Ni, 2014</xref>; <xref ref-type="bibr" rid="B15">Lou et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B48">Yang and Nepf, 2019</xref>), In this study, empirical equations proposed by <xref ref-type="bibr" rid="B37">Tanino and Nepf (2008)</xref>; <xref ref-type="bibr" rid="B32">Sonnenwald et&#xa0;al. (2019)</xref> and <xref ref-type="bibr" rid="B6">Ghisalberti and Nepf (2004)</xref> were adopted to estimate this coefficient.</p>
<sec id="s3_1">
<label>3.1</label>
<title>Calibrations</title>
<p>Datasets in experiments in our previous studies (<xref ref-type="bibr" rid="B3">Chen et al., 2020</xref>; <xref ref-type="bibr" rid="B14">Lou et al., 2022</xref>) were adopted to calibrate the model, including conditions of unidirectional flow (Test 1) and combined wave-current flows (Test 2-4).</p>
<p>In these simulations, the range of <italic>C<sub>D</sub>
</italic> was calculated as 0.3&#x2013;1.5 using the empirical equations (<xref ref-type="bibr" rid="B6">Ghisalberti and Nepf, 2004</xref>; <xref ref-type="bibr" rid="B37">Tanino and Nepf, 2008</xref>; <xref ref-type="bibr" rid="B32">Sonnenwald et&#xa0;al., 2019</xref>). For sediment simulations, the median grain size of 0.16&#xa0;mm was used, and the settling velocity <italic>&#x3c9;<sub>s</sub>
</italic> was selected as 0.002&#xa0;m/s computed by the equation of <xref ref-type="bibr" rid="B34">Soulsby (1997)</xref>. The critical Shields parameter was taken as 0.03 for flows without cylinders, calculated using the equation in <xref ref-type="bibr" rid="B7">Guo (2002)</xref>. The modified critical Shields parameters for flows with cylinders were calculated using Eqs. 15 and 16, and the values of 0.02, 0.0145, and 0.02&#x2013;0.026 were obtained for the sparse (Case-a with uniform cylinder height), dense (Case-b with uniform cylinder height), and vertically varying cylinder density (Case-c with different cylinder heights) cases.</p>
<sec id="s3_1_1">
<label>3.1.1</label>
<title>Unidirectional flow (Test 1)</title>
<p>The measured mean flow, <italic>TKE</italic>, and suspended sediment concentration under a unidirectional flow (Test 1) were adopted to validate the established model. <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref> shows the comparisons of measured and simulated mean flow velocities affected by cylinders. The data were measured within cylinders by <xref ref-type="bibr" rid="B3">Chen et&#xa0;al. (2020)</xref>. The simulated results generally agreed with the measurements. The results showed better agreements in Cases-a and b with vertically uniform cylinder density (uniform cylinder heights), with both the model skills and <italic>R<sup>2</sup>
</italic> values larger than 0.95. The simulated mean velocity in Case-a was greater than that in Case-b, in accordance with the measurements. The computed mean velocity in Case-c was relatively lower than the measurement, indicating that the cylinder-induced drag force was larger in the simulation than that in the measurements. In Case-c, the average model skills and <italic>R<sup>2</sup>
</italic> values were relatively low with the values smaller than 0.5, while the average <italic>RMSE</italic> value as an additional assessment was 0.06, which indicated acceptable model accuracy. The simulated results of mean velocities proved that the gradient of mean velocity was much smaller in Case-c with vertically varying cylinder density than that in case with uniform density.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Comparisons of measured and simulated mean velocities in Test 1: <bold>(A)</bold> Case-a, <bold>(B)</bold> Case-b, <bold>(C)</bold> Case-c.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g002.tif"/>
</fig>
<p>Under a unidirectional flow, the mean velocity was reduced significantly in the cylinder region. Cylinder density significantly affected the mean velocity. The <italic>U</italic> in Case-b with a higher cylinder density was slightly smaller than that in Case-a; the <italic>U</italic> in Case-c with vertically varying cylinder density was higher near the bottom and lower in the upper layer than that in Cases-a and b with uniform density. The gradient of mean velocity in Case-c was much smaller than that in Cases-a and b. In Case-c, some cylinders had a height of the water depth, thereby resulting in a more uniform vertical distribution of the mean velocity. The distributions of mean velocity in all the cases have shown a turbulent boundary-layer profile, which was due to the small canopy drag compared with the bed drag with the cylinder mainly contributing to the bed roughness.</p>
<p>The simulated <italic>TKE</italic> normalized by friction velocity <italic>u*</italic> are compared with the measured data in <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>, assuming that the lowest grid and measurement points near the bottom fall in the logarithmic velocity profile. Generally, the agreements between simulations and measurements were acceptable in Cases-a and b with uniform cylinder density (both the model skills and <italic>R<sup>2</sup>
</italic> values were larger than 0.90). The model could effectively simulate the <italic>TKE</italic> peak around the cylinder top in the cases with uniform cylinder density. Under the unidirectional condition, the gradually increasing of <italic>TKE</italic> from the cylinder top to the bottom was found in both measurement and simulation. The largest values of <italic>TKE</italic> occurred near the cylinder top due to the shear stress caused by the great velocity gradient. Relatively larger deviations occurred in Case-b near the bottom, where the simulations were greater than the measurements. The results implied that cylinder-induced turbulence production was lower in the measurements than that in the simulations in Case-b. In the laboratory experiments, the equipment (ADV) in the measurement was set up among the cylinders, which cannot catch all of the small stem-scale turbulence in the case with dense cylinders. The <italic>TKE</italic> values in Case-c were much lower than those in Cases-a and b. Affected by the vertically varying cylinder density, the turbulence was relatively higher around the layer interfaces under unidirectional currents, i.e. z =0.2, 0.4 and 0.6&#xa0;m, where cylinder density varies abruptly and strong shear was produced. The simulated vertical profile of <italic>TKE</italic> showed an extremely small gradient in Case-c, consistent with the measurements. However, the simulated <italic>TKE</italic> value was relatively lower than the measured data in Case-c with the model skills and <italic>R<sup>2</sup>
</italic> values lower than 0.3 and <italic>RMSE</italic> of 0.02. In the simulation, the peaks of <italic>TKE</italic> value at z =0.4 and 0.6&#xa0;m were reproduced, while that at z =0.2 m was absent compared with measured data. In Case-c, turbulence was produced by both the multi-shear layers due to the different cylinder heights and the stem-scale turbulence. The <italic>TKE</italic> value affected by the vertically varying cylinder density was complicate especially near the bottom (z=0.2 m), which could not be well resolved in the simulation.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Comparisons of measured and simulated TKEs in Test 1: <bold>(A)</bold> Case-a, <bold>(B)</bold> Case-b, <bold>(C)</bold> Case-c.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g003.tif"/>
</fig>
<p>Under the unidirectional flow (Test 1), the largest values of turbulent kinetic energy (<italic>TKE</italic>) occurred around the cylinder top in cases a and b owing to strong shear turbulence. Affected by the submerged cylinders, the discontinuity in drag appearing at the top of the cylinders generated a region of shear resembling a free shear layer near the top of the canopy, which resulted in higher <italic>TKE</italic> values around the canopy top and relatively lower near the bottom. The maximum value of <italic>TKE</italic> was higher in Case-b than that in Case-a, which was consistent with the results of (<xref ref-type="bibr" rid="B23">Nepf, 1999</xref>), that dense cylinders generate stronger shear and greater array-scale turbulence near the cylinder top, but large-scale shear eddies are constrained in the upper layer. In Case-c with vertically varying cylinder density, the vertical distribution of the <italic>TKE</italic> showed more than one reflection point. The varying cylinder heights in Case-c yielded multiple shear layers throughout the water column (0.2&#x2013;0.6 m), resulting in a higher turbulence near the bottom than that in case with uniform density.</p>
<p>The comparisons of the vertical distributions of measured and simulated sediment concentrations for Cases-a, b and c in Test 1 are shown in <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>. The simulation results from the improved and the traditional sediment transport model were both compared with the measured data. The computed suspended sediment concentrations using the improved sediment transport model matched relatively well with the measurements, especially near the bottom layer (z&lt;0.2m). The model skills were 0.85&#x2013;0.88 and <italic>R<sup>2</sup>
</italic> was 0.8&#x2013;0.88 in the three cases. The near-bed suspended sediment concentration can directly reflect the level of sediment suspension. The simulated sediment concentrations near bottom were in good agreement with the measurements, which proved that the proposed model can reproduce sediment suspension affected by cylinders based on the modified critical Shield number. However, the suspended sediment concentrations in the upper layer of water column were underestimated in Case-c. In Case-c, the simulated values of mean velocity and <italic>TKE</italic> were relatively smaller than the measurements, resulting in smaller bottom shear stress and diffusivity in the simulation. To better simulate the sediment concentration in the upper layer, the vegetation-induced turbulence including shear and stem-scale wake turbulent kinetic energies should be both considered in the future simulation. Sediment concentrations from the traditional sediment transport model were much lower than the measurements, and there was even no sediment suspension in Case-b and Case-c. Due to the higher cylinder densities in Case-b and Case-c, the flow velocities were damped severely, which resulted in small bed shear stress and the underestimation of sediment suspension using the traditional sediment transport model.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Comparisons of measured and simulated sediment concentrations in Test 1: <bold>(A)</bold> Case-a, <bold>(B)</bold> Case-b, <bold>(C)</bold> Case-c.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g004.tif"/>
</fig>
<p>Comparing the sediment concentrations in Case-a with sparse cylinders and Case-b with dense cylinders, it was found that the sediment concentrations increased with the growing cylinder density, especially near the bottom, which agreed well with the results in <xref ref-type="bibr" rid="B39">Tinoco and Coco, 2018</xref>; <xref ref-type="bibr" rid="B38">Tinoco and Coco, 2016</xref>). It was also proved that the bottom sediment suspension was mainly contributed by the turbulence generated around cylinder stems. In other words, the incipient motion velocity decreased with the increasing cylinder density in flows with cylinders. The sediment concentrations were higher in Case-b and Case-c than that in Case-a under a unidirectional flow (Test 1). Although the solid volume fractions were the same near the bottom in these cases, the vertically varying cylinder density promoted sediment suspension due to higher turbulence near the bottom than that in case with uniform density.</p>
</sec>
<sec id="s3_1_2">
<label>3.1.2</label>
<title>Combined wave-current flows (Tests 2-4)</title>
<p>To further validate the sediment model, experimental data from Tests 2, 3, and 4 were used to simulate sediment suspension in combined wave-current flows. <xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5</bold>
</xref>&#x2013;<xref ref-type="fig" rid="f7">
<bold>7</bold>
</xref> show that the computed suspended sediment concentrations using the improved sediment transport model agreed with the measurements near the bottom, while some deviations appeared in the upper layer, where the simulated sediment concentrations were smaller than the measured data. The model skills were 0.83&#x2013;0.85, 0.66&#x2013;0.88, and 0.54&#x2013;0.73 in Tests 2, 3, and 4, respectively (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). The <italic>R<sup>2</sup>
</italic> values were 0.84&#x2013;0.88, 0.84&#x2013;0.88, and 0.84&#x2013;0.89 in Tests 2, 3, and 4, respectively (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). The sediment concentration increased gradually with the growing wave height. Higher wave height induced larger bed shear stress, which enhanced bottom sediment flux and sediment suspension. More suspended sediment was kept stay in the upper layer by the relatively larger turbulent kinetic energy due to higher wave height. Using the modified Shields number given in Eqs. 15 and 16, the model could predict sediment suspension near bottom under combined wave&#x2013;current flows with a certain accuracy. Due to the similar reason as in Test 1, the suspended sediment concentrations in the upper layer of water column were also underestimated in the combined wave-current flows, especially in Case-c. Affected by cylinders with vertically varying density, the sediment simulation was more difficult due to the complicate turbulence in both canopy- and stem- scale. On the contrary, the sediment suspension was significantly underestimated using the traditional sediment transport model. It is necessary to consider the cylinder effects on the threshold condition of sediment suspension in the flows with cylinders.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Measured and simulated sediment concentrations in Test 2: <bold>(A)</bold> Case-a, <bold>(B)</bold> Case-b, <bold>(C)</bold> Case-c.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g005.tif"/>
</fig>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Measured and simulated sediment concentrations in Test 3: <bold>(A)</bold> Case-a, <bold>(B)</bold> Case-b, <bold>(C)</bold> Case-c.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g006.tif"/>
</fig>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Measured and simulated sediment concentrations in Test 4: <bold>(A)</bold> Case-a, <bold>(B)</bold> Case-b.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g007.tif"/>
</fig>
<table-wrap-group id="T4" position="float">
<label>Table&#xa0;4</label>
<caption>
<p>Accuracy of model calibrations and validations.</p>
</caption>
<table-wrap>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="2" align="center">Calibrations</th>
<th valign="middle" colspan="2" align="center">With uniform cylinder density</th>
<th valign="top" colspan="3" align="center">With vertically varying cylinder density</th>
</tr>
<tr>
<th valign="middle" align="center">
<italic>R<sup>2</sup>
</italic>
</th>
<th valign="top" align="center">
<italic>M<sub>S</sub>
</italic>
</th>
<th valign="middle" align="center">
<italic>R<sup>2</sup>
</italic>
</th>
<th valign="top" align="center">
<italic>M<sub>S</sub>
</italic>
</th>
<th valign="top" align="center">
<italic>RMSE</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">Test 1</td>
<td valign="middle" align="center">0.82-0.97</td>
<td valign="top" align="center">0.88-0.97</td>
<td valign="top" align="center">0.30-0.80</td>
<td valign="middle" align="center">0.20-0.88</td>
<td valign="top" align="center">0.02-0.08</td>
</tr>
<tr>
<td valign="middle" align="center">Test 2</td>
<td valign="middle" align="center">0.84-0.87</td>
<td valign="top" align="center">0.83-0.96</td>
<td valign="top" align="center">0.88</td>
<td valign="middle" align="center">0.88</td>
<td valign="top" align="center">0.15</td>
</tr>
<tr>
<td valign="middle" align="center">Test 3</td>
<td valign="middle" align="center">0.84-0.89</td>
<td valign="top" align="center">0.65-0.92</td>
<td valign="top" align="center">0.88</td>
<td valign="middle" align="center">0.87</td>
<td valign="top" align="center">0.17</td>
</tr>
<tr>
<td valign="middle" align="center">Test 4</td>
<td valign="middle" align="center">0.84-0.89</td>
<td valign="top" align="center">0.54-0.73</td>
<td valign="top" align="center">&#x2013;</td>
<td valign="middle" align="center">&#x2013;</td>
<td valign="top" align="center">&#x2013;</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="2" align="center">Validations</th>
<th valign="middle" colspan="5" align="center">With uniform cylinder density</th>
</tr>
<tr>
<th valign="middle" colspan="2" align="center">
<italic>R<sup>2</sup>
</italic>
</th>
<th valign="top" colspan="3" align="center">
<italic>M<sub>S</sub>
</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">C15-1</td>
<td valign="middle" colspan="2" align="center">0.93-0.97</td>
<td valign="top" colspan="3" align="center">0.93-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C15-2</td>
<td valign="middle" colspan="2" align="center">0.92-0.98</td>
<td valign="top" colspan="3" align="center">0.94-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C18-1</td>
<td valign="middle" colspan="2" align="center">0.93-0.98</td>
<td valign="top" colspan="3" align="center">0.97-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C18-2</td>
<td valign="middle" colspan="2" align="center">0.95-0.99</td>
<td valign="top" colspan="3" align="center">0.95-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C24-1</td>
<td valign="middle" colspan="2" align="center">0.97-0.98</td>
<td valign="top" colspan="3" align="center">0.95-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C24-2</td>
<td valign="middle" colspan="2" align="center">0.97-0.99</td>
<td valign="top" colspan="3" align="center">0.88-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C30-1</td>
<td valign="middle" colspan="2" align="center">0.96-0.98</td>
<td valign="top" colspan="3" align="center">0.91-0.98</td>
</tr>
<tr>
<td valign="middle" align="center">C30-2</td>
<td valign="middle" colspan="2" align="center">0.96-0.99</td>
<td valign="top" colspan="3" align="center">0.68-0.98</td>
</tr>
</tbody>
</table>
</table-wrap>
</table-wrap-group>
<p>Comparing the results in Cases-a and b, it was clear that the suspended sediment concentration increased with cylinder density, indicative of the effects of cylinders on suspended sediment transport in stem scale. Although the stem numbers were the same in Case-b and Case-c, the sediment distributions were slightly different. In Case-c, which was affected by vertically varying cylinder density, the array-scale turbulence caused by both the multi-shear layers and stem-scale turbulence affected sediment suspension. The imposed waves increased sediment suspension compared to the results under different hydrodynamic conditions. The co-occurring waves greatly enhanced bottom shear stress and promoted sediment suspension, and the sediment concentration was found to be increased with the growing wave height, especially near the bottom. In the upper layer (<italic>z</italic>&gt;0.2m), there were clear deviations between the measured and simulated data. The sediment concentration near the bottom was mainly controlled by the bottom sediment flux, while it was determined by turbulent diffusivity in the upper layer.</p>
</sec>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Validations</title>
<p>The proposed model was further verified using another experimental dataset by <xref ref-type="bibr" rid="B17">Lu (2008)</xref>, which has also been applied to validate the model in <xref ref-type="bibr" rid="B12">Li et&#xa0;al. (2020)</xref>. The material of sediment particles used in the experiments (<xref ref-type="bibr" rid="B17">Lu, 2008</xref>) was plastic with the median grain size of 0.217&#xa0;mm and sediment density of 1.082 g/cm<sup>3</sup>. The <italic>C<sub>D</sub>
</italic> values were 0.6&#x2013;0.7 in the simulations according to the empirical equations (<xref ref-type="bibr" rid="B16">Lou et&#xa0;al., 2021</xref>). The settling velocity <italic>&#x3c9;<sub>s</sub>
</italic> was 0.0047&#xa0;m/s based on the experimental measurements (<xref ref-type="bibr" rid="B17">Lu, 2008</xref>; <xref ref-type="bibr" rid="B12">Li et&#xa0;al., 2020</xref>). The critical Shields parameter was calculated as 0.035 for flows without cylinder using the equation in <xref ref-type="bibr" rid="B7">Guo (2002)</xref>. The modified critical Shields parameters for flows with cylinders were in the range of 0.02-0.03 calculated using Eq. (15).</p>
<p>To validate the proposed model, model results were compared with experimental datasets in 8 tests (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). <xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref> shows the comparisons of measured and simulated mean velocities. There was a good correlation between measured and predicted values with both the model skills and <italic>R<sup>2</sup>
</italic> values larger than 0.93 (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). Deviations between model results and experimental data were found at <italic>z/h</italic>&lt;0.25 in tests C24-2 and at <italic>z/h</italic>&lt;0.2 in C30-2. Larger values of flow velocity in the simulations were caused by the underestimation of <italic>C<sub>D</sub>
</italic>.</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Comparisons of measured and simulated mean velocities in tests carried out by <xref ref-type="bibr" rid="B17">Lu (2008)</xref>: <bold>(A)</bold> C15-1, <bold>(B)</bold> C15-2, <bold>(C)</bold> C18-1, <bold>(D)</bold> C18-2, <bold>(E)</bold> C24-1, <bold>(F)</bold> C24-2, <bold>(G)</bold> C30-1, <bold>(H)</bold> C30-2.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g008.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9</bold>
</xref> presents the comparisons of measured and simulated sediment concentrations normalized by the concentration at the cylinder height (<italic>C<sub>h</sub>
</italic>). The simulated results matched with the measured data in tests C15-1, C15-2, C18-1 and C18-2 with the averaged model skill and <italic>R<sup>2</sup>
</italic> value of 0.96 and 0.94. Relatively larger deviations occurred in tests C24-1 (model skill of 0.95) and C30-1 (model skill of 0.91), while obvious disagreements were found in C24-2 (model skill of 0.88) and C30-2 (model skill of 0.68). The simulated concentration was greater than the measured data especially near the bottom, which was greatly affect by the sediment suspension processes. In tests C24-1, C24-2, C30-1 and C30-2, the cylinder density was relatively low. The cylinder densities (<italic>&#x3d5;</italic>) in tests C24-2 and C30-2 were both less than 0.01 (0.006 and 0.007, respectively). It was found that the proposed sediment transport model was less accurate in cases with small cylinder density (<italic>&#x3d5;&lt;</italic>0.01).</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>Comparisons of measured and simulated sediment concentrations normalized by the concentrations at cylinder top (C<sub>h</sub>) in tests carried out by <xref ref-type="bibr" rid="B17">Lu (2008)</xref>: <bold>(A)</bold> C15-1, <bold>(B)</bold> C15-2, <bold>(C)</bold> C18-1, <bold>(D)</bold> C18-2, <bold>(E)</bold> C24-1, <bold>(F)</bold> C24-2, <bold>(G)</bold> C30-1, <bold>(H)</bold> C30-2.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1105112-g009.tif"/>
</fig>
</sec>
</sec>
<sec id="s4">
<label>4</label>
<title>Rationale and limitations of the proposed model</title>
<p>The traditional sediment transport model for bare bed (Eq. 14) without the consideration of cylinders was modified to reproduce the processes of sediment suspension affected by cylinders. Sediment flux near bed was calculated based on a modified critical Shields number due to the presence of cylinders in the improved model (Eqs. 15 and 16). This improvement allowed to capture the strong influence of cylinders on sediment suspension due to the cylinder induced turbulence. This proposed model can be simply used to simulate sediment suspension affected by cylinders in both unidirectional and the combined wave-current flows under similar conditions as the Eqs. 15 and 16 deduced. However, it was found from the model verifications that better agreements between simulations and measurements were observed in cases with uniform density cylinders. Sediment suspension affected by vertically varying density cylinders should be further studied considering the combined influences of shear and stem-scale wake turbulent kinetic energies.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusions</title>
<p>Numerical simulations on sediment suspension affected by cylinders under unidirectional and combined wave&#x2013;current flows were presented. A sediment transport model with a modified critical Shields number considering cylinder density (solid volume fraction <italic>&#x3d5;</italic>) proposed by <xref ref-type="bibr" rid="B38">Tinoco and Coco (2016)</xref> was employed to simulate sediment suspension in the flows with cylinders. Results demonstrated that the improved model could reasonably reproduce sediment suspension in the flows with cylinders, while the sediment transport model with the traditional critical Shields number failed to capture sediment suspension in the cylinder canopy. Numerical simulations showed that denser cylinders and larger wave heights resulted in stronger sediment suspensions. The stem-scale turbulence induced by cylinders and enhanced bed shear stress by waves greatly promoted sediment suspension. Flow&#x2013;cylinder&#x2013;sediment interactions are complex and require more investigations, and the proposed sediment transport model should be further validated using field observations.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>SL: Conceptualization, Methodology, Writing &#x2013; original draft, Funding acquisition. XC: Conceptualization, Methodology, Data analysis. SZ: Conceptualization, Methodology, Editing and analysis. GM: Writing &#x2013; review and editing. SLi: Conceptualization, Project administration. LR: Formula derivation and theoretical analysis. EN: Formula derivation and theoretical analysis, Formal analysis. IF: Conceptualization, Writing &#x2013; review and editing. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>This work was sponsored by the National Natural Science Foundation of China (42072281), Shanghai Science and Technology Innovation Action Plan (20230742500, 22ZR1464200, 22230712900), the Fundamental Research Funds for the Central Universities (22120210576), and Top Discipline Plan of Shanghai Universities-Class I (2022-3-YB-03).</p>
</sec>
<sec id="s9" 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="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Augustin</surname> <given-names>L. N.</given-names>
</name>
<name>
<surname>Irish</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Lynett</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Laboratory and numerical studies of wave damping by emergent and near-emergent wetland vegetation</article-title>. <source>Coast. Eng.</source> <volume>56</volume> (<issue>3</issue>), <fpage>332</fpage>&#x2013;<lpage>340</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2008.09.004</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Balke</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Klaassen</surname> <given-names>P. C.</given-names>
</name>
<name>
<surname>Garbutt</surname> <given-names>A.</given-names>
</name>
<name>
<surname>van der Wal</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Herman</surname> <given-names>P. M. J.</given-names>
</name>
<name>
<surname>Bouma</surname> <given-names>T. J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Conditional outcome of ecosystem engineering: A case study on tussocks of the salt marsh pioneer spartina anglica</article-title>. <source>Geomorphology</source> <volume>153</volume>, <fpage>232</fpage>&#x2013;<lpage>238</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.geomorph.2012.03.002</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lou</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S. G.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>G. F.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>H. Z.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>G. H.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Velocity and turbulence affected by submerged rigid vegetation under waves, currents and combined wave-current flows</article-title>. <source>Coast. Eng.</source> <volume>159</volume>, <fpage>103727</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2020.103727</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Pang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A time-efficient CFD approach for hydrodynamic coefficient determination and model simplification of submarine</article-title>. <source>Ocean Eng.</source> <volume>154</volume>, <fpage>16</fpage>&#x2013;<lpage>26</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.oceaneng.2018.02.003</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghisalberti</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Mixing layers and coherent structures in vegetated aquatic flows</article-title>. <source>J. Geophysical Research-Oceans</source> <volume>107</volume> (<issue>C2</issue>), <page-range>3&#x2013;1&#x2013;3&#x2013;11</page-range>. doi: <pub-id pub-id-type="doi">10.1029/2001JC000871</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghisalberti</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>The limited growth of vegetated shear layers</article-title>. <source>Water Resour. Res.</source> <volume>40</volume>, (<issue>7</issue>), <fpage>3-1-3-11</fpage>. doi: <pub-id pub-id-type="doi">10.1029/2003WR002776</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname> <given-names>J. K.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Hunter rouse and shields diagram</article-title>. <source>Adv. Hydraulic Water Eng.</source> <volume>1-2</volume>, <fpage>1096</fpage>&#x2013;<lpage>1098</lpage>. doi: <pub-id pub-id-type="doi">10.1142/9789812776969_0200</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Weng</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Numerical simulation of dam-break flow and bed change considering the vegetation effects</article-title>. <source>Int. J. Sediment Res.</source> <volume>32</volume> (<issue>1</issue>), <fpage>105</fpage>&#x2013;<lpage>120</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijsrc.2015.04.004</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Horppila</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Nurminen</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Effects of submerged macrophytes on sediment resuspension and internal phosphorus loading in lake hiidenvesi (southern Finland)</article-title>. <source>Water Res.</source> <volume>37</volume> (<issue>18</issue>), <fpage>4468</fpage>&#x2013;<lpage>4474</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0043-1354(03)00405-6</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Hartig</surname> <given-names>E. K.</given-names>
</name>
<name>
<surname>Orton</surname> <given-names>P. M.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Numerical modeling of salt marsh morphological change induced by hurricane sandy</article-title>. <source>Coast. Eng.</source> <volume>132</volume>, <fpage>63</fpage>&#x2013;<lpage>81</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2017.11.001</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>Y. M.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Flow structure and sediment transport with impacts of aquatic vegetation</article-title>. <source>J. Hydrodynamics</source> <volume>20</volume> (<issue>4</issue>), <fpage>461</fpage>&#x2013;<lpage>468</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S1001-6058(08)60081-5</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Estimating the distribution of suspended sediment concentration in submerged vegetation flow based on gravitational theory</article-title>. <source>J. Hydrology</source> <volume>587</volume>, <fpage>124921</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jhydrol.2020.124921</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lopez</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Garcia</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Open-channel flow through simulated vegetation: Suspended sediment transport modeling</article-title>. <source>Water Resour. Res.</source> <volume>34</volume> (<issue>9</issue>), <fpage>2341</fpage>&#x2013;<lpage>2352</lpage>. doi: <pub-id pub-id-type="doi">10.1029/98WR01922</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lou</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Sediment suspension affected by submerged rigid vegetation under waves, currents and combined wave&#x2013;current flows</article-title>. <source>Coast. Eng.</source> <volume>173</volume>, <fpage>104082</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2022.104082</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lou</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Laboratory study of the effect of vertically varying vegetation density on waves, currents and wave-current interactions</article-title>. <source>Appl. Ocean Res.</source> <volume>79</volume>, <fpage>74</fpage>&#x2013;<lpage>87</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.apor.2018.07.012</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lou</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Modelling of stem-scale turbulence and sediment suspension in vegetated flow</article-title>. <source>J. Hydraulic Res</source>. <volume>59</volume>, <fpage>355</fpage>&#x2013;<lpage>377</lpage> doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00221686.2020.1780491</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Lu</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2008</year>). <source>Experimental study on distribution law of suspended sediment in water flow of rigid plants (in Chinese)</source> (<publisher-loc>China</publisher-loc>: <publisher-name>Hohai University</publisher-name>). doi:&#xa0;<pub-id pub-id-type="doi">10.7666/d.y1268411</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Niroomandi</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Lou</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Numerical study of sediment transport on a tidal flat with a patch of vegetation</article-title>. <source>Ocean Dynamics</source> <volume>65</volume> (<issue>2</issue>), <fpage>203</fpage>&#x2013;<lpage>222</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s10236-014-0804-8</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Kirby</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Su</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Figlus</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Numerical study of turbulence and wave damping induced by vegetation canopies</article-title>. <source>Coast. Eng.</source> <volume>80</volume>, <fpage>68</fpage>&#x2013;<lpage>78</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2013.05.007</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Kirby</surname> <given-names>J. T.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Shock-capturing non-hydrostatic model for fully dispersive surface wave processes</article-title>. <source>Ocean Model.</source> <volume>43-44</volume>, <fpage>22</fpage>&#x2013;<lpage>35</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ocemod.2011.12.002</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mendez</surname> <given-names>F. J.</given-names>
</name>
<name>
<surname>Losada</surname> <given-names>I. J.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>An empirical model to estimate the propagation of random breaking and nonbreaking waves over vegetation fields</article-title>. <source>Coast. Eng.</source> <volume>51</volume> (<issue>2</issue>), <fpage>103</fpage>&#x2013;<lpage>118</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2003.11.003</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moller</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Spencer</surname> <given-names>T.</given-names>
</name>
<name>
<surname>French</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Leggett</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Dixon</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Wave transformation over salt marshes: A field and numerical modelling study from north Norfolk, England</article-title>. <source>Estuar. Coast. Shelf Sci.</source> <volume>49</volume> (<issue>3</issue>), <fpage>411</fpage>&#x2013;<lpage>426</lpage>. doi: <pub-id pub-id-type="doi">10.1006/ecss.1999.0509</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Drag, turbulence, and diffusion in flow through emergent vegetation</article-title>. <source>Water Resour. Res.</source> <volume>35</volume> (<issue>2</issue>), <fpage>479</fpage>&#x2013;<lpage>489</lpage>. doi: <pub-id pub-id-type="doi">10.1029/1998WR900069</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Flow and transport in regions with aquatic vegetation</article-title>. <source>Annu. Rev. Fluid Mechanics</source> <volume>44</volume>, <fpage>123</fpage>&#x2013;<lpage>142</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev-fluid-120710-101048</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
<name>
<surname>Vivoni</surname> <given-names>E. R.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Flow structure in depth-limited, vegetated flow</article-title>. <source>J. Geophysical Research-Oceans</source> <volume>105</volume> (<issue>C12</issue>), <fpage>28547</fpage>&#x2013;<lpage>28557</lpage>. doi: <pub-id pub-id-type="doi">10.1029/2000JC900145</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Neumeier</surname> <given-names>U.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Velocity and turbulence variations at the edge of saltmarshes</article-title>. <source>Continental Shelf Res.</source> <volume>27</volume> (<issue>8</issue>), <fpage>1046</fpage>&#x2013;<lpage>1059</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.csr.2005.07.009</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Ni</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2014</year>). <source>Laboratory investigations on the drag coefficient of cylinder subjected to wave motion</source> (<publisher-loc>The Netherlands</publisher-loc>: <publisher-name>Master thesis from Delft University of Technology</publisher-name>). 10.2001.</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodi</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>1987</year>). <article-title>Examples of calculation methods for flow and mixing in stratified fluids</article-title>. <source>J. Geophysical Research-Oceans</source> <volume>92</volume> (<issue>C5</issue>), <fpage>5305</fpage>&#x2013;<lpage>5328</lpage>. doi: <pub-id pub-id-type="doi">10.1029/JC092iC05p05305</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ros</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Colomer</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Serra</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Pujol</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Soler</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Casamitjana</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Experimental observations on sediment resuspension within submerged model canopies under oscillatory flow</article-title>. <source>Continental Shelf Res.</source> <volume>91</volume>, <fpage>220</fpage>&#x2013;<lpage>231</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.csr.2014.10.004</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sheng</surname> <given-names>Y. P.</given-names>
</name>
<name>
<surname>Lapetina</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The reduction of storm surge by vegetation canopies: Three-dimensional simulations</article-title>. <source>Geophysical Res. Lett.</source> <volume>39</volume>, <fpage>L20601</fpage>. doi: <pub-id pub-id-type="doi">10.1029/2012GL053577</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Snyder</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Hsu</surname> <given-names>T. J.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A numerical investigation of convective sedimentation</article-title>. <source>J. Geophysical Research-Oceans</source> <volume>116</volume>, <fpage>C09024</fpage>. doi: <pub-id pub-id-type="doi">10.1029/2010JC006792</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sonnenwald</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Stovin</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Guymer</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Estimating drag coefficient for arrays of rigid cylinders representing emergent vegetation</article-title>. <source>J. Hydraulic Res.</source> <volume>57</volume> (<issue>4</issue>), <fpage>591</fpage>&#x2013;<lpage>597</lpage>. doi: <pub-id pub-id-type="doi">10.1080/00221686.2018.1494050</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Soulsby</surname> <given-names>R. L.</given-names>
</name>
</person-group> (<year>1995</year>). <source>Bed shear-stresses due to combined waves and currents, section 4.5 in advances in coastal, morphodynamics</source>. Ed. <person-group person-group-type="editor">
<name>
<surname>Stive</surname> <given-names>M. J. F.</given-names>
</name>
</person-group> (<publisher-loc>Delft</publisher-loc>: <publisher-name>Delft Hydraulics</publisher-name>).</citation>
</ref>
<ref id="B34">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Soulsby</surname> <given-names>R. L.</given-names>
</name>
</person-group> (<year>1997</year>). <source>Dynamics of marine sands</source> (<publisher-loc>Thomas Telford, London</publisher-loc>: <publisher-name>A manual for practical applications</publisher-name>).</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Lei</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Impact of vegetation-generated turbulence on the critical, near-bed, wave-velocity for sediment resuspension</article-title>. <source>Water Resour. Res.</source> <volume>55</volume> (<issue>7</issue>), <fpage>5904</fpage>&#x2013;<lpage>5917</lpage>. doi: <pub-id pub-id-type="doi">10.1029/2018WR024335</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Liang</surname> <given-names>D. F.</given-names>
</name>
<name>
<surname>Lv</surname> <given-names>S. Q.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Incipient motion of sediment in the presence of emergent rigid vegetation</article-title>. <source>J. Hydro-environment Res.</source> <volume>7</volume>, <fpage>202</fpage>&#x2013;<lpage>208</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jher.2012.11.002</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tanino</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Laboratory investigation of mean drag in a random array of rigid, emergent cylinders</article-title>. <source>J. Hydraulic Eng.</source> <volume>134</volume> (<issue>1</issue>), <fpage>34</fpage>&#x2013;<lpage>41</lpage>. doi: <pub-id pub-id-type="doi">10.1061/(ASCE)0733-9429(2008)134:1(34)</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tinoco</surname> <given-names>R. O.</given-names>
</name>
<name>
<surname>Coco</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A laboratory study on sediment resuspension within arrays of rigid cylinders</article-title>. <source>Adv. Water Resour.</source> <volume>92</volume>, <fpage>1</fpage>&#x2013;<lpage>9</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.advwatres.2016.04.003</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tinoco</surname> <given-names>R. O.</given-names>
</name>
<name>
<surname>Coco</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Turbulence as the main driver of resuspension in oscillatory flow through vegetation</article-title>. <source>J. Geophysical Research-Earth Surface</source> <volume>123</volume> (<issue>5</issue>), <fpage>891</fpage>&#x2013;<lpage>904</lpage>. doi: <pub-id pub-id-type="doi">10.1002/2017JF004504</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tinoco</surname> <given-names>R. O.</given-names>
</name>
<name>
<surname>Cowen</surname> <given-names>E. A.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The direct and indirect measurement of boundary stress and drag on individual and complex arrays of elements</article-title>. <source>Experiments fluids</source> <volume>54</volume> (<issue>4</issue>), <fpage>1509</fpage>&#x2013;<lpage>1522</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00348-013-1509-3</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Van Hoften</surname> <given-names>J. D. A.</given-names>
</name>
<name>
<surname>Karaki</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1976</year>). <article-title>Interaction of waves and a turbulent current</article-title>. <source>Proc. Int. Conf. Coast. Eng. Proceedings of the 15th International Conference on Coastal Engineering ASCE New York</source>, <fpage>404</fpage>&#x2013;<lpage>422</lpage>. doi: <pub-id pub-id-type="doi">10.9753/icce.v15.22</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Van Rijn</surname> <given-names>L. C.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Sediment transport, part II: Suspended load transport</article-title>. <source>J. Hydraulic Engineering-ASCE</source> <volume>110</volume> (<issue>11</issue>), <fpage>1613</fpage>&#x2013;<lpage>1641</lpage>. doi: <pub-id pub-id-type="doi">10.1061/(ASCE)0733-9429(1984)110:11(1613)</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Van Rijn</surname> <given-names>L. C.</given-names>
</name>
</person-group> (<year>1993</year>). <source>Principles of sediment transport in rivers, estuaries, and coastal seas</source> (<publisher-loc>The Netherlands</publisher-loc>: <publisher-name>Aqua Publications</publisher-name>).</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wilmott</surname> <given-names>C. J.</given-names>
</name>
</person-group> (<year>1981</year>). <article-title>On the validation of models</article-title>. <source>Phys. Geogr.</source> <volume>2</volume> (<issue>2</issue>), <fpage>184</fpage>&#x2013;<lpage>194</lpage>. doi: <pub-id pub-id-type="doi">10.1080/02723646.1981.10642213</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Cox</surname> <given-names>D. T.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Modeling wave attenuation induced by the vertical density variations of vegetation</article-title>. <source>Coast. Eng.</source> <volume>112</volume>, <fpage>17</fpage>&#x2013;<lpage>27</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.coastaleng.2016.02.004</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>S. S. Y.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>One-dimensional modeling of dam-break flow over movable beds</article-title>. <source>J. Hydraulic Eng.</source> <volume>133</volume> (<issue>1</issue>), <fpage>48</fpage>&#x2013;<lpage>58</lpage>. doi: <pub-id pub-id-type="doi">10.1061/(ASCE)0733-9429(2007)133:1(48)</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>J. Q.</given-names>
</name>
<name>
<surname>Chung</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The onset of sediment transport in vegetated channels predicted by turbulent kinetic energy</article-title>. <source>Geophys. Res. Lett.</source> <volume>43</volume>, <page-range>11,261&#x2013;11,268</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/2016GL071092</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>J. Q.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H. M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Impact of vegetation on bed load transport rate and bedform characteristics</article-title>. <source>Water Resour. Res.</source> <volume>55</volume> (<issue>7</issue>), <fpage>6109</fpage>&#x2013;<lpage>6124</lpage>. doi: <pub-id pub-id-type="doi">10.1029/2018WR024404</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Nepf</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Wave-driven sediment resuspension within a model eelgrass meadow</article-title>. <source>J. Geophysical Research: Earth Surface</source> <volume>124</volume>, <page-range>1035&#x2013;1053</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2018JF004984</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Nurminen</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>The influence of macrophytes on sediment resuspension and the effect of associated nutrients in a shallow and large lake</article-title>. (<publisher-name>Lake Taihu</publisher-name>, <publisher-loc>China</publisher-loc>) <source>PloS One</source> <volume>10</volume> (<issue>6</issue>), <fpage>e0127915</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0127915</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>