Bedload and Suspended Load Transport in the 140km Reach Downstream of the Mississippi River Avulsion to the Atchafalaya River

The Mississippi River Delta has been continuously losing land since the 1930s due to several factors, chief of which is the reduced sediment supply. A few recent studies have estimated individual components of short-term sediment transport, i.e., bedload and suspended load, separately for some locations along the Lowermost Mississippi River (LmMR, commonly considered as the last 500-km reach of the Mississippi River before entering the Gulf of Mexico). However, the combined effects of both components on the long-term sediment supply along the river reach are still unclear. One of the major obstacles here hindering our understanding is that it is difficult and impractical to accurately measure bedloads in large alluvial rivers, such as the Mississippi. In this study, we estimated bedloads of three medium grain sizes (D50 = 0.125, 0.25 and 0.5 mm) for three locations along the uppermost 140-km reach of the LmMR: Tarbert Landing (TBL) at river kilometer (rk) 493, St Francisville (St F) at rk 419, and Baton Rouge (BTR) at rk 367.5 during 2004–2015. We also estimated suspended loads at St F during 1978–2015 and at BTR during 2004–2015 in order to discern the interactive relationship between bedloads and suspended loads. We found gradually increasing bedloads for all the three medium grain sizes from TBL (83, 41.5 and 20.75 million tons (MT), respectively) to BTR (96, 48 and 24 MT). We also found that suspended loads at TBL (reported previously) were significantly higher than those at St F and BTR during corresponding overlapping periods. Bedloads increased almost linearly with suspended loads, river discharge and river stage at the upstream locations (TBL and St F); however, such a relationship was not evident downstream at BTR. The total sediment load (bedload + suspended sediment load) was substantially higher at TBL (931 MT), while lower and nearly equal at other two downstream locations (550 MT at St F and 544 MT at BTR) during 2004–2010 (the matching period of availability of both loads). These findings indicate that the uppermost 20–25 km LmMR reach (covering TBL) has potentially entrapped substantial suspended load over the last three to four decades, while bedload transport prevails in the lower reach (covering St F and BTR). We suggest that future sediment management in the river should seek engineering solutions for moving trapped coarse sediments downstream towards the coast for the Mississippi River Delta restoration.


Introduction
Studies have found that sediment supply from many rivers in the world to the oceans has decreased substantially in the past several decades [1][2][3].The reduction estimates are mostly reported in suspended sediment load; for instance, Meade and Parker [4] reported 3.5 times decline in suspended sediment load of the Colorado River from the 1930s to the 2000s; Yang et al. [5] reported 100 megatons (MT)/year decline in sediment load of the Yangtze River from the 1950s to the 1990s; and Meade and Moody [6] noted a 3.5 times decrease in sediment load of the Mississippi River from the early to the late 20th century.These studies have improved our understanding of the trend of riverine suspended sediment loads to the world's coasts; however, our knowledge of bed sediment transport in these rivers is very limited.One of the major obstacles hindering our understanding is that it is difficult and impractical to accurately measure bed sediment transport, especially in large alluvial rivers.However, such knowledge is important and can greatly aid in assessing riverine sediment transport, channel morphodynamics and long-term future of the river deltas.
Bed sediment loads can vary largely from 1% to 33% of the suspended loads in alluvial rivers, depending on sediment particle size [7][8][9].Direct measurements of bed sediment transport are difficult and scarce in most alluvial rivers globally [7,10].Bed sediments are generally intermixed with sand particles (coarser sediment particles >0.063 mm in diameter) of different sizes [7].Coarser sediments varying from 0.063 to 1.2 mm in sizes are difficult to be distinguished as true bedloads or suspended loads because sometimes they move by traction (rolling, sliding) and/or saltation and sometime they move by suspension [7,11].Leopold and Emmett [12] emphasized the problem of intractability between bedloads and suspended loads because they did not find any sampling device that provided reliable bedload measurements moving along the riverbed.Therefore, bedloads are generally reported in terms of their transport rates which can be defined as the mass of bed sediment moving across a wetted riverbed for a given discharge regime in unit time [13].Bedload transport rates, suspended sediment loads, and suspended sand loads have mostly been studied separately.Their combined effect on the total sediment supply in rivers, however, still needs to be investigated.Such an analysis can provide relevant information for river engineering and sediment management.
The total sediment supply constitutes of both bedload and suspended load and contributes substantially to spatiotemporal changes in channel adjustment and morphology in alluvial rivers [14][15][16].However, many previous studies on sediment supply have either ignored bedload [17][18][19] or applied it as a fixed fraction of the sediment load [20][21][22].In all latter studies, bedload was given a fixed percentage based on tables proposed by Maddock and Borland [23] and further expanded by Lane and Borland [24].Both Maddock and Borland [23] and Lane and Borland [24] put the tables as best estimates of bedload percentage in different river systems.However, their estimates could possibly be subjected to considerable errors as they did not measure and/or quantify either the suspended loads or bedloads in any of the river types.In addition to these limitations in suspended load and bedload partitioning, interactive relationship between the two components also remains poorly understood.Quantifying bedload transport and investigating the relationship between bedload and suspended load can be especially beneficial for sediment management along regulated alluvial rivers with reduced sediment supply.
The final 500 km stretch of the Mississippi River from its diversion to the Atchafalaya River at the Old River Control Structure to the Gulf of Mexico in the southern USA is a renowned example of alluvial rivers with reduced sediment supply over time.The reduced sediment supply has contributed to land loss in the Mississippi River Delta Plain (MRDP) over the last several decades [25][26][27].Meade and Moody [6] concluded that mean annual sediment load for the Lowermost Mississippi River (LmMR) at its uppermost location, Tarbert Landing (TBL) [at 493 river kilometers (rk) above the river's Gulf outlet] reduced by about 3.5 times [from ~400 million tons (MT)/year to 115 MT/year] from pre-to post-human interference period (early to late 20th century).Rosen and Xu [28] estimated annual sediment loads at Tarbert Landing for 1980-2010 and noticed a 12 MT/year decline in the loads from 1993.Horowitz [29] estimated annual sediment loads at Tarbert Landing and the next downstream location of the LmMR, St Francisville (St F) at rk 419 (about 71 rk downstream of Tarbert Landing) for 1980-2007.The study reported that the mean annual sediment load at Tarbert Landing and at St Francisville declined by 21 and 41 MT, respectively after 1993.Long-term trends for suspended sand load have also been reported previously [30,31], but only at Tarbert Landing.The critical role of suspended sand loads in building new lands for the MRDP has been previously recognized [32][33][34].
Water 2017, 9, 716 3 of 28 Thus, quantification of sand loads at locations downstream from Tarbert Landing can help in clarifying their potential sand storage and delta building capacity.
Recently, Wang and Xu [35] assessed channel bars and Joshi and Xu [36] assessed river stages and cross-sections along the uppermost 10 and 140 km reaches of the LmMR below its diversion to the Atchafalaya River (near the Old River Control Structure, or ORCS), respectively.Both studies found that the reaches had potentially trapped substantial amount of sediments over the last three decades.They argued about a need for quantifying long term bed and suspended loads at the locations in which riverbed adjustment was found, i.e., Tarbert Landing, Bayou Sara (BS) (at rk 427) (about 8 rk upstream of St Francisville) and Baton Rouge (BTR) (at rk 367.5).The United States Geological Survey (USGS) data for suspended sediment concentrations upon which sediment and sand load calculations are based can be found for at least four decades  at Tarbert Landing and one decade (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) at Baton Rouge, but no such data is available for Bayou Sara.Nearest long term data  which can analogue sediment loads at Bayou Sara are available at St Francisville, about 8 rk upstream from Bayou Sara.
Majority of the long term suspended sediment and sand load budgets for the LmMR have only been carried out at Tarbert Landing.There are limited and discontinuous studies of flow and sediment transport between Tarbert Landing and Baton Rouge.Allison and Meselhe [37] and Horowitz [29] analyzed sediment budgets at St Francisville during 1981-2004 and 1981-2007, respectively.Further, studies quantifying suspended loads for the other locations downstream from St. Francisville are scarce and are limited to very short periods (three years in maximum).For example, Allison et al. [38] carried out a short-term suspended load budget investigation at four sites covering approximately 372 km of the LmMR reach during 2008-2010.With respect to bedload, Nittrouer et al. [11] estimated bedform transport rates in the LmMR but only for its lowermost reach between rk 167 (near New Orleans) and 0 (near the river's gulf outlet).In addition, Nittrouer et al. [11] did not consider water discharge, river velocity and sediment transport equations in their analysis.Pereira et al. [39] applied a 1-dimensional mobile-bed model (HECRAS 4.0) along the LmMR reach from Tarbert Landing to Venice (about 17 km above the river's gulf outlet).However, their study compared between several bedload equations only and did not quantify bedload trends for any period.Nittrouer et al. [40] estimated the bedload transport rates along the LmMR reach, however, their analysis only considered backwater flow velocities with no consideration of size of bed sediments.Recently, Knox and Latrubesse [41] quantified bedload at Tarbert Landing for 2003-2011.To our best knowledge, no study exists on quantifying multi-decadal suspended loads synchronously at multiple locations along the uppermost 140 km of the LmMR.Similarly, comprehensive analysis of decadal trends in bedloads based on cross-sectional area, river discharge and stage, river velocity, and grain size are also unavailable along the uppermost LmMR.
Hence, this study aims to assess sediment transport by quantifying bedload and suspended load along the first 140 km of the Lowermost Mississippi River below its diversion at the Old River Control Structure to the Atchafalaya River.The assessment can help in comparing spatiotemporal trends between bedloads and suspended loads, understanding their interactive relation, and allocating sites with maximum sediment trapping capacity, which can have relevant implications for the ongoing sediment diversion projects (both executed and planned) [42,43] and deltaic land restoration efforts.The specific objectives of this study include: (1) quantifying bedload transport at three locations of the uppermost 140 km reach of the LmMR (including Tarbert Landing); (2) analyzing long term daily, seasonal and annual suspended sediment and sand loads at two locations downstream of Tarbert Landing; (3) assessing suspended load transport from Tarbert Landing to Baton Rouge; and (4) investigating the relationship between bedload and suspended load along the studied river reach.

Study Site
This study focused on three locations covering the first 140 km of the LmMR: Tarbert Landing, St Francisville and Baton Rouge (Figure 1).The first location is approximately 10-15 km downstream of the river's diversion to the Atchafalaya River.The diversion is facilitated by an engineering complex, the Old River Control Structure (ORCS), which was constructed in 1963 for preventing the Atchafalaya River (AR) from capturing the bulk of Mississippi River water (>30%).Post-ORCS construction, about 25% of the Mississippi River water is diverted into the Atchafalaya River artificially under normal flow conditions.Atchafalaya River (AR) from capturing the bulk of Mississippi River water (>30%).Post-ORCS construction, about 25% of the Mississippi River water is diverted into the Atchafalaya River artificially under normal flow conditions.The three locations consist of monitoring stations jointly managed by the United States Geological Survey (USGS) and United States Army Corps of Engineers (USACE).The duration and availability of daily discharge, daily river stage and sediment records useful for this study vary according to the locations specifically (see Sections 2.2 and 2.3).Since long-term sediment and sandload analyses have been carried about at Tarbert Landing previously, this site was used only for bedload estimation.However, sediment load, sand load and bedload were analyzed at St Francisville and Baton Rouge.

River Stage, Discharge and Bathymetric Data
Daily River Stage records (RSd in m) for Tarbert Landing were collected from the USACE gauge station at Red River Landing (about 4 rk downstream from Tarbert Landing) during 1973-2015, while RSds for St Francisville were collected at Bayou Sara (about 8 rk upstream of St Francisville) during 1978-2015.In addition, RSds were also collected at Baton Rouge during 2004-2015 from USACE.Similarly, daily discharge records (Qd in cms) were collected from 1 January 1978 to 31 December 2015 and 1 January 2004 to 31 December 2015 at Tarbert Landing (from USACE) and Baton Rouge (from USGS) respectively.Both USACE and USGS did not have daily discharge records for St Francisville.However, the Qds at Tarbert Landing were used for daily sediment and sand load estimations at St Francisville because the river discharge is approximately same for corresponding river stages at Tarbert Landing and St Francisville [36].River discharge measurements were used for The three locations consist of monitoring stations jointly managed by the United States Geological Survey (USGS) and United States Army Corps of Engineers (USACE).The duration and availability of daily discharge, daily river stage and sediment records useful for this study vary according to the locations specifically (see Sections 2.2 and 2.3).Since long-term sediment and sand-load analyses have been carried about at Tarbert Landing previously, this site was used only for bedload estimation.However, sediment load, sand load and bedload were analyzed at St Francisville and Baton Rouge.

River Stage, Discharge and Bathymetric Data
Daily River Stage records (RS d in m) for Tarbert Landing were collected from the USACE gauge station at Red River Landing (about 4 rk downstream from Tarbert Landing) during 1973-2015, while RS d s for St Francisville were collected at Bayou Sara (about 8 rk upstream of St Francisville) during 1978-2015.In addition, RS d s were also collected at Baton Rouge during 2004-2015 from USACE.Similarly, daily discharge records (Q d in cms) were collected from 1 January 1978 to 31 December 2015 and 1 January 2004 to 31 December 2015 at Tarbert Landing (from USACE) and Baton Rouge (from USGS) respectively.Both USACE and USGS did not have daily discharge records for St Francisville.However, the Q d s at Tarbert Landing were used for daily sediment and sand load estimations at St Francisville because the river discharge is approximately same for corresponding river stages at Tarbert Landing and St Francisville [36].River discharge measurements were used for suspended load analysis specifically and both river stage and discharge measurements were used for bedload analysis.
For this study, bathymetric data were collected at Tarbert Landing, St Franciville and Baton Rouge during 2004 and 2013 from USACE.USACE conducted hydrographic surveys along the LmMR reach from Black Hawk, Louisiana (rk 521.4,just above the ORCS) to the river's Gulf Outlet at Head of Passes (rk 0) in 1992, 2004 and 2013.The surveys consist of river channel cross-sections (CS) at 400 m intervals collected using single beam fathometer and multibeam side scan sonar.Each CS runs across the river (perpendicular to the river flow) at a given point and includes river bank and river bed elevation at given widths.The annual bed elevation was deduced from the bathymetric measurements of three cross-sections, which were further used to quantify bed load transport rates at the sites.

Sediment and Sand Concentration Data
Records on suspended sediment concentration (SSC) (mg/L) and percentage of silt clay fractions in SSC (particles with diameter <0.0625 mm) were collected from USGS at St Francisville (during 1978-2015) and Baton Rouge (during 2004-2015).USGS uses depth-integrated suspended sampling measures with the help of several isokinetic point samplers for collecting SSC records in the sites.A total of 437 SSC samples were collected at St Francisville, while 149 SSC samples were collected at Baton Rouge during their corresponding sampling periods, respectively.The sampling records are assumed to be sufficient and unbiased across all years, seasons and flow regimes because USGS had carried out 1 to 3 sediment measurements for each month during sampling periods at the corresponding sites.More details about SSC collection, processing techniques, error adjustments and reliability of the records can be found in several reports published by USGS [44][45][46].

Bedload Transport Estimation
In this study, we used the Engelund-Hansen [47] sediment transport equation to estimate daily bedload transport rates at Tarbert Landing, St Francisville and Baton Rouge of the uppermost 140 km of the LmMR reach during 2004-2015.This equation ignores the use of critical shear stress and has been reported by other studies (e.g., [48][49][50]) to be more accurate than other sediment transport equations.Engelund-Hansen used the stream power concept [51] and the similarity principle to develop the following relationship between friction factor (f ), dimensionless bedload discharge (Φ) and dimensionless bed shear stress (θ) for alluvial channels (as noted in its simplest and most direct form by Stevens and Yang [52]): The components f, Φ, and θ in Equation ( 2) can be given by the following formulae: where g is the acceleration due to gravity (9.81 m/s 2 ); S is the slope of the channel; d is the channel depth in meters; V is the water velocity along the channel in m/s; q s is the bedload transport rate in kilogram/sec per meter of width (kg/sec-m); Υ s is the specific weight of the sediment (2650 kg/m 3 ); D 50 is the grain size (i.e., average diameter for corresponding sediment particles in meter); and S g is the specific gravity of sediment (2.68).Furthermore, the bedload transport rate can be calculated by substituting for f, Φ, θ, and Υ s in Equation (1) as follows: Finally, the daily bedload transport rate (in tons/day) at the three locations along the LmMR was calculated by: Q s = 95.24× RW × q s (6) where 95.24 is a unit conversion factor for converting bedload in kg/s-m to tons/day-m and RW is the daily river width at the location.
The methods for estimating channel depth, slope, flow velocity, river width, and grain-size in Equations ( 5) and ( 6) are explained in the following Sections 2.4.1-2.4.5.

River Depth (d)
The daily river depth at Tarbert Landing, St Francisville and Baton Rouge during 2004-2015 was obtained by subtracting the average annual riverbed elevation (RE bed ) at a given location from the daily river stage in the location.The average annual RE bed at each location was obtained from USACE bathymetric measurements for the years 2004

Riverbed Slope (S)
Riverbed slope between two consecutive cross-sections at a given location is the difference between riverbed elevations of both cross-sections (higher-lower) divided by the distance between the cross-sections (400 m-see Section 2.2).For this study, riverbed slopes were calculated annually for each location during 2004-2015.

Average Flow Velocity (V)
The velocity of water discharge (in m/s) at Tarbert Landing, St Francisville and Baton Rouge for each day during the period 2004-2015 was calculated as: where Q d is the daily water discharge in cms and A is the daily cross-sectional area of the river channel at each location in m 2 .The daily cross-sectional area at each location during the period was further calculated as the product of daily river depth (d in meter-see Section 2.4.1) and daily river width (RW in meter-see Section 2.4.4 below).

River Width (RW)
River width measurements are not available for any of the three sites along the LmMR.Therefore, to estimate the daily river width, Landsat 7 Satellite images with geographic references (in .tiffformat) were first collected during 2004-2015 using the USGS Earth Explorer tool (Figure 2).USGS Earth Explorer had 2-4 images for each month during 2004-2015, such that each image was captured on separate days.For this study, 55 images of LmMR at Tarbert Landing, 37 images at St Francisville and 42 images at Baton Rouge were collected during 2004-2015.The images were collected for days representing the lowest to highest river stages at Tarbert Landing, St Francisville and Baton Rouge along the LmMR.All images were opened in ArcMap 10.2 and the corresponding river widths at all locations were measured using the 'Measure' tool.Now, a linear-trend line was fitted between the river widths (dependent: y) and their corresponding river stages (independent: x) at each location.Thus obtained were trend-line equations at Tarbert Landing, St Francisville and Baton Rouge (Table 1).These equations gave robust R 2 values (more than 0.8) and considerably less root mean square errors (RMSE) between measured and predicted river widths at all locations (Table 1); therefore, we did not consider using other trend-line functions (such as polynomial, power, exponential and/or moving average).The equations were further used to estimate daily river widths at these locations during 2004-2015 respectively.
Water 2017, 9, 716 7 of 30 collected for days representing the lowest to highest river stages at Tarbert Landing, St Francisville and Baton Rouge along the LmMR.All images were opened in ArcMap 10.2 and the corresponding river widths at all locations were measured using the 'Measure' tool.Now, a linear-trend line was fitted between the river widths (dependent: y) and their corresponding river stages (independent: x) at each location.Thus obtained were trend-line equations at Tarbert Landing, St Francisville and Baton Rouge (Table 1).These equations gave robust R 2 values (more than 0.8) and considerably less root mean square errors (RMSE) between measured and predicted river widths at all locations (Table 1); therefore, we did not consider using other trend-line functions (such as polynomial, power, exponential and/or moving average).The equations were further used to estimate daily river widths at these locations during 2004-2015 respectively.The image was directly opened in Arc Map 10.2 and then the river width was measured using the 'Measure' tool.The same procedure was followed for measuring river width at Tarbert Landing, St Francisville and Baton Rouge for, respectively, 55, 37 and 42 days during 2004-2015.
Table 1.Trend-line equations between river width (y: dependent) and river stage (x: independent) for Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the LmMR during 2004-2015.These equations were further used for estimating daily river width from daily river stages at the three locations.

Representative Grain Sizes
The representative grain sizes (D 50 ) for each site for the calculation of bedload transport rates were deduced from the information table (Table 1) provided by Nordin and Queen [53].According to their data, 1.7% coarser sediments (>0.063 mm) in the LmMR were between 0.063 and 0.125 mm in size, 96.2% coarser sediments were between 0.125 and 0.25 mm, and 99.9% of the coarser sediments were between 0.25 and 0.5 mm.Therefore in this study, we used three sediment sizes of 0.125, 0.25 and 0.5 mm for all study sites to represent 99.9% of the bedload sediments in the river reach.

Development of Discharge-Sediment and Discharge-Sand Rating Curves
We followed the procedure used by Joshi and Xu [31] to develop discharge-sediment and discharge-sand rating curves at St Francisville and Baton Rouge, respectively.Daily sediment load (DSL in tonnes/day) was computed by multiplying SSC with the corresponding Q d for all the sampling dates at both sites as: where 0.0864 is a unit conversion factor for converting the sand mass to tons per day.Similarly, sand concentration (SSC s in mg/L) was calculated by subtracting silt/clay concentration [SSC f (mg/L) = % silt/clay × SSC] from SSC.Then SSC in Equation ( 8) was replaced by SSC s for calculating daily sand loads (DSL s in tonnes/day) for all sampling dates at both sites.
The dependent variables, DSL or DSL s (y) and the independent variable, Q d (x) were transformed to their natural logarithm (ln), and the relationship between them was evaluated by subjecting them to various linear and polynomial rating curves.Then, the "log transformed (ln) Q d s" were fitted in the rating curve equations to get "predicted ln DSL and DSL s ".Then, the "predicted ln DSL and DSL s " were transformed back (by taking their exponential values) to attain the "predicted DSL and DSL s ".Correction factor (CF) given by Gray et al. [54] as a modification of Duan [55] was applied to check the potential log-biasing in the retransformation procedure.Finally, all rating curves were evaluated according to the following four criteria: regression coefficient of the curves (R 2 ≥ 0.8), root mean square errors of the predicted DSL and DSL s (RMSE) (the lower the better), standard error (SE) of the curves (in ln units) (also, the lower the better) and a graphical assessment (good visual agreement between corresponding calibrated and predicted DSLs and DSL s s respectively) [31,56,57].
At first, single linear and polynomial sediment and sand rating curves were applied for the corresponding whole periods at both stations (i.e., 1978-2015 at St Francisville and 2004-2015 at Baton Rouge).At St Francisville, the four evaluation criteria did not meet for single sediment and sand rating curves: low R 2 (0.7 and 0.75 for linear and polynomial sediment rating curves; 0.74 and 0.78 for linear and polynomial sand rating curves) (Tables 2 and 3); comparatively higher RMSE (186,620 and 174,238 for linear and polynomial sediment rating curves; 66,486 and 59,781 for linear and polynomial sediment rating curves) and SE (0.805 for sediment and 0.496 for sand rating curves) (Tables 4 and 5); and a poor visual agreement between corresponding measured and predicted DSLs and DSL s s (Figures 3 and 5).Therefore, approximately decadal linear and polynomial rating curves were further applied at this station because decadal intervals have been previously reported to minimize both intra-and inter-annual variability in sediment samples [30,31].At Baton Rouge, interestingly, single polynomial rating curves met all four evaluation criteria (Tables 2-5, Figures 3 and 5).
Finally, all the evaluation criteria were met at St Francisville using polynomial sediment rating curves during the periods: 1978-1987, 1988-1997, 1998-2007 and 2008-2015, and polynomial sand rating curves during the periods: 1978-1987, 1988-1996, 1998-2007, 2008-2015.They gave the predicted DSLs and DSL s s most approximate to their corresponding calibrated measurements from USGS (Tables 4 and 5, Figures 4 and 6).We also found that the use of a correction factor improved the RMSEs and decreased differences between predicted and calibrated DSLs, while overestimated DSL s s (Tables 4 and 5).Hence, we used correction factors for sediment rating curves only at this Water 2017, 9, 716 9 of 28 station.The year 1997 had a very low R 2 for sand rating curves (Lin: 0.58 and Poly: 0.61) which was reducing R 2 s and increasing RMSEs and SEs when combined with its preceding or succeeding decade (i.e., 1988-1996 and 1998-2007).Therefore, DSL s for 1997 at St Francisville was calculated using a non-rating curve approach (described in the following Section 2.6).Similarly, at Baton Rouge a single polynomial rating curve with no correction factor gave the best approximation between predicted and calibrated DSLs and DSL s s, respectively (Tables 4 and 5, Figures 3 and 5).2).

Non-Rating Curve Approach for Sand Load Calculation
A non-rating curve approach was used for prediction DSLs at St Francisville during 1997 because of the same reasons mentioned in Section 2.5.In this method, the sand concentration for the first sampling date of the year was supposed to be equal to all consecutive days of missing concentrations until the next sampling date.In addition, the concentration of the first sampling date was used for all non-available dates earlier that year.Finally, all DSLss during 1997 at St Francisville were calculated using Equation (8).

Annual and Seasonal Bedloads and Suspended Loads
We calculated annual bedload transport rates (Qs) at Tarbert Landing, St Francisville and Baton Rouge as the sum of daily Qs between 1 January and 31 December for the three locations.The error range value for annual bedload transport rates was considered as ±24% based on the accuracy of Engelund-Hansen equation noted by Van Rijn [49].Van Rijn [49] reported that bedloads measured in three field sites and one lab condition were 61 to 87% (average 76%) of the corresponding bedloads computed from the Engelund-Hansen equation.
We also calculated annual suspended sediment and sand loads at St Francisville and Baton Rouge as the sum of DSL and DSLs between 1 January and 31 December for both sites.The error range value of ±18% given by Joshi and Xu [31] was used to adjust for errors in annual sediment and sand loads.These loads were used to analyze daily and annual sediment and sand load trends at both locations.Similarly, monthly sediment and sand loads were used to analyze their seasonal trends.3. It is to note that both predicted and measured SL s during 1997 were eliminated in this comparison because of the low R 2 value of both rating curves during this year (please see Table 2).

Non-Rating Curve Approach for Sand Load Calculation
A non-rating curve approach was used for prediction DSL s at St Francisville during 1997 because of the same reasons mentioned in Section 2.5.In this method, the sand concentration for the first sampling date of the year was supposed to be equal to all consecutive days of missing concentrations until the next sampling date.In addition, the concentration of the first sampling date was used for all non-available dates earlier that year.Finally, all DSL s s during 1997 at St Francisville were calculated using Equation (8).

Annual and Seasonal Bedloads and Suspended Loads
We calculated annual bedload transport rates (Q s ) at Tarbert Landing, St Francisville and Baton Rouge as the sum of daily Q s between 1 January and 31 December for the three locations.The error range value for annual bedload transport rates was considered as ±24% based on the accuracy of Engelund-Hansen equation noted by Van Rijn [49].Van Rijn [49] reported that bedloads measured in three field sites and one lab condition were 61 to 87% (average 76%) of the corresponding bedloads computed from the Engelund-Hansen equation.
We also calculated annual suspended sediment and sand loads at St Francisville and Baton Rouge as the sum of DSL and DSL s between 1 January and 31 December for both sites.The error range value of ±18% given by Joshi and Xu [31] was used to adjust for errors in annual sediment and sand loads.These loads were used to analyze daily and annual sediment and sand load trends at both locations.Similarly, monthly sediment and sand loads were used to analyze their seasonal trends.The annual and monthly sediment and sand loads were also compared with pre-existing sediment and sand loads at Tarbert Landing for overlapping periods.

Bedload Transport Rates
Total and annual bedload rates (Q s ) increased gradually from Tarbert Landing to Baton Rouge along the LmMR during 2004-2015 for all representative grain sizes (Figures 7-9).The total Q s for grain size of 0.125 mm were found to be 82.8,83.6 and 96 MT at Tarbert Landing, St Francisville and Baton Rouge respectively during the 12 years (Figures 7-10).Similarly, annual Q s for the same grain size of 0.125 mm averaged 6.9 million tons (MT)/year, ranging from 2.  9 and 10).When the representative grain size decreased by half, the annual Q s for all three locations along the LmMR also decreased by half during 2004-2015 (Figures 7-9).Therefore, the annual Q s for the grain size of 0.25 mm were half of the annual Q s for the grain size of 0.125 mm at Tarbert Landing, St Francisville and Baton Rouge in each of the 12 years during 2004-2015 (Figures 7-9).Further, the annual Q s for the grain size of 0.5 mm were half of the annual Q s for the grain size of 0.25 mm in each year (Figures 7-9).
and sand loads at Tarbert Landing for overlapping periods.

Bedload Transport Rates
Total and annual bedload rates (Qs) increased gradually from Tarbert Landing to Baton Rouge along the LmMR during 2004-2015 for all representative grain sizes (Figures 7-9).The total Qs for grain size of 0.125 mm were found to be 82.8,83.6 and 96 MT at Tarbert Landing, St Francisville and Baton Rouge respectively during the 12 years (Figures 7-10).Similarly, annual Qs for the same grain size of 0.125 mm averaged 6.9 million tons (MT)/year, ranging from 2.  9 and 10).When the representative grain size decreased by half, the annual Qs for all three locations along the LmMR also decreased by half during 2004-2015 (Figures 7-9).Therefore, the annual Qs for the grain size of 0.25 mm were half of the annual Qs for the grain size of 0.125 mm at Tarbert Landing, St Francisville and Baton Rouge in each of the 12 years during 2004-2015 (Figures 7-9).Further, the annual Qs for the grain size of 0.5 mm were half of the annual Qs for the grain size of 0.25 mm in each year (Figures 7-9 7).At Baton Rouge, the trends in bed load transport rates were exactly the opposite.In this location, annual Qs for all grain sizes experienced a decreasing trend throughout the study period (19.4 MT in 2004 decreased to 2.3 MT in 2015 for grain size of 0.125 mm; 9.7 MT in 2004 decreased to 1.2 MT in 2015 for grain size of 0.25 mm; and 4.8 MT in 2004 decreased to 0.6 MT in 2015 for grain size of 0.5 mm) (Figure 9).At St Francisville, the location between Tarbert Landing and Baton Rouge, annual Qs for all grain sizes did not experience any significant change during 2004-2015 (Figure 8).On a few occasions, annual Qs for all grain sizes increased or decreased abruptly for two to three years in continuum at all three locations (Figures 7-10).for grain size of 0.5 mm) (Figure 7).At Baton Rouge, the trends in bed load transport rates were exactly the opposite.In this location, annual Q s for all grain sizes experienced a decreasing trend throughout the study period (19.4 MT in 2004 decreased to 2.3 MT in 2015 for grain size of 0.125 mm; 9.7 MT in 2004 decreased to 1.2 MT in 2015 for grain size of 0.25 mm; and 4.8 MT in 2004 decreased to 0.6 MT in 2015 for grain size of 0.5 mm) (Figure 9).At St Francisville, the location between Tarbert Landing and Baton Rouge, annual Q s for all grain sizes did not experience any significant change during 2004-2015 (Figure 8).On a few occasions, annual Q s for all grain sizes increased or decreased abruptly for two to three years in continuum at all three locations (Figures 7-10).
Monthly averages of daily bedload rates (Q s ) along the LmMR during 2004-2015 showed linearly increasing trends with monthly averages of daily river discharge (Q d ) and river stage (RS d ) at Tarbert Landing and St Francisville respectively, however, the increasing trends with both variables were not eminent at Baton Rouge (Figures 11 and 12).The rate of increase of Q s s with both variables was highest at St Francisville (R 2 = 0.95 for Q s vs. Q d linear trend line and 0.84 for Q s vs. RS d trend lines) (Figures 11 and 12) and second highest at Tarbert Landing (R 2 = 0.83 for Q s vs. Q d linear trend line and 0.77 for Q s vs. RS d trend lines) (Figures 11 and 12).At Baton Rouge, however, the relationship of monthly Q s with monthly Q d and RS d appeared to be more scattered during the 12-year study period (R 2 = 0.33 for Q s vs. Q d linear trend line and 0.32 for Q s vs. RS d trend lines) (Figures 11 and 12).These relationships were noted down for representative grain size of 0.125 mm.Exactly same relationships of monthly Q s with monthly Q d and RS d were found for the other two grain sizes (0.25 and 0.5 mm) at each location along the LmMR during 2004-2015 (figures for these two grain sizes were not kept here to avoid repetition).
Water 2017, 9, 716 16 of 30 Monthly averages of daily bedload rates (Qs) along the LmMR during 2004-2015 showed linearly increasing trends with monthly averages of daily river discharge (Qd) and river stage (RSd) at Tarbert Landing and St Francisville respectively, however, the increasing trends with both variables were not eminent at Baton Rouge (Figures 11 and 12).The rate of increase of Qss with both variables was highest at St Francisville (R 2 = 0.95 for Qs vs. Qd linear trend line and 0.84 for Qs vs. RSd trend lines) (Figures 11 and 12) and second highest at Tarbert Landing (R 2 = 0.83 for Qs vs. Qd linear trend line and 0.77 for Qs vs. RSd trend lines) (Figures 11 and 12).At Baton Rouge, however, the relationship of monthly Qs with monthly Qd and RSd appeared to be more scattered during the 12-year study period (R 2 = 0.33 for Qs vs. Qd linear trend line and 0.32 for Qs vs. RSd trend lines) (Figures 11 and 12).These relationships were noted down for representative grain size of 0.125 mm.Exactly same relationships of monthly Qs with monthly Qd and RSd were found for the other two grain sizes (0.25 and 0.5 mm) at each location along the LmMR during 2004-2015 (figures for these two grain sizes were not kept here to avoid repetition).

Suspended Sediment Loads
Annual sediment loads (SL) at St Francisville averaged 97.1 million tons (MT) from 1978 to 2015, ranging from 40 MT in 2012 to 178.7 MT in 1979 (Figure 13).The total sediment amount for the entire 38-year study period was found to be 3688.3MT (Figure 13).Annual SL was lower than 50 MT for 3 of the 41 years (2000, 2006, and 2012), 50 to 100 MT for 19 years (including the Mississippi flood year of 2011), 100 to 150 MT for 12 years and more than 150 MT for the remaining 4 years (Figure 13).Annual SL was higher for approximately the first-half of the study period (117.3MT from 1978 to 1999), however, it decreased substantially towards the latter-half (69 MT from 2000 to 2015) (Figure 13).Also, annual SL was more than 100 MT only for one year (103.4MT in 2004) during this period.Long-term continuous increasing or decreasing trend (even for 3/4 years) in annual SL was not found at St Francisville despite a few instances of abrupt decrease during for 2 years in continuum (Figure    At Baton Rouge, annual SL averaged 67.2 MT from 2004 to 2015 and ranged 28.3 MT (in 2006) to 92.9 MT (in 2009) (Figure 13).During this period, the total sediment amount was found to be 806.1 MT (Figure 13).Annual SL was lower than 50 MT for 2 out of 12 years (2006 and 2012), while either >50 or <100 MT for the remaining 10 years (Figure 13).The long-term continuous increasing or decreasing trends in annual SL were not found at this station (Figure 13).This result was similar to that of St Francisville.
Both sites had slightly different seasonal trends of average DSLs.At St Francisville, average DSL during 1978-2015 increased each year from January to its maximum in April (i.e., from 313,955 to 427,400 tonnes/day), then decreased from May to its minimum in September (from 401,143 to 99,043 tonnes/day), showing the highest sediment transport during early spring months (Figure 14).Average DSL increased again for the remaining three months of the year i.e., 110,000 tonnes/day in October to 287,000 tons/day in December (Figure 14).At Baton Rouge, however, average DSL during 2004-2015 decreased slightly from January to February (i.e., from 229,317 to 228,956 tonnes/day) (Figure 14).It further increased during the three spring months from 268,017 tonnes/day in March to 300,080 tonnes/day in May (Figure 14), followed by a drop to its minimum in mid-fall (i.e., from 252,035 tonnes/day in June to 75,735 tonnes/day in September) (Figure 14).The final three months had a continuous increase in average DSL like that found at St Francisville i.e., 79,789 tonnes/day in October to 179,541 tonnes/day in December (Figure 14).

Suspended Sand Loads
St Francisville had an annual sand load (SL s ) averaging 19.9 million tons (MT) during 1978-2015, varying from 5 MT in 2006 to 43 MT in 1979 (Figure 15).The total SL s for the entire 38-year study period was 754 MT (Figure 15).Annual SL s was lower than 10 MT for 3 of 38 years (1997, 2000, and 2006), 10 to 25 MT for 26 years (including the Mississippi flood year of 2011), 25 to 40 MT for 6 years and more than 40 MT for the remaining 3 years (1979, 1983, and 1993) (Figure 15).Average annual SL s was higher for approximately the first-half of the study period (24 MT/year from 1978 to 1999), however, like annual SL, it decreased substantially towards the later-half (13.1 MT from 2000 to 2015) (Figure 15).In fact, annual SL s were less than 20 MT for all 16 years during the later-half as compared to 5 out of 22 years for the first-half.Like annual SL, annual SL s at St Francisville did not have any long-term increasing or decreasing trend (even for 3/4 years) (Figure 15).At Baton Rouge, annual SL s averaged 16.1 MT during 2004-2015, ranging from 2.9 MT in 2006 to 28.3 MT in 2009 with the total sand amount of 193 MT (Figure 15).Annual SL was lower than 10 MT for 4 out of 12 years (2006, 2007, 2012, and 2014), 10 to 25 MT for 6 years and higher than 25 MT for the remaining 2 years (Figure 15).The continuous increasing or decreasing trend (even for 3/4 years) were not found in annual SL s s at Baton Rouge as well (like its annual sediment loads and St Francisville's annual sediment and sand loads) (Figure 15).
The seasonal trends of average DSL s s varied little during April and May.At St Francisville, average monthly DSL s followed the same trend as average monthly DSL during 1978-2015: it increased from January to its maximum in April (i.e., 64,329 to 102,934 tonnes/day), then decreased from May to its minimum in September (97,724 to 9687 tonnes/day) (Figure 16).During the remaining   We found gradual increase in estimated bedload transport rates (Qs) from Tarbert Landing to Baton Rouge during 2004-2015.In addition, the 12 years showed significantly higher Qs at Baton Gomez et al. [65] examined the river reach between St. Louis and Thebes (approximately from the confluence of the Missouri and Mississippi Rivers to the confluence of Ohio and Mississippi Rivers) and concluded that several levee breaks along this reach could have been caused by sediment and circumjacent sand deposits from the 1993 flood.The overbank sedimentation could have removed substantial stored bed/suspended sediments in the reach and eliminated several SSC sources for the LmMR reach [66].There were no sediment records available to find the effects of 1993 flood on Baton Rouge suspended loads; however, we did not observe any significant trend in annual SL and SL s at this site during the last decade.Even the LmMR flood in 2011 did not seem to make a significant impact on the suspended loads at Baton Rouge, since there was negligible difference in average annual SL and SL s during the two periods: 2004-2011 (pre-2011 flooding) (mean annual SL and SL s = 67.7 and 17 MT respectively) and 2012-2015 (post-2011 flooding) (mean annual SL and SL s = 66.2 and 14.1 MT respectively).These intra-location pre-and post-flooding suspended load trends at St Francisville and Baton Rouge indicate that floods probably have substantial impacts on suspended loads for longer periods, while negligible impacts for shorter periods.

Relationship between Bedload and Suspended Load in the LmMR
Our monthly averages of daily bedload rates (Q s ) along the LmMR showed linearly increasing trends with monthly averages of daily suspended sediment loads (DSL) and daily suspended sand loads (DSL s ) at Tarbert Landing (during the corresponding matching periods: 2004-2010 for Q s -DSL trends and 2004-2013 for Q s -DSL s ) and St Francisville (during 2004-2015), respectively (Figures 17 and 18).However, monthly Q s showed more scattered trends with monthly DSL and DSL s respectively at Baton Rouge (during 2004-2015) (Figures 17 and 18).These relationships were noted down for representative grain size of 0.125 mm and were exactly same for other grain sizes of 0.25 and 0.5 mm.The trends of Q s with DSL, DSL s , river discharge and river stage (see Section 3.1 for the latter two) possibly imply that bedload is more dependent on suspended load (specially suspended sand load), river discharge and river stage in upstream river locations with shallower beds, while the dependence gradually decreases as the river moves downstream and the bed elevation increases continuously.
Total and annual bedload rates (Q s ) for grain sizes 0.125, 0.25 and 0.5 mm at Tarbert Landing were about 4, 2 and 1% of the suspended sediment load and 18, 9 and 4.5% of the suspended sand load respectively during 2004-2010 (matching period for both loads at Tarbert Landing) (Figure 19).Furthermore, total and annual Q s at St Francisville for the aforementioned grain sizes were respectively about 10, 5 and 2.5% of its suspended sediment load and 53, 26.5 and 13% of its suspended sand load during 2004-2015 (matching period for both loads at St Francisville) (Figure 19).Finally, the total and average annual Q s at Baton Rouge for these grain sizes were approximately 12, 6 and 3% of its suspended sediment load and 50, 25 and 12.5% of its suspended sand load respectively during 2004-2015 (Figure 19).Most of these values were within the ranges of percentage of bedload in suspended sediment (1% to 10%) and suspended sand load (10% to 33%) noted by previous studies ([ [7][8][9]; see Section 1).However, the only noticeable outliers in this regard were substantially higher percentage of Q s in suspended sand load of St Francisville and Baton Rouge for grain size of 0.125 mm.These higher ranges were likely because of higher bed elevation at these two locations, gradually decreasing grain sizes from upstream to downstream in the LmMR (as explained in first paragraph of Section 4.1) and substantially lower suspended sand loads in both locations as compared to the upstream location at Tarbert Landing (see Section 4.1).
respectively during 2004-2015 (Figure 19).Most of these values were within the ranges of percentage of bedload in suspended sediment (1% to 10%) and suspended sand load (10% to 33%) noted by previous studies ([ [7][8][9]; see Section 1).However, the only noticeable outliers in this regard were substantially higher percentage of Qs in suspended sand load of St Francisville and Baton Rouge for grain size of 0.125 mm.These higher ranges were likely because of higher bed elevation at these two locations, gradually decreasing grain sizes from upstream to downstream in the LmMR (as explained in first paragraph of Section 4.1) and substantially lower suspended sand loads in both locations as compared to the upstream location at Tarbert Landing (see Section 4.1).The findings for bedload transport in this study are opposite to the trends found for total suspended load along the LmMR.Apart from the probable reasons listed for corresponding trends in both, the contrasting differences could also be possible because suspended loads were estimated from sediment/sand-water discharge rating curves, while bedload transport rates were quantified using the concept of boundary and critical shear stresses and pre-studied sediment transport equations.Nittrouer et al. [40] found opposite spatial trends in Qs than the findings of this study along the LmMR reach (i.e., gradual increase from Baton Rouge to Tarbert Landing) but only in  Difference between total bedload transport rates (Qs) for representative grain size of 0.125 mm (A), suspended sediment and sand loads, and total supply at Tarbert Landing [TBL, at river kilometer (rk) 490)], St Francisville (St F, rk 419), and Baton Rouge (BTR, rk 367.5) along the Lowermost Mississippi River during 2004-2010 (the matching period for these components at all three locations).X-axis distances are distance of the sites from the Head of the Passes (rk 0) near the LmMR's Gulf of Mexico outlet.Note: the suspended sand load is a coarser subset of suspended sediment load, hence, the total sediment supply at a given location along the river was calculated as the sum of total bedload and total suspended sediment load at the location.

Sediment Trapping in the LmMR
From this analysis, it seems likely that the first 20-25 km reach of the LmMR just below its diversion to the Atchafalaya River (covering Tarbert Landing) has been a major sediment trapping location over the last four decades.Lower bedload transport along this reach during the last decade suggests that bedloads possibly acted as sources of sand loads near the riverbed because both loads can be intermixed especially in laterally migrating sand dunes under low and medium flows [40,67,68].However, the sand particles possibly get trapped in the sediment channel bars present along this reach.Wang and Xu [35] noted high three-decadal sediment accumulation in the three sediment bars present at 18, 24 and 26 km downstream of the Old River Control Structure respectively (total sediment load of 36 million tons during 1985-2013).This may have been the primary reason of lower proportion of sand outflow between Tarbert Landing and St Francisville as compared to sediment outflow over the last three and half decades (mean annual SL at Tarbert Landing was 27% higher than St Francisville during 1980-2010, while mean annual SLs was about 20% higher during the same period).The total sediment supply at Tarbert Landing was higher than St Francisville and Baton Rouge during 2004-2010 despite its lower Qs during the same period (Figure 19; Section 4.2).This observation can probably imply that the sediment trapping capacity gradually reduces downstream along the LmMR reach particularly below its uppermost 20-25 km stretch.
The importance of suspended loads for restoration of the MRDP has been previously well documented [32][33][34]69].Continuous movement of suspended load is critical in long-term delta building and restoration processes [34,70].Therefore, sediment management in the LmMR should focus on increasing outflow of suspended loads trapped in the uppermost 20-25 km reach to the downstream locations (especially to St Francisville) which can further help in planning and executing sediment diversion projects along the river subsequently.In addition, the bed and suspended load estimates in this study also indicate that significant spatiotemporal variability can occur in sediment supplies from upstream to downstream reaches of highly engineered alluvial rivers on a global scale.X-axis distances are distance of the sites from the Head of the Passes (rk 0) near the LmMR's Gulf of Mexico outlet.Note: the suspended sand load is a coarser subset of suspended sediment load, hence, the total sediment supply at a given location along the river was calculated as the sum of total bedload and total suspended sediment load at the location.
The findings for bedload transport in this study are opposite to the trends found for total suspended load along the LmMR.Apart from the probable reasons listed for corresponding trends in both, the contrasting differences could also be possible because suspended loads were estimated from sediment/sand-water discharge rating curves, while bedload transport rates were quantified using the concept of boundary and critical shear stresses and pre-studied sediment transport equations.Nittrouer et al. [40] found opposite spatial trends in Q s than the findings of this study along the LmMR reach (i.e., gradual increase from Baton Rouge to Tarbert Landing) but only in specific river discharges limited to 11,750 and 17,500 cms.They estimated negligible differences in Q s during high flows throughout the reach (from Tarbert Landing to Baton Rouge).The differences in the Q s estimates of the two studies could be because they considered back water flows in the LmMR which were ignored in this study.
The total sediment supply (bedload + suspended sediment load) for all grain sizes during the overlapping period of 2004-2010 at all locations was highest at Tarbert Landing (931, 911 and 893 MT for grain sizes of 0.125, 0.25, 0.5 respectively), and decreased gradually at two downstream locations along the LmMR (550, 525 and 500 MT at St Francisville and 544, 504 and 485 MT at Baton Rouge for the three aforementioned grain sizes respectively) (Figure 19).Furthermore, for the two locations downstream of Tarbert Landing, the total sediment loads were a little higher at Baton Rouge than St Francisville for the grain sizes of 0.125 and 0.25 mm respectively (904 MT for a grain size of 0.125 mm and 855 MT for 0.25 mm at Baton Rouge as compared to 893 for 0.125 mm and 850 MT for 0.25 mm at St Francisville) during 2004-2015 (the matching period for both locations) (Figure 19).For the remaining grain size of 0.5 mm the total suspended loads were almost equal at both locations during the 12-year period (~830 MT) (Figure 19).

Sediment Trapping in the LmMR
From this analysis, it seems likely that the first 20-25 km reach of the LmMR just below its diversion to the Atchafalaya River (covering Tarbert Landing) has been a major sediment trapping location over the last four decades.Lower bedload transport along this reach during the last decade suggests that bedloads possibly acted as sources of sand loads near the riverbed because both loads can be intermixed especially in laterally migrating sand dunes under low and medium flows [40,67,68].However, the sand particles possibly get trapped in the sediment channel bars present along this reach.Wang and Xu [35] noted high three-decadal sediment accumulation in the three sediment bars present at 18, 24 and 26 km downstream of the Old River Control Structure respectively (total sediment load of 36 million tons during 1985-2013).This may have been the primary reason of lower proportion of sand outflow between Tarbert Landing and St Francisville as compared to sediment outflow over the last three and half decades (mean annual SL at Tarbert Landing was 27% higher than St Francisville during 1980-2010, while mean annual SL s was about 20% higher during the same period).The total sediment supply at Tarbert Landing was higher than St Francisville and Baton Rouge during 2004-2010 despite its lower Q s during the same period (Figure 19; Section 4.2).This observation can probably imply that the sediment trapping capacity gradually reduces downstream along the LmMR reach particularly below its uppermost 20-25 km stretch.
The importance of suspended loads for restoration of the MRDP has been previously well documented [32][33][34]69].Continuous movement of suspended load is critical in long-term delta building and restoration processes [34,70].Therefore, sediment management in the LmMR should focus on increasing outflow of suspended loads trapped in the uppermost 20-25 km reach to the downstream locations (especially to St Francisville) which can further help in planning and executing sediment diversion projects along the river subsequently.In addition, the bed and suspended load estimates in this study also indicate that significant spatiotemporal variability can occur in sediment supplies from upstream to downstream reaches of highly engineered alluvial rivers on a global scale.Information on long-term sediment supplies along several different sediment-starved rivers can benefit the river management practices.

Uncertainty and Constraints in Bedload Estimation
The estimation of bedloads (Q s ) in this study is one of the few attempts in quantifying daily Q s at multiple sites along a large regulated alluvial river for more than a decade.Most previous studies quantifying Q s along the LmMR have applied measurements of bedform heights using multibeam echosounder surveys in their two or three dimensional models [11,41,71,72].Since, we did not have such measurements; we only used Engelund-Hansen equation primarily developed for and tested in laboratory conditions.Furthermore, we did not have site-specific bedload measurements like suspended sediment concentrations to validate the estimates of bedload transport rates along the LmMR.This is a main uncertainty in the bedload estimates of this study.
There was a need to estimate annual riverbed elevation and other daily variables such as flow velocity, river width and cross-sectional area at the three sites for quantifying daily Q s from Equation (1).The calculations for all these parameters were subjected to potential errors.Riverbed elevations (Re bed ) were only available for 2004 and 2013 at all locations; therefore, Re bed had to be interpolated for the period of 2005-2012 and of 2014-2015.In interpolation, the difference in Re bed for 2004-2013 was only considered as overall increase or decrease, but any specific increase or decrease in Re bed within the years between 2004 and 2013 was ignored.The uncertainty for average flow velocity was dependent on the errors for cross-sectional area.The calculation of daily cross-sectional area was subjected to the potential errors in estimating river width through earth explorer and arc map tools.Also, the daily cross-sectional areas at all locations were probably overestimated because the river channel was considered as rectangle and the cross-sectional areas were measured as a product of river depth and river width.This study could not find a method to measure the cross-sectional area such that it exactly coincided with the half oval shaped riverbed cross-section.
Despite the above-mentioned uncertainties, this study provides plausible estimates of bedload rates that are within previously reported Q s estimation ranges, i.e., 1% to 10% of suspended sediment load, or 10% to 33% of suspended sand load.A ±24% margin of error is given for all Q s estimates (see Section 2.7) to properly document estimation performance.

Conclusions
This study analyzed multi-decadal sediment transport along the first 140 kilometers of the Lowermost Mississippi River downstream of its diversion to the Atchafalaya River.The study found that in the past decade, bedload transport rates at Tarbert Landing upstream increased slightly, while the rates at Baton Rouge downstream declined largely.Over the past two decades, both suspended sediment and sand loads were consistently higher at Tarbert Landing than those at St Francisville and Baton Rouge.These results indicate an accumulation of coarser sediments in the uppermost river reach, i.e., 20-25 kilometers below the Old River Control Structure, covering Tarbert Landing.Bedloads seemed to increase with the increasing suspended loads, river discharge and river stage at upstream locations, i.e., at Tarbert Landing and St Francisville; however, such a relationship was not noticeable downstream at Baton Rouge.A clear decreasing trend in annual suspended loads was identified at St Francisville after the 1993 major flood, showing that floods play an important role in longer-term sediment transport.This study implies that substantial spatiotemporal variability can exist in sediment transport of a large alluvial river, and that future sediment management in the Lowermost Mississippi River should find solutions for mobilizing coarse sediments trapped in the studied river reach (i.e., the uppermost 140 kilometers of LmMR from Tarbert Landing to Baton Rouge) to the river reach below Baton Rouge for coastal restoration and new land building.

Figure 1 .
Figure 1.The Lowermost Mississippi River (LmMR) and three locations used in this study to analyze bedload and suspended load below the Old River Control Structure (ORCS): (1) Tarbert Landing (TBL) at river kilometer (or rk) 493, (2) St Francisville (St F) at rk 419, and (3) Baton Rouge (BTR) at rk 367.5.The Red River Landing (RRL) (at rk 486.5) just below TBL consists of long term river stage records from USGS gauge present at the site.Also, Bayou Sara (BS) (at rk 427) (about 9 rk above St F) is the nearest reference station for long term river stage at St F. The Bonnet Carre Spillway (BCS) is about 19 rk upstream of New Orleans.

Figure 1 .
Figure 1.The Lowermost Mississippi River (LmMR) and three locations used in this study to analyze bedload and suspended load below the Old River Control Structure (ORCS): (1) Tarbert Landing (TBL) at river kilometer (or rk) 493, (2) St Francisville (St F) at rk 419, and (3) Baton Rouge (BTR) at rk 367.5.The Red River Landing (RRL) (at rk 486.5) just below TBL consists of long term river stage records from USGS gauge present at the site.Also, Bayou Sara (BS) (at rk 427) (about 9 rk above St F) is the nearest reference station for long term river stage at St F. The Bonnet Carre Spillway (BCS) is about 19 rk upstream of New Orleans.

Figure 2 .
Figure 2. Screenshot of the Landsat 7 Satellite image for Tarbert Landing of the Lowermost Mississippi River.The image was collected in .tiffformat using the USGS Earth Explorer tool.The red horizontal line represents the river width of the LmMR at Tarbert Landing on the date of imagery.The image was directly opened in Arc Map 10.2 and then the river width was measured using the 'Measure' tool.The same procedure was followed for measuring river width at Tarbert Landing, St Francisville and Baton Rouge for, respectively, 55, 37 and 42 days during 2004-2015.

Figure 2 .
Figure 2. Screenshot of the Landsat 7 Satellite image for Tarbert Landing of the Lowermost Mississippi River.The image was collected in .tiffformat using the USGS Earth Explorer tool.The red horizontal line represents the river width of the LmMR at Tarbert Landing on the date of imagery.The image was directly opened in Arc Map 10.2 and then the river width was measured using the 'Measure' tool.The same procedure was followed for measuring river width at Tarbert Landing, St Francisville and Baton Rouge for, respectively, 55, 37 and 42 days during 2004-2015.

Figure 3 .
Figure 3.Comparison of sediment loads calculated from SSC records measured, processed, and calibrated by USGS (Measured SL) with those predicted from single sediment rating curves [either linear (Lin) or polynomial (poly)] (Predicted SL) at St Francisville during 1978-2015 (a-d) and Baton Rouge along the Lowermost Mississippi River during 2004-2015 (e-h).It is to note that Duan correction factors were applied for the curves denoted by "CF" in the figure (b, d, f, h), while the remaining curves denoted by "No CF" in the figure were analyzed without correction factors (a, c, e, g).

Figure 3 .
Figure 3.Comparison of sediment loads calculated from SSC records measured, processed, and calibrated by USGS (Measured SL) with those predicted from single sediment rating curves [either linear (Lin) or polynomial (poly)] (Predicted SL) at St Francisville during 1978-2015 (a-d) and Baton Rouge along the Lowermost Mississippi River during 2004-2015 (e-h).It is to note that Duan correction factors were applied for the curves denoted by "CF" in the figure (b, d, f, h), while the remaining curves denoted by "No CF" in the figure were analyzed without correction factors (a, c, e, g).

Figure 4 .
Figure 4. Comparison of Measured SL (please see caption of Figure 3 for definition) with those predicted from several sediment rating curves (Predicted SL) at St Francisville of the Lowermost Mississippi River during 1978-2015.Specific terminologies found in this figure such as Lin, Poly, CF, and No CF are same as explained in Figure 3.

Figure 4 . 30 Figure 5 .
Figure 4. Comparison of Measured SL (please see caption of Figure 3 for definition) with those predicted from several sediment rating curves (Predicted SL) at St Francisville of the Lowermost Mississippi River during 1978-2015.Specific terminologies found in this figure such as Lin, Poly, CF, and No CF are same as explained in Figure 3.Water 2017, 9, 716 12 of 30

Figure 5 .
Figure 5.Comparison of sand loads calculated from SSC records measured, processed and calibrated by USGS (Measured SL s ) with those predicted from single sediment rating curves [either linear (Lin) or polynomial (poly)] (Predicted SL s ) at St Francisville of the Lowermost Mississippi River during 1978-2015 (a-d) and BTR during 2004-2015 (e-h).Specific terminologies found in this figure i.e., 'CF' and 'No CF' are same as explained in Figure 3.

Figure 6 .
Figure 6.Comparison of Measured SLs and Predicted SLs at St Francisville of the Lowermost Mississippi River during 1978-2015.Specific terminologies found in this figure such as Lin, Poly, CF, and No CF are same as explained in Figure3.It is to note that both predicted and measured SLs during 1997 were eliminated in this comparison because of the low R 2 value of both rating curves during this year (please see Table2).

Figure 6 .
Figure 6.Comparison of Measured SL s and Predicted SL s at St Francisville of the Lowermost Mississippi River during 1978-2015.Specific terminologies found in this figure such as Lin, Poly, CF, and No CF are same as explained in Figure3.It is to note that both predicted and measured SL s during 1997 were eliminated in this comparison because of the low R 2 value of both rating curves during this year (please see Table2).
3 (in 2006) to 12.8 MT/year (in 2015) at Tarbert Landing (Figures 7 and 10); 7 MT/year, ranging from 3.5 (in 2006) to 10.3 MT/year (in 2008) at St Francisville (Figures 8 and 10); and 8.1 MT/year, ranging from 2.2 (in 2014) to 19.4 MT/year (in 2004) at Baton Rouge (Figures 3 (in 2006) to 12.8 MT/year (in 2015) at Tarbert Landing (Figures 7 and 10); 7 MT/year, ranging from 3.5 (in 2006) to 10.3 MT/year (in 2008) at St Francisville (Figures 8 and 10); and 8.1 MT/year, ranging from 2.2 (in 2014) to 19.4 MT/year (in 2004) at Baton Rouge (Figures ).Specifically, annual Qs for all grain sizes experienced an increasing trend at Tarbert Landing from 2006 to 2015 (2.3 MT in 2006 increased to 12.8 MT in 2015 for grain size of 0.125 mm; 1.2 MT in 2006 increased to 6.4 MT in 2015 for grain size of 0.25 mm; and 0.6 MT in 2006 increased to 3.2 MT in 2015 for grain size of 0.5 mm) (Figure

Figure 7 .
Figure 7. Annual bedload transport rates (Qs) at Tarbert Landing of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D50) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 7 . 30 Figure 8 .
Figure 7. Annual bedload transport rates (Q s ) at Tarbert Landing of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D 50 ) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).Water 2017, 9, 716 15 of 30

Figure 8 .
Figure 8. Annual bedload transport rates (Q s ) at St Francisville of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D 50 ) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 8 .
Figure 8. Annual bedload transport rates (Qs) at St Francisville of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D50) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 9 .
Figure 9. Annual bedload transport rates (Qs) at Baton Rouge of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D50) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 10 .
Figure 10.Annual bedload transport rates (Qs) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River for a grain size (D50) of 0.125 mm.Note: Annual Qs for grain size of 0.25 mm were exactly half of annual Qs for grain size of 0.125 mm, while annual Qs for grain size of 0.5 mm were exactly half of annual Qs for grain size of 0.25 mm at all stations.Therefore, annual Qs for grain sizes of 0.25 and 0.5 mm had similar trends as annual Qs for grain size of 0.125 mm at all locations.

Figure 9 .
Figure 9. Annual bedload transport rates (Q s ) at Baton Rouge of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D 50 ) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 8 .
Figure 8. Annual bedload transport rates (Qs) at St Francisville of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D50) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 9 .
Figure 9. Annual bedload transport rates (Qs) at Baton Rouge of the Lowermost Mississippi River.Note: A, B and C represent the transport rates for representative grain size (D50) of 0.125, 0.25 and 0.5 mm respectively as selected for this analysis (please see Section 2.4.5).

Figure 10 .
Figure 10.Annual bedload transport rates (Qs) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River for a grain size (D50) of 0.125 mm.Note: Annual Qs for grain size of 0.25 mm were exactly half of annual Qs for grain size of 0.125 mm, while annual Qs for grain size of 0.5 mm were exactly half of annual Qs for grain size of 0.25 mm at all stations.Therefore, annual Qs for grain sizes of 0.25 and 0.5 mm had similar trends as annual Qs for grain size of 0.125 mm at all locations.

Figure 10 .
Figure 10.Annual bedload transport rates (Q s ) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River for a grain size (D 50 ) of 0.125 mm.Note:Annual Q s for grain size of 0.25 mm were exactly half of annual Q s for grain size of 0.125 mm, while annual Q s for grain size of 0.5 mm were exactly half of annual Q s for grain size of 0.25 mm at all stations.Therefore, annual Q s for grain sizes of 0.25 and 0.5 mm had similar trends as annual Q s for grain size of 0.125 mm at all locations.

Figure 11 .
Figure 11.Relationship of monthly average of daily discharge (Qd in cms) with the monthly average of daily bedload transport rates (Qs in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the Lowermost Mississippi River during 2004-2015.

Figure 11 . 30 Figure 12 .
Figure 11.Relationship of monthly average of daily discharge (Q d in cms) with the monthly average of daily bedload transport rates (Q s in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the Lowermost Mississippi River during 2004-2015.Water 2017, 9, 716 17 of 30

Figure 12 .
Figure 12.Relationship of monthly average of daily river stages (RS d in m) with the monthly average of daily bedload transport rates (Q s in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River during 2004-2015.

Figure 13 .
Figure 13.Annual Sediment Loads (SL) at Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the Lowermost Mississippi River.Note: The annual SLs at Tarbert Landing were taken from Rosen and Xu [28].

Figure 13 .
Figure 13.Annual Sediment Loads (SL) at Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the Lowermost Mississippi River.Note: The annual SLs at Tarbert Landing were taken from Rosen and Xu [28].

Figure 13 .
Figure 13.Annual Sediment Loads (SL) at Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the Lowermost Mississippi River.Note: The annual SLs at Tarbert Landing were taken from Rosen and Xu [28].

Figure 15 .
Figure 15.Annual Sand Loads (SLs) at Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the Lowermost Mississippi River.Note: The annual SLss at TBL were taken from Joshi and Xu [31].

4 . Discussion 4 . 1 .
Spatial and Temporal Variation in Bedload and Suspended Load along the LmMR

Figure 16 .
Figure 16.Monthly averages of daily sand load at Tarbert Landing (TBL, during 1973-2013), St Francisville (St F, during 1978-2015) and Baton Rouge (BTR, during 2004-2015) of the Lowermost Mississippi River.Note: The seasonal trend for DSL s at Tarbert Landing was taken from Joshi and Xu [31].

Figure 17 .
Figure 17.Relationship of monthly average of daily suspended sediment load (DSL in tons/day) with the monthly average of daily bedload transport rates (Qs in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River during.The relationship for monthly DSL and monthly Qs was during 2004-2010 at TBL (taken from Rosen and Xu [28]) and during 2004-2015 at St F and BTR.

Figure 17 .
Figure 17.Relationship of monthly average of daily suspended sediment load (DSL in tons/day) with the monthly average of daily bedload transport rates (Qs in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River during.The relationship for monthly DSL and monthly Qs was during 2004-2010 at TBL (taken from Rosen and Xu [28]) and during 2004-2015 at St F and BTR.Water 2017, 9, 716 24 of 30

Figure 18 .
Figure 18.Relationship of monthly average of daily suspended sand load (DSLs in tons/day) with the monthly average of daily bedload transport rates (Qs in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River during 2004-2015.The relationship for monthly DSLs and monthly Qs was during 2004-2013 at TBL (taken from Joshi and Xu [31]) (2015) and during 2004-2015 at St F and BTR.

Figure 18 .
Figure 18.Relationship of monthly average of daily suspended sand load (DSLs in tons/day) with the monthly average of daily bedload transport rates (Qs in tons/day) for representative grain-size of 0.125 mm (A) at Tarbert Landing (TBL), St Francisville (St F), and Baton Rouge (BTR) of the Lowermost Mississippi River during 2004-2015.The relationship for monthly DSLs and monthly Qs was during 2004-2013 at TBL (taken from Joshi and Xu [31]) (2015) and during 2004-2015 at St F and BTR.

Figure 19 .
Figure19.Difference between total bedload transport rates (Qs) for representative grain size of 0.125 mm (A), suspended sediment and sand loads, and total supply at Tarbert Landing [TBL, at river kilometer (rk) 490)], St Francisville (St F, rk 419), and Baton Rouge (BTR, rk 367.5) along the Lowermost Mississippi River during 2004-2010 (the matching period for these components at all three locations).X-axis distances are distance of the sites from the Head of the Passes (rk 0) near the LmMR's Gulf of Mexico outlet.Note: the suspended sand load is a coarser subset of suspended sediment load, hence, the total sediment supply at a given location along the river was calculated as the sum of total bedload and total suspended sediment load at the location.

Figure 19 .
Figure19.Difference between total bedload transport rates (Qs) for representative grain size of 0.125 mm (A), suspended sediment and sand loads, and total supply at Tarbert Landing [TBL, at river kilometer (rk) 490)], St Francisville (St F, rk 419), and Baton Rouge (BTR, rk 367.5) along the Lowermost Mississippi River during 2004-2010 (the matching period for these components at all three locations).X-axis distances are distance of the sites from the Head of the Passes (rk 0) near the LmMR's Gulf of Mexico outlet.Note: the suspended sand load is a coarser subset of suspended sediment load, hence, the total sediment supply at a given location along the river was calculated as the sum of total bedload and total suspended sediment load at the location.
If RE bed in 2004 <RE bed in 2013, RE bed was increasing.(b) If RE bed in 2004 >RE bed in 2013, RE bed was decreasing.RE bed in 2005 = RE bed in 2004 + ROI, RE bed in 2006 = RE bed in 2005 + ROI and so on), while the rate of decrease (ROD) in RE bed was subtracted consecutively in decreasing riverbeds (i.e., RE bed in 2005 = RE bed in 2004 − ROD, RE bed in 2006 = RE bed in 2005 − ROD and so on).

Table 1 .
Trend-line equations between river width (y: dependent) and river stage (x: independent) for Tarbert Landing (TBL), St Francisville (St F) and Baton Rouge (BTR) of the LmMR during 2004-2015.These equations were further used for estimating daily river width from daily river stages at the three locations.

Table 2 .
Discharge-suspended sediment load rating curves developed for St Francisville (St F) and Baton Rouge (BTR) of the LmMR.Here, ln (Q d ) is the independent variable and ln (DSL) is the dependent variable.

Table 3 .
Discharge -suspended sand load rating curves developed for St Francisville (St F) and Baton Rouge (BTR) of the LmMR.Here, ln (Q d ) is the independent variable and ln (DSL s ) is the dependent variable.

Table 4 .
Root mean square errors (RMSEs) (tonnes) of Daily Sediment Loads (DSL) predicted through discharge-sediment load rating curves for each period at St Francisville (St F) and Baton Rouge (BTR) as shown in Table1.Here, SE is the standard error and CF-Poly is the Duan correction factor used in polynomial rating curves, while CF-Lin is the Duan correction factor used in linear rating curves."No CF" represents DSLs calculated without applying correction factors during their retransformation from predicted ln DSLs while "CF" represents DSLs calculated by applying the correction factors during the retransformation procedure.

Table 5 .
Root mean square errors (RMSEs) (tonnes) of Daily Sand Loads (DSLs) predicted through discharge-sand load rating curves for each period at St Francisville (St F) and Baton Rouge (BTR) as shown in Table1.Specific terminologies in the table i.e., No CF, SE, CF-Poly, CF-Lin, are same as explained in Table4.

Table 5 .
Root mean square errors (RMSEs) (tonnes) of Daily Sand Loads (DSLs) predicted through discharge-sand load rating curves for each period at St Francisville (St F) and Baton Rouge (BTR) as shown in Table1.Specific terminologies in the table i.e., No CF, SE, CF-Poly, CF-Lin, are same as explained in Table4.