We present measurements which show that the tidal Ems River in Germ any is extremely muddy over a 30 km + turbid zone, with fluid mud o f 1-2 m thickness covering the bed with suspended sediment concentrations (SSC) o f greater than 50 kg.m-3. Moreover, we show that these elevated SSC contain large quantities of organic material which deplete dissolved oxygen (DO) and produce summertime hypoxic zones. Using mathematical modeling, we develop simplified representations o f the estuary physics that reproduce the tidally-averaged circulation, SSC distribution, and oxygen depletion. These models show that SSC and oxygen concentrations are extremely sensitive to factors such as the mean depth, the mixing due to bottom friction (turbulence), and river flow. The observed increase in SSC and decrease in DO over the past 25 years is linked to the progressive deepening o f the tidal Em s from 4-5 m to 7 m between 1985- 1994, which moved the turbid zone upstream and decreased mixing. A review of scientific literature and data from the Em s suggests that hum an intervention (dyking,channel modification) combines with more gradual natural changes (sea level rise, climate variation) to continually modify sediment transport

de Swart, Huib E.

Talke, Stefan A.

PDXScholar (Portland State University)

Portland State UniversityPDXScholarCivil and Environmental Engineering FacultyPublications and Presentations Civil and Environmental Engineering2011Circulation, Sediment Concentration and Oxygen Depletion inthe Tidal Ems RiverStefan A. TalkePortland State University, talke@pdx.eduHuib E. de SwartUtrecht UniversityLet us know how access to this document benefits you.Follow this and additional works at: http://pdxscholar.library.pdx.edu/cengin_facPart of the Civil Engineering Commons, and the Environmental Engineering CommonsThis Article is brought to you for free and open access. It has been accepted for inclusion in Civil and Environmental Engineering Faculty Publicationsand Presentations by an authorized administrator of PDXScholar. For more information, please contact pdxscholar@pdx.edu.Citation DetailsTalke, S.A.; de Swart, H.E. (2011). Circulation, sediment concentration and oxygen depletion in the tidal Ems River, in: Heip, C. et al.(Ed.) (2011). Aspects of coastal research in contribution to LOICZ in the Netherlands and Flanders (2002-2010). LOICZ Research& Studies, 38: pp. 42-483. Circulation, sedim ent concentration and oxygen depletion in the tidal Em s River Stefan A. Talke1 and H uib E. de Swart21 Civil and Environmental Engineering Department, Portland State University, 202T EngineeringBuilding, 1930 SW Fourth Avenue, Portland, OR, 97201, USA,e-mail: s.a.talke@pdx.edu2 Institute for Marine and Atmospheric Research, Utrecht University, Princetonplein 5, 3584 CCUtrecht, the Netherlands, e-mail: h.e.deswart@uu.nlAbstractWe present measurements which show that the tidal Em s River in Germany is extremely muddy over a 30 km  + turbid zone, with fluid m ud o f  1-2 m  thickness covering the bed w ith suspended sedim ent concentrations (SSC) o f  greater than 50 kg.mT M oreover, we show that these elevated SSC contain large quantities o f organic material which deplete dissolved oxygen (DO) and produce summertime hypoxic zones. Using mathematical modeling, we develop simplified representations o f the estuary physics that reproduce the tidally-averaged circulation, SSC distribution, and oxygen depletion. These models show that SSC and oxygen concentrations are extremely sensitive to factors such as the mean depth, the mixing due to bottom  friction (turbulence), and river flow. The observed increase in SSC and decrease in D O  over the past 25 years is linked to the progressive deepening o f  the tidal Em s from  4-5 m  to 7 m  between 1985- 1994, which m oved the turbid zone upstream  and decreased mixing. A  review o f scientific literature and data from  the Em s suggests that hum an intervention (dyking,channel modification) combines with more gradual natural changes (sea level rise, climate variation) to continually modify sedim ent transport.i. OverviewThe tidal Em s River is part o f  the Em s-D ollard estuary, which is located at the border o f Germ any and the Netherlands (Fig. A-3.1). The estuary extends approximately 100 km  from  the barrier island o f Borkum  (km 0) to a tidal weir (gate) in H erbrum  (km 100). Tidal forcing is significant, and the mean tidal range o f up to 3.5 m  produces currents o f 1-2 m .s 1 during peak flood, with smaller currents occurring during the ebb.42GermanyT e rb o rgB o rk u m ,’OgumDelfzi|iHabrumNetherlands -100 km(a)0 km 50 kmR iver(b )N um ber o f days below  a threshold  dissolved oxygen  concentration, Terborg1990 1995 2000 2005YearFig. A-3.1 : Location o f Ems Estuary (a), and increasing oxygen depletion over time (b); the crosses (X) denote location o f monitoring stations, solid line denotes location o f experimental transects.A gradual transition between salt water (S =30) and freshwater (S < 1) occurs in the middle o f  the estuary, with a peak gradient on the order o f  S = 1 km“1. Channel depths are maintained at a depth o f  8 m  in the outer estuary (maximum is 30 m), and are approximately 7 m  from  Em den (km 47) to Papenburg (km 87). River discharge is primarily from  the Em s River, and varies from  an average o f  35 m bs“1 during the warm  summer m onths to peak flows o f  up to 600 m3.s_1 during w inter m onths. The yearly average flow is 80 mbs"1, with a median flow o f  ~65 mbs"1. O ver the past 25 years, the num ber o f  days with severe oxygen depletion has greatly increased (Fig. A-3.1). D uring the same period, the river between Em den and Papenburg was deepened from  4-5 m  to 7 m  (1985-1994), and an increase in tidal range o f  up to 1.5 m  was observed near the tidal weir (Talke & de Swart, 2006).W e undertook extensive hydrodynamic and water quality measurements in m onthly increments from  Feb. 2005 to Dec. 2007 along the longitudinal axis o f  the Em s estuary. Selected cross- sectional and fixed channel m easurem ents occurred over a tidal period in M arch 2005, February 2006, and Septem ber 2006 in the area near Pogum  (km 53) and Terborg (km 63). As explained in m ore detail in Talke et al. (2009a), all these experiments included measurem ents o f  water quality parameters such as nutrient concentrations, SSC, and chlorophyll. Additionally, many o f  these field experiments included m easurem ents o f  flow velocity, salinity, tem perature, dissolved oxygen (DO) concentration, and optical backscatter (for turbidity and SSC). We focus here on these latter measurements.1 0 0(/)>v(0(D5 mg/L< 2 mg/L43ii. Field MeasurementsDuring low river discharge conditions during the summer, salinity is well mixed vertically and transitions from  marine conditions (S > 20) to freshwater conditions (S < 1) over nearly 30 km (Fig. A-3.2a). As occurs in many estuaries, SSCs vary spatially and collect in the so-called estuary turbidity maximum (ETM; Fig. A-3.2b). Unlike m ost other estuaries, this ETM  extends outside o f  the brackish water region (S= 1-2) and covers the entire freshwater zone up to the tidal weir at km  100. SSCs in the ETM  are enorm ous and form  a fluid m ud layer (SSC > 1 0  kg.rn3) o f  1-2 m depth over a 35 km  region. The density gradient in the vertical suppresses mixing. The SSC varies sharply in the brackish water region, and produces a longitudinal density gradient that is up to an order o f  magnitude larger than the longitudinal salinity gradient (on the order o f 10 kg.rn3 near the bed).Fig. A-3.2 : Measurements o f salinity (a), suspended sediment concentration (b), and dissolved oxygen concentration (c) along the axis o f the Ems River Estuary on August 2nd, 2006 (Reprinted from Talke et al., 2009b)..A t the water surface, turbidity at S = 4 salinity attenuates light intensity by 3 orders o f  magnitude within the first 15 cm o f the water column (not shown). Under these conditions, phytoplankton cannot grow. The effect o f  turbidity on the onset o f  phytoplankton blooms was theoretically investigated by De Swart et al. (2009), and shows that m uch lower SSC concentrations than measured in the Em s ETM  either move the phytoplankton bloom  seaward or prevent one from  occurring.44The oxygen demand from  organic material (approximately 10 % o f  SSC os organic carbon) in the ETM  is sufficient to produce a large zone o f hypoxia (< 2 m g.f1) over m uch o f the freshw ater/brackish water region (Fig. A-3.2c; see Talke et al., 2009b). To first order, we find that oxygen demand is proportional to SSC, with the smallest Dissolved Oxygen (DO) occurring in the region o f  largest SSC (Talke et al., 2009b). Fish populations are stressed and cannot live at D O  < 5 m g.f1 and D O  < 2 m g.f1, respectively. Flence, understanding the physical conditions and processes that produce large SSC, fluid mud, and low D O  are essential to devising solutions for the stressed ecology o f the Ems.iii. M odellingTo understand why SSC is elevated and oxygen is depleted, we developed mathematical (analytical) models that identify im portant controlling parameters in a simplified, idealized domain. While such models neccessarily exclude detail such as complex bathymetry, its simplicity allows for param eter sensitivity studies that provide insight into the effect o f  changing conditions. The equations can be found in Talke et al. (2008, 2009a,b); here, we qualitatively describe our methodology.As a first step, we assume that the distribution o f  SSC and depleted oxygen is determined primarily by tidally averaged circulation and by mixing, which can be captured by constant vertical and horizontal mixing coefficients. M oreover, we assume a constant depth and an exponentially decreasing w idth (funnel shape), w ith an e-folding scale that is based on the natural environm ent. We assume a rigid lid and make the Boussinesq assum ption, following standard practice. River discharge is assumed to be constant. Salinity is observed to be well mixed vertically and varies like a hyperbolic tangent in the longitudinal direction (see Fig. A-3.2).The assumptions listed above were classically applied by Flansen & Rattray (1965) to the tidally averaged shallow water equations, which describe how pressure— either from  a surface slope or from  a horizontal density gradient-produces circulation. We modify the classical relationship by recognizing that dynamically significant density gradients occur no t only due to salinity, bu t also due to the longitudinal distribution o f SSC (Fig. A-3.2). M oreover, we note that SSC is distributed approximately exponentially in the vertical, with the largest concentrations near the bed (Talke et al., 2009a). This is modelled by assuming a balance vertically between diffusion and a constant settling velocity. These insights are then mathematically incorporated into the shallow w ater equations. The resulting novel expression describes circulation patterns that arise due the salt field, river discharge, and the distribution o f SSC. As shown in Figure A-3.3a and A-3.3b, surface flow is seaward everywhere. N ear the bed, a cell o f  upstream  circulation occurs due the salinity gradient. Turbidity induced circulation downstream  o f the ETM  opposes, bu t does not reverse, this upstream  flow. A  large dependence on depth is observed; increasing the depth from  5 m  to 7 m, which occurred in the Em s between 1985 and 1994, results in a factor o f  nearly 3x greater circulation.The distribution o f SSC presented in Figure A-3.3c and A-3.3d is solved for by assuming equilibrium conditions. Mathematically, we assume that the local deposition o f  sediment equals local erosion and that the quantity o f  suspended sediment in the model domain is fixed. Consequently, the total sedim ent flux through a cross-section vanishes, i.e., the area-integrated sedim ent fluxes due to a combination o f advection and diffusion are zero.45o-2-4H = 5 m  (a) u (x ,z )^J   ^ [AGOiO 20 40 60 80 100SN0(c) C(x,z)■2■40 20 40 60 80 100H =  7 m(b) u ( x \ z f ~20 40 60 80 100(d) C(x,z)20 en -4 20 en0-2-4(e) 0 2fx,z;IR0 20 40 60 80 100distance (km)86 V4  en E200 20 40 60 80 1000-2-4-6(f) O /x .z )in-sr0 20 40 60 80 100distance (km)86 V 4  en20Fig. A-3.3 : Modeled circulation, SSC, and dissolved oxygen for different depths (Reprinted from Talke et al., 2009b).The result is an analytical expression for SSC that depends on the salinity field, the freshwater discharge, depth, the horizontal and vertical mixing coefficient, the settling velocity o f sediment, and the prescribed total am ount o f  suspended sediment in the system. Com paring Figures A-3.2b and A-3.3d, we observe that the modelled turbidity zone (Fig. A-3.3d) qualitatively captures the features o f  the measured distribution o f  SSC (Fig. A-3.2b). Elevated SSC o f greater than 10 kg.rn3 extends over the freshwater and brackish water zone and covers the bed by 1-2 m  in both  the m easurem ent and model. M oreover, the SSC gradient in the brackish water transition zone is steep in both  cases.Interestingly, increasing the depth from  5 m  (Fig. A-3.3c) to 7 m  (Fig. A-3.3d) results in a larger turbidity zone w ith nearly an order o f magnitude greater SSC that is centered farther upstream. Similar changes have been observed in the Em s River, in which channel deepening from  4-5 m  to 7 m  between 1985 and 1994 resulted in an upstream  m ovem ent o f the ETM  and greatly increased SSC (e.g., Talke & de Swart, 2006). O ur model suggests that contributors to the historically observed variation include: (1) decreased river flow (because o f greater cross-sectional area), (2) increased gravitational circulation (cubic dependence on depth), and (3) a changed vertical distribution o f SSC (which is now proportionally lower in the water column). O ther factors include decreased mixing and increased sediment supply (not shown). M ore recent models also point out the im portance o f tidal asymmetries, increased tidal range, settling lag46effects, reduced bottom  friction, and flocculation processes in producing an upstream  m ovem ent o f  SSC (Chemetsky et al., 2010; W interwerp 2010).The distribution o f circulation and SSC is next used to model oxygen depletion. A  mass-balance approach to dissolved oxygen (DO) yields an advection-diffusion equation with a sink term  that is proportional to organic material, and hence SSC. Reaeration at the surface is held proportional to the local oxygen deficit, while sedim ent oxygen demand at the bed produces an additional source o f  oxygen depletion. Additional details and model verification are given in Talke et al. (2009b). As shown in Figure A-3.3e and A-3.3f, the modeled oxygen depletion is closely related to the distribution o f SSC. The changes described earlier for SSC between the H = 5  m  case and the H = 7  m  case cause a switch from  well aerated conditions (DO > 5 mg.1-1) to conditions that are hypoxic (DO < 2 m g.f1) over a large section o f  the estuary. From  these simulations we conclude that it is the dynamics o f  turbidity maxima— rather than the residence time o f  water or the input o f nutrients— that determines the depletion o f oxygen in highly turbid estuaries such as the Ems.A review o f  scientific literature and data from  the Em s suggests that hum an intervention and natural variation continually modify sedim ent transport (Talke & de Swart, 2006). Until approximately 1950, the largest change to estuarine physics (as m easured for example by tidal range) occurred due to land reclamation. Since 1950, channel m odification and deepening has amplified tidal range, increased the tidal asymmetry in currents, magnified SSC and turbidity, and produced hypoxia. O ther sources o f  variation include rising sea level and increasing tidal range in the N orth  Sea, as well as variations to wind, waves, and river discharge caused by the N orthern Atlantic Oscillation cycle (Talke & de Swart, 2006). Although natural changes are occurring m ore gradually than anthropogenic modifications, both factors alter depth, flow and mixing characteristics (Talke & de Swart, 2006). The models and measurements presented here suggest how these factors affect estuarine sedim ent transport patterns and water quality.iv. ConclusionsIn this paper we summarize recent measurements and models which show that the Em s estuary is a highly stressed environm ent m arked by extremely turbid conditions, fluid mud, and hypoxia. A  tidally-averaged circulation model is used to solve for the distribution o f  suspended sediment concentrations and dissolved oxygen. The model results suggest that deepening o f the Em s River is responsible for the increased SSC and degraded w ater quality over the past 25 years.ReferencesChemetsky, A.S., Schuttelaars, H.M. & Talke, S.A. (2010). The effect o f tidal asymmetry and temporal settling lag on sediment trapping in tidal estuaries. Ocean Dynamics 60: 1219-1241, doi: 10.1007/sl0236-010-0329-8.Hansen, D.Y. & Rattray, M. Jr. (1965). Gravitational circulation in straits and estuaries. Journal o f Marine Research 23:104—122.De Swart, H.E., Schuttelaars, H.M. & Talke, S.A. (2009). Phytoplankton growth in turbid estuaries: a simple model. Continental Shelf Research 29:136—147, doi:10.1016/j.csr.2007.09.006.47Talke, S.A. & De Swart, H.E. (2006). Hydrodynamics and morphology in the Ems/Dollard Esutary: review of models, measurements, scientific literature, and the effects of changing conditions. Report for Rijkswaterstaat (Ministry of Public Works, Netherlands), IMAU Report 06-01, 78 pp.Talke, S.A., De Swart, H.E. & Schuttelaars, H.M. (2008). An analytical model for the equilibrium distribution of sediment in an estuary. In: River, coastal and estuarine morphodynamics, ed. C.M. Dohmen-Janssen and S.J.M.H. Hulscher, pp. 403M-12. Tondon: Taylor & Francis.Talke, S.A., De Swart, H.E. & Schuttelaars, H.M. (2009a). Feedback between residual circulations and sediment distribution in highly turbid estuaries: an analytical model. Continental Shelf Research 29:119-135, doi:10.1016/j.csr.2007.09.002.Talke, S.A., De Swart, H.E. & De Jonge, Y.N. (2009b). An idealized model and systematic process study of oxygen depletion in highly turbid estuaries. Estuaries and Coasts 32:602—620, doi: 10.1007/sl2237- 009-9171-y.http://depts.washington.edu/uwefm/publications/Talke_deSwart_deJonge_Idealized_Model_oxygen_depletion_EC_2009.pdfWinterwerp, J.C. (2010). Fine sediment transport by tidal asymmetry in the high-concentrated Ems River: indications for a regime shift in response to channel deepening. Ocean Dynamics , doi: 10.1007/sl0236-010-0332-0.48

Circulation, Sediment Concentration and Oxygen Depletion in the Tidal Ems River

http://pdxscholar.library.pdx.edu/cgi/viewcontent.cgi?article=1080&context=cengin_fac

Circulation, Sediment Concentration and Oxygen Depletion in the Tidal Ems River

Abstract

Similar works

Full text

Available Versions

PDXScholar (Portland State University)