A gridded sea surface salinity data set for the tropical Pacific with sample applications (1950–2008)

We present a gridded data set of Sea Surface Salinity (SSS) for the tropical Pacific (120°E–70°W; 30°N–30°S), with a grid resolution of 1° longitude, 1° latitude and 1 month, from 1950 to 2008. The product, together with its associated error field, is derived from an objective analysis of about 10 million validated SSS records, with most of the data originating from Voluntary Observing Ships, TAO/TRITON moorings and Argo profilers (during the most recent period). We expect this product to benefit studies in oceanography, meteorology and paleoceanography. As examples of applications, we analyse: (a) the seasonal and ENSO (El Niño Southern Oscillation) modes of observed SSS variability, (b) the ability of 23 coupled models used in the Intergovernmental Panel for Climate Change 4th Assessment Report (IPCC AR4) to simulate the mean SSS and these two time varying modes, and (c) the usefulness of the SSS product and of its associated error field in calibrating and validating the paleo-salinity time series. We anticipate improvements and regular updates to our product, as more SSS data become available from in situ networks and from the ongoing and near-future satellite-derived observations by SMOS (Soil Moisture and Ocean Salinity) and Aquarius.     

Corrections to ocean surface forcing in the California Current System using 4D variational data assimilation

The option for surface forcing correction, recently developed in the 4D-variational (4DVAR) data assimilation systems of the Regional Ocean Model System (ROMS), is presented. Assimilation of remotely-sensed (satellite sea surface height anomaly and sea surface temperature) and in situ (from mechanical and expendable bathythermographs, Argo floats and CTD profiles) oceanic observations has been applied in a realistic, high resolution configuration of the California Current System (CCS) to sequentially correct model initial conditions and surface forcing, using the Incremental Strong constraint version of ROMS-4DVAR (ROMS-IS4DVAR). Results from both twin and real data experiments are presented where it is demonstrated that ROMS-IS4DVAR always reduces the difference between the model and the observations that are assimilated. However, without corrections to the surface forcing, the assimilation of surface data can degrade the temperature structure at depth. When using surface forcing adjustment in ROMS-IS4DVAR the system does not degrade the temperature structure at depth, because differences between the model and surface observations can be reduced through corrections to surface forcing rather than to temperature at depth. However, corrections to surface forcing can generate abnormal spatial and temporal variability in the structure of the wind stress or surface heat flux fields if not properly constrained. This behavior can be partially controlled via the choice of decorrelation length scales that are assumed for the forcing errors. Abnormal forcing corrections may also arise due to the effects of model error which are not accounted for in IS4DVAR. In particular, data assimilation tends to weaken the alongshore wind stress in an attempt to reduce the rate of coastal upwelling, which seems to be too strong due to other sources of error. However, corrections to wind stress and surface heat flux improve systematically the ocean state analyses. Trends in the correction of surface heat fluxes indicate that, given the ocean model used and its potential limitations, the heat flux data from the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) used to impose surface conditions in the model are generally too low except in spring-summer, in the upwelling region, where they are too high. Comparisons with independent data provide confidence in the resulting forecast ocean circulation on timescales ～14 days, with less than 1.5 °C, 0.3 psu, and 9 cm RMS error in temperature, salinity and sea surface height anomaly, respectively, compared to observations.     

Long-term (1998–2007) trends on the spatial distribution of heterotrophic ciliates in the East China Sea in summer: effect of the Three Gorges Dam construction

The spatial distribution of heterotrophic ciliates, environmental factors and potential food items (bacteria, Synechococcus spp. and nanoflagellates) were measured in the East China Sea to examine which variables contributed importantly to the long-term distribution of ciliates between 1998 and 2007. In July 1998 and June 2003, heterotrophic ciliates were found to be abundant (1,000–2,000 × 10³ cells m⁻³) in regions where surface salinity <32 but extremely low (<500 × 10³ cells m⁻³) in shelf waters of surface salinity >32. After August 2003, shortly after the completion of the Three Gorges Dam, we found no significant areal differences in the abundance of heterotrophic ciliates (HC). However, we found a significantly negative correlation between temperature and HC abundance of surface water after the completion of the dam, suggesting that temperature had a greater influence on HC abundance, once the original saline state had changed. For the long-term trends on the vertical distribution of HC, their abundance was significantly higher in the upper 50 m of the water column than at either 75 or 100 m. Abundance of Synechococcus spp. at these levels varied significantly in regions of surface salinity <32, suggesting that ciliates and picophytoplankton contribute greatly to mediating the transfer of organic matter to higher trophic levels in this marine ecosystem.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part II - Performance and application to the California Current System

The Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation systems have been systematically applied to the mesoscale circulation environment of the California Current to demonstrate the performance and practical utility of the various components of ROMS 4D-Var. In particular, we present a comparison of three approaches to 4D-Var, namely: the primal formulation of the incremental strong constraint approach; the dual formulation “physical-space statistical analysis system”; and the dual formulation indirect representer approach. In agreement with theoretical considerations all three approaches converge to the same ocean circulation estimate when using the same observations and prior information. However, the rate of convergence of the dual formulation was found to be inferior to that of the primal formulation. Other aspects of the 4D-Var performance that relate to the use of multiple outer-loops, preconditioning, and the weak constraint are also explored. A systematic evaluation of the impact of the various components of the 4D-Var control vector (i.e. the initial conditions, surface forcing and open boundary conditions) is also presented. It is shown that correcting for uncertainties in the model initial conditions exerts the largest influence on the ability of the model to fit the available observations. Various important diagnostics of 4D-Var are also examined, including estimates of the posterior error, the information content of the observation array, and innovation-based consistency checks on the prior error assumptions. Using these diagnostic tools, we find that more than 90% of the observations assimilated into the model provide redundant information. This is a symptom of the large percentage of satellite data that are used and to some extent the nature of the data processing employed. This is the second in a series of three papers describing the ROMS 4D-Var systems.     

Grazing suppression of dimethylsulphoniopropionate (DMSP) accumulation in iron-fertilised, sub-Antarctic waters

The impact of in situ iron fertilisation on the production of particulate dimethylsulphoniopropionate (DMSPp) and its breakdown product dimethyl sulphide (DMS) was monitored during the SOLAS air-sea gas exchange experiment (SAGE). The experiment was conducted in the high nitrate, low chlorophyll (HNLC) waters of the sub-Antarctic Southern Ocean (46.7°S 172.5°E) to the south-east of New Zealand, during March-April, 2004. In addition to monitoring net changes in the standing stocks of DMSPp and DMS, a series of dilution experiments were used to determine the DMSPp production and consumption rates in relation to increased iron availability. In contrast to previous experiments in the Southern Ocean, DMS concentrations decreased over the course of the 15-d iron-fertilisation experiment, from an integrated volume-specific concentration in the mixed layer on day 0 of 0.78 nM (measured values 0.65-0.91 nM) to 0.46 nM (measured values 0.42-0.47 nM) by day 15, in parallel with the surrounding waters. DMSPp, chlorophyll a and the abundance of photosynthetic picoeukaryotes exhibited indiscernible or only moderate increases in response to the raised iron availability, despite an obvious physiological response by the phytoplankton. High specific growth rates of DMSPp, equivalent to 0.8-1.2 doublings d⁻¹, occurred at the simulated 60% light level of the dilution experiments. Despite the high production rates, DMSPp accumulation was suppressed in part by microzooplankton grazers who consumed between 61% d⁻¹ and 126% d⁻¹ of the DMSPp production. Temporal trends in the rates of production and consumption illustrated a close coupling between the DMSP-producing phytoplankton and their microzooplankton grazers. Similar grazing and production rates were observed for the eukaryotic picophytoplankton that dominated the phytoplankton biomass, partial evidence that picoeukaryotes contributed a substantial proportion of the DMSP synthesis. These rates for DMSPp and picoeukaryotes were considerably higher than for chlorophyll a, indicating higher cycling rates of the DMSP-producing taxa than for the bulk phytoplankton community. When compared to the total phytoplankton community, there was no evidence of selection against the DMSP-containing phytoplankton by the microzooplankton grazers; the opposite appeared to be the case. SAGE demonstrated that increased iron availability in the HNLC waters of the Southern Ocean does not invariably lead to enhanced DMS sea-air flux. The potential suppression of DMSPp accumulation by grazers needs to be taken into account in future attempts to elevate DMS emission through in situ iron fertilisation and in understanding the hypothesised link between levels of Aeolian iron deposition in the Southern Ocean, DMS emission and global albedo.     

Altered feeding habits and strategies of a benthic forage fish (Fundulus heteroclitus) in chronically polluted tidal salt marshes

Responses in feeding ecology of a benthic forage fish, mummichogs (Fundulus heteroclitus), to altered prey resources were investigated in chronically polluted salt marshes (the Arthur Kill-AK, New York, USA). The diet niche breadth of the AK populations of mummichogs was significantly lower than that of the reference population, reflecting reduced benthic macroinfaunal species diversity. Most of the AK populations also had 2-3 times less food in their gut than the reference population. This disparity in gut fullness among the populations appeared to be partly due to ingested prey size shifts; some of the AK populations ingested fewer large prey than the reference population. Furthermore, benthic assemblages were strongly associated with sediment-associated mercury; gut fullness of the AK populations also significantly decreased with increasing mercury body burdens. These results indicate that chronic pollution may have directly (chemical bioaccumulation) and indirectly (reduced prey availability) altered the feeding ecology of mummichogs.     

Numerical simulation of groundwater flowing to horizontal seepage wells under a river

A horizontal seepage well, consisting of an interconnected vertical well, galleries, chambers and small-diameter radiating bores, is used to acquire relatively clean water that has been filtered through natural alluvial deposits in a riverbed. It has wide application, especially in arid and semi-arid areas. The lack of calculation formulae or models for horizontal seepage wells, up until now, has resulted in several false applications. Based on the analysis of groundwater flow characteristics, it has been concluded that several flow regimes coexist and hydraulic head loss exists in the horizontal seepage well. To avoid the difficulty of confirming the flux or head distribution in such a complex system, the model boundary of the whole horizontal seepage well has been moved to that of just the vertical well, and the well-aquifer system was treated as a heterogeneous medium, where the horizontal seepage well itself is a highly permeability medium. A mathematical model has been developed, based on the coupled seepage-pipe flow, by the introduction of equivalent hydraulic conductivity according to different flow regimes. Then a three-dimensional finite difference numerical model, based on the mathematical model, was developed and applied to a horizontal seepage well in China. The numerical model verified the groundwater flow characteristics of the horizontal seepage well. An erratum to this article can be found at     

Autochthonous inheritance of zircon through Cretaceous partial melting of Carboniferous plutons: the Arthur River Complex,Fiordland, New Zealand

TIMS and SHRIMP U–Pb analyses of zircons from Milford Orthogneiss metadiorite (P = 1–1.4 GPa; T ≥ 750°C) of the Arthur River Complex of northern Fiordland reveal a bimodal age pattern. Zircons are predominantly either Paleozoic (357.0 ± 4.2 Ma) and prismatic with oscillatory zoning, or Cretaceous (133.9 ± 1.8 Ma) and ovoid with sector or patchy zoning. The younger age component is not observed overgrowing older grains. Most grains of both ages are overgrown by younger Cretaceous (~120 Ma) metamorphic zircon with very low U and Th/U (0.01). We interpret the bimodal ages as indicating initial igneous emplacement and crystallisation of a dioritic protolith pluton at ~357 Ma, followed by Early Cretaceous granulite-facies metamorphism at ~134 Ma, during which a significant fraction (~60%) of the zircon grains dissolved, and subsequently reprecipitated, effectively in situ, in partial melt pockets. The remaining ~40% of original Paleozoic grains were apparently not in contact with the partial melt, remained intact, and show only slight degrees of Pb loss. Sector zoning of the Cretaceous grains discounts their origin by solid state recrystallisation of Paleozoic grains. The alternative explanation—that the Paleozoic component represents a 40% inherited component in an Early Cretaceous transgressive dioritic magma—is considered less likely given the relatively high solubility of zircon in magma of this composition, the absence of 134 Ma overgrowths, the single discrete age of the older component, equivalent time-integrated ¹⁷⁷Hf/¹⁷⁶Hf compositions of both age groups, and the absence of the Cambrian-Proterozoic detrital zircon that dominates regional Cambro-Ordovician metasedimentary populations. Similar bimodal Carboniferous-Early Cretaceous age distributions are characteristic of the wider Arthur River Complex; 8 of 12 previously dated dioritic samples have a Paleozoic component averaging 51%. Furthermore, the age and chemical suite affinity of these and several more felsic rocks can be matched with those of the relatively unmetamorphosed Carboniferous plutonic terrane along the strike of the Mesozoic margin in southern Fiordland, also supporting the in situ derivation of the Carboniferous “inherited” component.     

Seasonal Dynamics of Microbial Processes in Bottom Sediments of the Sevastopol Coastal Area

Oceanology - A seasonal study of microbial methane oxidation (MO) and sulfate reduction (SR) processes in bottom sediments was carried out at five stations in the estuarian Sevastopol coastal area,...