The Columbia River plume as cross-shelf exporter and along-coast barrier

An intensive Lagrangian particle-tracking analysis of the July 2004 upwelling period was conducted in a hindcast model of the US Pacific Northwest coast, in order to determine the effect of the Columbia River plume on the fate of upwelled water. The model, implemented using Regional Ocean Modeling System (ROMS), includes variable wind and atmospheric forcing, variable Columbia river flow, realistic boundary conditions from Navy Coastal Ocean Model (NCOM), and 10 tidal constituents. Model skill has been demonstrated in detail elsewhere [MacCready, P., Banas, N.S., Hickey, B.M., Dever, E.P., Liu, Y., 2008. A model study of tide- and wind-induced mixing in the Columbia River estuary and plume. Continental Shelf Research, this issue, doi:10.1016/j.csr.2008.03.015]. Particles were released in the Columbia estuary, along the Washington coastal wall, and along the model's northern boundary at 48°N. Particles were tracked in three dimensions, using both velocities from ROMS and a vertical random displacement representing turbulent mixing. When 25 h of upwelling flow is looped and particles tracked for 12 d, their trajectories highlight a field of transient eddies and recirculations on scales from 5 to 50 km both north and south of the Columbia. Not all of these features are caused by plume dynamics, but the presence of the plume increases the entrainment of inner-shelf water into them. The cumulative effect of the plume's interaction with these transient features is to increase cross-shelf dispersion: 25% more water is transported laterally past the 100 m isobath when river and estuarine effects are included than when they are omitted. This cross-shelf dispersion also disrupts the southward transport of water along the inner shelf that occurs in the model when the Columbia River is omitted. This second effect—increased retention of upwelled water on the Washington shelf—may be partly responsible for the regional-scale alongcoast gradient in chlorophyll biomass, although variations in shelf width, the Juan de Fuca Eddy to the north, and the intermittency of upwelling-favorable winds are likely also to play important roles.     

Dynamic characteristics and horizontal transports of internal solitons generated at the Columbia River plume front

River plume front-generated internal solitons play an important role in the interaction between the plume and coastal waters. The internal solitons drive a non-harmonic velocity field, resulting in a horizontal transport that carries plume water seaward and redistributes nutrients and sediments. In this study, we present observations of internal solitons generated at the Columbia River plume front that separates the new, tidal plume, older plume and coastal waters. Scale analyses suggest that the plume front-generated internal solitons are highly non-linear waves, and their dynamic properties do not conform to any weakly non-linear theory. Thus, a high-order Korteweg–de Vries (KdV) theory is used to analyze the internal solitons. The comparison between theoretical values and cruise data shows that the high-order KdV model is much better than the weakly non-linear theories for prediction of the soliton dynamic parameters. Based on the model, we develop theoretical and numerical solutions of the soliton-induced upper layer horizontal transport and Lagrangian water parcel transport distance, which shows that the water particle drift, during the internal soliton passage, is as far as 1 km, and demonstrates the role of the internal solitons on the exchange between the plume and ambient coastal water. Energy fluxes caused by the internal solitons are estimated using the high-order KdV theory. The leading soliton fluxes 2.0×10³ W m⁻¹ per unit crest length, and carries energy of 4.2×10⁵ J m⁻¹. The total energy carried by the eight internal solitons is 1.6×10⁶ J m⁻¹, about 70% of the total frontal energy.     

On the sensitivity of tidal residuals off the west coast of Britain to mesh resolution

An irregular mesh model of the west coast of Britain is used to examine the sensitivity of tidal residuals to mesh resolution in the region. Computed residuals are compared with earlier published results determined with a high resolution (1 km grid) finite difference model of the eastern Irish Sea. Initial calculations show that tidal residuals are largest in nearshore regions particularly in the vicinity of headlands. Local refinement of the mesh in these regions leads to a more detailed picture of the flow field, particularly adjacent to the coast. Although large scale offshore features of the flow can be resolved using the high resolution finite difference model, such an approach leads to a “stair case” representation of the coastal boundary with an adjacent near coastal region of spurious tidal residuals. By using an irregular mesh that follows the coast, this effect is removed. In the Mersey river region the tidal residual is resolved with a mesh resolution of 120 m, although calculations show that its distribution is particularly sensitive to small scale features of the topography. A variable mesh that can accurately represent the lateral variations in river width and details of topography in both the nearshore and estuarine environment appears essential in modelling the coastal spread of freshwater plumes from rivers and pollutants discharged into the near coastal environment.     

Patos Lagoon outflow within the Río de la Plata plume using an airborne salinity mapper: Observing an embedded plume

Major river systems discharging into continental shelf waters frequently form buoyant coastal currents that propagate along the continental shelf in the direction of coastal trapped wave propagation (with the coast on the right/left, in the northern/southern hemisphere). The combined flow of the Uruguay and Paraná Rivers, which discharges freshwater into the Río de la Plata estuary (Lat. ∼36°S), often gives rise to a buoyant coastal current (the ‘Plata plume’) that extends northward along the continental shelf off Uruguay and Southern Brazil. Depending upon the prevailing rainfall, wind and tidal conditions, the Patos/Mirim Lagoon complex (Lat. ∼32°S) may also produce a freshwater outflow plume that expands across the inner continental shelf. Under these circumstances the Patos outflow plume can be embedded in temperature, salinity and current fields that are strongly influenced by the larger Plata plume. The purpose of this paper is to present observations of such an embedded plume structure and to determine the dynamical characteristics of the ambient and embedded plumes.     

Numerical modeling of the Amazon River plume

Marine circulation above the northern Brazilian continental shelf is subject to energetic forcing factors of various origins: high water buoyancy fluxes induced by the Amazon River freshwater discharge, a strong coastal current associated with a mesoscale current (North Brazil Current (NBC)), a forcing by semidiurnal tide and by Northeast or Southeast trade winds according to the season. Using a three-dimensional (3-D) hydrodynamic numerical model (MOBEEHDYCS), and realistic bathymetry and coastline of the northern Brazilian shelf, this paper aims at studying the influence of some specific physical processes on the morphology of the Amazon plume. The very large volume discharge (180 000 m³/s on average) and the weak effect of Coriolis force are additional characteristics of the studied system, which induce a particular dynamics. The various forcing factors are successively introduced into the model in order to simulate and to determine their respective influences upon the plume extent and the hydrodynamics at the shelf scale. Simulation reveal that the coastal current is at the origin of the permanent northwestward Amazon plume extension while wind effect can either reinforce or moderate this situation. The tide intervenes also to modify the position of the salinity front: a horizontal migration of salinity front is observed under its action.     

Model of the Long Island Sound outflow: Comparison with year-long HF radar and Doppler current observations

A three-dimensional primitive-equation model is used to simulate the Long Island Sound (LIS) outflow for a 1-year (2001) period. The model domain includes LIS and New York Bight (NYB). Tidal and wind forcing are included, and seasonal salinity and temperature variations are assimilated. The model results are validated with the HF radar, moored acoustic Doppler current profiler (ADCP), and ferry-based ADCP observations. The agreement between simulated and observed flow patterns generally is very good. The difference in seasonal mean currents between the model and moored ADCP is about 0.01 m/s; the correlation of dominant velocity fluctuations between the model and HF radar is 0.83; and the difference in mean LIS transport between the model and shipboard ADCP is about 5%. However, the model predicts a prominent tidally generated headland eddy not supported by the HF radar observation. The model sensitivity study indicates that the tides, winds, and ambient coastal front all have important impact on the buoyant outflow. The tides and winds cause stronger vertical mixing, which reduces the surface plume strength. The ambient coastal front, on the other hand, tends to enhance the plume.     

A model study of tide- and wind-induced mixing in the Columbia River Estuary and plume

A numerical simulation of circulation in the Columbia River estuary and plume during the summer of 2004 is used to explore the mixing involved as river water is transformed into shelf water. The model is forced with realistic river flow, tides, wind stress, surface heat flux, and ocean boundary conditions. Simulated currents and water properties on the shelf near the mouth are compared with records from three moorings (all in 72 m of water) and five CTD sections. The model is found to have reasonable skill; statistically significant correlations between observed and modeled surface currents, temperature, and salinity are all 0.42–0.72 for the mooring records. Equations for the tidally averaged, volume-integrated mechanical energy budget (kinetic and potential) are derived, with attention to the effects of: (i) Reynolds averaging, (ii) a time varying volume due to the free surface, and (iii) dissipation very close to the bottom. It is found that convergence of tidal pressure work is the most important forcing term in the estuary. In the far field plume (which has a volume 15 times greater than that of the estuary), the net forcing is weaker than that in the estuary, and may be due to either tidal currents or wind stress depending on the time period considered. These forcings lead to irreversible mixing of the stratification (buoyancy flux) that turns river water into shelf water. This occurs in both the plume and estuary, but appears to be more efficient (17% vs. 5%), and somewhat greater (4.2 MW vs. 3.3 MW), in plume vs. estuary. This demonstrates the importance of both wind and tidal forcing to watermass transformation, and the need to consider the estuary and plume as part of a single system.     

Determination of transport rates in the Yellow River–Bohai Sea mixing zone via natural geochemical tracers

In light of the current problems facing the Yellow River and surrounding areas (e.g., periods of zero river discharge, increasing nitrate concentrations of the Bohai Sea), we examined the coastal mixing dynamics around the mouth of the Yellow River. Naturally occurring radium isotopes (²²³Ra, ²²⁴Ra, ²²⁶Ra, and ²²⁸Ra) and other geochemical tracers (Ba, Si, and salinity) were employed to determine river plume transport scales and rates. Barium and radium exhibit elevated concentrations within the salinity gradient where they are desorbed from particles via ion-exchange. Once they are added to the system, they decrease offshore from dilution with lower concentration Bohai Sea water, and in the case of ²²⁴Ra and ²²³Ra, by radioactive decay. Using radium “ages” to assess the dissolved material transport scales and rates proved to be a useful tool in this environment. The ages based on the ²²⁴Ra/²²⁸Ra activity ratio increased gradually until salinities reached ∼25 when they rapidly increased due to decreased mixing at higher salinities. Integrated net transport rates through the salinity front ranged from 1.4 to 1.6 cm/s and did not vary significantly with river discharge. Thus, tidal mixing appears to dominate in this system, at least over the range of discharges investigated (80–600 m³/s). Determining the temporal scale of flow across the coastal zone in this region is a valuable first step toward examining whether the Yellow River is contributing to the increasing inorganic nitrogen concentrations in the central Bohai Sea.     

Numerical study of generation of the tidal shear front off the Yellow River mouth

Tidal shear front off the Yellow River mouth has been observed and modeled in the previous studies. However, a detailed investigation of the front generation has not been conducted. The aim of this paper is to use a three-dimensional tidal model coupled to a sediment transport module to examine the front formation. The model predicted a tidal shear front that propagated offshore and lasted 1–2 h at both flood and ebb phase off the Yellow River mouth. The sensitivity numerical experiments showed that the topography with a strong slope off the Yellow River mouth was a determining factor for the front generation, and a parallel orientation between the major axes of ellipses and co-tidal lines of maximum tidal current was a necessary condition. While the bottom friction and the river runoff had no effect on the front location but affected the front intensity, the front generation was not sensitive to the coastline variation. The study concluded that the bottom slope off the river mouth induces a strong variation in the bottom stress in a cross-shore direction, which produces both maximum phase gradient and sediment concentration variability across the tidal shear front. With the extending Yellow River delta, the tidal shear front under the new bathymetry of year 2003 has been strengthened and pushed further offshore due to an increased bottom slope.     

Observation of the supercritical Pearl River plume front under the downwelling-favorable winds

For the Pearl River plume, the supercritical, distinct plume front appears in downwelling-favorable winds, which is easily observed due to the distinct boundary between the plume water and the ambient water. In this paper, in situ and satellite observations of a plume front are utilized to explore the Pearl River plume front properties under the downwelling-favorable winds. Field observations clearly show frontal structure, especially the two-layer structure in the plume water and the downward-motion of water in the frontal region. The Advanced Synthetic Aperture Radar(ASAR) images are also analyzed to unveil the plume front: there is a white stripe on the west side out of the river mouth under downwelling-favorable winds, which is identified as a supercritical plume front, and the width of the plume front is about 250 m. The normalized velocity gradient shows the intense velocity convergence in the front region. Also, analyses of ASAR images imply that the river discharge plays an important role in controlling the location and shape of the front.     

River plume patterns and dynamics within the Southern California Bight

Stormwater river plumes are important vectors of marine contaminants and pathogens in the Southern California Bight. Here we report the results of a multi-institution investigation of the river plumes across eight major river systems of southern California. We use in situ water samples from multi-day cruises in combination with MODIS satellite remote sensing, buoy meteorological observations, drifters, and HF radar current measurements to evaluate the dispersal patterns and dynamics of the freshwater plumes. River discharge was exceptionally episodic, and the majority of storm discharge occurred in a few hours. The combined plume observing techniques revealed that plumes commonly detach from the coast and turn to the left, which is the opposite direction of Coriolis influence. Although initial offshore velocity of the buoyant plumes was ∼50 cm/s and was influenced by river discharge inertia (i.e., the direct momentum of the river flux) and buoyancy, subsequent advection of the plumes was largely observed in an alongshore direction and dominated by local winds. Due to the multiple day upwelling wind conditions that commonly follow discharge events, plumes were observed to flow from their respective river mouths to down-coast waters at rates of 20–40 km/d. Lastly, we note that suspended-sediment concentration and beam-attenuation were poorly correlated with plume salinity across and within the sampled plumes (mean r²=0.12 and 0.25, respectively), while colored dissolved organic matter (CDOM) fluorescence was well correlated (mean r²=0.56), suggesting that CDOM may serve as a good tracer of the discharged freshwater in subsequent remote sensing and monitoring efforts of plumes.     

Response of salinity distribution around the Yellow River mouth to abrupt changes in river discharge

To investigate how salinity changes with abrupt increases and decreases in river discharge, three surveys were conducted along six sections around the Yellow River mouth before, during and after a water regulation event during which the river discharge was increased from ∼200 to >3000 m³ s⁻¹ for the first 3 days, was maintained at >3000 m³ s⁻¹ for the next 9 days and was decreased to <1000 m³ s⁻¹ for the final 4 days. The mean salinity in the Yellow River estuary area during the event varied ∼1.21, which is much larger than its seasonal variation (∼0.50) and interannual variation (∼0.05). Before the event, a small plume was observed near the river mouth. During the event, the plume extended over 24 km offshore in the surface layer in the direction of river water outflow. After the event, the plume diminished in size but remained larger than before the event. The downstream propagation of the plume (as in a Kelvin wave sense) was apparent in the bottom layer during the second survey and in both the surface and bottom layers during the third survey. The plume sizes predicted by the formulas from theoretical studies are larger than those we observed, indicating that factors neglected by theoretical studies such as the temporal variation in river discharge and vertical mixing in the sea could be very important for plume evolution. In addition to the horizontal variation of the plume, we also observed the penetration of freshwater from the surface layer into the bottom layer. A comparison of two vertical processes, wind mixing and tidal mixing, suggests that the impact of wind mixing may be comparable with that of tidal mixing in the area close to the river mouth and may be dominant over offshore areas. The change in Kelvin number indicates an alteration of plume dynamics due to the abrupt change in river discharge during the water regulation event.     

Factors influencing the spreading of a low-discharge river plume

Coastal plumes resulting from the continuous discharge of brackish or fresh river water are common features of continental and shelf seas. They are important for several aspects of the coastal environment, and can influence the local socio-economy to some degree. It is known from many studies that the evolution of plumes depends on various factors, such as the local bathymetry, hydrodynamics and meteorological conditions; most of these works; however, have focused on medium to large-scale rivers, while the smaller-scale discharges commonly found in the microtidal environments of the Mediterranean Sea have been less studied. This paper is centred on the behaviour of a freshwater plume arising from one of such outflows, in terms of both the physical configuration of the waterbody and the characteristics of the main driving mechanisms (discharge rate and wind stress). The modelled cases correspond to an open shallow bay, limited at one end by a large headland, and into which a typical Mediterranean waterway discharges. This particular setup is representative of a number of different bays existing on the Eastern Spanish coast. The numerical results highlight the large influence of the bay's topography on the river plume's extension and inner structure.     

Interaction of a river plume with coastal upwelling in the northeastern South China Sea

Observational and modeling studies were conducted to investigate the Pearl River plume and its interaction with the southwesterly driven upwelling circulation in the northern South China Sea during the summer. After exiting the Pearl River Estuary, the discharged freshwater generates a nearly stationary bulge of freshwater near the entrance of the estuary. Forced by the wind-driven coastal upwelling current, the freshwater in the outer part of the bulge flows downstream at the speed of the current and forms a widening and deepening buoyant plume over the shelf. The plume axis gradually shifts offshore of the current maximum as a result of currents induced by the contrasting density at the nose of plume and by the intensified Ekman drift in the plume. In this plume–current system, the fraction of the discharged freshwater volume accumulated in the bulge reaches a steady state and the volume of newly discharged freshwater is transported downstream by the upwelling current. Enhancement of stratification by the plume thins the surface frictional layer and enhances the cross-shelf circulation in the upper water column such that the surface Ekman current and compensating flow beneath the plume are amplified while the shoaling of the deeper dense water in the upwelling region changes minimally. The pressure gradient generated between the buoyant plume and ambient seawater accelerates the wind-driven current along the inshore edge of the plume but retards it along the offshore edge. Along the plume, downward momentum advection is strong near the highly nonlinear source region and a weaker upward momentum advection occurs in the far field over the shelf. Typically, the plume is shaped by the current over the shelf while the current itself is adjusting to a new dynamic balance invoked by the plume-induced changes of vertical viscosity and the horizontal pressure gradient. The spatial variation of this new balance leads to a coherent change in the cross-isobath transport in the upper water column during upwelling.     

Rapid formation of hyperpycnal sediment gravity currents offshore of a semi-arid California river

Observations of sediment dispersal from the Santa Clara River of southern California during two moderately sized river discharge events suggest that river sediment rapidly formed a negatively buoyant (hyperpycnal) bottom plume along the seabed within hours of peak discharge. An array of acoustic and optical sensors were placed at three stations 1 km from the Santa Clara River mouth in 10-m water depth during January–February 2004. These combined observations suggest that fluid mud concentrations of suspended sediment (>10 g/l) and across-shore gravity currents (∼5 cm/s) were observed in the lower 20–40 cm of the water column 4–6 h after discharge events. Gravity currents were wave dominated, rather than auto-suspending, and appeared to consist of silt-to-clay sized sediment from the river. Sediment mass balances suggest that 25–50% of the discharged river sediment was transported by these hyperpycnal currents. Sediment settling purely by flocs (∼1 mm/s) cannot explain the formation of the observed hyperpycnal plumes, therefore we suggest that some enhanced sediment settling from mixing, convective instabilities, or diverging plumes occurred that would explain the formation of the gravity currents. These combined results provide field evidence that high suspended-sediment concentrations from rivers (>1 g/l) may rapidly form hyperpycnal sediment gravity currents immediately offshore of river mouths, and these pathways can explain a significant portion of the river-margin sediment budget. The fate of this sediment will be strongly influenced by bathymetry, whereas the fate of the remaining sediment will be much more influenced by ocean currents.     

Quantifying the suspended sediment discharge to the ocean from the Markham River,Papua New Guinea

The Markham River is a small river draining a tropical mountain range with altitudes between 1000 and 3000 m and discharges directly into a submarine canyon, the head of which is at 30 m depth and reaches depths of 500 m only 4 km from the shore. As such, the Markham discharge system serves as a possible analogue for rivers discharging onto margins during low stands of sea-level. Located in a tectonically active area and with high rainfall, sediment supply is high and episodic and is sometimes related to catastrophic mountain landslides. The river has an estimated sediment load of 12 Mt yr⁻¹. Occasionally, high energy flows are generated at the river mouth which is evident from the channel morphology and sediment distribution. Profiles of salinity and suspended sediment concentrations (SSC) show that sediment is dispersed via a plume with components at both the surface, intermediate depth along isopycnal surfaces and near the sea bed. The dispersal pattern of the surface freshwater plume is largely determined by the buoyancy force. The surface plume is very thin with salinity gradients 15 ppt m⁻¹ while a Richardson number greater than unity suggested that the mixing zone is highly stratified. Estimates of the horizontal sediment flux gradient of the surface plume along the estuary axis suggest that about 80% of the sediment discharged is lost from the plume within a distance of 2 km from the river mouth. Particle fall velocities estimated from the vertical flux indicate values less than those of flocculated material. Layers of sediment with SSCs between 500 and 1000 mg l⁻¹ were observed at intermediate depths and near the seabed during periods of both high and intermediate discharge. The mass of sediment in a SSC layer at intermediate depths between 150 and 250 m within the canyon channel was estimated to be equivalent to an average of 2 to 3 days of Markham sediment discharge. SSCs near the seabed of between 250 and 750 mg l⁻¹ suggest that layers of significantly elevated density exist near the seabed, moving under the influence of gravity down steep seabed slopes of the Markham canyon.     

The formation of headland/island sandbanks

There are four extensive sandbanks in the vicinity of the Isle of Portland, a headland in the English Channel. The formation and maintenance of the two most prominent of these sandbanks (one on either side of the headland) can largely be explained by net bedload convergence, driven by instantaneous headland eddies generated by tidal flow past the headland. However, there are also two less prominent sandbanks (again, one on either side of the headland), which are not located in zones of bedload convergence. It is suggested here that these latter two sandbanks were formed when the Isle of Portland was isolated from the mainland by a tidal strait. Relative sea-level data and radiocarbon dates indicate that this would have occurred ca. 9–7 ka BP, prior to the closure of the strait by sedimentation. Tidal flow through this strait generated eddy systems in addition to the headland eddies, leading to the formation of associated headland/island sandbanks. At 7 ka BP, sedimentation resulted in closure of the strait, leading to the present-day headland configuration, and subsequent reworking of these now moribund sandbanks formed by the strait. A series of idealised morphological model experiments, parameterised using bedrock depths and glacial isostatic adjustment model output of relative sea level, are here used to simulate this hypothesised sequence of sandbank evolution over the Holocene. The results of the model experiments are corroborated by in situ observations of bedforms and sediment characteristics, and by acoustic Doppler current profiler (ADCP) data applied to predictions of bedload transport over the sandbanks. In addition to demonstrating the mechanism which leads to the formation of sandbanks by tidal flow through a strait, the model results show that upon subsequent closure of such a strait, these sandbanks will no longer be actively maintained, in contrast to sandbanks which are continuously maintained by headland eddies.     

Space-time variability of the Amazon River plume based on satellite ocean color

Satellite ocean color images were used to determine the space-time variability of the Amazon River plume from 2000–2004. The relationship between sea-surface salinity (SSS) and the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) absorption coefficient for dissolved and detrital material (a_dg) (r²=0.76, n=30, rmse=0.4) was used to identify the Amazon River plume low-salinity waters (<34 psu). The plume's spatial information was extracted from satellite bi-weekly time series using two metrics: plume area and plume shape. These metrics identified the seasonal variability of plume dimensions and dispersion patterns. During the study period, the plume showed the largest areas from July to August and the smallest from December to January. The mean annual amplitude and the mean, maximum and minimum plume areas were 1020×10³ km², 680×10³ km², 1506×10³ km² and 268×10³ km², respectively. Three main shapes and dispersion pattern periods were identified: (1) flow to the northeastern South American coast, in a narrow band adjacent to the continental shelf, from January to April; (2) flow to the Caribbean region, from April to July; and (3) flow to the Central Equatorial Atlantic Ocean, from August to December. Cross-correlation techniques were used to quantify the relationship between the plume's spatial variability and environmental forcing factors, including Amazon River discharge, wind field and ocean currents. The results showed that (1) river discharge is the main factor influencing plume area variability, (2) the wind field regulates the plume's northwestward flow velocity and residence time near the river mouth, and (3) surface currents have a strong influence over river plume dispersion patterns.     

Connectivity of the Apalachicola River flow variability and the physical and bio-optical oceanic properties of the northern West Florida Shelf

Maps of satellite-derived estimates of monthly averaged chlorophyll a concentration over the northern West Florida Shelf show interannual variations concentrated near the coastline, but also extending offshore over the shelf in a tongue-like pattern from the Apalachicola River during the late winter and early spring. These anomalies are significantly correlated with interannual variability in the flow rate of the Apalachicola River, which is linked to the precipitation anomalies over the watershed, over a region extending 150–200 km offshore out to roughly the 100 m isobath. This study examines the variability of the Apalachicola River and its impacts on the variability of water properties over the northern West Florida Shelf. A series of numerical model experiments show that episodic wind-driven offshore transport of the Apalachicola River plume is a likely physical mechanism for connecting the variability of the river discharge with oceanic variability over the middle and outer shelf.