Calculating ecological carrying capacity of shellfish aquaculture using mass-balance modeling: Narragansett Bay, Rhode Island

Increasing growth in the aquaculture industry demands ecosystem-based techniques for management if that growth is to be ecologically sustainable and promote equity among users of the ecosystems in which it occurs. Models of carrying capacity can be used to responsibly limit the growth of aquaculture in increasingly crowded coastal areas. Narragansett Bay, Rhode Island, USA is one such crowded coastal region experiencing a rapid increase in bivalve aquaculture. An ecosystem mass-balance model was used to calculate the ecological carrying capacity of bivalve aquaculture. Cultured oyster biomass is currently at 0.47 t km⁻² and could be increased 625 times without exceeding the ecological carrying capacity of 297 t km⁻². This translates to approximately 38,950 t of harvested cultured oysters annually which is 4 times the total estimated annual harvest of finfish. This potential for growth is due to the high primary productivity and large energy throughput to detritus of this ecosystem. Shellfish aquaculture has potential for continued growth and is unlikely to become food limited due, in part, to the large detritus pool.     

Feeding behaviour of the sea cucumber <Emphasis Type="Italic">Cucumaria frondosa</Emphasis> (Echinodermata: Holothuroidea) in the laboratory and the field: relationships between tentacle insertion rate,flow speed,and ingestion

The sea cucumber, Cucumaria frondosa, is a benthic suspension feeder that captures food particles on its tentacles and then inserts them into its mouth one at a time. Previous studies have suggested that tentacle insertion rate (TIR) could be a useful indicator of food intake. The present study determined whether flow velocity affects TIR and whether TIR is a good indicator of ingestion. Video observations of sea cucumbers in Passamaquoddy Bay (45°01.70N, 66°55.74W) in August 1995 showed that TIRs increased with velocities up to 55 cm s⁻¹ and decreased steadily at flows above that up to 130 cm s⁻¹. In October 2006, laboratory flume studies were carried out on specimens collected from the same site in the previous August. Temperature and salinity (12°C and 32) in the flume were the same as in the field at the time of collection. There was high individual variation in feeding behavior at free-stream velocities of 4–40 cm s⁻¹ and TIR was independent of flow. As the number of tentacle insertions increased in the flume experiments, the amount of chloropigments in the digestive tracts of the sea cucumbers also increased. This suggests that TIR, which can be measured non-intrusively using remote video techniques, could be a good indicator of feeding behavior and ingestion in C. frondosa.     

Rookery Bay and Naples Bay circulation simulations: Applications to tides and fresh water inflow regulation

A three-dimensional, baroclinic numerical circulation model is applied to the combined Rookery Bay and Naples Bay estuary complex. We focus on two aspects of the model simulations: the tidal circulation, for which we have observations of sea level and currents at several stations spanning Rookery Bay, and the salinity variations, for which we have observations at one station. After establishing model veracity through quantitative comparisons with observations, we use the model to address salinity changes that could arise due to regulatory fresh water flow rate modifications. Two different sets of experiments are performed. The first considers the salinity changes in Rookery Bay by increasing the fresh water flow rates into Henderson Creek through the addition of either 1.4 m³/s or 2.8 m³/s to the ambient flow rates. The second considers the salinity changes in Naples Bay by decreasing the fresh water flow rates through the subtraction of 9.9 m³/s from the Golden Gate Canals’ ambient flow rates. All of these scenarios result in significant, quantifiable salinity changes within the Rookery Bay and Naples Bay estuary complexes, demonstrating the environmental management utility of using high resolution, three-dimensional, numerical circulation models, as applied herein, for assessing the water properties of complex, multi-connected estuarine systems.     

桑沟湾溶解态汞的生物地球化学行为

作为一种毒性极强的重金属元素,汞在海洋中可以通过海洋生物的呼吸、摄食和吸附等过程被利用,并在沿食物链传递的过程中不断富集,最终危害人类的健康。利用冷原子荧光光度法(CV-AFS)对2011年4月、8月、10月和2012年1月桑沟湾溶解态无机汞(DIHg)和溶解态有机汞(DOHg)的含量进行了测定。结果表明,桑沟湾四个季节DIHg的浓度范围分别为52~865、131~359、31~134和22~119 pmol·L-1,DOHg的浓度范围分别为37~214、52~635、21~98和未检出~51 pmol·L-1。桑沟湾DIHg和DOHg具有相似的分布特征,均呈现出从近岸向外海逐渐降低的趋势,有明显的季节变化。桑沟湾DIHg和DOHg的周日变化与潮汐呈现出较好的负相关关系,且存在明显的昼夜差异。影响桑沟湾DIHg和DOHg分布的主要因素包括河流和地下水的输入、大气的干湿沉降、与黄海的交换、活性气态汞(Hg0)在海-气界面的交换、表层水体发生的光化学还原反应以及养殖生物的清除等。通过初步计算,桑沟湾溶解态汞(TDHg)的存留时间约为(1.27±0.53)年,远远低于大洋。根据美国国家环境保护局汞的质量标准和我国地表水环境质量标准,桑沟湾没有明显的汞污染。但海产品体内富集的汞可能会带来潜在的生态危机和食品安全问题,需要相关部门加以重视,确保桑沟湾养殖产业的平衡发展。     

A general process-based mass-balance model for phosphorus/eutrophication as a tool to estimate historical reference values for key bioindicators,as exemplified using data for the Gulf of Riga

This work describes how a general, process-based mass-balance model (CoastMab) for phosphorus for coastal areas may be used as a tool to estimate realistic values of “natural” or preindustrial reference levels of key bioindicators in coastal science, including the Secchi depth, a standard measure of water clarity, the chlorophyll-a concentration, an operational measure of phytoplankton biomass and the concentration of cyanobacteria, a measure of the concentration of harmful algae. The CoastMab-model is an ecosystem model giving monthly predictions to achieve seasonal variations of basin-wide properties. The selected case-study area, the Gulf of Riga, is sensitive to nutrient loading because of its shallowness and low openness towards the Baltic Proper. The morphometry of any coastal area, as given by the size and form parameters, influences all internal processes, such as sedimentation, resuspension, diffusion in water and from sediments to water, biouptake and retention in biota, stratification, mixing and outflow. There has been no mass-balance modeling for nitrogen (N) in this work because empirical data (from the HELCOM database) clearly indicate that the monthly primary production in the Gulf of Riga is regulated by phosphorus (P) – the mean monthly total-N to total-P ratios are well over 7.2 (the Redfield-ratio) and generally higher than 15 for the data used in this study (from 1992 to 2005). At present anthropogenic loads, the average modeled monthly values for Secchi depth, chlorophyll (Chl), cyanobacteria (CB) and total-P (TP) are 3.2 m, 3.8 μg/l, 78 μg/l and 31.3 μg/l, respectively. If 50% of all anthropogenic sources to the Gulf of Riga via rivers, point sources and diffuse sources were to be removed, these values would be 3.6 m, 3.4 μg Chl/l, 63 μg CB/l and 29.1 μg TP/l. If 60% of the anthropogenic phosphorus fluxes to the Baltic Proper were to be omitted and as well as 75% of all direct anthropogenic sources to the Gulf of Riga, the values would be 4.6 m, 2.7 μg Chl/l, 45 μg CB/l and 25.4 μg TP/l. These values represent the “natural” reference levels and it is not realistic to expect that remedial measures would improve the conditions more than that. Using the CoastWeb-model, similar calculations can be made for any given coastal area and the data necessary for such calculations are discussed in this work.     

Interactions between two introduced species: Zostera japonica (dwarf eelgrass) facilitates itself and reduces condition of Ruditapes philippinarum (Manila clam) on intertidal flats

Dwarf eelgrass (duckgrass; Zostera japonica) and Manila clams (Ruditapes philippinarum) are two introduced species that co-occur on intertidal flats of the northeast Pacific. Through factorial manipulation of clam (0, 62.5, 125 clams m⁻²) and eelgrass density (present, removed by hand, harrowed), we examined intra- and interspecific effects on performance, as well as modification of the physical environment. The presence of eelgrass reduced water flow by up to 40% and was also observed to retain water at low tide, which may ameliorate desiccation and explain why eelgrass grew faster in the presence of conspecifics (positive feedback). Although shell growth of small (20–50 mm) clams was not consistently affected by either treatment in this 2-month experiment, clam condition improved when eelgrass was removed. Reciprocally, clams at aquaculture densities had no effect on eelgrass growth, clam growth and condition, or porewater nutrients. Overall, only Z. japonica demonstrated strong population-level interactions. Interspecific results support an emerging paradigm that invasive marine ecosystem engineers often negatively affect infauna. Positive feedbacks for Z. japonica may characterize its intraspecific effects particularly at the stressful intertidal elevation of this study (+1 m above mean lower low water).     

Assessing biomass gains from marsh restoration in Delaware Bay using Ecopath with Ecosim

The Delaware Bay ecosystem has been the focus of extensive habitat restoration efforts to offset finfish losses due to mortality associated with power plant water intake. As a result, a 45 km² or a 3% increase in total marsh area was achieved by 1996-1997 through the restoration efforts of the Public Service Enterprise Group (PSEG). To quantify the impact of restoration efforts on system productivity, an Ecopath with Ecosim model was constructed that represented all major components of the ecosystem. The model consisted of 47 functional groups including: 27 fish species, 5 invertebrate groups, 4 multi-species benthic groups, 6 multi-species fish groups, 3 plankton groups, 1 shorebird group and 1 marine mammal group. Biomass, abundance, catch, and demographic data were obtained from the literature or from individual stock assessments conducted for principal ecosystem components. A base Ecosim model was fitted to time series of key species in the Bay representing the period 1966-2003. To access the gains from marsh restoration, model simulations reflecting no restoration were conducted to estimate the productivity that would have been lost if restoration efforts had not occurred. The results indicated that restoration increased total ecosystem biomass by 47.7 t km⁻² year⁻¹. Simulations indicated increased biomasses across a wide range of species including important forage and commercially important species. The marsh restoration also significantly impacted ecosystem structure increasing the ratio of production-to-respiration, increasing system path length and decreasing the ratio of production-to-biomass.     

Modelling the effects of ‘coastal’ acidification on copper speciation

We present here a copper speciation model that accounts for the long-term (‘coastal-acidification’) and short-term (daily and seasonal variation) variability in water pH and water temperature. The developed model is applied to a sub-tropical estuary (Moreton Bay, Australia) at a one hundred year time scale so that outputs are consistent with climate change projections. The model predicts that the mean cupric ion concentration (Cu²⁺) in the estuary will increase by 115% over the next 100 years as a result of the projected decrease in pH and increase in water temperature. Through calibration, the estimated concentration of copper-complexing dissolved organic matter (DOM) in the estuary is found to be 22.5 nM. An increase in the concentration of Cu²⁺, which is the most toxic and bioavailable form of copper, has implications for ecosystem health and may have a negative effect on the detoxifying capacity of DOM. Models that provide a framework for coupling biological, chemical and physical processes are important for providing a holistic perspective of coastal systems, especially for better understanding a system within the context of climatic and non-climatic drivers.     

Calculation of the unit emergy value of water in an Italian watershed

Emergy is a thermodynamics-based entity that enables the implementation of a holistic environmental accounting system. It contributes to identify and measure all the inputs (energy and matter) supporting a given system, expressed in a common unit, namely solar emergy joule (sej). The emergy per unit product (called unit emergy value, UEV), is a measure of the environmental cost of a given resource. It is specific of the system/process and gives information on the dynamics, components and functioning of it. This paper presents the emergy evaluation of water resources within the watershed of the river Sieve, located in the Province of Florence (Italy). Along the river, an artificial basin has been created by means of a dam to preserve water quantity and quality, and to protect the Florentine area from dangerous floods and inundations. Different UEVs of water can be identified along the course of the river, especially upstream and downstream of the dam. These values quantify both the environmental and human efforts made to ensure and regulate the presence of water at different points of the river. The UEV of water flowing in the river increases from 1.35 × 10⁵ sej/g upstream, to 5.80 × 10⁵ sej/g downstream of the dam, depending mainly on man-made infrastructure. Along the watershed, three different systems of extraction, purification and distribution of water have been chosen on the basis of their dimension, type and location. UEVs of water distributed and the emergy investment necessary to implement different water management strategies are presented. The value of water purified and distributed decreases from 2.00 × 10⁶ sej/g for the smallest plant in the mountainous area, to 1.72 × 10⁶ sej/g for the largest plant, in the city of Florence, depending on production efficiency.     

Modeling estuarine-shelf exchanges in a deltaic estuary: Implications for coastal carbon budgets and hypoxia

The export of wetland-derived materials to the coastal ocean (i.e., the “Outwelling” hypothesis) has received considerable attention over the past several decades. While a number of studies have shown that estuaries export appreciable amounts of nutrients and carbon, few studies have attempted to estimate the importance of estuarine sources for the coastal carbon budgets in river-dominated coastal ecosystems. A novel tidal prism model was developed to examine estuarine-shelf exchanges in the Barataria estuary, a deltaic estuary located in the north-central Gulf of Mexico. This estuary has been the site of a massive wetland loss, and it has been hypothesized that carbon export from the eroding coastal wetlands supports the development of a large hypoxic zone in the coastal Gulf of Mexico. The model results show that the Barataria estuary receives nitrogen through the tidal passes and releases carbon to the coastal ocean. The mean calculated tidal water discharge of 6930 m³ s⁻¹ is equivalent to about 43% of the lower Mississippi River discharge. The annual total organic carbon (TOC) export is 109 million kg, or 57 gC m² yr⁻¹ when prorated to the total water area of the estuary. This carbon export is equivalent to a loss of 0.5 m of wetland soil horizon over an area of 8.4 km², and accounts for about 34% of the observed annual wetland loss in the estuary between 1978 and 2000. Compared to the lower Mississippi River, the Barataria estuary appears to be a very small source of TOC for the northern Gulf of Mexico (2.7% of riverine TOC), and is unlikely to have a significant influence on the development of the Gulf's hypoxia.     

A system dynamics model for the management of the Manila clam, Ruditapes philippinarum (Adams and Reeve, 1850) in the Bay of Arcachon (France)

The Manila clam Ruditapes philippinarum (Adams and Reeve, 1850) is one of the mollusc species that, driven mainly by the shellfish market industry, has extended throughout the world, far beyond the limits of its original habitat. The Manila clam was introduced into France for aquaculture purposes, between 1972 and 1975. In France, this venerid culture became increasingly widespread and, since 1988, this species has colonised most of the embayments along the French Atlantic coast. In 2004, this development resulted in a fishery of ca. 520 t in Arcachon Bay.     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

Estimation of annual potential evapotranspiration at regional scale based on the effect of moisture on soil respiration

Potential evapotranspiration (PET) is an important component of water cycle. For traditional models derived from the principle of aerodynamics and the surface energy balance, its calculation always includes many parameters, such as net radiation, water vapor pressure, air temperature and wind speed. We found that it can be acquired in an easier way in specific regions. In this study, a new PET model (PETP model) derived from two empirical models of soil respiration was evaluated using the Penman-Monteith equation as a standard method. The results indicate that the PETP model estimation concur with the Penman-Monteith equation in sites where annual precipitation ranges from 717.71 mm to 1727.37 mm (R² = 0.68, p = 0.0002), but show large discrepancies in all sites (R² = 0.07, p = 0.1280). Then we applied our PETP model at the global scale to the regions with precipitation higher than 700 mm using 2.5° CMAP data to obtain the annual PET for 2006. As expected, the spatial pattern is satisfactory overall, with the highest PET values distributed in the lower latitudes or coastal regions, and with an average of 1292.60 ± 540.15 mm year⁻¹. This PETP model provides a convenient approach to estimate PET at regional scales.     

Local calibration of coliforms parameters of water quality problem at Igapó I Lake, Londrina, Paraná, Brazil

This article presents results concerning the local calibration of the transport parameters (longitudinal and transversal diffusions and decay coefficient) for a two-dimensional problem of water quality at Igapó I Lake, located in Londrina, Paraná, Brazil, using fecal coliforms as an indicator of water quality. The simulation of fecal coliforms concentrations all over the water body is conducted by means of a structured discretization of the geometry of Igapó I Lake, together with the finite difference and finite element methods. By using the velocity field, modeled by the Navier-Stokes and Poisson equations, the flow of fecal coliforms is described by means of a transport model, which considers advective and diffusive processes, as well as a process of fecal coliforms decay. In the checkpoint, the longitudinal and transversal diffusion coefficients and the coliforms decay coefficient that best fitted the value of the fecal coliforms concentration were D_x = D_y = 0.001 m²/h and k = 0.5 d⁻¹ = 0.02083 h⁻¹. A qualitative and quantitative analysis of the numerical simulations conducted in function of the diffusion coefficients and of the coliforms decay parameter provided a better understanding of the local water quality at Igapó I Lake.     

Potential of bivalves' secondary settlement differs with species: a comparison between cockle (Cerastoderma edule) and clam (Ruditapes philippinarum) juvenile resuspension

Juvenile bivalves may be dispersed by entering a bysso-pelagic phase where they drift through the water mass aided by a long thread. The ability to resuspend and control the specific weight in two bivalve species, the cockle Cerastoderma edule (L.) and the Japanese clam Ruditapes philippinarum (Adams and Reeves), was documented with juveniles through flume and still-water experiments. Cockle juveniles initially placed on an unsuitable substratum were exposed to two shear velocities (u _*). At the end of the experiment, 42 (±15)% (for u _*=0.51 cm s⁻¹) and 79 (±9)% (for u _*=0.99 cm s⁻¹) of individuals were retrieved from the sand area which represents only 8% of the total flume surface. Most juveniles (70.5%) with shell lengths <2.5 mm migrated from the unsuitable Plexiglas substratum to the sand array by resuspension in the water column. The percentage was lower (21.5%) for larger individuals. The same experimental design was applied to clams, which immediately adhered to the Plexiglas substratum and remained attach to it. Sinking rates of live and dead specimens of both species were measured in a 1 m long transparent PVC tube. Cockle fall velocities showed severe deceleration, probably due to byssus secretion (up to 15-fold slower than dead cockles), sometimes interrupted by brutal acceleration probably due to byssal rupture. Cockles were able to reduce their sinking rate for shell lengths up to 4.25 mm. By contrast, clam sinking rates were constant, and similar to dead clam sinking rates. Specific weights of all experimental juveniles were calculated in relation to their lengths, and their passive motion into the boundary layer was theoretically assessed with Shields curve. In short, C. edule and R. philippinarum can both exhibit dense populations in the field with a good capacity to colonize, although juveniles display different abilities to resuspend in the water column. Received: 27 January 1997 / Accepted: 13 February 1997     

Development of TRIPLEX-Management model for simulating the response of forest growth to pre-commercial thinning

In order to simulate forest growth response to pre-commercial thinning (PCT), TRIPLEX1.0 - a process-based model designed to predict forest growth as well as carbon (C) and nitrogen (N) dynamics - was modified and improved to also simulate managed forest ecosystem thinning practices. A three-parameter Weibull distribution model was integrated to simulate thinning treatments within the newly developed TRIPLEX-Management model. The thinning intensity component within the model allows users to simulate thinning treatments by applying basal area, stand density and volume to quantify thinning intensity. Natural mortality decreased following thinning due to an increase in growing space for residual stems. Predicted litterfall pools also increased after thinning events took place. The TRIPLEX-Management model was tested against published observational data for Jack Pine (Pinus banksiana Lamb.) stands subjected to PCT in Northwestern Ontario, Canada. The coefficients of determination (R²) between the predicted and observed variables including stand density, mean DBH (diameter at breast height), the quadratic mean DBH, total volume and merchantable volume as well as belowground, aboveground, and total biomass ranged from 0.50 to 0.88 (n = 20, P < 0.001) with the exception of mean tree height (R² = 0.25, n = 20, P < 0.05). Overall, the Willmott index of agreement between predicted and observed variables ranged from 0.97 to 1.00. Results show that the TRIPLEX-Management model is generally capable of simulating growth response to PCT for Jack Pine stands.     

Individual-based modeling of flooding and salinity effects on a coastal swamp forest

Coastal swamps are among the rapidly vanishing wetland habitats in Louisiana. Increased flooding, nutrient and sediment deprivation, and salt-water intrusion have been implicated as probable causes of the decline of coastal swamps. We developed a two-species individual-based forest succession model to compare the growth and composition of a cypress-tupelo swamp under various combinations of flooding intensity and salinity levels, using historical time-series of stage and salinity data as inputs. Our model simulates forest succession over 500 years by representing the growth, mortality, and reproduction of individual Taxodium distichum (baldcypress) and Nyssa aquatica (water tupelo) trees in a 1-km² spatial grid of 10 m × 10 m cells that vary in water levels and salinity through differences in elevation. We independently adjusted the elevations of each cell to obtain different grid-wide mean elevations and standard deviations of elevation; this affected the temporal and spatial pattern of flooding. We calibrated the model by adjusting selected parameters until averaged basal area, stem density and wood production rates under two different mean elevations (partially versus highly flooded) were qualitatively similar to comparable values reported for swamps in the literature. Corroboration involved comparing model predictions to four well-monitored contrasting habitat sites within the Maurepas Basin, Louisiana, USA. Model predictions of both species combined showed the same patterns among sites as the data, but the model overestimated wood production and the dominance of T. distichum. Exploratory simulations predicted that increased flooding leads to swamps with reduced basal areas and stem densities, while increased salinity resulted in lower basal areas at low salinity concentration (∼1-3 psu) and complete tree mortality at higher salinity concentrations (∼2-6 psu). Our model can provide insight into the succession dynamics of coastal swamps and information for the effective design of restoration actions.     

Assessing the carrying capacity of tilapia in an intertidal mangrove-based polyculture system of Pearl River Delta, China

A trophic model of an intertidal mangrove-based polyculture system in Pearl River Delta, China, was constructed using the Ecopath with Ecosim software. This polyculture system was chosen since it is the first integrated multi-trophic aquaculture (IMTA) system that was constructed on the basis of mangrove planting in China. The energy flows, ecosystem property, and carrying capacity of tilapia in the polyculture system were analyzed and evaluated. The results show the trophic level of 1.00 for primary producers and detritus to 2.85 for grass carp. The geometric mean of the trophic transfer efficiencies was 7.0%, with 7.2% from detritus and 6.8% from primary producers within the system. The ecosystem property indices show that this polyculture system has a high value of total primary production/total respiration (TPP/TR) and total primary production/total biomass (TPP/TB), together with low Finn's cycled index (FCI), Finn's mean path length (FML), and connectance index (CI), indicating that this system is at a development stage according to Odum's theory. The principal fish cultured in the system is tilapia, and mixed trophic impacts (MTI) show that tilapia has a marked impact on most compartments in this system, and the carrying capacity was found to be a tilapia culture biomass of 5.8 t ha⁻¹ in the system.     

Carbon, nitrogen, oxygen and sulfide budgets in the Black Sea: A biogeochemical model of the whole water column coupling the oxic and anoxic parts

Carbon, nitrogen, oxygen and sulfide budgets are derived for the Black Sea water column from a coupled physical-biogeochemical model. The model is applied in the deep part of the sea and simulates processes over the whole water column including the anoxic layer that extends from ?115 m to the bottom (?2000 m). The biogeochemical model involves a refined representation of the Black Sea foodweb from bacteria to gelatinous carnivores. It includes notably a series of biogeochemical processes typical for oxygen deficient conditions with, for instance, bacterial respiration using different types of oxidants (i.e denitrification, sulfate reduction), the lower efficiency of detritus degradation, the ANAMMOX (ANaerobic AMMonium OXidation) process and the occurrence of particular redox reactions. The model has been calibrated and validated against all available data gathered in the Black Sea TU Ocean Base and this exercise is described in Gregoire et al. (2008). In the present paper, we focus on the biogeochemical flows produced by the model and we compare model estimations with the measurements performed during the R.V. KNORR expedition conducted in the Black Sea from April to July 1988 (Murray and the Black Sea Knorr Expedition, 1991). Model estimations of hydrogen sulfide oxidation, metal sulfide precipitation, hydrogen sulfide formation in the sediments and water column, export flux to the anoxic layer and to the sediments, denitrification, primary and bacterial production are in the range of field observations.With a simulated Gross Primary Production (GPP) of 7.9 mol C m⁻² year⁻¹ and a Community Respiration (CR) of 6.3 mol C m⁻² year⁻¹, the system is net autotrophic with a Net Community Production (NCP) of 1.6 mol C m⁻² year⁻¹. This NCP corresponds to 20% of the GPP and is exported to the anoxic layer. In order to model Particulate Organic Matter (POM) fluxes to the bottom and hydrogen sulfide profiles in agreement with in situ observations, we have to consider that the degradation of POM in anoxic conditions is less efficient that in oxygenated waters as it has often been observed (see discussion in Hedges et al., 1999). The vertical POM profile produced by the model can be fitted to the classic power function describing the oceanic carbon rate (C_R=Z^−α) using an attenuation coefficient α of 0.36 which is the value proposed for another anoxic environment (i.e. the Mexico Margin) by Devol and Hartnett (2001). Due to the lower efficiency of detritus degradation in anoxic conditions and to the aggregation of particles that enhanced the sinking, an important part of the export to the anoxic layer (i.e. 33%, 0.52 mol C m⁻² year⁻¹) escapes remineralization in the water column and reaches the sediments. Therefore, sediments are active sites of sulfide production contributing to 26% of the total sulfide production.In the upper layer, the oxygen dynamics is mainly governed by photosynthesis and respiration processes as well as by air-sea exchanges. ?71% of the oxygen produced by phytoplankton (photosynthesis+nitrate reduction) is lost through respiration, ?21% by outgasing to the atmosphere, ?5% through nitrification and only ?2% in the oxidation of reduced components (e.g. Mn²⁺, Fe²⁺, H₂S).The model estimates the amount of nitrogen lost through denitrification at 307 mmol N m⁻² year⁻¹ that can be partitioned into a loss of ?55% through the use of nitrate for the oxidation of detritus in low oxygen conditions, ?40% in the ANAMMOX process and the remaining ?5% in the oxidation of reduced substances by nitrate.In agreement with data analysis performed on long time series collected since the 1960s (Konovalov and Murray, 2001), the sulfide and nitrogen budgets established for the anoxic layer are not balanced in response to the enhanced particle fluxes induced by eutrophication: the NH₄ and H₂S concentrations increase.     

Settlement of <Emphasis Type="Italic">Macoma balthica</Emphasis> larvae in response to benthic diatom films

The role of multi-species benthic diatom films (BDF) in the settlement of late pediveliger larvae of the bivalve Macoma balthica was investigated in still-water bioassays and multiple choice flume experiments. Axenic diatom cultures that were isolated from a tidal mudflat inhabited by M. balthica were selected to develop BDF sediment treatments characterized by a different community structure, biomass, and amount of extracellular polymeric substances (EPS). Control sediments had no added diatoms. Although all larvae settled and initiated burrowing within the first minute after their addition in still water, regardless of treatment, only 48–52% had completely penetrated the high diatom biomass treatments after 5 min, while on average 80 and 69% of the larvae had settled and burrowed into the control sediments and BDF with a low diatom biomass (<3.5 μg Chl a g⁻¹ dry sediment), respectively. The percentage of larvae settling and burrowing into the sediment was negatively correlated with the concentration of Chl a and EPS of the BDF. This suggests higher physical resistance to bivalve penetration by the BDF with higher diatom biomass and more associated sugar and protein compounds. The larval settlement rate in annular flume experiments at flow velocities of 5 and 15 cm s⁻¹ was distinctly lower compared to the still-water assays. Only 4.6–5.8% of the larvae were recovered from BDF and control sediments after 3 h. Nonetheless, a clear settlement preference was observed for BDF in the flume experiments; i.e., larvae settled significantly more in BDF compared to control sediments irrespective of flow speed. Comparison with the settlement of polystyrene mimics and freeze-killed larvae led to the conclusion that active selection, active secondary dispersal and, at low flow velocities (5 cm s⁻¹), passive adhesion to the sediment are important mechanisms determining the settlement of M. balthica larvae in estuarine biofilms.