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.     

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.     

Estimation of ecologically significant circulation features of the Columbia River estuary and plume using a reduced-dimension Kalman filter

A data assimilation method was used to estimate the variability of three ecologically significant features of the Columbia River estuary and plume: the size of the plume, the orientation of the plume, and the length of the salinity intrusion in the estuary. Our data assimilation method was based on a reduced-dimension Kalman filter that enables fast data assimilation of nonlinear dynamics in the estuary and plume. The assimilated data included measurements of salinity, temperature, and water levels at 13 stations in the estuary and at five moorings in the plume.     

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.     

A buoyant plume adjacent to a headland—Observations of the Elwha River plume

Small rivers commonly discharge into coastal settings with topographic complexities - such as headlands and islands - but these settings are underrepresented in river plume studies compared to more simplified, straight coasts. The Elwha River provides a unique opportunity to study the effects of coastal topography on a buoyant plume, because it discharges into the Strait of Juan de Fuca on the western side of its deltaic headland. Here we show that this headland induces flow separation and transient eddies in the tidally dominated currents (O(100 cm/s)), consistent with other headlands in oscillatory flow. These flow conditions are observed to strongly influence the buoyant river plume, as predicted by the “small-scale” or “narrow” dynamical classification using Garvine's (1995) system. Because of the transient eddies and the location of the river mouth on the headland, flow immediately offshore of the river mouth is directed eastward twice as frequently as it is westward. This results in a buoyant plume that is much more frequently “bent over” toward the east than the west. During bent over plume conditions, the plume was attached to the eastern shoreline while having a distinct, cuspate front along its westernmost boundary. The location of the front was found to be related to the magnitude and direction of local flow during the preceding O(1 h), and increases in alongshore flow resulted in deeper freshwater mixing, stronger baroclinic anomalies, and stronger hugging of the coast. During bent over plume conditions, we observed significant convergence of river plume water toward the frontal boundary within 1 km of the river mouth. These results show how coastal topography can strongly influence buoyant plume behavior, and they should assist with understanding of initial coastal sediment dispersal pathways from the Elwha River during a pending dam removal project.     

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.     

Helium and neon isotopes in the Imnaha Basalt, Columbia River Basalt Group: Evidence for a Yellowstone plume source

Helium, neon, and argon isotopic compositions were measured in two flows of the Columbia River flood basalt. The Imnaha Basalt has a ³He/⁴He ratio of 11.4 times atmospheric and ²⁰Ne/²²Ne and ²¹Ne/²²Ne ratios characteristic of a plume component. The measured ³He/⁴He is a lower limit, due to possible preferential ³He loss and/or addition of radiogenic ⁴He. A Wanapum Basalt flow, erupted approximately 2 Ma later in the waning stages of volcanism, has more MORB-like noble gases. The He, Nd and Sr isotopic compositions of these lavas suggest that the Columbia River basalts were derived from the Yellowstone plume head which contained both ‘high-helium’ plume material and entrained depleted mantle. As the eruptions progressed the plume component in the melting region was gradually diluted or replaced.     

Some aspects of the hydroclimatology of the Quesnel River Basin,British Columbia,Canada

In conjunction with available climate data, surface runoff is investigated at 12 gauges in the Quesnel watershed of British Columbia to develop its long‐term (1926–2004) hydroclimatology. At Quesnel itself, annual mean values of air temperature, precipitation and runoff are 4·6 °C, 517 and 648 mm, respectively. Climate data reveal increases in precipitation, no significant trend in mean annual air temperature, but an increasing trend in mean minimum temperatures that is greatest in winter. There is some evidence of decreases in winter snow depth. On the water year scale (October–September), a strong positive correlation is found between discharge and precipitation (r = 0·70, p < 0·01) and a weak negative correlation is found between precipitation and temperature (r = ? 0·36, p < 0·01). Long‐term trends using the Mann‐Kendall test indicate increasing annual discharge amounts that vary from 8 to 14% (12% for the Quesnel River, p = 0·03), and also a tendency toward an earlier spring freshet. River runoff increases at a rate of 1·26 mm yr^?1 m^?1 of elevation from west to east along the strong elevation gradient in the basin. Discharge, temperature and precipitation are correlated with the large‐scale climate indices of the Pacific Decadal Oscillation (PDO) and El‐Niño Southern Oscillation (ENSO). Copyright © 2009 John Wiley & Sons, Ltd.     

Modeling the response of primary production and sedimentation to variable nitrate loading in the Mississippi River plume

Increases in nitrate loading to the Mississippi River watershed during the last 50 years are considered responsible for the increase in hypoxic zone size in Louisiana-Texas shelf bottom waters. There is currently a national mandate to decrease the size of the hypoxic zone to 5000 km² by 2015, mostly by a 30% reduction in annual nitrogen discharge into the Gulf of Mexico. We developed an ecosystem model for the Mississippi River plume to investigate the response of organic matter production and sedimentation to variable nitrate loading. The nitrogen-based model consisted of nine compartments (nitrate, ammonium, labile dissolved organic nitrogen, bacteria, small phytoplankton, diatoms, micro- and mesozooplankton, and detritus), and was developed for the spring season, when sedimentation of organic matter from plume surface waters is considered important in the development of shelf hypoxia. The model was forced by physical parameters specified along the river-ocean salinity gradient, including residence time, light attenuation by dissolved and particulate matter, mixed layer depth, and dilution. The model was developed using measurements of biological biomasses and nutrient concentrations across the salinity gradient, and model validation was performed with an independent dataset of primary production measurements for different riverine NO₃ loads. Based on simulations over the range of observed springtime NO₃ loads, small phytoplankton contributed on average 80% to primary production for intermediate to high salinities (>15), and the main contributors to modeled sedimentation at these salinities were diatom sinking, microzooplankton egestion, and small phytoplankton mortality. We investigated the impact of limiting factors on the relationship between NO₃ loading and ecosystem rates. Model results showed that primary production was primarily limited by physical dilution of NO₃, followed by abiotic light attenuation, light attenuation due to mixing, and diatom sinking. Sedimentation was mainly limited by the first three of these factors. Neither zooplankton grazing or plume residence times acted as limiting factors of ecosystem rates. Regarding nutrient reductions to the watershed, simulations showed that about half of the percent decrease in NO₃ load was reflected in decreased plume sedimentation. For example, a 30% decrease in NO₃ load resulted in a 19% decrease in average plume primary production and a 14% decrease in sedimentation. Finally, our model results indicated that the fraction of primary production exported from surface waters is highly variable with salinity (7–87%), a finding which has important implications for predictive models of hypoxic zone size that assume a constant value for this ratio.     

Estimation of Columbia River virgin flow: 1879 to 1928

The Columbia River is a major source of and conduit for Pacific Northwest economic activity, and is one of the more heavily modified rivers in North America. Understanding human and climate‐induced changes in its hydrologic properties is, therefore, vital. Long streamflow records are essential to determining how runoff has changed over time, and Columbia River daily streamflow record at The Dalles began in 1878. To understand and separate anthropogenic and climate effects, however, it is also necessary to have a basin‐scale estimate of virgin or naturalized flow. The United States Geological Survey has calculated a monthly averaged adjusted river flow at The Dalles for 1879–1999 that accounts for the effects of flow regulation. The Bonneville Power Administration has estimated the monthly averaged virgin flow at The Dalles, i.e. the flow in the absence of both flow regulation and irrigation depletion for 1929–89. We have estimated the monthly virgin flow of the Columbia River at The Dalles from records of irrigated area for the missing early years, i.e. for the period 1879–1928. In addition, to allow hindcasting of a virgin flow sediment transport for the system, a daily virgin flow index with realistic higher moments and spectral properties has been calculated. Examination of the virgin flow record shows that climate change since the late 19th century has decreased annual average flow volume by > 7%; irrigation depletion has reduced the flow by another ∼7%. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Estimating turbidity and total suspended matter in the Adour River plume (South Bay of Biscay) using MODIS 250-m imagery

The Basque coastal waters (South Bay of Biscay) are directly influenced by the Adour River freshwater plume. The Adour outflow leads to important variations of suspended matter concentrations and turbidity, which in turn may affect biological productivity and water quality. This study aims at both developing specific algorithms and testing the efficiency of atmospherically corrected MODIS-Aqua 250-m surface reflectance product (MYD09) to map total suspended matter concentrations and turbidity within the Adour coastal region. First, regional empirical algorithms based on in-situ data were tested to retrieve the concentration of total suspended matter and turbidity from the remote sensing reflectance. Then, the respective sensitivity of MODIS surface reflectance bands 1 and 2 for water quality application was investigated as well as the quality of atmospheric corrections. Finally, selected algorithms were applied to the MYD09 product. The resulting 250-m resolution maps were then compared to 1000-m maps produced by IFREMER and comparisons between satellite measurements and in-situ sampling points were performed. Results show that MODIS-Aqua band 1 (620–670 nm) is appropriate for predicting turbidity and total suspended matter concentrations using polynomial regression models, whilst band 2 is unadapted. Comparison between total suspended matter concentration 250-m resolution maps and mineral suspended matter 1000-m maps (generated by IFREMER) produced consistent results. A high correlation was obtained between turbidity measured in-situ and turbidity retrieved from MODIS-Aqua satellite data.     

水平螺旋度与沙尘暴的动力学关系研究

应用2002～2006年高空流场和地面观测资料,计算近地面至500 hPa水平螺旋度的大小.结果表明：水平螺旋度负值中心值越大,500 hPa到近地面风速越大,西风增强,风速垂直切变越大,辐合上升运动越强,形成沙尘暴的强度就越强.水平螺旋度负值中心常常在河西走廊附近最强,导致其下游东南方发生沙尘暴.水平螺旋度负值中心与其下游沙尘暴发生强度有一致的对应关系：当水平螺旋度负值中心≤-200 m2/s2时,未来24小时内该区下游将有沙尘天气出现,当≤-600 m2/s2时,6小时内该区下游将有能见度小于500 m的强沙尘暴天气出现,当≤-1000 m2/s2时,6小时内该区下游将有能见度小于50 m的特强沙尘暴天气出现.     

Distinguishing human and climate influences on hydrological disturbance processes in the Columbia River,USA

Abstract

This paper distinguishes human and climate influences on the Columbia River streamflow disturbance regime, examines how this disturbance regime has changed over the last 150 years, and discusses downstream impacts. Flow management and withdrawal have greatly curtailed exceedence of the natural bankfull level of ～20 000 m³ s^-1. The frequency distribution of Columbia River flow has also changed. Sediment transport is positively correlated with streamflow standard deviation, and has been greatly reduced by flow regulation. Three kinds of spring freshet style have been identified; there are also three kinds of winter freshet. Flow regulation and regional climate warming have changed freshet styles and reduced maximum flows during the spring season. Downstream effects of hydrological alterations include increased salinity intrusion length, loss of shallow water habitat area during the freshet season, increased tides throughout most of the year, and a decrease in area of the Columbia River plume during spring and summer. Although climate changes and variations have played a substantial role in changing the hydrological disturbance regime, their influence is still less than that of human manipulation of the flow cycle.

Citation Jay, D. A. & Naik, P. K. (2011) Distinguishing human and climate influences on hydrological disturbance processes in the Columbia River, USA. Hydrol. Sci. J. 56(7), 1186–1209.     

A multi‐species dendroclimatic reconstruction of Chilko River streamflow,British Columbia,Canada

Dendroclimatological data were used to reconstruct the discharge history of Chilko River, which drains a glacierized watershed in the Coast Mountains of British Columbia. We correlated ring‐width records from Engelmann spruce (ES) (Picea engelmanni) and mountain hemlock (MH) (Tsuga mertensiana) trees to historical hydroclimate data. Over the period of record, spruce and hemlock radial growth correlates significantly with temperature and snow depth, respectively. We found that a multi‐species approach provided a better model fit and reconstructive power. Using these relationships, we developed generalized linear models for mean June, July, and June‐July discharge. The proxy records provide insights into streamflow variability of a typical Coast Mountains river over the past 240 years and confirm the long‐term influence of the Pacific Decadal Oscillation (PDO) on hydroclimatic regimes in the region. A relationship also exists between the reconstructed June‐July discharge record and the North Pacific (NP) Index, suggesting that winter atmospheric patterns over the North Pacific influence the hydrology of coastal British Columbia. Copyright © 2010 John Wiley & Sons, Ltd.     

Tidal effects on estuarine circulation and outflow plume in the Chesapeake Bay

The response of the Chesapeake Bay to river discharge under the influence and absence of tide is simulated with a numerical model. Four numerical experiments are examined: (1) response to river discharge only; (2) response to river discharge plus an ambient coastal current along the shelf outside the bay; (3) response to river discharge and tidal forcing; and (4) response to river discharge, tidal forcing, and ambient coastal current. The general salinity distribution in the four cases is similar to observations inside the bay. Observed features, such as low salinity in the western side of the bay, are consistent in model results. Also, a typical estuarine circulation with seaward current in the upper layer and landward current in the lower layer is obtained in the four cases. The two cases without tide produce stronger subtidal currents than the cases with tide owing to greater frictional effects in the cases with tide. Differences in salinity distributions among the four cases appear mostly outside the bay in terms of the outflow plume structure. The two cases without tide produce an upstream (as in a Kelvin wave sense) or northward branch of the outflow plume, while the cases with tide produce an expected downstream or southward plume. Increased friction in the cases with tide changes the vertical structure of outflow at the entrance to the bay and induces large horizontal variations in the exchange flow. Consequently, the outflow from the bay is more influenced by the bottom than in the cases without tide. Therefore, a tendency for a bottom-advected plume appears in the cases with tide, rather than a surface-advected plume, which develops in the cases without tide. Further analysis shows that the tidal current favors a salt balance between the horizontal and vertical advection of salinity around the plume and hinders the upstream expansion of the plume outside the bay.     

Case studies on mesoscale structures of heavy rainfall system in the Yangtze River generated by Meiyu front

The heavy rainfall generated by Meiyu front in the Yangtze River is popular topics in recent years. The heavy rainfall is formed in advantageous weather con-dition and interaction between different mesoscale processes. The previous studies included formation and maintenance of heavy rainfall, atmospheric circu-lation, relationship between the Meiyu frontal and mesoscale disturbance[1]. In recent years, the numer-ous studies and filed experiments are focus at the mesoscale structures of heavy…     

Heat transfer between a thermochemical plume channel and the surrounding mantle in the presence of horizontal mantle flow

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.     

Detection of anthropogenic gadolinium in the Brisbane River plume in Moreton Bay, Queensland, Australia

Wastewater effluent is known to contain macro and micropollutants, which may be deleterious to environmental health. One such class of micropollutants is chelated gadolinium, which are used as MRI contrast agents. As these MRI contrast agents can be assumed to behave conservatively during estuarine mixing, it is possible to calculate how much wastewater is represented in any particular sample. In this study, the percentage contribution of wastewater at specific locations in Moreton Bay, Qld, were determined by calculating the additional anthropogenic gadolinium contribution to the total rare earth element concentrations. Wastewater contributions were measured at concentrations as low as 0.2%, demonstrating the applicability of this technique for wastewater effluent plume mapping.