River plumes are generated by the flow of buoyant river water into the coastal ocean, where they significantly influence water properties and circulation. They comprise dynamically distinct regions spanning a large range of spatial and temporal scales, each contributing to the dilution and transport of freshwater as it is carried away from the source. River plume structure varies greatly among different plume systems, depending on the forcing and geometry of each system. Individual systems may also exhibit markedly different characteristics under varied forcing conditions. Research over the past decade, including a series of major observational efforts, has significantly improved our understanding of the dynamics and mixing processes in these regions. Although these studies have clarified many individual processes, a holistic description of the interaction and relative importance of different mixing and transport processes in river plumes has not yet been realized.

No abstract

No abstract

No abstract

[1] Every summer, a large area (15,000 km 2 on average) over the Texas-Louisiana shelf in the northern Gulf of Mexico turns hypoxic due to decay of organic matter that is primarily derived from nutrient inputs from the Mississippi/Atchafalaya River System. Interannual variability in the size of the hypoxic zone is large. The 2008 Action Plan put forth by the Mississippi River/Gulf of Mexico Watershed Nutrient Task Force, an alliance of multiple state and federal agencies and tribes, calls for a reduction of the size of the hypoxic zone through nutrient management in the watershed. Comprehensive models help build mechanistic understanding of the processes underlying hypoxia formation and variability and are thus indispensable tools for devising efficient nutrient reduction strategies and for building reasonable expectations as to what responses can be expected for a given nutrient reduction.Here we present such a model, evaluate its hypoxia simulations against monitoring observations, and assess the sensitivity of the hypoxia simulations to model resolution, variations in sediment oxygen consumption, and choice of physical horizontal boundary conditions. We find that hypoxia simulations on the shelf are very sensitive to the parameterization of sediment oxygen consumption, a result of the fact that hypoxic conditions are restricted to a relatively thin layer above the bottom over most of the shelf. We show that the strength of vertical stratification is an important predictor of dissolved oxygen concentration in bottom waters and that modification of physical horizontal boundary conditions can have a large effect on hypoxia simulations because it can affect stratification strength.Citation: Fennel, K., J. Hu, A. Laurent, M. Marta-Almeida, and R. Hetland (2013), Sensitivity of hypoxia predictions for the northern Gulf of Mexico to sediment oxygen consumption and model nesting,

Abstract. The Texas-Louisiana shelf in the Northern Gulf of Mexico receives large inputs of nutrients and freshwater from the Mississippi/Atchafalaya River system. The nutrients stimulate high rates of primary production in the river plume, which contributes to the development of a large and recurring hypoxic area in summer, but the mechanistic links between hypoxia and river discharge of freshwater and nutrients are complex as the accumulation and vertical export of organic matter, the establishment and maintenance of vertical stratification, and the microbial degradation of organic matter are controlled by a non-linear interplay of factors. Unraveling these interactions will have to rely on a combination of observations and models. Here we present results from a realistic, 3-dimensional, physical-biological model with focus on a quantification of nutrient-stimulated phytoplankton growth, its variability and the fate of this organic matter. We demonstrate that the model realistically reproduces many features of observed nitrate and phytoplankton dynamics including observed property distributions and rates. We then contrast the environmental factors and phytoplankton source and sink terms characteristic of three model subregions that represent an ecological gradient from eutrophic to oligotrophic conditions. We analyze specifically the reasons behind the counterintuitive observation that primary production in the light-limited plume region near the Mississippi River delta is positively correlated with river nutrient input, and find that, while primary production and phytoplanktonCorrespondence to: K. Fennel (katja.fennel@dal.ca) biomass are positively correlated with nutrient load, phytoplankton growth rate is not. This suggests that accumulation of biomass in this region is not primarily controlled bottom up by nutrient-stimulation, but top down by systematic differences in the loss processes.

The structure of a river plume is related to the vertical mixing using an isohaline-based coordinate system. Salinity coordinates offer the advantage of translating with the plume as it moves or expanding as the plume grows. This coordinate system is used to compare the relative importance of different dynamical processes acting within the plume and to describe the effect each process has on the structure of the plume. Vertical mixing due to inertial shear in the outflow of a narrow estuary and wind mixing are examined using a numerical model of a wind-forced river plume. Vertical mixing, and the corresponding entrainment of background waters, is greatest near the estuary mouth where inertial shear mixing is large. This region is defined as the near field, with the more saline, far-field plume beyond. Wind mixing increases the mixing throughout the plume but has the greatest effect on plume structure at salinity ranges just beyond the near field. Wind mixing is weaker at high salinity classes that have already been mixed to a critical thickness, a point where turbulent mixing of the upper layer by the wind is reduced, protecting these portions of the plume from further wind mixing. The work done by mixing on the plume is of similar magnitude in both the near and far fields.

This numerical modeling study quantifies for the first time the contribution of various processes to estuarine circulation in periodically stratified tidal flow under the impact of a constant horizontal buoyancy gradient. The one-dimensional water column equations with periodic forcing are first cast into nondimensional form, resulting in a multidimensional parameter space spanned by the modified inverse Strouhal number and the modified horizontal Richardson number, as well as relative wind speed and wind direction and the residual runoff. The along-tide momentum equation is then solved for the tidal-mean velocity profile in such a way that it is equated to the sum of the contributions of tidal straining (resulting from the temporal correlation between eddy viscosity and vertical shear), gravitational circulation (resulting from the depth-varying forcing by a constant horizontal buoyancy gradient), wind straining, and depth-mean residual flow (resulting from net freshwater runoff). This definition of tidal straining does not only account for tidal asymmetries resulting from horizontal buoyancy gradients but also from wind straining and residual runoff. For constant eddy viscosity, the well-known estuarine circulation analytical solution with polynomial residual profiles is directly obtained. For vertically parabolic and constant-in-time eddy viscosity, a new analytic solution with logarithmic residual profiles is found, showing that the intensity of the gravitational circulation scales with the horizontal Richardson number. For scenarios with realistic spatially and temporally varying eddy viscosity, a numerical water column model equipped with a state-of-the-art two-equation turbulence closure model is applied to quantify the individual contributions of the various processes to estuarine circulation. The fundamental outcome of this study is that, for irrotational flow with periodic stratification and without wind forcing and residual runoff, the tidal straining is responsible for about two-thirds and gravitational circulation is responsible for about one-third of the estuarine circulation, proportionally dependent on the horizontal Richardson number, and weakly dependent on the Strouhal number. This new and robust result confirms earlier estimates by H. Burchard and H. Baumert, who suggested that tidal straining is the major generation mechanism for estuarine turbidity maxima. However, a sensitivity analysis of the model results to details of the turbulence closure model shows some uncertainty with respect to the parameterization of sheared convection during flood. Increasing down-estuary wind straining and residual runoff reduce the quantitative contribution of tidal straining. For relatively small horizontal Richardson numbers, the tidal straining contribution to estuarine circulation may even be reversed by down-estuary wind straining.