Oceanic CO2 uptake and biogeochemical variability during the formation of the Eastern North Atlantic Central water under two contrasting NAO scenarios

Long-term variability of the biogeochemical properties during the formation of central waters in the Eastern North Atlantic were analyzed between 42–47°N and 10–20°W from the dataset gathered during the Galicia VII (GVII) and C. Darwin 58/59 (CD58/59) cruises. These cruises that showed important changes in the thermohaline properties and the nutrient abundance of the upper layers were carried out under contrasting conditions of the North Atlantic Oscillation (NAO) index. The different climate forcing led a meridional shift of the transition zone between the formation regions of subpolar and subtropical Eastern North Atlantic Central Water (ENACWp and ENACWt, respectively). This displacement conditioned the presence of each ENACW in the study region and so the thermohaline and biogeochemical properties. The effect of the observed variability at decadal scale on the air–sea CO₂ gradient (ΔfCO₂) and exchange (F_CO2) was analyzed using 1D model approach throughout 11 weekly-steps that simulated the development of a spring bloom during the shoaling of the mixed layer. The outputs of the model showed an intensification of the ocean CO₂ uptake due to higher biological CO₂ drawdown, during positive NAO conditions and its weakening under negative NAO influence.     

Evaluation of two atmospheric models for wind–wave modelling in the NW Mediterranean

The NW Mediterranean experiences, as illustrated by the last decade, strong and rapidly varying storms with severe waves and winds. This has motivated a continuous validation of models and the efforts to improve wave and wind predictions. In this paper we use two atmospherics models, MASS (from SMC-Meteorological Office of Catalunya) and ARPEGE (from Météo-France), to force two third generation wave models: WAM and SWAN. The evaluation and comparison has been carried out for two severe storms registered in November 2001 and March–April 2002.The ARPEGE and MASS models predicted higher 10 m wind speeds than coastal meteorological stations, a fact attributed to local land influences. Regarding the 10 m wind direction, models do not present large differences, although considerable deviations from recorded data were found during some dates. ARPEGE presents less scatter and lower errors than MASS when compared with QuikSCAT data.The 10m wind fields from both atmospheric models were used to force the two selected wave models and analyse the errors and sensitivities when predicting severe wave storms. The wave model simulations show some interesting results; during the storm, the spatial wave pattern using ARPEGE showed a higher maximum, although the values of significant wave height at the buoys were lower than the ones forced by MASS (with both WAM and SWAN). The SWAN simulations show a better agreement in predicting the growing and waning of the storm peaks. The prediction of mean period was improved when using the ARPEGE wind field. However the underestimation by SWAN due to the large energy at high frequencies was evident. Validation of spectral shape predictions showed that it still has considerable error when predicting the full frequency spectra. The storms showed bimodal spectral features which were not always reproduced by wave models and are likely to be responsible for part of the discrepancies.     

Changes in size at maturity of European hake Atlantic populations in relation with stock structure and environmental regimes

European hake (Merluccius merluccius) female size at maturity is estimated on an annual basis for Bay of Biscay and Galician coast, which are parts of the distribution of the Northern and Southern stocks, respectively. Clear trends in this reproductive parameter are observed along the time series and the potential factors affecting these trends have been investigated. Total biomass, different indexes of SSB, age diversity index, fishing mortality at age, NAO winter index, upwelling index and temperature were included in multiple regressions models to assess the relative importance of each of them on shifts in size at maturity.Bay of Biscay and Galician coast hake have followed different evolution in patterns of changes in size at maturity. In Bay of Biscay, a steadily decline of 15 cm has been observed from 1987 to 2004, which is well predicted by fishing mortality and age diversity, but also the environment may have played an important role. However, on the Galician coast a drastic decline of 16 cm from 1980 to 1988 was followed by a rapid increase in size at maturity during the next 10 years to original values and a stable period in the last 6 years. Decreasing biomass may explain the decline in size at maturity in the first period. However, total biomass and spawning biomass declined even during the period when size at maturity increased, which is contrary to compensatory theory. Shifts in environmental regime, NAO and upwelling, may have contributed to a decelerated growth during this period that might explain the later maturation.     

Environmental causes of the fluctuations of red shrimp (Aristeus antennatus) landings in the Catalan Sea

The fluctuations of catches of the deep-sea shrimp Aristeus antennatus (Risso, 1816) were analysed in the 6 ports of Catalonia where production is concentrated, based on monthly landings from 1988 to 2004. The 6 selected ports produced 300 t of red shrimp in 2004, or 80% of the total production in Catalonia. The series for each port showed clear interannual variability, with peaks of production in the early 1990s and more recently from 2001 to 2003. The time periods of the monthly data series, estimated by frequency analysis based on Fourier transform, varied around 7–8 years in the four central ports and 12–13 years in the two northern ports. Additionally, the different curves were not in phase: even in nearby ports, the maximum production is observed in different years. Since the North Atlantic Oscillation (NAO) index is an excellent proxy for long-term series of environmental variables, we aimed to explore relationships between the fluctuation of the NAO index and A. antennatus landings in the Catalan Sea. The correlation between the mean annual NAO index and the annual catches in each port was positive and significant with some time lags (from 1 to 3 years). The existence of clear patterns linking the NAO with marine ecological processes has been demonstrated in many studies, but the underlying ecological mechanisms are far from being well understood. The variations in environmental parameters linked to the NAO may act on biological organisms at different levels (individual, population) through physiology (metabolic and reproductive processes) or through trophic relationships, including ecological cascade effects. We propose that NAO-induced environmental variability may enhance food supply to A. antennatus and hence strengthen the reproductive potential of particular year classes, which result in increased catches 1 to 3 years later, although other possible effects of environmental variability on the population dynamics of this species are worth investigating.     

The contribution of total suspended solids to the Bay of Biscay by Cantabrian Rivers (northern coast of the Iberian Peninsula)

Information of suspended sediments fluxes of small rivers to the coastal zone is sparse, and this is particularly so for the Iberian Rivers. To help address this shortage of information, the relationship between fluvial discharge and total suspended solids (TSS) for the main 28 Cantabrian Rivers using data from 22 years monitoring by the COCA network has been analysed, and their particulate material fluxes to the Bay of Biscay coasts have been quantified. The Cantabrian Fluvial System (drainage basin area of 20,333 km²) may be considered as a quasi-homogeneous fluvial system with an average discharge of 561 m^− 3 s^− 1 and average loads of 35 kg_TSS s^− 1 with rivers showing similar average yields of 56 t km^− 2 a^− 1. The average TSS contribution is 1.2 ± 0.2 10⁹ kg a^− 1. This seaward flux of sediment is dispersed along the entire North Iberian coast and is rather modest (25% of the total supply) in comparison with the output from the French Rivers to the Bay of Biscay. The TSS loads of Cantabrian Rivers indicate they are similar to world upland rivers and those of other parts of Northern Europe according to Milliman and Syvistki [Milliman and Syvistki, 1992. Geomorphic/tectonic control of sediment discharge to the ocean: the importance of small mountainous rivers. Journal of Geology, 100: 525–544] and Milliman [Milliman, 2001. Delivery and fate of fluvial water and sediment to the sea: a marine geologist's view of European rivers. Scientia Marina, 65: 121–132]. Although their TSS flux is practically negligible (13,000 times lower) when compared to the world average flux, they provide a good example of the role of small Atlantic temperate rivers.     

Ocean wave spectrum estimation using measured vessel motions from an in-service container ship

This article is about the use of measured wave-induced vessel motions for estimation of ocean wave spectra by application of the wave buoy analogy. In the study, data from a larger, in-service container ship is considered. The estimation of wave spectra, equivalently sea state parameters, is based on measurements from, respectively, a gyro and two accelerometers leading to the simultaneous use of the pitching motion together with the horizontal and vertical accelerations in a position close to the forward perpendicular. The study of in-service data leads to contemplations about the vessel's advance speed, as the possible existence of sea current means that speed-over-ground (SOG) and speed-through-water (STW) will be different. The article discusses aspects related to advance speed in the context of the wave buoy analogy, and a smaller sensitivity study is conducted. Preceding to the sensitivity study, a comparison is made between sea state estimates by the wave buoy analogy and estimates obtained from a hindcast study. The article shows an acceptable agreement between the two sets of estimates. Following, the main conclusion from the sensitivity study on advance speed is that errors and uncertainties in the speed log have an effect on the estimates of the wave buoy analogy. In fact, the effect can be severe if reliable STW measurements are not available. In the final part, the article includes a few discussions about (non)stationary conditions in the context of the wave buoy analogy, and, although the effect on results is not necessarily detrimental, care must be shown when the wave buoy analogy is applied during in-service conditions.     

The influence of oceanographic scenarios on the population dynamics of demersal resources in the western Mediterranean: Hypothesis for hake and red shrimp off Balearic Islands

The aim of the present paper is to study the relationships between some climatic indices and parental stock, recruitment and accessibility to trawl fishery of hake (Merluccius merluccius) and red shrimp (Aristeus antennatus) off Balearic Islands (western Mediterranean). Available annual catch per unit effort, recruitment and spawning stock biomass have been used as biological data. As environmental data, the meso-scale IDEA index and the large-scale North Atlantic Oscillation (NAO) and Mediterranean Oscillation (MO) indices have been used. To analyze possible links between these indices with the population dynamics of demersal resources, two non-linear approaches have been applied: (i) stock–recruitment relationships from Ricker and Beverton–Holt models, by sequentially incorporating environment factors; (ii) generalized additive modelling, both classical general and threshold non-additive models were considered. The latter simulate an abrupt change in explicative variables across different phases (time periods or climatic index values). The results have shown that two oceanographic scenarios around the Balearic Islands, associated with macro and meso-scale climate regimes, can influence the population dynamics of hake and red shrimp. This is especially true for recruitment, which seems to be enhanced during low NAO and IDEA indices periods. During these periods, colder-than-normal winters generate high amounts of cold Western Mediterranean Intermediate Waters (WIW) in the Gulf of Lions, which flow southwards and reach the Balearic Islands channels in spring, increasing the productivity in the area. This oceanographic scenario could also be favourable to the distribution of hake on the fishing grounds where the trawl fleet targets this species, increasing its accessibility to the fishery. Both spawning stock and abundance of red shrimp seems to be also enhanced by high MO index periods, which could reflect the increased presence of the saline and warm Levantine Intermediate Waters (LIW) in the study area, extending over the fishing grounds of this species. The proposed interactions can be useful to assess and manage these important demersal resources.     

基于马尔可夫链的极值波高预测

本文研究了由日最大波高系列估算设计极值波高时，相邻日最大波高间的相关性对极值预测的影响，从日最大波高系列遵从马尔可夫链的假定出发，考虑到国内外经常采用对数一正态分布的韦布尔分布拟合波高长期分布的现实，本文用解析法求解了对数一正态分布情况下的极值预测。同时，对解析法难以求解的非正态随机变量情况（如韦布尔分布），用计算机随机模拟方法求解其极值预测，用上述两种方法对北大西洋和北海有关日最大波高系列的预测     

Operational Oceanography System applied to the Prestige oil-spillage event

This contribution describes the procedure used during the Prestige oil-spillage event, by means of an Operational Oceanography System, and the behaviour of the present prediction tools (hydrodynamic and dispersion models) applied to it. The accuracy of these tools is estimated by a reanalysis of field data transmitted by a sea surface drifting buoy, released at the time of the oil spill. The numerical models applied were the Regional Ocean Modeling System (ROMS), fed by the available six-hourly NCEP atmospheric information, together with a Lagrangian Particle-Tracking Model (LPTM). ROMS has been used to estimate the current fields for the Bay of Biscay, whilst the LPTM has provided the oil spill trajectories. The results demonstrate that the accuracy of the numerical models depends upon the quality of the meteorological input data. In this case, the current fields at the sea surface, derived by ROMS, have been underestimated by the wind fields of the NCEP reanalysis data. An efficient calibration of these wind fields, with data provided by the Gascony buoy (fixed oceanic and atmospheric station), achieves more realistic looking results; this is reflected in the comparison between the buoy trajectory predicted numerically and the tracked movements of the drifting buoy.     

比斯开湾风浪形成与航法

比斯开湾是世界著名大风浪区之一。强大的西风环流、冷空气和气旋活动频繁以及独特的地形,形成了季节性狂风恶浪,为来往于直布罗陀海峡—多佛尔海峡航线的船舶带来威胁。只有准确预报该海域的天气状况,采取正确航法,才能安全驶过这一险恶海域。     

Filling up gaps in wave data with genetic programming

A given time series of significant wave heights invariably contains smaller or larger gaps or missing values due to a variety of reasons ranging from instrument failures to loss of recorders following human interference. In-filling of missing information is widely reported and well documented for variables like rainfall and river flow, but not for the wave height observations made by rider buoys. This paper attempts to tackle this problem through one of the latest soft computing tools, namely, genetic programming (GP). The missing information in hourly significant wave height observations at one of the data buoy stations maintained by the US National Data Buoy Center is filled up by developing GP models through spatial correlations. The gap lengths of different orders are artificially created and filled up by appropriate GP programs. The results are also compared with those derived using artificial neural networks (ANN). In general, it is found that the in-filling done by GP rivals that by ANN and many times becomes more satisfactory, especially when the gap lengths are smaller. Although the accuracy involved reduces as the amount of gap increases, the missing values for a long duration of a month or so can be filled up with a maximum average error up to 0.21 m in the high seas.     

Assimilation of buoy and satellite data in wave forecasts with integral control variables

The quality of numerical wave forecasts can be improved significantly by assimilating wave observations prior to the forecast. In the present study a technique for such assimilation is developed that exploits (a) the efficiency of a limited number of integral control variables, and (b) the effectiveness of variational (model-consistent) assimilation. The formal procedure is independent of the type of control variables and of the wave model (moreover, no adjoint wave model is required). In the present study, integral control variables are chosen to represent large-scale errors in the driving wind fields and uncertainties in the wave model. The assimilation technique is validated with observations of the ERS-1 satellite altimeter and two waverider buoys in two consecutive storms in the Norwegian Sea. The assimilation of the observations reduced the errors in the forecasted significant wave height at the buoy locations typically from 25% to 12%. For low-frequency waves the effect of the assimilation is similarly significant at one buoy location but marginal at the other buoy location.     

Decadal sea level trends in the Bay of Biscay from tide gauges, GPS and TOPEX

Sea level time series derived from TOPEX altimeter and from tide gauge measurements in the Bay of Biscay (Eastern North Atlantic Ocean) are used to investigate the regional sea level rise. Altimetry sea level anomalies are computed from TOPEX measurements and resampled into fixed along track bins in order to obtain consistent records and to approach as close as possible to the coast. Tidal corrections are critical in the estimation of sea level trends; therefore an additional analysis has been carried out in order to identify and correct for residual tidal signals. The obtained mean sea level rise in the region is + 3.09 ± 0.21 mm/year over the period 1993–2002. Tide gauges have been corrected for vertical land motions by means of collocated GPS measurements. GPS-corrected tide gauges and nearby altimetry trends are found to be statistically consistent. The relationship with the atmospheric forcing is also investigated revealing that 15% of this sea level rise is attributed to the atmospheric pressure effects.     

Monitoring marine plankton ecosystems: Identification of the most relevant indicators of the state of an ecosystem

The number of variables involved in the monitoring of an ecosystem can be high and often one of the first stages in the analysis is to reduce the number of variables. We describe a method developed for geological purposes, using the information theory, that enables selection of the most relevant variables. This technique also allows the examination of the asymmetrical relationships between variables. Applied to a set of physical and biological variables (plankton assemblages in four areas of the North Sea), the method shows that biological variables are more informative than physical variables although the controlling factors are mainly physical (sea surface temperature in winter and spring). Among biological variables, diversity measures and warm-water species assemblages are informative for the state of the North Sea pelagic ecosystems while among physical variables sea surface temperature in late winter and early spring are highly informative. Although often used in bioclimatology, the utilisation of the North Atlantic Oscillation (NAO) index does not seem to provide a lot of information. The method reveals that only the extreme states of this index has an influence on North Sea pelagic ecosystems. The substantial and persistent changes that were detected in the dynamic regime of the North Sea ecosystems and called regime shift are detected by the method and corresponds to the timing of other shifts described in the literature for some European Systems such as the Baltic and the Mediterranean Sea when both physical and biological variables are considered.     

Estimation of extreme significant wave heights and the associated uncertainties: A case study using NORA10 hindcast data for the Barents Sea

The reliability of an offshore structure is dependent on its response to the extreme wave climate; therefore, an adequate knowledge of the wave climate at a location is a prerequisite during design as well as in marine operations. This study aims to contribute to the knowledge of the extreme wave climate in the Norwegian Sector of the Barents Sea, using wave hindcast datasets from the Norwegian Reanalysis 10 km (NORA10) database for four locations.We have considered three commonly used methods for the estimation of extreme wave heights, that is, the initial distribution method, the peak over threshold method, and the annual maxima method. The parametric bootstrap concept is considered in the estimation of the epistemic uncertainty related to sample size. The estimated 100-year significant wave heights obtained from the three methods differ, and the degree of variability in the estimates varies, depending on the dataset. The epistemic uncertainty due to sample size is wider when considering the annual maxima method.The estimates obtained from the three methods show the importance of considering different methods and their associated uncertainties when estimating extreme wave values for design. While it is difficult to single out the best method among the three, the estimated values give knowledge of the possible range of the extreme significant wave heights at the locations. Generally, the datasets considered in this study suggest that the wave climate is less harsh further north compared to the southern region of the Barents Sea. The datasets do not suggest any temporal trends in the historical significant wave heights at any of the locations.     

Decreasing oxygen concentrations in the Gullmar Fjord, Sweden, as confirmed by benthic foraminifera, and the possible association with NAO

The deep Gullmar Fjord on the Swedish west coast may be subject to a trend of decreasing oxygen content of the bottom-water. Due to a non-uniform frequency of oxygen measurements before and after 1980, and small differences in the minimum concentrations, it is difficult to assess with certainty if this is true or not. We have used benthic foraminifera from three ultra-high resolution sediment records to try to find out if such a trend exists. The sediment records comprise the time between 1930 and 1996, with an almost year-based temporal resolution. The foraminiferal records show a significant change in faunal composition during the mid-1970's when it altered from a common Skagerrak fauna to an assemblage characterized by an opportunistic, low oxygen tolerant species; Stainforthia fusiformis. This species is an indicator of low oxygen environments of many Scandinavian fjords. The timing of this faunal change coincides with one of the first recorded severe low oxygen events in the Gullmar Fjord in the winter of 1973/1974. The change suggests that the Gullmar Fjord has actually been subject to a decrease in bottom-water oxygen concentrations. The reason for this change is not clear but here, for the first time, we present sediment records that suggest that there is a causal relation between variations in the North Atlantic Oscillation (NAO) indices and the marine environment. In the early 1970s the NAO-indices altered from a mainly negative phase towards a mainly positive period which is still ongoing. Positive NAO-indices (winter values) in the Skagerrak region are generally marked by westerly winds which normally prevent the exchange of bottom-water in the fjords, i.e., the oxygenation of the sea floor and the supply of oxygen to benthic life.     

Assessment of wind quality for oceanographic modelling in semi-enclosed basins

The quality of surface winds derived from four meteorological models is assessed in the semi-enclosed Adriatic Sea over a 2-month period: a global hydrostatic model ECMWF T511 (40 km resolution), a hydrostatic limited area model LAMBO (20 km), and two non-hydrostatic limited area models: LAMI (7 km) and COAMPS™ (4 km). These wind models are used to drive a 2 km resolution wave model (SWAN) of the Adriatic, and wind and wave results are compared with observations at the ISMAR oceanographic tower off Venice. Waves are also compared at buoy locations near Ancona and Ortona. Consistently with earlier studies, the ECMWF fields underestimate the wind magnitude and do not reproduce the known spatial structure of strong wind events. The results show that the higher-resolution, limited area models LAMI and COAMPS exhibit better amplitude response than the coarser ECMWF: there is a 3- to 4-fold reduction of the wind underestimation at the platform (from 36% to 8–11%). The wave response is also improved with LAMI and COAMPS: there is a 2-fold reduction in the underestimation of wave heights at the platform. These non-hydrostatic models also produce wind fields with more realistic small-scale, spatial structure during strong wind events. The temporal correlation between observed and modelled wind, however, is highest with the global ECMWF model due to the fact that large-scale features can be predicted deterministically, whereas small-scale features can only be predicted stochastically. Models with less small-scale structure have better correlation because they have less “noise.” This explanation is supported by increased correlation between modelled and observed waves, the waves representing a smoothing of the wind over fetch and duration. Although there is room for improvement, the high-resolution, non-hydrostatic models (LAMI and COAMPS) offer significant advantages for driving oceanographic simulations in semi-enclosed basins such as the Adriatic Sea.     

Wave spectrum estimation conditioned on machine learning-based output using the wave buoy analogy

In this work, a hybrid approach for wave spectrum estimation is proposed. Fundamentally, the approach is based on the wave buoy analogy, processing ship response measurements, via a framework combining machine learning and a physics-based method dependent on available transfer functions. Specifically, a non-parametric (Bayesian) estimate is obtained of the directional wave spectrum conditioned on integral wave parameters established by a convolutional neural network. The developed method is assessed in a case study considering about two years of data obtained from an in-service container ship. The method produces good results, significantly improved when compared to the initial estimate made without constraints.     

Prediction of wave parameters by using fuzzy inference system and the parametric models along the south coasts of the Black Sea

Forecasting of sea-state characteristics has a great importance in coastal and ocean engineering studies. Therefore, the purpose of this study was to investigate performances of Adaptive-Network-Based Fuzzy Inference System (ANFIS) and several parametric methods in the Black Sea. For this purpose, different fuzzy models with different input combinations were developed for two different wind data sources (TSMS and ECMWF) at two offshore buoy stations. It also aimed to apply several approaches to event-based data sets for wave predictions. Generally, in literature the tendency is to use time series data for wave predictions. In this kind of prediction approach, lagged time series data are taken as inputs and current or future variables are taken as output. In this study, event-based data for each independent storm were extracted from time series data. Simultaneous or concurrent data of wind speed, blowing duration, fetch length and wave characteristics were detected for each single storm. These event data were then used to set up models. The hindcast results were validated with significant wave height and mean wave period data recorded in Hopa and Sinop buoy stations. The performance of developed fuzzy models were also compared with that of four different parametric methods (Wilson, SPM, Jonswap, and CEM methods) applied for two wind data sources at both buoy stations. Finally, it was determined that in the prediction of both wave parameters (H _s and T _z) the ANFIS models (R = 0.66, squared correlation coefficient, and MAE = 0.37 m, mean absolute error, for the best model in prediction of H _s) were more accurate than the parametric methods (R = 0.63 and MAE = 0.75 m for the best model in prediction of H _s).     

Influence of model domain size,wind directions and Ekman transport on storm surge development inside the Chesapeake Bay: A case study of extratropical cyclone Ernesto, 2006

The Chesapeake Bay is vulnerable to severe flooding caused by hurricanes and strong Northeasters. A 3D storm surge model of the Chesapeake Bay is developed for studying the impact of model domain size, wind directions and Ekman transport on the storm surge in the Chesapeake Bay. The model encompasses the Chesapeake Bay and the US East Coast shelf to reduce the influence of model domain size on surge prediction inside the Chesapeake Bay and to account for both local and remote wind effects. This study used 3D model experiments, with respect to different wind directions, to diagnose the relative influences of the local and remote wind effects and Ekman transport on spatial surge distribution during storm events. The model results confirmed that spatial surge distribution can be well explained by the superposition of two distinct physically driven mechanisms during a storm event: incoming surge wave caused by remote effects and local wind forcings. A large model domain is a necessity for predicting storm surge accurately inside the Chesapeake Bay. The model results suggest that the interactions of the incoming surge propagating into the Bay and the local wind forcing from N and NE directions result in an enhanced setup in the lower to middle portions of the Bay, whereas the superposition of incoming surge and the local wind forcing from S and SE directions enhance the surge in the upper Bay region. A combined northwesterly wind over the middle to upper portions of the Bay and southwesterly wind over the lower Bay can cause a large setdown throughout the entire Bay. The Ekman setup along the coast contributes significantly to the water level variations during storm events. It enhances (reduces) surge inside the Bay under the wind forcings from N and NE (SW, S, and SE) directions.