Metal exposure and biological responses in resident and transplanted blue mussels (Mytilus edulis) from the Scheldt estuary

The Western Scheldt river and estuary is known to be highly polluted as it receives industrial, agricultural and domestic effluents from one of the most populated and industrialised zones in Europe. Aquatic organisms are exposed to pollutants, specifically metals that are present in different environmental phases, e.g. dissolved, suspended material or sediment phases. The objective of this study was to study the relationship that exists between environmental metal levels, the degree of metal uptake by aquatic organisms with the concomitant biological responses. For this purpose the bivalve mollusk, Mytilus edulis, was selected as bioaccumulation indicator organism. Environmental (water and sediment) and mussel samples were collected during the late winter (March 2000) from sampling sites in the Scheldt estuary. Sites were selected to represent a salinity gradient from the mouth of the estuary to the furthest distribution area of mussels upstream in the system. Together with standard water quality parameters (e.g. salinity, dissolved oxygen, dissolved organic carbon, etc.) concentrations of twelve metals were analysed in the water (dissolved and suspended matter) and sediments. Levels of these metals were also measured in the soft tissue of M. edulis, together with concomitant biomarker responses in resident mussel populations at each site. The biomarkers that were included in this study were condition index, scope for growth, survival in air, cell membrane stability, DNA damage, and metallothioneins. Data were subjected to multivariate statistical analysis. The physico-chemical parameters and metals in the environmental samples clustered the sites to reflect the distribution based on the salinity gradient. Bioaccumulation results revealed increased metal uptake along a pollution gradient with the highest metal bioaccumulation occurring at the upstream most sites and therefore closest to the in the industrial activities. However, the biomarker responses clustered the sites in a manner that reflected the influence of combination of internal exposure (bioaccumulation) and external exposure (physico-chemical conditions). These differences in biomarker responses clearly demonstrated were attributed to abiotic factors other than metal pollution alone e.g. localized short-term increases in increased suspended sediment concentrations and decreased dissolved oxygen concentrations.     

Double-averaged velocity profiles over fixed dune shapes

Spatially averaged profiles of time averaged velocity, using integrals over thin horizontal slabs (Cartesian double average), are employed in characterizing the flow over fixed dune shapes. For comparison the spatial averaging method of Smith and McLean (1977) that averages along lines at constant distance from the local bed elevation is also investigated. The Cartesian double averaged profiles of the inverse of the velocity shear are nearly constant below the crest elevation, but increase rapidly above the crest level. This results in velocity profiles that increase linearly with distance from the bed below the crest. Above the crest it can be argued that the velocity increases logarithmically, but a power law profile can also be argued. Spatially averaged eddy viscosity profiles are calculated by multiplying the average Reynolds stress by the inverse shear. The resulting profile is more complex than the uniform flow counterpart.     

Power-law index for velocity profiles in open channel flows

Turbulence theory has demonstrated that the log law is one of the established theoretical results for describing velocity profiles, which is in principle applicable for the near-bed overlap region, being less than about 20% of the flow depth. In comparison, the power law that is often presented in an empirical fashion could apply to larger fraction of the flow domain. However, limited information is available for evaluating the power-law exponent or index. This paper attempts to show that the power law can be derived as a first-order approximation to the log law, and its power-law index is computed as a function of the Reynolds number as well as the relative roughness height. The result obtained also coincides with the fact that the one-sixth power included in the Manning equation is of prevalent acceptance, while higher indexes would be required for flows over very rough boundaries.     

Levels and profiles of PCBs and OCPs in marine benthic species from the Belgian North Sea and the Western Scheldt Estuary

Various benthic invertebrates (flying crab, common shrimp, and red starfish), small fish (sand goby), benthic flatfish (dab, plaice, and sole) and gadoids (bib and whiting) were collected in the Belgian North Sea and along the Scheldt Estuary, both representing areas impacted by various contaminants to different degrees. The levels of 25 polychlorinated biphenyls (PCBs) and 15 organochlorine pesticides (OCPs), which included penta- and hexachlorobenzene, alpha-, beta-, and gamma-hexachlorocyclohexane isomers, chlordanes, and DDT and metabolites, were determined. Sum of PCBs and OCPs in benthic invertebrates and goby ranged from 1.5 to 280 ng/g wet weight (ww) and from 0.27 to 23 ng/g ww, respectively. The fish livers revealed total PCB and OCP levels ranging from 20 to 3200 ng/g ww and from 6.0 to 410 ng/g ww, respectively. Levels of both contaminant groups were significantly higher in samples from the Scheldt Estuary compared to the Belgian North Sea. For most species a highly inverse correlation was found between the concentration of contaminants and the distance to Antwerp (r between 0.812 and 0.901, p < 0.05), pointing to a higher degree of exposure further upstream. PCB and OCP exposures are highly correlated (r between 0.836 and 1.000, p < 0.05), which suggests that the pollution can be classified as historical. However, because urban and industrial centres may still be emitting these compounds, more recent point and non-point sources cannot be ruled out.     

Estimation of groundwater velocity in localized fracture zones from well temperature profiles

In fractured rocks, well temperature logs often exhibit abrupt temperature changes over a small distance, typically surging or dropping several degrees over several tens of meters. The abrupt temperature changes usually occur in localized fracture zones or small faults. Away from these zones, temperature quickly restores itself to the background linear conductive profile. A theoretical model is presented herein to interpret the pattern of such abrupt temperature changes. I illustrate that this type of temperature profile originates from a different physical process than the physical processes found in some other typical temperature patterns in large scales. I also demonstrate that, using the field data, the theoretical model can be effectively used to estimate fluid flow velocity in fracture zones or local faults. Using the temperature profiles from two boreholes, fluid flow velocity in fracture zones is estimated to be as high as 1.1×10⁻⁶ m/s. Although fracture flow is a highly localized feature, it can reach three orders of magnitude higher than fluid flow in the hosting rock matrices.     

Estimation of shear wave velocity profiles by the inversion of spatial autocorrelation coefficients

The subsurface shear-wave velocity (Vs) is considered to be a key parameter for site characterization and assessment of earthquake hazard because of its great influence on local ground-motion amplification. Array microtremor measurements are widely used for the estimation of shear-wave velocities. Compared to other methods such as frequency-wavenumber (f-k) methods, the spatial autocorrelation (SPAC) method requires fewer sensors and thus is relatively easier to implement and gives robust estimations of shear-wave velocity profiles for depths down to a few hundred meters. The quantity derived from observed data is the SPAC coefficient, which is a function of correlation distance, frequency and phase velocity. Generally, estimation of Vs profiles is a two stage process: Estimation of the dispersion data from the SPAC coefficients and inversion of the dispersion data for shear-wave velocity structure. In this study, instead of inverting dispersion curves, a more practical approach is used; that is, observed SPAC coefficients are directly inverted for the S-wave velocities. A synthetic case and a field data application are presented to test the potential of the inversion algorithm. We obtain an iterative damped least-squares solution with differential smoothing. The differential smoothing approach constrains the change in shear-wave velocities of the adjacent layers and thus stabilizes the inversion.     

Flow time estimation with spatially variable hillslope velocity in ungauged basins

The automated spatial estimation of the hillslope runoff dynamics is used as a valuable tool for the estimation of the travel time distribution (flow time), a major factor for the hydrologic prediction in ungauged basins. In fact, while the flow time function is usually obtained by rescaling the flow paths with constant channel and hillslope velocities, in this work a spatially distributed kinematic component, as a function of terrain features and in particular slope and land use, is implemented and its influence on the hydrologic response is tested by means of the Width Function Instantaneous Unit Hydrograph (WFIUH) framework. Hillslope surface flow velocities are evaluated by applying different uniform flow formulas within an automated DEM-based (terrain analysis) algorithm. A comparison test of the performances of the Manning, Darcy, Maidment and Soil Conservation Service uniform flow equations is performed for several case studies in Italy pertaining to different climatic and geomorphic conditions. Results provide new insights for a better understanding of the flow time function also introducing a more parsimonious and physically-based calibration scheme of the WFIUH.     

Flow velocity and primary production influences carbon utilization in nascent epilithic stream biofilms

In small streams, the majority of carbon turnover is due to microbial activity in biofilms. Flow velocity is a key factor influencing biofilm function, and nascent biofilms with high energy need for growth might be especially sensitive to hydrodynamics. The major part of carbon supply is allochthonous but algae can provide easily available exudates for biofilm bacteria. In this study, epilithic biofilms were grown for 2 weeks in a third order stream in Thuringia, Germany, and then incubated in replicate flow channels in climate-controlled chambers. Glucose and arbinose were added immediately to all channels, and the effects of flow velocity and light availability on rates of sugar removal were examined. Phosphate addition did not influence sugar decrease rates. Flow velocities of either 0.3 m s⁻¹ or 0.7 m s⁻¹ resulted in 1.3 to 3.1 times higher decrease rates under the higher flow velocity. Light exclusion resulted in a 2.2 to 2.6 times faster sugar decrease but only a 0.5 times dissolved organic carbon increase compared to channels with light input, suggesting a strong internal coupling of primary producers and heterotrophs. Our results indicate that carbon uptake from the water column is fostered at higher flow velocities and that primary production is an important internal carbon source in nascent epilithic biofilms.     

Barium build-up in the Teign estuary

Accumulations of barium occur in the River Teign, in south-west England, and its estuary, largely in the form of barytes. These concentrations are derived from weathering in the Teign Valley orefield, an area of mineralization scattered with mine dumps. It appears that although barium values may increase in the upper confines of the estuary the net effect on lower reaches will be negligible.     

The influence of changes in tidal asymmetry on residual sediment transport in the Western Scheldt

This study fits into a wider research program from RWS RIKZ concerning the exchange of sediment between the coast and the tidal basins, Western Scheldt and Wadden Sea, which are the largest basins along the Dutch Coast, over different time-scales. For both basins, questions about the evolution of the import/export at the mouth recently arose. In case of the Western Scheldt, which is the subject of this study, mainly the uncertainty about the future developments after the change from import to export at the mouth was noticed in the 1990s, which necessitated a more detailed study of this area.     

Inversion of ambient seismic noise HVSR to evaluate velocity and structural models of the Lower Tagus Basin,Portugal

During its history, several significant earthquakes have shaken the Lower Tagus Valley (Portugal). These earthquakes were destructive; some strong earthquakes were produced by large ruptures in offshore structures located southwest of the Portuguese coastline, and other moderate earthquakes were produced by local faults. In recent years, several studies have successfully obtained strong-ground motion syntheses for the Lower Tagus Valley using the finite difference method. To confirm the velocity model of this sedimentary basin obtained from geophysical and geological data, we analysed the ambient seismic noise measurements by applying the horizontal to vertical spectral ratio (HVSR) method. This study reveals the dependence of the frequency and amplitude of the low-frequency (HVSR) peaks (0.2–2 Hz) on the sediment thickness. We have obtained the depth of the Cenozoic basement along a profile transversal to the basin by the inversion of these ratios, imposing constraints from seismic reflection, boreholes, seismic sounding and gravimetric and magnetic potentials. This technique enables us to improve the existing three-dimensional model of the Lower Tagus Valley structure. The improved model will be decisive for the improvement of strong motion predictions in the earthquake hazard analysis of this highly populated basin. The methodology discussed can be applied to any other sedimentary basin.     

Circulation and mixing in estuary outflows

The outflow of brackish waters from the Tees estuary on an ebb tide forms a plume which is initially spread by the action of gravity. Under quiescent conditions, the rate of vertical mixing of the brackish field is relatively slow and an area of low salinity may form in the surface waters as a result of the superimposition of the outflow plumes from successive tides. A field area of some 70 km² has been observed to be enclosed by the 34.2 × 10^?3 contour.Averaged over a tidal period, the estuary outflow from the Tees forms a broad plume of low-salinity water which extends to seaward of the 20-m contour. The movement of the field thereafter depends on the residual coastal drift which runs either to the northwest or to the southeast with a typical speed of 0.015 m s^?1. The tidally averaged vertical mixing coefficient in this offshore field was estimated to be about 10^?4 m^?2 s^?1, which is similar to that off the Wear estuary and an order of magnitude smaller than that calculated for the waters off the Humber.     

Polychlorinated biphenyls in the Brisbane River estuary, Australia

The concentrations of PCBs in water, sediments, fish (6 species), crustacea (3 species), molluscs (1 species), polychaetes (mixed species) and birds (1 species) from the Brisbane River estuary were determined. Levels were comparable to other similar situations elsewhere. No relationship was found between trophic class and PCB concentration although pelicans at the highest trophic level examined exhibited the highest muscle tissue concentration (i.e., 15.7 mg kg⁻¹ on a wet weight basis). The pattern of distribution of PCBs within body tissues of fish was closely related to lipid content.     

Intra-annual variability in nutrients in the Godavari estuary,India

Daily variations in nutrients were monitored for 15 months (September 2007–November 2008) in the Godavari estuary, Andhra Pradesh, India, at two fixed locations. River discharge has significant influence on nutrients loading to the estuary, which peaks during June–August (peak discharge period; monsoon) whereas exchanges at the sediment–water interface, groundwater and rainwater contribute significantly during other period. Despite significant amount of nutrients brought by discharge to the study region, phytoplankton biomass, in terms of chlorophyll-a (Chl a), did not increase significantly due to high suspended load and shallow photic depth. Nutrients showed downward gradient towards downstream of the estuary from upstream due to dilution by nutrient poor seawater and biological uptake. The N:P ratios were higher than Redfield ratio in both upstream and downstream of the estuary during no discharge period suggesting PO₄ to be a limiting nutrient for phytoplankton production, at levels <0.10 μmol L⁻¹. On the other hand, Si:N ratios were always more than unity during entire study period at both the stations indicating that Si(OH)₄ is not a limiting nutrient. Our results suggest that suspended matter limits phytoplankton biomass during peak discharge period whereas PO₄ during no discharge period.     

Modeling flocculation in a hypertidal estuary

When fine particles are involved, cohesive properties of sediment can result in flocculation and significantly complicate sediment process studies. We combine data from field observations and state-of-the-art modeling to investigate and predict flocculation processes within a hypertidal estuary. The study site is the Welsh Channel located at the entrance of the Dee Estuary in Liverpool Bay. Field data consist of measurements from a fixed site deployment during 12–22 February 2008. Grain size, suspended sediment volume concentration, and current velocity were obtained hourly from moored instruments at 1.5 m above bed. Near-bottom water samples taken every hour from a research vessel are used to convert volume concentrations to mass concentrations for the moored measurements. We use the hydrodynamic model Proudman Oceanographic Laboratory Coastal Ocean Modelling System (POLCOMS) coupled with the turbulence model General Ocean Turbulence Model (GOTM) and a sediment module to obtain three-dimensional distributions of suspended particulate matter (SPM). Flocculation is identified by changes in grain size. Small flocs were found during flood and ebb periods—and correlate with strong currents—due to breakup, while coarse flocs were present during slack waters because of aggregation. A fractal number of 2.4 is found for the study site. Turbulent stresses and particle settling velocities are estimated and are found to be related via an exponential function. The result is a simple semiempirical formulation for the fall velocity of the particles solely depending on turbulent stresses. The formula is implemented in the full three-dimensional model to represent changes in particle size due to flocculation processes. Predictions from the model are in agreement with observations for both settling velocity and SPM. The SPM fortnight variability was reproduced by the model and the concentration peaks are almost in phase with those from field data.