A generalized model of pelagic biogeochemistry for the global ocean ecosystem. Part II: Numerical simulations

This paper presents a global ocean implementation of a multi-component model of marine pelagic biogeochemistry coupled on-line with an ocean general circulation model forced with climatological surface fields (PELAgic biogeochemistry for Global Ocean Simulations, PELAGOS). The final objective is the inclusion of this model as a component in an Earth System model for climate studies. The pelagic model is based on a functional stoichiometric representation of marine biogeochemical cycles and allows simulating the dynamics of C, N, P, Si, O and Fe taking into account the variation of their elemental ratios in the functional groups. The model also includes a parameterization of variable chlorophyll/carbon ratio in phytoplankton, carrying chl as a prognostic variable. The first part of the paper analyzes the contribution of non-local advective–diffusive terms and local vertical processes to the simulated chl distributions. The comparison of the three experiments shows that the mean chl distribution at higher latitudes is largely determined by mixing processes, while vertical advection controls the distribution in the equatorial upwelling regions. Horizontal advective and diffusive processes are necessary mechanisms for the shape of chl distribution in the sub-tropical Pacific. In the second part, the results have been compared with existing datasets of satellite-derived chlorophyll, surface nutrients, estimates of phytoplankton community composition and primary production data. The agreement is reasonable both in terms of the spatial distribution of annual means and of the seasonal variability in different dynamical oceanographic regions. Results indicate that some of the model biases in chl and surface nutrients distributions can be related to deficiencies in the simulation of physical processes such as advection and mixing. Other discrepancies are attributed to inadequate parameterizations of phytoplankton functional groups. The model has skill in reproducing the overall distribution of large and small phytoplankton but tends to underestimate diatoms in the northern higher latitudes and overestimate nanophytoplankton with respect to picoautotrophs in oligotrophic regions. The performance of the model is discussed in the context of its use in climate studies and an approach for improving the parameterization of functional groups in deterministic models is outlined.     

Surface CO2 measurements in the English Channel and Southern Bight of North Sea using voluntary observing ships

Ships of opportunity have been used to investigate ocean–atmosphere CO₂ fluxes in the English Channel and Southern Bight of the North Sea. Continuous underway measurements of the fugacity of seawater carbon dioxide (fCO₂^sw), chlorophyll, temperature and salinity have been performed along 26 transects during the spring and autumn periods. The spatial fCO₂^sw distribution along the Channel and Southern Bight is modulated by the photosynthetic activity, temperature changes and water mixing between inputs from the North Atlantic Ocean and riverine discharges. The seasonal variability of fCO₂^sw is assessed and discussed in terms of the biology and temperature effects, these having similar impacts. The variation of fCO₂^sw shows similar interannual patterns, with lower values in spring. The annual average of air–sea CO₂ fluxes places the English Channel as neutral area of CO₂ uptake. The spring and autumn data allow differentiating between distal and proximal continental areas. The Southern Bight shows a tendency towards net CO₂ uptake on the distal continental shelf, whereas the Scheldt and Thames Plumes show a CO₂ source behaviour on the proximal continental shelves.     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

An approach to estimating the hull girder response of a ship due to springing

This article describes an estimation method for the hull girder response of a ship due to springing. The linear and nonlinear springing effects on the hull girder are evaluated. Previous studies on the springing response focused mainly on the symmetric response, or vertical response. In this article, however, the springing analysis is extended to asymmetric responses, or horizontal and torsional responses. The Timoshenko beam model was used to calculate the hull girder response and the quadratic strip method was employed to calculate hydrodynamic forces and moments on the hull. To remove irregular frequencies, a rigid lid was adopted on the hull free surface level and hydrodynamic coefficients were interpolated for asymptotic values. Applications to two ships for the symmetric and asymmetric responses were carried out and the effect of springing responses is also discussed.     

Fluxes of dissolved inorganic carbon in three estuarine systems of the Cantabrian Sea (north of Spain)

The diffusive and in situ fluxes of dissolved inorganic carbon (DIC) and total alkalinity (TA) have been measured and an estimation has been made of the water–atmosphere fluxes of CO₂ in three estuarine systems of the Cantabrian Sea during the spring of 1998. Each of these systems undergoes a different anthropogenic influence. The diffusive fluxes of dissolved inorganic carbon and total alkalinity obtained present values ranging between 0.54–2.65 and 0.0–2.4 mmol m⁻² day⁻¹, respectively. These ranges are in agreement with those of other coastal systems. The in situ fluxes are high and extremely variable (35–284 mmol TA m⁻² day⁻¹, 43–554 mmol DIC m⁻² day⁻¹ and 22–261 mmol dissolved oxygen (DO) m⁻² day⁻¹), because the systems studied are very heterogeneous. The values of the ratio of the in situ fluxes of TA and DIC show on average that the rate of dissolution of CaCO₃ is 0.37 times that of organic carbon oxidation. Equally, the interval of variation of the relationship between the benthic fluxes of inorganic carbon and oxygen (F_DIC/F_DO) is very wide (0.3–13.9), which demonstrates the different contributions made by the processes of aerobic and anaerobic degradation of the organic matter, as well as by the dissolution–precipitation of CaCO₃. The water–atmosphere fluxes of CO₂ present a clear dependence on the salinity. The brackish water of these systems (salinity<20), where maximum fluxes of 989 mmol m⁻² day⁻¹ have been estimated, act as a source of CO₂ to the atmosphere. The more saline zones of the estuary (salinity>30) act as a sink of CO₂, with fluxes between −5 and −10 mmol m⁻² day⁻¹.     

Hydrographic cruises off northwest Africa: the Canary Current and the Cape Ghir region

We present hydrographic data for several sections located along the African coastline and off Cape Ghir, carried out at times of weak surface winds (October 1995 and September 1997). The main sections are near the continental slope, at mean distances between 100 and 150 km from the coastline. North of Cape Ghir (31°N) the geostrophic transport (referenced to 650 m) of North Atlantic Central Water through these sections is 3.7 and 2.0 Sv for 1995 and 1997, respectively. This confirms that a major fraction of the water transport by the Canary Current flows east, into the continental slope off northwest Africa, at latitudes above Cape Ghir. Most of this flow continues south past Cape Ghir, along the coast and probably through the eastern passages of the Canary Archipelago. A significant fraction, however, may escape through surface Ekman transport (0.3–0.5 Sv during the early fall season) and by offshore flow at Cape Ghir (1.1 Sv in September 1997, referenced to 650 m). Despite the weak winds the Cape Ghir filament was clearly visible, characterized by localized coastal upwelling associated to a cyclonic shallow structure and cold (and fresh) waters stretching offshore as a very shallow feature (50–100 m deep). The satellite images show that the surface temperature field is highly variable, in rapid response to the surface winds, always with a core region of relatively cold water and commonly with one or two associated eddies. Our results support the existence of two recirculation cells in the area: a horizontal one that connects the interior eastern boundary currents with the coastal region and a vertical one related to both wind-induced and filament upwelling. The data also show a salinity subsurface maximum at the root of the filament, linked to water inflow from northern latitudes, and a subsurface anticyclonic eddy over the Agadir canyon, likely related to the poleward slope undercurrent.     

Coupled transport-reaction pathways and distribution patterns between siliciclastic-carbonate sediments at the Ria de Vigo

This paper examines the linkages between the space-distribution of grain sizes and the relative percentage of the amount of mineral species that result from the mixing process of siliciclastic and carbonate sediments at the Ria de Vigo (NW of Spain).The space-distribution of minerals was initially determined, starting from a detailed mineralogical study based on XRD-Rietveld analysis of the superficial sediments. Correlations between the maps obtained for grain sizes, average fractions of either siliciclastic or carbonates, as well as for individual-minerals, were further stabilised. From this analysis, spatially organized patterns were found between carbonates and several minerals involved in the siliciclastic fraction. In particular, a coupled behaviour is observed between plagioclases and carbonates, in terms of their relative percentage amounts and the grain size distribution.In order to explain these results a conceptual model is proposed, based on the interplay between chemical processes at the seawater–sediment interface and hydrodynamical factors. This model suggests the existence of chemical control mechanisms that, by selective processes of dissolution-crystallization, constrain the mixed environment's long-term evolution, inducing the formation of self-organized sedimentary patterns.     

Patterns in diversity and community structure of epipelagic copepods from the Brazil–Malvinas Confluence area, south-western Atlantic

Diversity and structure of epipelagic copepod communities were investigated using 70 zooplankton samples collected from the top 50 m of the Brazil–Malvinas Confluence area between 35° and 55°S. Biogeographic patterns were investigated using multivariate analysis. Biodiversity patterns were examined using different univariate indexes. Representatives of 35 species of copepods from 23 genera and 13 families were found. Two zones were delimited based on their copepod composition, one in the subtropical waters of the Brazil current, and the other in the subantarctic waters of the Malvinas current. Both environments displayed contrasting patterns of biodiversity depending upon which element is measured (regional species richness as well as mean point species richness were significantly higher in the subtropical group of stations, whereas taxonomic distinctness was significantly higher in the subantarctic community). Based on these contrasting results, we suggest the use of both kind of indexes when defining priority areas for conservation.     

Management of coastal eutrophication: Integration of field data, ecosystem-scale simulations and screening models

A hybrid approach for eutrophication assessment in estuarine and coastal ecosystems is presented. The ASSETS screening model (http://www.eutro.org) classifies eutrophication status into five classes: High (better), Good, Moderate, Poor and Bad (worse). This model was applied to a dataset from a shallow coastal barrier island system in southwest Europe (Ria Formosa), with a resulting score of Good. A detailed dynamic model was developed for this ecosystem, and the outputs were used to drive the screening model. Four scenarios were run on the research model: pristine, standard (simulates present loading), half and double the current nutrient loading. The Ria Formosa has a short water residence time and eutrophication symptoms are not apparent in the water column. However, benthic symptoms are expressed as excessive macroalgal growth and strong dissolved oxygen fluctuations in the tide pools. The standard simulation results showed an ASSETS grade identical to the field data application. The application of the screening model to the other scenario outputs showed the responsiveness of ASSETS to changes in pressure, state and response, scoring a grade of High under pristine conditions, Good for half the standard scenario and Moderate for double the present loadings. The use of this hybrid approach allows managers to test the outcome of measures against a set of well-defined metrics for the evaluation of state. It additionally provides a way of testing and improving the pressure component of ASSETS. Sensitivity analysis revealed that sub-sampling the output of the research model at a monthly scale, typical for the acquisition of field data, may significantly affect the outcome of the screening model, by overlooking extreme events such as occasional night-time anoxia in tide pools.     

Second-order perturbation added mass fluctuation on vertical vibration of tension leg platforms

Total motion in waves can be considered as a superposition of the motion of the body in still water and the forces on the restrained body. In this study the effect of added mass fluctuation on vertical vibration of tension leg platforms (TLP) in the case of vibration in still water for both free and forced vibration subjected to axial load at the top of the leg is presented. This effect is more important when the amplitude of vibration is large. Also this is important in fatigue life study of tethers. The structural model used here is very simple. Perturbation method is used to formulate and solve the problem. First- and second-order perturbations are used to solve the free and forced vibrations, respectively.