Assessment of Suspended Sand Availability under Different Flow Conditions of the Lowermost Mississippi River at Tarbert Landing during 1973–2013

Rapid land loss in the Mississippi River Delta Plain has led to intensive efforts by state and federal agencies for finding solutions in coastal land restoration in the past decade. One of the proposed solutions includes diversion of the Mississippi River water into drowning wetland areas. Although a few recent studies have investigated flow-sediment relationships in the Lowermost Mississippi River (LmMR, defined as the 500 km reach from the Old River Control Structure to the river's Gulf outlet), it is unclear how individual sediment fractions behave under varying flow conditions of the river. The information can be especially pertinent because the quantity of coarse sands plays a critical role for the Mississippi-Atchafalaya River deltaic development. In this study, we utilized long-term (1973–2013) records on discharge and sediments at Tarbert Landing of the LmMR to assess sand behavior and availability under different river flow regimes, and extreme sand transport events and their recurrence. We found an average annual sand load (SL) of 27.2 megatonnes (MT) during 1973 and 2013, varying largely from 3.37 to 52.30 MT. For the entire 41-year study period, a total of approximately 1115 MT sand were discharged at Tarbert Landing, half of which occurred during the peak 20% flow events. A combination of intermediate, high and peak flow stages (i.e., river discharge was ě18,000 cubic meter per second) produced about 71% of the total annual SL within approximately 120 days of a year. Based on the long-term sediment assessment, we predict that the LmMR has a high likelihood to transport 4 to 446 thousand tonnes of sand every day over the next 40 years, during which annual sand loads could reach a maximum of 51.68 MT. Currently, no effective plan is in place to utilize this considerably high sand quantity and we suggest that river engineering and sediment management in the LmMR consider practices of hydrograph-based approach for maximally capturing riverine sediments.


Introduction
River deltas comprise approximately five percent of the Earth's total land area and over 500 million people reside in them [1,2].They are important regions in both economic and environmental perspectives as they act as commercial centers and also provide a plethora of natural resources [3,4].However, many river deltas around the world are facing land loss as the consequence of human and natural factors including river engineering [5,6], accelerated subsidence [7,8], reduced riverine sediment supply [6,9,10], disconnection of the river with its floodplains [11], coastal land erosion [12], and relative sea level rise [13].
One renowned example of river deltas facing land loss problems is the Mississippi River Delta Plain (MRDP) in the southern USA.The MRDP has been losing a substantial amount of land since the early 20th century [6,7,14,15].A recent study on MRDP land loss [16] reported a disappearing rate of 43 km 2 /year since 1985.Potential of the MRDP land loss due to river engineering has long been recognized.More than a century ago, Corthell [17] already warned: "If certain levee structures were placed in a manner that fresh water and sediments, along with vital nutrients, were laid to waste off the mouth of the Mississippi River, their deltaic regenerative properties would be lost and unrecoverable."However, the land loss issue captured major public attention only from the late 1970s (e.g., [7,14,15,18]).During the last decade, in the wake of Hurricane Katrina, there has been an increasing concern over the issue in general public and scientific communities, which led to intensive efforts by the state and federal governments for finding solutions in offsetting coastal land loss [19,20].One such solution is diversion of the Mississippi River for outsourcing the river water and sediments to coastal marshes for stabilizing the deltaic system.The United States Army Corps of Engineers (USACE) has constructed three notable river diversions in the lowermost Mississippi River reach, namely, Caernarvon (river kilometer, or rk 131), Davis Pond (rk 190) and West Bay (rk 8), of which only West Bay focuses solely on sediment retention and capture [21].In addition, the State of Louisiana's Master Plan for coastal restoration (2012) proposed six large to small water (discharge from 141 to 7079 cubic meter per second, or cms) and sediment diversions [20,22], which are still in different planning phases.
Recent reports suggested that the three executed diversions have not gained significant steps towards their objectives despite careful planning and several years of operation [23,24].The Caenarvon and Davis Pond freshwater diversions have not induced significant salinity reduction [24,25] and have been subjected to more vegetation loss and nutrient overloading post hurricanes Katrina and Rita [24,26].Similarly, the planned discharge of the West Bay diversion has been increased from 396 cms at its inception in 2003 to 765 cms currently at the cost of adverse effects in navigational route; however, it has not produced desired land growth in the surrounding area [23,27,28].
A primary goal of the river diversion focusing entirely on sediment retention and capture is to divert flow carrying the maximum amount of sediments into adjoining drowned areas for delta restoration without hampering the ecological, structural, hydrological and functional integrity of the river at all [20,22,29].In the context of this goal, Rosen and Xu [21] analyzed the flow-sediment relationship for Tarbert Landing of the Lowermost Mississippi River (abbreviated hereafter as LmMR and defined as the 500-km reach from the Old River Control Structure to the river's Gulf outlet) to quantify sediment loads carried by varying flows during 1980-2010.They found that about half of the total annual sediment yield was produced within about 120 days every year when the river was at intermediate (18,000-25,000 cms) and high (25,000-32,000 cms) flows; hence, they recommended these flows to be diverted according to their natural cycle of occurrence during the year.Allison et al. [30] quantified short-term sediment budgets for different locations along the LmMR and suggested that years with high annual flow yielded high sediment input in the system and vice-versa, highlighting the advantage of sediment diversions during high flows.
Suspended solids in the Mississippi River have been found to be composed of a high proportion of fine clay/silt particles (<0.0625 mm) and a low proportion of coarser particles (>0.0625 mm) [30,31].Studies by Nepf [32] and Nittrouer et al. [33] postulated that sand may play a much more critical role in new land building in the Mississippi River delta than fine clay/silt.The importance of sand transport for the Mississippi River deltaic development has been increasingly recognized [33][34][35][36].Therefore, analysis of the relationship of total suspended sediment with river hydrology alone is not enough for developing effective sediment management plans.Although the sediment assessments by Rosen and Xu [21], Allison et al. [30] and Nittrouer and Viparelli [36] provide critical information for understanding sediment availability in the LmMR, they give little insights into the actual quantity of riverine sand under different flow regimes.The information about the actual sand availability can be crucial for developing management practices in maximizing sand capture in the LmMR.
The Tarbert Landing river gauge station (31 ˝00'30" N, 91 ˝7'25" W) is located at river kilometer 493 (river mile 306.3) of the Mississippi River.The station is below the Old River Control Structure (ORCS) [38] (Figure 1) that diverts approximately 25% (under the normal flow conditions) of the Mississippi River's water into the Atchafalaya River.The site provides the most updated and most comprehensive sediment records for the LmMR where both the United States Geological Survey (USGS) and United States Army Corps of Engineers (USACE) have a monitoring station (USGS Station ID: 07295100 and USACE Gauge ID: 01100).ORCS was built in 1963 with the primary goal of preventing a large amount of Mississippi River water (>30%) from entering the Atchafalaya River [39,40].Discharge at Tarbert Landing is, therefore, manipulated by ORCS based on specific river flow conditions.

Flow and Sediment Concentration Data
Records on mean daily discharge (Q d in cms) were collected for Tarbert Landing from USACE for the period from 1 January 1973 to 31 December 2013.For the same period, measurements on suspended sediment concentrations (SSC) in milligram per liter (mg/L) and corresponding percentage of silt/clay (fine sediment) fractions in SSC (i.e., diameter < 0.0625 mm) were collected from USGS.USGS carries out depth-integrated suspended sediment sampling every 12 to 26 days using several isokinetic point samplers (i.e., P-61, P-63, D-96, D-99) ranging from four to eight verticals and each vertical consisting of two to five samples ( [10,[41][42][43]-these studies also cover in depth analysis of SSC collection and processing techniques and error adjustments).From 1973 to 2013, a total of 1043 SSC samples were collected, processed and documented for Tarbert Landing.During Water 2015, 7, 7022-7044 these 41 years, each month had 1 to 3 sampling dates; therefore, it is assumed that the SSC data have a sufficient, unbiased representation across all seasons and flow regimes.The discharge and sediment concentration data were used to compute sand loads as described in section 2.5 below.

Flow and Sediment Concentration Data
Records on mean daily discharge (Qd in cms) were collected for Tarbert Landing from USACE for the period from 1 January 1973 to 31 December 2013.For the same period, measurements on suspended sediment concentrations (SSC) in milligram per liter (mg/L) and corresponding percentage of silt/clay (fine sediment) fractions in SSC (i.e., diameter < 0.0625 mm) were collected from USGS.USGS carries out depth-integrated suspended sediment sampling every 12 to 26 days using several isokinetic point samplers (i.e., P-61, P-63, D-96, D-99) ranging from four to eight verticals and each vertical consisting of two to five samples ( [10,[41][42][43]-these studies also cover in depth analysis of SSC collection and processing techniques and error adjustments).From 1973 to 2013, a total of 1043 SSC samples were collected, processed and documented for Tarbert Landing.During these 41 years, each month had 1 to 3 sampling dates; therefore, it is assumed that the SSC data have a sufficient, unbiased representation across all seasons and flow regimes.The discharge and sediment concentration data were used to compute sand loads as described in section 2.5 below.

NOAA's River Stages and their Corresponding Flow Regimes
For its flood warning prediction, NOAA defined five river stages at Red River Landing.In their study on long-term suspended sediment transport at Tarbert Landing, Rosen and Xu [21] identified corresponding discharge for these stages: discharge <13,000 cms for Low Flow Stage, 13,000-18,000 cms for Action Flow Stage, 18,000-25,000 cms for Intermediate Flow Stage, 25,000-32,000 cms for High Flow stage, and >32,000 cms for Peak Flow Stage.These regimes were further

NOAA's River Stages and their Corresponding Flow Regimes
For its flood warning prediction, NOAA defined five river stages at Red River Landing.In their study on long-term suspended sediment transport at Tarbert Landing, Rosen and Xu [21] identified corresponding discharge for these stages: discharge <13,000 cms for Low Flow Stage, 13,000-18,000 cms for Action Flow Stage, 18,000-25,000 cms for Intermediate Flow Stage, 25,000-32,000 cms for High Flow stage, and >32,000 cms for Peak Flow Stage.These regimes were further used several times in our study, e.g., in frequency analysis, duration curves, and sand distribution and transport trends across these regimes.

Sand Concentration in River Discharge
Using the percentage of silt/clay fractions in suspended sediment concentration, we first calculated silt/clay concentration (SSC f ) by multiplying the percentage with SSC.Sand concentration (SSC s in mg/L) for each sampling event was then quantified by subtracting the SSC f (mg/L) from SSC (mg/L).The distribution of SSC s across a daily river discharge (Q d ) range was then analyzed by building two types of SSC s (y-axis)-Q d (x-axis) plots: P-1 and P-2.Average SSC s and their percentage changes within pre-selected Q d intervals (every 3000 cms) were plotted against those Q d intervals in P-1.Similarly, individual SSC s were fitted against their corresponding Q d s in P-2.The upper limit of the Q d range in P-1 after which average SSC s began to decrease following a continuous increase Water 2015, 7, 7022-7044 (27,000 cms) gave the point for separating increasing and decreasing SSC s in P-2.Hence, all SSC s values for Q d ď 27,000 cms were defined as increasing sand, while all SSC s values for Q d > 27,000 cms were defined as decreasing sand in P-2.

Development of Discharge-Sand Load Rating Curves
Daily sand load (DSL in tonnes/day) was computed by multiplying SSC s with the corresponding daily discharge (Q d in cms) for all the sampling dates during 1 January 1973 and 31 December 2013 as: where 0.0864 is a unit conversion factor for converting the sand mass to tonnes per day.There were nine outliers out of a total of 1043 sediment sampling events for the entire 41-year study period.These ~1% outliers were four to six times higher than the long-term standard deviation of sand concentrations; hence, we decided to remove them from further analysis.Now, a natural logarithm (ln) was taken for the two variables, DSL (dependent; y) and Q d (independent; x), and both linear and polynomial rating curves were applied for the relation between them.The evaluation of all applied rating curves were based on four criteria: regression coefficient of the curves (R 2 must be ě0.8),root mean square errors of the predicted (or calculated) DSLs (RMSE) (the lower the better), standard error (SE) of the rating curves (in ln units) (also, the lower the better) and a graphical assessment (good visual agreement between corresponding calibrated and predicted DSLs) [44][45][46].
To achieve the "predicted DSLs", we fitted "log transformed (ln) Q d s" in the rating curve equations to get "predicted ln DSLs" at first and then transformed back thus obtained "predicted ln DSLs" by taking their exponential values.We also checked potential log-biasing in this retransformation procedure using the following correction factor (CF) given by Duan [47] and modified by Gray et al. [48] because, firstly, it does not require normality of residuals and, secondly, residuals for a few rating curves in our analyses were not normally distributed (p-values < 0.05 in Shapiro-Wilk tests).
CF " where e i is the difference between i th observations of "measured log DSLs" and "predicted log DSLs" and n is the total number of samples used in the given rating curve.Single linear and polynomial rating curves were applied for the whole period at first, however, all four criteria to evaluate rating curves approach were not met here: lower R 2 (0.69 for linear and 0.7 for polynomial rating curve) (Table 1), comparatively higher RMSE (71,067 for linear and 67,950 for polynomial rating curve) (Table 2), comparatively higher SE (0.823) (Table 2), and poor visual agreement between corresponding measured and predicted DSLs (Figure 2).Low sample size of sediment concentrations during each year stopped us from applying sand rating curves annually during 1973-2013.In addition, rating curves in decadal intervals can minimize year to year variability in sediment samples and give robust average-annual predictions over decadal periods (supplementary information in [36]).Therefore, linear and polynomial rating curves were further applied for the following approximately decadal intervals in continuum: 1973-1985 (n = 463), 1986-1995 (n = 242), 1996-2005 (n = 187), and 2006-2013 (n = 142).The prerequisite of R 2 ě 0.8 was met for three of the four periods (1973-85: linear R 2 = 0.8, polynomial R 2 = 0.84; 1996-2005: linear R 2 = 0.81, polynomial R 2 = 0.83; 2006-13: linear R 2 = 0.82, polynomial R 2 = 0.87) (Table 1), so corresponding rating curves for these periods were subjected to further evaluation using other three criteria.However, the period 1986-95 had R 2 s (R-squares) < 0.8 (0.57) for both rating curves (Table 1).Thus, each year was checked with annual linear and polynomial rating curves to find the years responsible for lowering the combined linear and polynomial R 2 s in this period (Table A1 in Appendix).We found all R 2 s during 1986-1990 (0.15-0.51) in one cluster, substantially lower than Water 2015, 7, 7022-7044 all R 2 s during 1991-1995 (0.69-0.92) in another cluster.Hence, based on approximation of individual R 2 s of annual rating curves, we combined the two periods 1986-90 (n = 118) and 1991-95 (n = 124) for further evaluation of their corresponding rating curves.Finally, we found that polynomial discharge-sand load rating curves during the four durations: 1973-1985, 1991-1995, 1996-2005 and 2006-2013 met all the four criteria and provided DSL estimates most approximate to the measured DSLs (Tables 1 and 2; Figure 3).We also found that the use of correction factors overestimated DSLs slightly (for polynomial curves) as well as substantially (for linear curves) as compared to their corresponding calibrated measurements (Table 2; Figures 2 and 3).Hence, based on evaluation of these overestimations and previous arguments regarding unreliability of the correction factors [49,50], we decided to use polynomial sand rating curves categorized into aforementioned four periods without correction factor to calculate sand loads for each day from 1973 to 2013 except for the period 1986-1990.The reason for excluding rating curve analysis from 1986 to 1990 and the procedure followed to calculate daily sand loads during this period have been explained in Section 2.6 further down.Finally, we found that polynomial discharge-sand load rating curves during the four durations: 1973-1985, 1991-1995, 1996-2005 and 2006-2013 met all the four criteria and provided DSL estimates most approximate to the measured DSLs (Tables 1 and 2; Figure 3).We also found that the use of correction factors overestimated DSLs slightly (for polynomial curves) as well as substantially (for linear curves) as compared to their corresponding calibrated measurements (Table 2; Figures 2 and 3).Hence, based on evaluation of these overestimations and previous arguments regarding unreliability of the correction factors [49,50], we decided to use polynomial sand rating curves categorized into aforementioned four periods without correction factor to calculate sand loads for each day from 1973 to 2013 except for the period 1986-1990.The reason for excluding rating curve analysis from 1986 to 1990 and the procedure followed to calculate daily sand loads during this period have been explained in Section 2.6 further down.2. It is noted that both predicted and measured SLs during the period 1986-1990 were eliminated in this comparison because of the low R 2 value of both rating curves during this period (please see Table 1).

Non-Rating Curve Approach for Sand Load Calculation
For the period of 1986-1990, three of the four criteria to evaluate rating curves approach were not met: lower R 2 (0.36 for both polynomial and linear rating curves) (Table 1), comparatively higher SE (1.132) (Table 2) and poor visual agreement between corresponding measured and predicted DSLs (Figure 4).Therefore, calibrated sand concentration measurements (117 samples) and daily discharge were used to calculate sand loads for each day during this period.Here, starting from 1986, the earliest available sand concentration of the year was assumed to be equal to all consecutive days of missing concentration until the next value was available.In addition, last available concentration of the earlier year was used for filling values of missing days of the current year if the earliest concentrations did not start from the first day of the year.Finally, DSLs for non-  2. It is noted that both predicted and measured SLs during the period 1986-1990 were eliminated in this comparison because of the low R 2 value of both rating curves during this period (please see Table 1).

Non-Rating Curve Approach for Sand Load Calculation
For the period of 1986-1990, three of the four criteria to evaluate rating curves approach were not met: lower R 2 (0.36 for both polynomial and linear rating curves) (Table 1), comparatively higher SE (1.132) (Table 2) and poor visual agreement between corresponding measured and predicted DSLs (Figure 4).Therefore, calibrated sand concentration measurements (117 samples) and daily discharge were used to calculate sand loads for each day during this period.Here, starting from 1986, the earliest available sand concentration of the year was assumed to be equal to all consecutive days of missing concentration until the next value was available.In addition, last available concentration of the earlier year was used for filling values of missing days of the current year if the earliest concentrations did not start from the first day of the year.Finally, DSLs for non-sand rating curve years were calculated using the formula in Equation ( 1).

Range of Errors Associated with Predicted Sand Loads
We considered two types of errors (E-1 and E-2) in the SL estimates (it must be noted that the standard errors discussed earlier in Section 2.5 accounted for the entire models rather than individual estimates).E-1 is associated with the methods used by USGS for depth-integrated sampling and calibration of SSCs.It has previously been reported to be approximately ±10% of the total calibrated SSCs, SSCss and SSCfs [51][52][53].E-2 is associated with dependent variables (ln SL) in the rating curves.The confidence intervals (CI) for each "ln predicted SL" at 95% level of significance in their rating curves were provided with the help of their corresponding E-2s in our analysis (Figure 5).We estimated an approximate E-2 of ±15% in all SLs predicted from rating curve approach (based on confidence interval plots, RMSEs and percentage difference between measured and predicted SLs which averaged −13.4% during the four periods).Thus, the total error in SL measurements and predictions (E-1 + E-2) during rating curve years was about ±25%.We only selected E-1 for all estimates during 1986-1990 because we did not use rating curve approach in this period.Therefore, error range for SL estimates during 1986-1990 was ±10%.For convenience and consistency in reporting, we used an error range of ±18% for all the SL estimates during 1973-2013 (approximately ~average of 25 and 10).

Range of Errors Associated with Predicted Sand Loads
We considered two types of errors (E-1 and E-2) in the SL estimates (it must be noted that the standard errors discussed earlier in Section 2.5 accounted for the entire models rather than individual estimates).E-1 is associated with the methods used by USGS for depth-integrated sampling and calibration of SSCs.It has previously been reported to be approximately ˘10% of the total calibrated SSCs, SSC s s and SSC f s [51][52][53].E-2 is associated with dependent variables (ln SL) in the rating curves.The confidence intervals (CI) for each "ln predicted SL" at 95% level of significance in their rating curves were provided with the help of their corresponding E-2s in our analysis (Figure 5).We estimated an approximate E-2 of ˘15% in all SLs predicted from rating curve approach (based on confidence interval plots, RMSEs and percentage difference between measured and predicted SLs which averaged ´13.4% during the four periods).Thus, the total error in SL measurements and predictions (E-1 + E-2) during rating curve years was about ˘25%.We only selected E-1 for all Water 2015, 7, 7022-7044 estimates during 1986-1990 because we did not use rating curve approach in this period.Therefore, error range for SL estimates during 1986-1990 was ˘10%.For convenience and consistency in reporting, we used an error range of ˘18% for all the SL estimates during 1973-2013 (approximately ~average of 25 and 10).

Daily, Annual and Seasonal Sand Load Trends
DSLs were calculated using aforementioned methods of rating curve (1973-1985, 1991-1995, 1996-2005, 2006-2013) and non-rating curve (1986-1990) approaches.The sum of DSLs from 1st January to 31st December during each year gave their corresponding annual SLs.Maximum, minimum and average DSLs and annual SLs were plotted against their corresponding years for information regarding daily and annual SL trends throughout the study period.Similarly, monthly SLs were calculated by averaging DSLs for each month separately from 1973 to 2013.Maximum, minimum and average monthly SLs were plotted against their corresponding months to analyze their seasonal trends.

Frequency Analysis of Sand Loads
In this study, we analyzed the amount of sand transported at Tarbert Landing during 1973-2013 under different river flow conditions, i.e., six frequencies on the flow duration curves (1%, 5%, 10%, 20%, 50%, and 75%) and five river stages (Low, Action, Intermediate, High, and Peak).The Gumbel distribution [54,55] was used for analyzing annual maximum and minimum DSLs, while the Weibull distribution [56,57] was used for analyzing total annual sand loads (SL) at Tarbert Landing.All annual maximum/minimum DSLs and total annual SLs during 1973 and 2013 were sorted in descending order separately at first.The non-exceedance probabilities {F(X)} for maximum and minimum DSLs were obtained with the Gumbel distribution (Equation (3)), while the nonexceedance probabilities for total annual SLs were obtained with the Weibull distribution (Equation ( 4)) as given below:

Daily, Annual and Seasonal Sand Load Trends
DSLs were calculated using aforementioned methods of rating curve (1973-1985, 1991-1995, 1996-2005, 2006-2013) and non-rating curve (1986)(1987)(1988)(1989)(1990) approaches.The sum of DSLs from 1st January to 31st December during each year gave their corresponding annual SLs.Maximum, minimum and average DSLs and annual SLs were plotted against their corresponding years for information regarding daily and annual SL trends throughout the study period.Similarly, monthly SLs were calculated by averaging DSLs for each month separately from 1973 to 2013.Maximum, minimum and average monthly SLs were plotted against their corresponding months to analyze their seasonal trends.

Frequency Analysis of Sand Loads
In this study, we analyzed the amount of sand transported at Tarbert Landing during 1973-2013 under different river flow conditions, i.e., six frequencies on the flow duration curves (1%, 5%, 10%, 20%, 50%, and 75%) and five river stages (Low, Action, Intermediate, High, and Peak).The Gumbel distribution [54,55] was used for analyzing annual maximum and minimum DSLs, while the Weibull distribution [56,57] was used for analyzing total annual sand loads (SL) at Tarbert Landing.All annual maximum/minimum DSLs and total annual SLs during 1973 and 2013 were sorted in descending order separately at first.The non-exceedance probabilities {F(X)} for maximum and minimum DSLs Water 2015, 7, 7022-7044 were obtained with the Gumbel distribution (Equation (3)), while the non-exceedance probabilities for total annual SLs were obtained with the Weibull distribution (Equation ( 4)) as given below: where X is annual maximum/minimum DSL (tonnes/day) or total annual SL (megatonnes), m is the rank of the annual SL, n is the total number of years in the distribution, and a and b are the Gumbel distribution parameters that were calculated through: where µ x is the average and S x is the standard deviation of the annual maximum and minimum DSLs.
Maximum and minimum DSLs (Q p ) for the return periods [T(X)] of 2-, 5-, 10-, 20-, and 40-years were calculated using the Gumbel distribution as: where the frequency factor K(T) is defined as: A frequency factor is computed for a certain return.In this study, we computed frequency factors of ´0.1643, 0.7195, 1.3046, 1.8658, and 2.4163 for the 2-, 5-, 10-, 20-, and 40-year return periods, respectively.Annual SLs for the same return periods were estimated using a linear interpolation from the Weibull distribution of annual SLs (i.e., 1/{1-F(X)}).
Seasonally, average Q d increased continuously from January to its maximum in April (16,550 to 22,468 cms), then decreased continuously from May to its minimum in October (21,696 to 8171 cms) inferring maximum river flow during the spring months (March, April and May) (Figure 8).For the remaining two months in the year, average Q d followed an increasing trend again (9702 cms in Water 2015, 7, 7022-7044 November and 14,581 cms in December) (Figure 8).In addition, the maximum Q d was observed in May (45,844 cms) while the minimum Q d was observed in July (3143 cms) (Figure 8).

Sand Concentrations under Different Flow Regimes
Average SSCss and their percentage changes at Tarbert Landing showed early increasing trend from the lowest Qd interval (3000-6000 cms) (Figure 9).Average SSCs for the lowest Qd interval (3000-6000 cms) was about 11 mg/L which increased up to about 91 mg/L between 24,000 and 27,000 cms Qd (715% increase) (Figure 9).Further, average SSCss fluctuated for higher Qd intervals, i.e., they first decreased up to approximately 70 mg/L for Qd interval between 30,000 and 36,000 cms (715% to 529%), then increased up to approximately 103 mg/L for the next interval between 36,000 and 39,000 cms (529% to 832%), then decreased up to about 72 mg/L (832% to 552%) for the next interval between 39,000 and 42,000 cms and finally increased up to about 82 mg/L (552% to 642%) for the highest Qd interval (42,000-45,000 cms) (Figure 9).Similar trends were also observed in individual SSCs values across the entire Qd range.SSCss showed an early increasing trend even with lower discharge levels (about 6000 cms) which continued until substantially high flows (27,000 cms with R 2 = 0.36) (Figure 10).The elevated concentrations remained almost constant (R 2 = 0.0061) for higher flows (>27,000 cms) (Figure 10).

Sand Concentrations under Different Flow Regimes
Average SSC s s and their percentage changes at Tarbert Landing showed early increasing trend from the lowest Q d interval (3000-6000 cms) (Figure 9).Average SSC s for the lowest Q d interval (3000-6000 cms) was about 11 mg/L which increased up to about 91 mg/L between 24,000 and 27,000 cms Q d (715% increase) (Figure 9).Further, average SSC s s fluctuated for higher Q d intervals, i.e., they first decreased up to approximately 70 mg/L for Q d interval between 30,000 and 36,000 cms (715% to 529%), then increased up to approximately 103 mg/L for the next interval between 36,000 and 39,000 cms (529% to 832%), then decreased up to about 72 mg/L (832% to 552%) for the next interval between 39,000 and 42,000 cms and finally increased up to about 82 mg/L (552% to 642%) for the highest Q d interval (42,000-45,000 cms) (Figure 9).Similar trends were also observed in individual SSC s values across the entire Q d range.SSC s s showed an early increasing trend even with lower discharge levels (about 6000 cms) which continued until substantially high flows (27,000 cms with R 2 = 0.36) (Figure 10).The elevated concentrations remained almost constant (R 2 = 0.0061) for higher flows (>27,000 cms) (Figure 10).
Water 2015, 7, 7022-7044 As with average Q d , years with higher average DSLs had higher minimum and maximum DSLs as compared to years with lower average DSLs (Figure 11).
Annual sand load from 1973 to 2013 averaged 27.2 MT, ranging from 3.37 to 52.30 MT and producing a total sand amount of 1114.82MT for the entire 41-year study period (Figure 11).Annual SL was lower than 10 MT for four of 41 years (1986, 1987, 1988 and 1989), higher than 40 MT for eight years including the Mississippi flood years of 1973, 1993 and 2011, and either > 10 MT or < 40 MT for the remaining 29 years (Figure 11).In addition, annual SL averaged 26.3 MT from 1973 to 1999 (28.8 MT from 1972to 1979, 19.3 MT from 1980to 1989, 32 MT from 1990 to 1999), which later increased to approximately 29 MT from 2000 to 2013 (calculations from table of Figure 11).Despite the low average annual SL between 1986 and 1989, there was no continuous increasing or decreasing trend (even for 2/3 years) in annual SL (Figure 11).
Annual sand load from 1973 to 2013 averaged 27.2 MT, ranging from 3.37 to 52.30 MT and producing a total sand amount of 1114.82MT for the entire 41-year study period (Figure 11).Annual SL was lower than 10 MT for four of 41 years (1986, 1987, 1988 and 1989) 11).Despite the low average annual SL between 1986 and 1989, there was no continuous increasing or decreasing trend (even for 2/3 years) in annual SL (Figure 11).Seasonal trend of the average DSL was similar to that of the average Qd.Average DSL increased continuously from January to its maximum in May (i.e., 86,315 to 137,387 tons/day), then decreased from June to its minimum in September (106,426 to 14,935 tons/day) inferring maximum sand transport during spring (March, April and May) (Figure 12).For the remaining three months in the year, average DSL followed an increasing trend again, i.e., 20,871 tons/day in October to 70,092 tons/day in December (Figure 12).In addition, maximum DSL was observed in June (444,626 tons/day) while minimum DSL was observed in August (258 tons/day) for the whole period (Figure 12).Seasonal trend of the average DSL was similar to that of the average Q d .Average DSL increased continuously from January to its maximum in May (i.e., 86,315 to 137,387 tons/day), then decreased from June to its minimum in September (106,426 to 14,935 tons/day) inferring maximum sand transport during spring (March, April and May) (Figure 12).For the remaining three months in the year, average DSL followed an increasing trend again, i.e., 20,871 tons/day in October to 70,092 tons/day in December (Figure 12).In addition, maximum DSL was observed in June (444,626 tons/day) while minimum DSL was observed in August (258 tons/day) for the whole period (Figure 12).

Sand Load Distribution with River Discharge
Hydrologically, with respect to the NOAA's river stages, Intermediate, High and Peak Flow Stages together carried majority of the total SL (793.4MT; 71%) within an average of 122 days per year (Table 4).Individually, Intermediate, High and Peak Flow Stages carried approximately 384 (34%), 266 (24%) and 143 (13%) MT of SLs for average durations of 76, 34 and 12 days/year respectively, while, Low and Action Flow Stages carried approximately 146 (13%) and 175 (16%) MT of SL respectively for relatively longer average durations of 181 and 61 days/year (Table 4).

Sand Load Distribution with River Discharge
Hydrologically, with respect to the NOAA's river stages, Intermediate, High and Peak Flow Stages together carried majority of the total SL (793.4MT; 71%) within an average of 122 days per year (Table 4).Individually, Intermediate, High and Peak Flow Stages carried approximately 384 (34%), 266 (24%) and 143 (13%) MT of SLs for average durations of 76, 34 and 12 days/year respectively, while, Low and Action Flow Stages carried approximately 146 (13%) and 175 (16%) MT of SL respectively for relatively longer average durations of 181 and 61 days/year (Table 4).Also, 50 and 75% of flow regimes produced majority of total SL in the area, i.e., 966 (87%) and 1082 MT (97%) respectively (Table 5).Similarly, 1%, 5%, 10% and 20% of flow regimes also produced few to slightly more than half of total SL, i.e., 46 (4%), 196 (18%), 340 (31%) and 571 MT (51%), respectively (Table 5).Linear relationship between all annual flow volumes (km 3 ) above the low and action flow stages (discharge ě 18,000 cms) and corresponding annual sand loads (MT) showed that intermediate, high and peak flow stages jointly accounted for about 66% variation in sand loads during each year from 1973 to 2013 (R 2 = 0.66) (Figure 13).The relationship further suggested that about 66% of sand loads were produced by these three stages jointly during each year of the study period.
Linear relationship between all annual flow volumes (km 3 ) above the low and action flow stages (discharge ≥ 18,000 cms) and corresponding annual sand loads (MT) showed that intermediate, high and peak flow stages jointly accounted for about 66% variation in sand loads during each year from 1973 to 2013 (R 2 = 0.66) (Figure 13).The relationship further suggested that about 66% of sand loads were produced by these three stages jointly during each year of the study period.

Maximum and Minimum Sand Loads for Different Return Periods
Averages of maximum and minimum DSLs for each year throughout the whole study period at Tarbert Landing were 226,981 and 4865 tonnes/day, while their standard deviations were 90,595 and 4015 tonnes/day respectively (Table 6).Based on these means and standard deviations, parameters a and b in Gumbel distribution were found to be 186,209 and 70,637 for highest DSLs and 3058 and 3131 for lowest DSLs during each year respectively (Table 6).Non-exceedance probabilities of maximum and minimum DSLs {Gumbel F(X)} and total SLs {Weibull F(X)} for each year from 1973 to 2013 are represented by Figure 14.For longer return periods of 20 and 40 years, we predicted that DSLs at Tarbert Landing can reach a maximum of 396 and 446 thousand tons and a minimum of 12 and 15 thousand tons, respectively (Table 7).Similarly, for shorter return periods of 2, 5, and 10 years, maximum DSL predictions were 212, 292, and 345 thousand tons respectively, while minimum DSL predictions were 4, 8, and 10 thousand tons,

Maximum and Minimum Sand Loads for Different Return Periods
Averages of maximum and minimum DSLs for each year throughout the whole study period at Tarbert Landing were 226,981 and 4865 tonnes/day, while their standard deviations were 90,595 and 4015 tonnes/day respectively (Table 6).Based on these means and standard deviations, parameters a and b in Gumbel distribution were found to be 186,209 and 70,637 for highest DSLs and 3058 and 3131 for lowest DSLs during each year respectively (Table 6).Non-exceedance probabilities of maximum and minimum DSLs {Gumbel F(X)} and total SLs {Weibull F(X)} for each year from 1973 to 2013 are represented by Figure 14.For longer return periods of 20 and 40 years, we predicted that DSLs at Tarbert Landing can reach a maximum of 396 and 446 thousand tons and a minimum of 12 and 15 thousand tons, respectively (Table 7).Similarly, for shorter return periods of 2, 5, and 10 years, maximum DSL predictions were 212, 292, and 345 thousand tons respectively, while minimum DSL predictions were 4, 8, and 10 thousand tons, respectively (Table 7).We also predicted that in the next 2-, 5-, 10-, 20-, and 40-years total annual sand load can reach as much as 28.2, 40.28, 44.78, 45.76 and 51.68 MT, respectively (Table 7).7).We also predicted that in the next 2-, 5-, 10-, 20-, and 40-years total annual sand load can reach as much as 28.2, 40.28, 44.78, 45.76 and 51.68 MT, respectively (Table 7).

Long-Term Trend of Sand Loads in the LmMR
In this study we present sand transport estimates for each individual year along with their error range (±18%), providing a comprehensive range of sand transport in the past four decades.Previously, Nittrouer and Viraparelli [36] reported an average annual sand transport of about 24 MT for Tarbert Landing during 1973 and 2012, which is about 12% lower than our estimate for average annual sand transport (27 ± 4.9 MT) during 1973 and 2013.This difference falls within the total error range of all SL estimates in this study (±18%).The relatively small difference between the two estimates may be caused by different estimation approaches: Nittrouer and Viraparelli used linear rating curves at decadal intervals for daily sand load estimation while this study used a

Long-Term Trend of Sand Loads in the LmMR
In this study we present sand transport estimates for each individual year along with their error range (˘18%), providing a comprehensive range of sand transport in the past four decades.Previously, Nittrouer and Viraparelli [36] reported an average annual sand transport of about 24 MT for Tarbert Landing during 1973 and 2012, which is about 12% lower than our estimate for average annual sand transport (27 ˘4.9 MT) during 1973 and 2013.This difference falls within the total error range of all SL estimates in this study (˘18%).The relatively small difference between the two estimates may be caused by different estimation approaches: Nittrouer and Viraparelli used linear rating curves at decadal intervals for daily sand load estimation while this study used a combination of five-yearly to approximately decadal polynomial rating curves and monthly sand concentration records.Our SL estimate for 2013 (31 ˘5.6 MT), which was not included in Nittrouer and Viraparelli's Water 2015, 7, 7022-7044 analysis, was about 14% higher than our 41-year long average SL estimate; however, this difference also falls within our error range (˘18%).
In their short-term sediment budget study, Allison et al. [30] reported an annual SL of 73.5, 62.2 and 78.9 MT for the water years of 2008, 2009 and 2010, respectively, resulting in a total SL of 214.6 MT.These estimates are 62%, 36% and 108% higher than our calendar-based estimates for these three years (45.5 ˘8.1, 45.8 ˘8.2 and 38.0 ˘6.8 MT, respectively), or nearly doubled of our estimate for the entire three years (129.3˘23.3 MT).These two sets of estimates cannot be compared directly because of the different time base; however, it seems that Allison et al. [30] [58].In addition, the rating curve which Allison et al. [30] used to calculate DSLs for each day during the three years had a comparatively lower R 2 (0.62) (information from the supplementary files of Allison et al. [30]) and no other criteria for model validation which could have resulted in greater variability between calibrated and estimated SLs.Furthermore, their estimates do not fall within the error range of our annual estimates during these three years.These three arguments provide essential support for questioning the reliability of their estimates.
In our study, we found a very low sand transport during 1986 and 1989.Only 19 ˘1.9 MT of sand was discharged during this four-year period, making an average annual SL of just 5 MT, i.e., less than one fifth of the long-term average annual SL of 27 ˘4.9MT.Nittrouer and Viraparelli [36] also reported low average SL for a period longer than these four years (1980 to 1989) (about 11 MT); however, they did not separate SLs for each year during the four-year period.The notable abrupt drop in sand transport during these four years may have been a result from the four-year severe drought in the Southeastern and Midwestern United States between 1986 and 1989, as reported by Cook et al. [59] and Trenberth and Guillemot [60].During this drought period, both discharge and sediment concentrations of the LmMR were considerably lower than those of other years in the past four decades.
The declining trend of total suspended sediment load input in the LmMR during the last several decades has been well documented [6,[61][62][63].However, Rosen and Xu [21] contradicted this trend by suggesting that suspended sediment input has slightly increased in the river for the recent two decades (from 1990 to 2010), although no statistical significance was found.Our findings for sand transport (the courser sediment fraction) are identical with Rosen and Xu [21] as we also found a stronger increasing trend in sand loads starting from 1990-the average annual sand load from 1990 to 2013 (30 ˘5.4 MT) was clearly higher than that during 1973 and 1989 (23 ˘4.1 MT).The higher sand transport in the past two decades has mainly resulted from the increased discharge during the same period of time.

Hydrologic Control for Sand Transport in the LmMR
Our findings suggest that notable sand load present in the LmMR can best be diverted during its intermediate, high and peak flow stages (discharge ě 18,000 cms) within only approximately four months each year.We found that maximum river flow and sand transport both occurred during the spring months (March, April and May); therefore, it is highly likely that the three river stages are prevalent during spring and scarce during other seasons of the year.Highest sand transport by these three stages can be linked to their rapid increase in discharge regimes from nearly 6000 to 27,000 cms and very slow and inconsistent decrease in regimes post 27,000 cms.Sand concentrations seem to have reached their peaks during the intermediate and early high flow stages (discharge: 18,000-27,000 cms), hence resulting in the highest sand loads.
There is no previous work available for comparison to these findings of long term sand load availability with river discharge.The work by Rosen and Xu [21] seems to support our findings, but they analyzed the availability of total suspended sediments with river discharge as compared to Water 2015, 7, 7022-7044 the total sand load in our study and emphasized that intermediate and high river stages combined carry highest sediment loads, however, peak stage contributes relatively little in sand transport.Biedenharn and Thorne [64] argued that discharge between 17,000 and 40,000 cms transport more suspended sediments than other discharge regimes.Differences in weights and volumes of sediment and sand concentrations may be the reason behind subtle differences in flow regimes reported to carry highest sediment and sand loads in these studies.We also found that almost the entire total sand load (97%) in the LmMR was transported by 75% of total water discharge throughout the study period.This discharge regime includes almost half of the lower flow stage along with all other Mississippi River stages (>8325 cms) at Tarbert Landing.The sand behaviors found in this study-(1) increasing rapidly with increasing discharge; however, decreasing slowly beyond a given discharge regime just after reaching its climax; and (2) maximum sand percentage transported by substantially less flow volume percentage-can be compared with other large river systems in the world with a sinking delta.Such information can help planning for deltaic land protection through effective sand management.

Future Likelihood of Sand Transport in the LmMR
Our frequency analysis reveals that the LmMR at Tarbert Landing has the potential to transport substantial amount of sand every day in the next 40 years (4 to 446 thousand tonnes).Based on the sand yields at different river stages, we argue that the intermediate, high and/or peak river stages will possibly transport higher DSLs in shorter periods while low and/or action stages will possibly transport lower DSLs in longer periods of time.Our findings further indicated that annual sand load has the potential to reach as much as 52 MT in the next 40 years.Based on our observations regarding the linear relationship between annual flow volumes in intermediate, high and peak flow stages and annual SLs, we also argue that years with high average daily discharge and annual flow would produce high annual sand loads within the given return periods and vice-versa.Previous studies have analyzed several year peaks for suspended sediment loads and concentrations [65,66] and even for river flows and stages [67,68] in different rivers around the world.However, to the best of our knowledge, maximum and minimum DSLs and annual SLs for short-and long-term return periods have not been analyzed for any river location to date.Thus, we could not compare these sand estimates in the LmMR with any other study.
The peak high and low DSLs and peak annual SLs vary between Tarbert Landing and other sites in the LmMR and the variation is based on sand-flow relationship and sand percentage in sediment load.Quantification of peak SLs at other sites is beyond the scope of this study.However, the analysis of daily and annual sand loads for several return periods can be helpful in speculating the importance of sediment diversion as per sand amount present at the site.It is also possible to incorporate the maximum/minimum DSLs and annual SLs for several return periods into any proposed land loss model for the MRDP.This can be done by quantifying the percentage loss of sand when given land area (km 2 ) was lost in "n" (where n = 2, 5, 10, 20, or 40) years and/or the amount of sand required to attain a goal of restoring certain land (km 2 ) in "n" years.

Conclusions
This study is the first comprehensive analysis of four-decade long sand transport under different flow conditions in the Lowermost Mississippi River.Our findings show that the majority of sand at Tarbert Landing are transported during the intermediate, high and peak river flow stages, and that their most effective capture can be achieved within 120 days of a year when discharge is greater than 18,000 cms.We also predict that the LmMR will most likely transport 4 to 446 thousand tons of sand every day over the next 40 years, during which annual sand load can reach as high as 52 million tons.Such considerably high sand loads are a precious resource for coastal Louisiana and should be effectively captured for offsetting land loss in the Mississippi River Delta before they are lost to deep waters of the Gulf of Mexico.To achieve this goal, river engineering and sediment Water 2015, 7, 7022-7044 management should consider applications using the hydrograph-based sand availability approach for maximum sediment capture.This may have implications for impeding coastal land loss in other sediment-starving deltas in the world.

Figure 1 .
Figure 1.Study area map (modified from Rosen and Xu [21]) showing the location of the LmMR at Tarbert Landing (TBL) (USGS Station ID 07295100 and USACE Gage ID 01100).MR and AR denote the courses of the Mississippi and Atchafalaya Rivers respectively; ORCS is the Old River Control Structure; "Sim" denotes Simmesport (USGS Station ID 07381490) site of the Atchafalaya River; RRL is Red River Landing, the gauging station for USGS just below TBL consisting of river stage records and CAR denotes Carrolton, New Orleans.The three Mississippi River diversions introduced earlier have also been shown in the figure: Davis Pond Freshwater Diversion (DPFD), Caernarvon Freshwater Diversion (CFD) and West Bay Sediment Diversion (WBSD).

Figure 1 .
Figure 1.Study area map (modified from Rosen and Xu [21]) showing the location of the LmMR at Tarbert Landing (TBL) (USGS Station ID 07295100 and USACE Gage ID 01100).MR and AR denote the courses of the Mississippi and Atchafalaya Rivers respectively; ORCS is the Old River Control Structure; "Sim" denotes Simmesport (USGS Station ID 07381490) site of the Atchafalaya River; RRL is Red River Landing, the gauging station for USGS just below TBL consisting of river stage records and CAR denotes Carrolton, New Orleans.The three Mississippi River diversions introduced earlier have also been shown in the figure: Davis Pond Freshwater Diversion (DPFD), Caernarvon Freshwater Diversion (CFD) and West Bay Sediment Diversion (WBSD).

Figure 2 .
Figure 2. Scatter plots showing comparison between sand loads calculated from sand concentrations measured, processed and calibrated by USGS (Measured SL) and those predicted from single sandrating curve (either linear or polynomial) (Predicted SL) at Tarbert Landing from 1973 to 2013.Here, linear rating curves were used for predicting SLs in (a) and (c)while polynomial rating curves were used for predicting SLs in (b) and (d).In addition, Duan correction factors were applied in predicted SLs of curves (c) and (d) (denoted by "CF" in the figure) while the SLs in curves (a) and (b) were predicted without correction factors (denoted by "No CF" in the figure).

Figure 2 .
Figure 2. Scatter plots showing comparison between sand loads calculated from sand concentrations measured, processed and calibrated by USGS (Measured SL) and those predicted from single sand-rating curve (either linear or polynomial) (Predicted SL) at Tarbert Landing from 1973 to 2013.Here, linear rating curves were used for predicting SLs in (a) and (c)while polynomial rating curves were used for predicting SLs in (b) and (d).In addition, Duan correction factors were applied in predicted SLs of curves (c) and (d) (denoted by "CF" in the figure) while the SLs in curves (a) and (b) were predicted without correction factors (denoted by "No CF" in the figure).

Figure 3 .
Figure 3. Scatter plots showing comparison between sand loads calculated from sand concentrations measured, processed and calibrated by USGS (Measured SL) and those predicted from several sandrating curves (Predicted SL) at Tarbert Landing from 1973 to 2013.Specific terminologies pertaining to parts (a), (b), (c) and (d) of this figure i.e., Linear, Poly, CF, and No CF are same as explained in Figure2.It is noted that both predicted and measured SLs during the period 1986-1990 were eliminated in this comparison because of the low R 2 value of both rating curves during this period (please see Table1).

Figure 3 .
Figure 3. Scatter plots showing comparison between sand loads calculated from sand concentrations measured, processed and calibrated by USGS (Measured SL) and those predicted from several sand-rating curves (Predicted SL) at Tarbert Landing from 1973 to 2013.Specific terminologies pertaining to parts (a), (b), (c) and (d) of this figure i.e., Linear, Poly, CF, and No CF are same as explained in Figure2.It is noted that both predicted and measured SLs during the period 1986-1990 were eliminated in this comparison because of the low R 2 value of both rating curves during this period (please see Table1).

Figure 4 .
Figure 4. Scatter plots showing comparison between sand loads calculated from sand concentrations measured, processed and calibrated by USGS (Measured SL) and those predicted from single sand-rating curve (either linear or polynomial) (Predicted SL) at Tarbert Landing from 1986 to 1990.Specific terminologies pertaining to parts (a), (b), (c) and (d) of this figure i.e., Linear, Poly, CF, and No CF are same as explained in Figure 2.

Figure 4 .
Figure 4. Scatter plots showing comparison between sand loads calculated from sand concentrations measured, processed and calibrated by USGS (Measured SL) and those predicted from single sand-rating curve (either linear or polynomial) (Predicted SL) at Tarbert Landing from 1986 to 1990.Specific terminologies pertaining to parts (a), (b), (c) and (d) of this figure i.e., Linear, Poly, CF, and No CF are same as explained in Figure 2.

Water 2015, 7 , 10 Figure 5 .
Figure 5. Confidence interval (CI) for the all "ln predicted SL" values at 95% confidence level in accordance to their corresponding SSCs-Qd rating curves during the four periods as shown at Tarbert Landing.It is noted that in all four periods (a), (b), (c) and (d), uppermost and lowermost curves represent the upper and lower limits of the CI, while the middle curve represents all individual "ln predicted SL" values as given in a.

Figure 5 .
Figure 5. Confidence interval (CI) for the all "ln predicted SL" values at 95% confidence level in accordance to their corresponding SSC s -Q d rating curves during the four periods as shown at Tarbert Landing.It is noted that in all four periods (a), (b), (c) and (d), uppermost and lowermost curves represent the upper and lower limits of the CI, while the middle curve represents all individual "ln predicted SL" values as given in a.

Water 2015, 7 , 12 Figure 6 .
Figure 6.Annual mean, maximum, and minimum of daily discharge at Tarbert Landing of the LmMR.

Figure 7 .
Figure 7. Flow duration curve for Tarbert Landing of the LmMR during 1973-2013.The vertical dash lines represent the exceedance probabilities for five river stages as defined by NOAA, i.e., Peak, High, Intermediate, Action and Low Flow Stages.

Figure 6 . 12 Figure 6 .
Figure 6.Annual mean, maximum, and minimum of daily discharge at Tarbert Landing of the LmMR.

Figure 7 .
Figure 7. Flow duration curve for Tarbert Landing of the LmMR during 1973-2013.The vertical dash lines represent the exceedance probabilities for five river stages as defined by NOAA, i.e., Peak, High, Intermediate, Action and Low Flow Stages.

Figure 7 .
Figure 7. Flow duration curve for Tarbert Landing of the LmMR during 1973-2013.The vertical dash lines represent the exceedance probabilities for five river stages as defined by NOAA, i.e., Peak, High, Intermediate, Action and Low Flow Stages.

Figure 8 .
Figure 8. Seasonal trend of monthly mean, maximum, and minimum of daily discharge at Tarbert Landing of the LmMR during 1973-2013.

Figure 8 .
Figure 8. Seasonal trend of monthly mean, maximum, and minimum of daily discharge at Tarbert Landing of the LmMR during 1973-2013.

Water 2015, 7 , 14 Figure 9 .
Figure 9. Average sand concentrations and percentage change with each 3000 cms increment of discharge at Tarbert Landing of the LmMR during 1973-2013.

Figure 10 .
Figure 10.Distribution of sand concentration with discharge at Tarbert Landing of the LmMR during 1973-2013.

Figure 9 . 14 Figure 9 .
Figure 9. Average sand concentrations and percentage change with each 3000 cms increment of discharge at Tarbert Landing of the LmMR during 1973-2013.

Figure 10 .
Figure 10.Distribution of sand concentration with discharge at Tarbert Landing of the LmMR during 1973-2013.

Figure 10 .
Figure 10.Distribution of sand concentration with discharge at Tarbert Landing of the LmMR during 1973-2013.
, higher than 40 MT for eight years including the Mississippi flood years of 1973, 1993 and 2011, and either > 10 MT or < 40 MT for the remaining 29 years (Figure 11).In addition, annual SL averaged 26.3 MT from 1973 to 1999 (28.8 MT from 1972 to 1979, 19.3 MT from 1980 to 1989, 32 MT from 1990 to 1999), which later increased to approximately 29 MT from 2000 to 2013 (calculations from table of

Figure 11 .
Figure 11.Annual mean, maximum, and minimum of daily sand loads and total annual sand loads at Tarbert Landing of the LmMR.

Figure 11 .
Figure 11.Annual mean, maximum, and minimum of daily sand loads and total annual sand loads at Tarbert Landing of the LmMR.

Figure 12 .
Figure 12.Seasonal trend of monthly mean, maximum, and minimum of daily sand load at Tarbert Landing of the LmMR during 1973-2013.

Figure 12 .
Figure 12.Seasonal trend of monthly mean, maximum, and minimum of daily sand load at Tarbert Landing of the LmMR during 1973-2013.

Figure 13 .
Figure 13.Annual flow volume above the Low and Action flow stages (discharge ≥ 18,000 cms) versus annual sand load at Tabert Landing of the LmMR during 1973-2013.

Table 6 .
Mean (µx), Standard Deviations (SD) (Sx) and respective parameters (a and b) in Gumbel distribution for maximum and minimum DSL (tonnes/day) at Tarbert Landing of the LmMR from 1973 to 2013.

Figure 13 .
Figure 13.Annual flow volume above the Low and Action flow stages (discharge ě 18,000 cms) versus annual sand load at Tabert Landing of the LmMR during 1973-2013.

Figure 14 .Table 7 .
Figure 14.Non-exceedance probabilities of the Gumbel distribution for maximum and minimum DSLs and Weibull distribution for total SLs each year at Tarbert Landing of the LmMR from 1973 to 2013.Table 7. Gumbel distribution based prediction of maximum and minimum DSLs and Weibull distribution based prediction of total annual SLs for the 2-, 5-, 10-, 20-, and 40-year returns at Tarbert Landing of the LmMR.

Figure 14 .Table 7 .
Figure 14.Non-exceedance probabilities of the Gumbel distribution for maximum and minimum DSLs and Weibull distribution for total SLs each year at Tarbert Landing of the LmMR from 1973 to 2013.

Table 1 .
Discharge-sand load rating curves developed for Tarbert Landing of the Lowermost Mississippi River (LmMR).Here, x = ln daily discharge (Q d ) (the independent variable) and y = ln daily sand load (DSL) (the dependent variable).

Table 2 .
Root mean square errors (RMSEs) of DSLs predicted through discharge-sand load rating curves for each period 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 3 .
Long-term flow conditions based on the U.S. National Oceanic and Atmospheric Administration (NOAA's) Mississippi River flow stages at Tarbert Landing of the LmMR from 1973 to 2013.

Table 3 .
Long-term flow conditions based on the U.S. National Oceanic and Atmospheric Administration (NOAA's) Mississippi River flow stages at Tarbert Landing of the LmMR from 1973 to 2013.
[21]: * The discharge ranges for Intermediate, High, and Peak Flow Stages are adopted from Rosen and Xu[21].

Table 4 .
Sand transport under five river stages at Tarbert Landing of the LmMR from 1973 to 2013.

Table 4 .
Sand transport under five river stages at Tarbert Landing of the LmMR from 1973 to 2013.

Table 6 .
Mean (µ x ), Standard Deviations (SD) (S x ) and respective parameters (a and b) in Gumbel distribution for maximum and minimum DSL (tonnes/day) at Tarbert Landing of the LmMR from 1973 to 2013.
overestimated annual SLs for 2008, 2009 and 2010 water years at Tarbert Landing.Allison et al.'s SL estimates for 2008, 2009 and 2010, respectively, were numerically higher than 38, 36 and 42 of 45 annual SL estimates provided by USACE for all water years from 1952 to 1996