Parameterization of iron and manganese cycling in the Black Sea suboxic and anoxic environment

New and published data on the distribution and speciation of manganese and iron in seawater are analyzed to identify and parameterize major biogeochemical processes of their cycling within the suboxic (15.6σ_t16.2) and anoxic layers (σ_t16.2) of the Black Sea. A steady-state transport-reaction model is applied to reveal layering and parameterize kinetics of redox and dissolution/precipitation processes. Previously published data on speciation of these elements in seawater are used to specify the nature of the transformations. Two particulate species of iron (Fe(III) hydroxide and Fe(II) sulfide) are necessary to adequately parameterize the vertical profile of suspended iron, while three particulate species (hydrous Mn(IV) oxide, Mn(II) sulfide, and Mn(II) carbonate) are necessary to describe the profile of suspended manganese. In addition to such processes as mixing and advection, precipitation, sinking, and dissolution of manganese carbonate are found to be essential in maintaining the observed vertical distribution of dissolved Mn(II). These results are used to interpret the observed difference in the form of vertical distribution for dissolved Mn(II) and Fe(II). Redox transformations of iron and manganese are coupled via oxidation of dissolved iron by sinking suspended manganese at σ_t16.2±0.2 kg m⁻³. The particulate manganese, necessary for this reaction, is supplied through oxidation of dissolved Mn(II). The best agreement with observations is achieved when nitrate, rather than oxygen, is set to oxidize dissolved Mn(II) in the lower part of the suboxic layer (15.90σ_t16.2). The results support the idea that, after sulfides of these metals are formed, they sink with particulate organic matter. The sinking rates of the particles and specific rates of individual redox and dissolved-particulate transformations have been estimated by fitting the vertical profile of the net rate.     

Benthic foraminifera living in Gulf of Mexico bathyal and abyssal sediments: Community analysis and comparison to metazoan meiofaunal biomass and density

Benthic foraminiferal biomass, density, and species composition were determined at 10 sites in the Gulf of Mexico. During June 2001 and 2002, sediment samples were collected with a GoMex box corer. A 7.5-cm diameter subcore was taken from a box core collected at each site and sliced into 1-cm or 2-cm sections to a depth of 2 or 3 cm; the >63-μm fraction was examined shipboard for benthic foraminifera. Individual foraminifers were extracted for adenosine triphosphate (ATP) using a luciferin–luciferase assay, which indicated the total ATP content per specimen; that data was converted to organic carbon. Foraminiferal biomass and density varied substantially (2–53 mg C m⁻²; 3600–44,500 individuals m⁻², respectively) and inconsistently with water depth: although two 1000-m deep sites were geographically separated by only 75 km, the foraminiferal biomass at one site was relatively low (9 mg C m⁻²) while the other site had the highest foraminiferal biomass (53 mg C m⁻²). Although most samples from Sigsbee Plain (>3000 m) had low biomass, one Sigsbee site had >20 mg foraminiferal C m⁻². The foraminiferal community from all sites (i.e. bathyal and abyssal locales) was dominated by agglutinated, rather than calcareous or tectinous, species. Foraminiferal density never exceeded that of metazoan meiofauna at any site. Foraminiferal biomass, however, exceeded metazoan meiofaunal biomass at 5 of the 10 sites, indicating that foraminifera constitute a major component of the Gulf's deep-water meiofaunal biomass.     

Diversity patterns in macrobenthos across a continental slope in the NE Atlantic

Different estimates were used to assess the diversity of the total macrofauna and its major taxonomic groups separately from a broad bathymetric range at a site in the NE Atlantic. In the Goban Spur region, a transect was sampled from the shelf to the abyssal plain over a depth range from 200 to 4500 m and in the Porcupine Sea Bight two stations were sampled (at 3670 m and 4115 m). Species diversity (the number of species per number of individuals) increased with increasing water depth, both when expressed as Hurlbert's E(S_n) and as Shannon's H′log e. The expected number of species in a 100-individual sample E(S₁₀₀) of total macrofauna increased from 30 on the shelf to 68 on the abyssal plain. Evenness (the proportional abundance of species), estimated with Shannon's J′, also increased with water depth from 0.66 to 0.91, whereas dominance (Simpson's D) decreased from 0.09 to 0.01. Species richness (the number of species per unit of area), however, showed a parabolic pattern with a peak at the upper slope. The largest number of species was found at the slope station at 1425 m (232 species within 0.66 m²). It is argued that species richness is not a synonym of species diversity, but that species richness depends both on species density (which decreases with increasing water depth) and on species diversity. Across the whole bathymetric range (200 to 4500 m) a total of 696 species within 8327 specimens in a total sampled area of 4.12 m² were counted, yielding mean values of 12 individuals per species and 169 species per m². Different communities were found to exist on the shelf, slope and abyss. It is suggested that this could have been caused by different selection processes. Differences in life-history strategies and organic-matter supply could (at least partly) explain the different community structures and diversity patterns found along the depth gradient.     

Variability of physical and chemical characteristics along the 70-m isobath of the southeastern Bering Sea

From observations of ice cover, temperature, salinity, currents and nitrate, it is evident that along-shelf variability was significant over the middle shelf of the eastern Bering Sea, but less distinct than that observed in the cross-shelf domains. Along the 70-m isobath, three zones were evident in the summer: the southeastern cold pool (centered at 57°N); an intermediate zone, consisting of warmer water, with weaker stratification; and the northern cold pool, extending northward from 58°N. Small-scale (20 km) horizontal features that persisted for months were common. Nutrient concentrations were related to salinity and were replenished more uniformly over the southern shelf, than north of the Pribilof Islands. Although mean currents were weak (1 cm s⁻¹), short energetic advective events impacted the temperature and salinity structure.     

Physical dynamics and biological implications of Cyclone Noah in the lee of Hawaii during E-Flux I

Quasi-synoptic observations of the horizontal and vertical structure of a cold-core cyclonic mesoscale eddy feature (Cyclone Noah) were conducted in the lee of Hawai’i from November 4–22, 2004 as part of the E-Flux interdisciplinary collaborative research program. Cyclone Noah appears to have spun up to the southwest of the ‘Alenuihaha Channel (between Maui and Hawai’i) as a result of strong and persistent northeasterly trade winds through the channel. Shipboard hydrographic surveys 2.5 months later suggest that Noah weakened and was in a hypothesized spin-down phase of its life cycle. Although the initial surface expression of Noah was limited in scale to 40 km in diameter and, as evidenced by surface temperatures, 2–3 °C cooler than the surrounding waters, depth profiles revealed a fully developed semi-elliptical shallow feature (200 m), 144 km long and 90 km wide (based on sigma-t=23 kg m⁻³) with tangential speeds of 40–80 cm s⁻¹, and substantial isopycnal doming. Potential vorticity distribution of Noah suggests that radial horizontal flow of the core water was inhibited from the surface to depths of 75 m, with high vorticity confined above the sigma-t=23.5 kg m⁻³ isopycnal surface. Upward displacements of isopycnal surfaces in the eddy's center (50 m) were congruent with enhanced pigment concentrations (0.50 mg m⁻³). Comparisons of the results obtained for E-Flux I (Noah) and E-Flux III (Opal) suggest that translation characteristics of cyclonic Hawaiian lee eddies may be important in establishing the biogeochemical and biological responses of the oligotrophic ocean to cyclonic eddies.     

Probability distributions of wave heights in bimodal seas in an offshore basin

This study aims at assessing the adequacy for describing bimodal sea states of different non-linear probability distributions that have been developed for single sea states. It is based on data collected at an offshore test basin. The measurements represent three bimodal sea states with individual unidirectional wave systems propagating at 60, 90 and 120 from each other. The wave spectra are separated into swell and wind sea components and the relative energy ratio between the areas under the associated spectral curves is estimated and is related with the statistics of the time series considered. Dependence is found between the normalized high order cumulants, which describe the non-Gaussian surface, and the predominant contribution of the wind sea energy. Furthermore, the probabilities of exceedance of the individual wave heights are estimated and compared with the Rayleigh model and with other models that take into account either the effect of spectral bandwidth or the effect of wave nonlinearities. The results are discussed with respect to three classes of sea states that reflect the relative contribution of swell and wind sea energy.     

Dissociation constants of protonated methionine species in seawater media

The dissociation constants (pK₁ and pK₂) for methionine have been measured in artificial seawater as a function of salinity (S = 5 to 35) and temperature (5 to 45 °C). The seawater pK₂ values were lower than the values in NaCl at the same ionic strength while the pK₁ values in seawater were lower only at S = 35. In an attempt to understand these differences, we have made measurements of the constants in Na–Mg–Cl solutions at 25 °C. The measured values have been used to determine the formation of Mg²⁺ complexes and Pitzer interaction parameters for Mg²⁺ with methionine. The seawater model with the interaction parameters accounts for the differences between the value of pK₁ and pK₂ between NaCl and seawater. This study demonstrates that it is important to consider all of the ionic interactions in natural waters when examining the dissociation of organic acids.     

Dissolved organic carbon, nitrogen and phosphorus in the N-E Atlantic and the N-W Mediterranean with particular reference to non-refractory fractions and degradation

Dissolved organic carbon (DOC), dissolved organic nitrogen (DON) and dissolved organic phosphorus (DOP) measured in deep profiles in the N-E Atlantic and in the N-W Mediterranean in the period 1984–2002 are described. After accurate validation, they show close agreement with those previously published.Classic profiles were obtained, with concentrations decreasing in deep waters. In the Mediterranean and in the Atlantic comparable concentrations were found in the 1500–2000 m waters, 44–46 μmol l⁻¹ DOC, 2.6–2.8 μmol l⁻¹ DON and 0.02–0.03 μmol l⁻¹ DOP. In the surface layers, DOC concentrations were higher, but DON and DOP concentrations lower, in the Mediterranean than in the Atlantic, leading to higher element ratios in the Mediterranean. In autumn, values were, respectively, DOC:DON 17 vs. 14, DOC:DOP 950 vs. 500 and DON:DOP 55 vs. 35. The data suggest an increase in DOC and DON in the North Atlantic Central Water over 15 years, which may be linked to the North Atlantic climatic oscillations.Refractory DOM found in the 1500–2000 m layer exhibited C:N:P ratios of 1570:100:1. The labile+semi-labile (=non-refractory) DOM (nrDOM) pool was computed as DOM in excess of the refractory pool. Its contribution to total DOM above the thermocline in the open sea amounted to 25–35% of DOC, 30–35% of DON, and 60–80% of DOP. Element ratios of the nrDOM varied among stations and were lower than those of refractory DOM, except for C:N in the Mediterranean: nrDOC:nrDON 10–19, nrDOC:nrDOP 160–530 and nrDON:nrDOP 15–38. The specific stoichiometry of DOM in the Mediterranean led us to postulate that overconsumption of carbon is probably a main process in that oligotrophic sea.By coupling non-refractory DOM stoichiometry and relationships between the main DOM elements in the water column, the relative mineralization of C, N and P from DOM was studied. Below the thermocline, the preferential removal of phosphorus with regard to carbon from the semi-labile DOM can be confirmed, but not the preferential removal of nitrogen. In the ocean surface layers, processes depend on the oceanic area and can differ from deep waters, so preferential carbon removal seems more frequent. Bacterial growth efficiency data indicate that bacteria are directly responsible for mineralization of a high proportion of DON and DOP in the deep water.     

Determination of picomolar iron in seawater by double Mg(OH)2 precipitation isotope dilution high-resolution ICPMS

A low-blank pre-concentration procedure is described for the analysis of picomolar iron (Fe) in seawater by isotope dilution high-resolution inductively coupled plasma mass-spectrometry (HR-ICPMS). The procedure uses a two-step Mg(OH)₂ co-precipitation procedure to extract Fe from a 50 ml seawater sample into a 100 μl 4% nitric acid (HNO₃) solution followed by HR-ICPMS measurement. The high pre-concentration ratio ( 500:1) achieved by the procedure minimizes the Fe blank due to ICPMS instrumental Fe background and results in a detection limit of  2 pM and a precision of  4% at the 50 pM Fe level. The measurement of a low-Fe seawater sample spiked with gravimetric Fe standard shows that the method can clearly distinguish 0.01 nM Fe from 0.02 nM Fe in seawater with high accuracy. The method is demonstrated by the analysis of dissolved Fe in the equatorial Pacific Ocean.     

Dissolved manganese and silicic acid in the Columbia River plume: A major source to the California current and coastal waters off Washington and Oregon

The spatial distributions of dissolved manganese and nutrients were examined in the Columbia River plume off Oregon and Washington during the summer of 2004 and 2005 as part of the River Influence on Shelf Ecosystems (RISE) program. Factors influencing the hydrochemical characteristics of the freshly formed and aged Columbia River plume were investigated. Hydrographic data and nutrient concentrations were used to delineate three distinct water sources for the Columbia River Plume: California Current surface water, coastal upwelled water, and Columbia River water. The warm, intermediate salinity, nutrient poor California Current water contains low levels of dissolved manganese (< 5 nM) and silicic acid (< 5 μM), and is depleted in nitrate. The cold, high salinity, nutrient rich, freshly upwelled water is highly variable (2–20 nM) in dissolved manganese and can be as high as  45 μM in silicic acid and  30 μM nitrate. The variable Columbia River has summer temperatures ranging from  13 to 24 °C, high silicic acid concentrations (ranging from  120 to 200 μM), and lower nitrate concentrations (ranging from  2 to 20 μM). During the summer, the concentrations of silicic acid and dissolved manganese can exceed 100 μM and 200 nM, respectively, in near-field Columbia River plumes. These values are markedly greater than those of surface coastal waters (even during upwelling conditions). As the plume advects and mixes, the concentrations of these two constituents remain relatively high within plume waters. The concentrations of dissolved manganese in the near-field plume vary with tidal amplitude, exhibiting much higher concentrations for a given salinity during spring tides than during neap tides. For example, the Columbia River plume at a salinity of 20 has a concentration of dissolved manganese of  240 nM during spring tides, as compared to only  60 nM during low amplitude tides. Silicic acid concentrations in the near-field plume remain relatively constant throughout the tidal month. Calculations indicate there is roughly an equivalent yearly delivery of dissolved manganese and silicic acid to the coastal waters off Oregon and Washington by upwelled waters and by the Columbia River plume.     

Relevance of peat draining rivers in central Sumatra for the riverine input of dissolved organic carbon into the ocean

Sources and discharges of dissolved organic carbon (DOC) from the central Sumatran river Siak were studied. DOC concentrations in the Siak ranged between 560 and 2594 μmol l⁻¹ and peak out after its confluence with the river Mandau. The Mandau drains part of the central Sumatran peatlands and can be characterized as a typical blackwater river due to its high DOC concentration, its dark brown-coloured, acidic water (pH 4.4–4.7) and its low concentration of total suspended matter (12–41 mg l⁻¹). The Mandau supplies about half of the DOC that enters the Siak Estuary where it mixes conservatively with ocean water. The DOC input from the Siak into the ocean was estimated to be 0.3 Tg C yr⁻¹. Extrapolated to entire Indonesia the data suggest a total Indonesian DOC export of 21 Tg yr⁻¹ representing 10% of the global riverine DOC input into the ocean.     

Dissociation constants of protonated cysteine species in seawater media

The dissociation constants (pK₁, pK₂ and pK₃) for cysteine have been measured in seawater as a function of temperature (5 to 45 °C) and salinity (S = 5 to 35). The seawater values were lower than the values in NaCl at the same ionic strength. In an attempt to understand these differences, we have made measurements of the constants in Na–Mg–Cl solutions at 25 °C. The measured values have been compared to those calculated from the Pitzer ionic interaction model. The lower values of pK₃ in the Na–Mg–Cl solutions have been attributed to the formation of Mg²⁺ complexes with Cys²⁻ anions

Mg²⁺ + Cys²⁻ = MgCys

The stability constants have been fitted to

after corrections are made for the interaction of Mg²⁺ with H⁺.The pK₁ seawater measurements indicate that H₃Cys⁺ interacts with SO₄²⁻. The Pitzer parameters β⁰(H₃CysSO₄), β¹(H₃CysSO₄) and C(H₃CysSO₄) have been determined for this interaction. The formation of CaCys as well as MgCys are needed to account for the values of pK₂ and pK₃ in seawater.The consideration of the formation of MgCys and CaCys in seawater yields model calculated values of pK₁, pK₂ and pK₃ that agree with the measured values to within the experimental error of the measurements. This study shows that it is important to consider all of the ionic interactions in natural waters when examining the dissociation of organic acids.     

The shipboard analysis of trace levels of sulfur hexafluoride, chlorofluorocarbon-11 and chlorofluorocarbon-12 in seawater

Methods are described for the rapid (11 min) automated shipboard analysis of dissolved sulfur hexafluoride (SF₆) in small volume (200 cm³) seawater samples. Estimated precision for the SF₆ measurements is 2% or 0.02 fmol kg⁻¹ (whichever is greater). The method also allows for the simultaneous measurement of chlorofluorocarbon-11 (CFC11) and chlorofluorocarbon-12 (CFC12) on the same water sample, with significantly improved sensitivity over previous analytical methods.     

A dynamic equation for the potential energy anomaly for analysing mixing and stratification in estuaries and coastal seas

In this paper, a time-dependent dynamic equation for the potential energy anomaly, , is rigorously derived from dynamic equations for potential temperature and salinity, the continuity equation and the equation of state for sea water. The terms locally changing are (A) the -advection, (B) the depth-mean straining, (C) the non-mean straining, (D) the vertical advection, (E) the vertical mixing, (F) surface and bottom density fluxes, (G) inner sources of density e.g. due to absorption of solar radiation and the non-linearity of the equation of state, and (H) horizontal divergence of horizontal turbulent density fluxes. In order to derive the equation in concise form, a vertical velocity (linearly varying with depth) with respect to depth-proportional vertical coordinates had to be defined. The evaluation of the terms in the -equation is then carried out for a one-dimensional tidal straining study and a two-dimensional estuarine circulation study. Comparisons to empirical estimates for these terms are made for the one-dimensional study. It is concluded that the -equation provides a general reference for empirical bulk parameterisations of stratification and mixing processes in estuaries and coastal seas and that it is a tool for complete analysis of the relevant terms from numerical models.     

Experimental study on the gap entrance profile affecting rudder gap cavitation

The unsteady cavity patterns around the gap of the conventional and newly developed semi-spade rudders for marine ships are visualized qualitatively using a high-speed CCD camera. Time-resolved PIV analysis is also performed with the bubble tracers to study the flow behavior over the rudder surface. In addition, pressure measurements are conducted on the rudder surface and inside the gap to find out the flow characteristics around the gap entrance of the rudder. Both the rudders are tested without a propeller wake at the various cavitation numbers and at the rudder deflection angle of −8°θ10°. The strong cavitation patterns around the conventional rudder gap are significantly reduced by adopting a newly developed entrance profile, and a time-resolved velocity field is found to be very effective in catching the vortical cavity flow around the rudder gap. The stagnation point near the gap entrance of the conventional rudder can cause unsteady cavity flow. However, the developed rudder has very stable pressure distribution along the horn surface and decreases the pressure inside the gap because of the modification of the gap entrance. The pressure distribution around the gap of the suction side is closely related to the variation of the rudder deflection angle. The cavitation inception speed is delayed by about 4 knots in the angle range of −5°θ5° by employing the developed profile of the gap entrance.     

Geometry of scour hole around, and the influence of the angle of attack on the burial of finite cylinders under combined flows

Experiments were conducted to investigate the geometry of the scour hole and flow structure around short cylinders under the action of waves alone (WA) and combined flows (CF). The study is aimed at better understanding the dynamics of isolated objects on a sandy floor under oscillatory flows as occurs in shallow water regions in coastal areas. Flow velocities within the fluid core were recorded and 3D mapping of the bottom was performed with sub-aquatic acoustic sensors. Experiments were conducted for cylinder Reynolds wave number and Keulegan-Carpenter number within the ranges 10⁴R_e1.7×10⁵ and 2KC71, respectively. The present experimental evidence shows that the geometric characteristics of the scour hole (length and width) depend primarily on the Keulegan-Carpenter number (KC) and the cylinder aspect ratio (a_r=L_c/D). The effect of variation in the angle of attack of the flow with respect to the cylinder main axis was also investigated. Initial orientations of zero and ninety degrees were found to be stable while cylinders with intermediate initial orientations tended to orientate their main axes perpendicular to the flow direction. The final angle of orientation was found to be primarily a function of the Shields parameter, θ, and the initial angle of attack, α_i.     

An inflection in the rate of early mid-Holocene eustatic sea-level rise: A new sea-level curve from Singapore

This study presents a sea-level curve from 9500 to 6500 cal BP for the farfield location of Singapore, on the Sunda Shelf in southeast Asia. The curve is based on more than 50 radiocarbon dates from elevations of +1.43 m to −15.09 m representing sea-level index points in intertidal mangrove and shallow marine sediments deposited by sea-level rise accompanying deglaciation. The results indicate that mean sea level rose rapidly from around −17 m at 9500 cal BP to around −3 m by 8000 cal BP. After this time, the data suggest (but do not unequivocally prove) that the rate of sea-rise slowed for a period of 300–500 years centred on 7700 cal BP, shortly after the cessation of meltwater input to the oceans from the northern hemisphere. Renewed sea-level rise amounting to 3–5 m began around 7400 cal BP and was complete by 7000 cal BP. The existence of an inflection in the rate of sea-level rise, with a slow-down centred on 7700 cal BP, is broadly consistent with other available sea-level curves over this interval and is supported by evidence of stable shorelines and delta initiation elsewhere at this time, as well as evidence of comparatively rapid retreat of the West Antarctic ice sheet beginning around 7500 cal BP. ‘Stepped’ sea-level rise occurring shortly after 7500 cal BP and also earlier during deglaciation may have served to focus significant post-glacial episodes of human maritime/coastal dispersal, into comparatively narrow time intervals.     

Seasonal cycles and long-term trends of plankton in shelf and oceanic habitats of the Norwegian Sea in relation to environmental variables

The plankton abundance data of the Continuous Plankton Recorder (CPR) route from Bergen or Rotterdam to Weather Station Mike (6444^′N, 2E) from 1949 to 1981 were analysed for long-term trends and seasonal production cycles, and were related to environmental data. The data were explored using the canonical correlation analysis and nonparametric techniques like the Nadaraya–Watson regression. While large copepods such as Calanus spp. and Metridia lucens did not show any temporal trends, a sharp decrease in the abundances of smaller copepods and phytoplankton was observed after 1960. The temporal trends were not related to the NAO, but did show a correlation with the wind direction. Seasonal abundance curves showed that production of both phytoplankton and zooplankton taxa started earlier in coastal water compared to Atlantic water. From the 1950s to the 1970s most taxa showed a delay in the start of the seasonal production cycles, indicating a reduction in the length of the productive cycle. This may to some extent explain the reduced abundance of smaller copepods, phytoplankton and other species during the 1960s and 1970s.     

Uncertainties in the preparation of Ra Mn fiber standards

Delayed coincidence counters (RaDeCC), used for measuring ²²³Ra and ²²⁴Ra preconcentrated from water onto MnO₂-impregnated acrylic fiber (“Mn-fiber”), require a standard Mn-fiber column that has a precisely known activity of ²²⁴Ra for calibration. This may be done by adding an aged ²²⁸Th standard solution to adsorb both ²²⁸Th and its daughter ²²⁴Ra quantitatively onto a Mn fiber. We used both seawater and deionized water (DIW) for testing the adsorption efficiency of Th and Ra onto Mn fibers. Our experimental results show that more than 50% of thorium (²³²Th and ²²⁸Th) breaks through the Mn-fiber column when DIW is used as a medium. However, near quantitative recoveries are obtained if filtered (0.45 μm) seawater is used to prepare the standard. In the case of pure DIW, the pH (initial pH  5.3) rises to > 10 after passing through the column while seawater (initial pH  7.8) changes to  7.2. Thus, the lack of thorium adsorption in DIW may be attributed to this huge increase of pH and the consequent formation of Th(OH)₄ and polyhydroxyl colloids. Based on these observations, we recommend that one should use either artificial seawater or natural seawater (which has negligible ²²⁴Ra and ²²⁸Th) as a loading solution after 0.45 μm filtration. In addition, the thorium adsorption efficiency should be confirmed either by thorium analysis of the effluent solution or long-term monitoring of the supported ²²⁴Ra on the Mn fiber using the RaDeCC. Similar cautions are likely necessary for making ²²³Ra standards by adsorption of ²²⁷Ac onto Mn fibers.     

Annual dynamics of pelagic primary production and respiration in a shallow coastal basin

The Wadden Sea (North Sea, Europe) is a shallow coastal sea with high benthic and pelagic primary production rates. To date, no studies have been carried out in the Wadden Sea that were specifically designed to study the relation between pelagic respiration and production by comparable methods. Because previous studies have suggested that the import of primary-produced pelagic organic matter is important for benthic Wadden Sea carbon budgets, we hypothesised that on an annual average the northern Wadden Sea water column is autotrophic. To test this hypothesis, we studied annual dynamics of primary production and respiration at a pelagic station in a shallow tidal basin (List Tidal Basin, northern Wadden Sea). Since water depth strongly influences production estimates, we calculated primary production rates per unit area in two ways: on the basis of the mean water depth (2.7 m) and on the basis of 1 m depth intervals and their respective spatial extent in the List Tidal Basin. The latter more precise estimate yielded an annual primary production of 146 g C m^− 2 y^− 1. Estimates based on the mean water depth resulted in a 40% higher annual rate of 204 g C m^− 2 y^− 1. The total annual pelagic respiration was 50 g C m^− 2 y^− 1. The P/R ratio varied between seasons: from February to October the water column was autotrophic, with the highest P/R ratio of 4–5 during the diatom spring bloom in April/May. In autumn and winter the water column was heterotrophic. On an annual average, the water column of the List Tidal Basin was autotrophic (P/R 3). We suggest that a large fraction of the pelagic produced organic matter was respired locally in the sediment.