Response of single benthic metrics and multi-metric methods to anthropogenic pressure gradients, in five distinct European coastal and transitional ecosystems

In recent times many benthic indices have been proposed to assess the ecological quality of marine waters worldwide. In this study we compared single metrics and multi-metric methods to assess coastal and transitional benthic status along human pressure gradients in five distinct environments across Europe: Varna bay and lake (Bulgaria), Lesina lagoon (Italy), Mondego estuary (Portugal), Basque coast (Spain) and Oslofjord (Norway). Hence, 13 single metrics (abundance, number of taxa, and several diversity and sensitivity indices) and eight of the most common indices used within the European Water Framework Directive (WFD) for benthic assessment were selected: index of size spectra (ISS), Benthic assessment tool (BAT), Norwegian quality index (NQI), Multivariate AMBI (M-AMBI), Benthic quality index (BQI), (Benthic ecosystem quality index (BEQI), Benthic index based on taxonomic sufficiency (BITS), and infaunal quality index (IQI). Within each system, sampling sites were ordered in an increasing pressure gradient according to a preliminary classification based on professional judgement. The different indices are largely consistent in their response to pressure gradient, except in some particular cases (i.e. BITS, in all cases, or ISS when a low number of individuals is present). Inconsistencies between indicator responses were most pronounced in transitional waters (i.e. IQI, BEQI), highlighting the difficulties of the generic application of indicators to all marine, estuarine and lagoonal environments. However, some of the single (i.e. ecological groups approach, diversity, richness) and multi-metric methods (i.e. BAT, M-AMBI, NQI) were able to detect such gradients both in transitional and coastal environments, being these multi-metric methods more consistent in the detection than single indices. This study highlights the importance of survey design and good reference conditions for some indicators. The agreement observed between different methodologies and their ability to detect quality trends across distinct environments constitutes a promising result for the implementation of the WFD’s monitoring plans. Moreover, these results have management implications, regarding the dangers of misclassification, uncertainty in the assessment, use of conflicting indices, and testing and validation of indices.     

Carbon fluxes through major phytoplankton groups during the spring bloom and post-bloom in the Northwestern Mediterranean Sea

The carbon flux through major phytoplankton groups, defined by their pigment markers, was estimated in two contrasting conditions of the Northwestern Mediterranean open ocean ecosystem: the spring bloom and post-bloom situations (hereafter Bloom and Post-bloom, respectively). During Bloom, surface chlorophyll a (Chl a) concentration was higher and dominated by diatoms (53% of Chl a), while during Post-bloom Synechococcus (42%) and Prymnesiophyceae (29%) became dominant. The seawater dilution technique, coupled to high pressure liquid chromatography (HPLC) analysis of pigments and flow cytometry (FCM), was used to estimate growth and grazing rates of major phytoplankton groups in surface waters. Estimated growth rates were corrected for photoacclimation based on FCM-detected changes in red fluorescence per cell. Given the 30% average decrease in the pigment content per cell between the beginning and the end of the incubations, overlooking photoacclimation would have resulted in a 0.40 d^?1 underestimation of phytoplankton growth rates. Corrected average growth rates (μ_o) were 0.90±0.20 (SD) and 0.40±0.14 d^?1 for Bloom and Post-bloom phytoplankton, respectively. Diatoms, Cryptophyceae and Synechococcus were identified as fast-growing groups and Prymnesiophyceae and Prasinophyceae as slow-growing groups across Bloom and Post-bloom conditions. The higher growth rate during Bloom was due to dominance of phytoplankton groups with higher growth rates than those dominating in Post-bloom. Average grazing rates (m) were 0.58±0.20 d^?1 (SD) and 0.31±0.07 d^?1. The proportion of phytoplankton growth consumed by microzooplankton grazing (m/μ_o) tended to be lower in Bloom (0.69±0.34) than in Post-bloom (0.80±0.08). The intensity of nutrient limitation experienced by phytoplankton indicated by μ_o/μ_n (where μ_n is the nutrient-amended growth rate), was similar during Bloom (0.78) and Post-bloom (0.73). Primary production from surface water (PP) was estimated with ¹⁴C incubations. A combination of PP and Chl a synthesis rate yielded C/Chl a ratios of 34±21 and 168±75 (g:g) for Bloom and Post-bloom, respectively. Transformation of group-specific Chl a fluxes into carbon equivalents confirmed the dominant role of diatoms during Bloom and Synechococcus and Prymnesiophyceae during Post-bloom.     

Coastal and estuarine habitat mapping, using LIDAR height and intensity and multi-spectral imagery

The airborne laser scanning LIDAR (LIght Detection And Ranging) provides high-resolution Digital Terrain Models (DTM) that have been applied recently to the characterization, quantification and monitoring of coastal environments. This study assesses the contribution of LIDAR altimetry and intensity data, topographically-derived features (slope and aspect), and multi-spectral imagery (three visible and a near-infrared band), to map coastal habitats in the Bidasoa estuary and its adjacent coastal area (Basque Country, northern Spain). The performance of high-resolution data sources was individually and jointly tested, with the maximum likelihood algorithm classifier in a rocky shore and a wetland zone; thus, including some of the most extended Cantabrian Sea littoral habitats, within the Bay of Biscay. The results show that reliability of coastal habitat classification was more enhanced with LIDAR-based DTM, compared with the other data sources: slope, aspect, intensity or near-infrared band. The addition of the DTM, to the three visible bands, produced gains of between 10% and 27% in the agreement measures, between the mapped and validation data (i.e. mean producer's and user's accuracy) for the two test sites. Raw LIDAR intensity images are only of limited value here, since they appeared heterogeneous and speckled. However, the enhanced Lee smoothing filter, applied to the LIDAR intensity, improved the overall accuracy measurements of the habitat classification, especially in the wetland zone; here, there were gains up to 7.9% in mean producer's and 11.6% in mean user's accuracy. This suggests that LIDAR can be useful for habitat mapping, when few data sources are available. The synergy between the LIDAR data, with multi-spectral bands, produced high accurate classifications (mean producer's accuracy: 92% for the 16 rocky habitats and 88% for the 11 wetland habitats). Fusion of the data enabled discrimination of intertidal communities, such as Corallina elongata, barnacles (Chthamalus spp.), and stands of Spartina alterniflora and Phragmites australis, which presented misclassification when conventional visible bands were used alone. All of these results were corroborated by the kappa coefficient of agreement. The high classification accuracy found here, selecting data sources, highlights the value of integrating LIDAR data with multi-spectral imagery for habitat mapping in the intertidal complex fringe.     

Capabilities of the bathymetric Hawk Eye LiDAR for coastal habitat mapping: A case study within a Basque estuary

The bathymetric LiDAR system is an airborne laser that detects sea bottom at high vertical and horizontal resolutions in shallow coastal waters. This study assesses the capabilities of the airborne bathymetric LiDAR sensor (Hawk Eye system) for coastal habitat mapping in the Oka estuary (within the Biosphere Reserve of Urdaibai, SE Bay of Biscay, northern Spain), where water conditions are moderately turbid. Three specific objectives were addressed: 1) to assess the data quality of the Hawk Eye LiDAR, both for terrestrial and subtidal zones, in terms of height measurement density, coverage, and vertical accuracy; 2) to compare bathymetric LiDAR with a ship-borne multibeam echosounder (MBES) for different bottom types and depth ranges; and 3) to test the discrimination potential of LiDAR height and reflectance information, together with multi-spectral imagery (three visible and near infrared bands), for the classification of 22 salt marsh and rocky shore habitats, covering supralittoral, intertidal and subtidal zones. The bathymetric LiDAR Hawk Eye data enabled the generation of a digital elevation model (DEM) of the Oka estuary, at 2 m of horizontal spatial resolution in the terrestrial zone (with a vertical accuracy of 0.15 m) and at 4 m within the subtidal, extending a water depth of 21 m. Data gaps occurred in 14.4% of the area surveyed with the LiDAR (13.69 km²). Comparison of the LiDAR system and the MBES showed no significant mean difference in depth. However, the Root Mean Square error of the former was high (0.84 m), especially concentrated upon rocky (0.55–1.77 m) rather than in sediment bottoms (0.38–0.62 m). The potential of LiDAR topographic variables and reflectance alone for discriminating 15 intertidal and submerged habitats was low (with overall classification accuracy between 52.4 and 65.4%). In particular, reflectance retrieved for this case study has been found to be not particularly useful for classification purposes. The combination of the LiDAR-based DEM and derived topographical features with the near infrared and visible bands has permitted the mapping of 22 supralittoral, intertidal and subtidal habitats of the Oka estuary, with high overall classification accuracies of between 84.5% and 92.1%, using the maximum likelihood algorithm. The airborne bathymetric Hawk Eye LiDAR, although somewhat limited by water turbidity and wave breaking, provides unique height information obscured from topographic LiDAR and acoustic systems, together with an improvement of the habitat mapping reliability in the complex and dynamic coastal fringe.     

Growth and grazing rate dynamics of major phytoplankton groups in an oligotrophic coastal site

There has been more attention to phytoplankton dynamics in nutrient-rich waters than in oligotrophic ones thus requiring the need to study the dynamics and responses in oligotrophic waters. Accordingly, phytoplankton community in Blanes Bay was overall dominated by Prymnesiophyceae, remarkably constant throughout the year (31 ± 13% Total chlorophyll a, Tchl a) and Bacillariophyta with a more episodic appearance (20 ± 23% Tchl a). Prasinophyceae and Synechococcus contribution became substantial in winter (Prasinophyceae = 30% Tchl a) and summer (Synechococcus = 35% Tchl a). Phytoplankton growth and grazing mortality rates for major groups were estimated by dilution experiments in combination with high pressure liquid chromatography and flow cytometry carried out monthly over two years. Growth rates of total phytoplankton (range = 0.30–1.91 d⁻¹) were significantly higher in spring and summer (μ > 1.3 d⁻¹) than in autumn and winter (μ ∼ 0.65 d⁻¹) and showed a weak dependence on temperature but a significant positive correlation with day length. Microzooplankton grazing (range = 0.03–1.4 d⁻¹) was closely coupled to phytoplankton growth. Grazing represented the main process for loss of phytoplankton, removing 60 ± 34% (±SD) of daily primary production and 70 ± 48% of Tchl a stock. Chla synthesis was highest during the Bacillarophyceae-dominated spring bloom (Chl a_synt = 2.3 ± 1.6 μg Chl a L⁻¹ d⁻¹) and lowest during the following post-bloom conditions dominated by Prymnesiophyceae (Chl a_synt = 0.23 ± 0.08 μg Chl a L⁻¹ d⁻¹). This variability was smoothed when expressed in carbon equivalents mainly due to the opposite dynamics of C:chl a (range = 11–135) and chl a concentration (range = 0.07–2.0 μg chl a L⁻¹). Bacillariophyta and Synechococcus contribution to C fluxes was higher than to biomass because of their fast-growth rate. The opposite was true for Prymnesiophyceae.     

Pulsating subsidence during a rift episode: stratigraphic and tectonic consequences (Aptian–Albian, northern Spain)

Backstripping analysis has been carried out on nine outcrop sections of the Basque-Cantabrian Aptian and Albian of northern Spain which consist of shallow-water marine deposits up to 7000 m thick. It reveals four main pulses of common accelerated subsidence: I (Early Aptian), II (Early Albian), III (lower Late Albian), and IV (uppermost Late Albian) separated by tectonic quiescent intervals. Tectonic subsidence rates of 70 m Myr^–1 for the Aptian and 101 m Myr^–1 for the Albian were recorded (Sopuerta section). The calculated stretching factor referred to a pre-rift continental crust of 35 km is (β_s=1.32) and represents approximately 8.5 km of crustal thinning. Points of change towards lesser subsidence rate in the curves are correlated with major unconformities in the sections, in the Basque-Cantabrian Basin and in Boreal Europe, suggesting geodynamic changes related with the opening of Bay of Biscay. Minor variations of spreading direction and intensity are proposed to explain the subsidence pulses.     

Statistical method of precipitation prediction

Precipitation records from Salamanca (Spain) are being analyzed to detect temporal oscillations. This paper shows a method to identify these signals in order to determine the behaviour of the precipitation time series in this region.     

Can light intensity cause shifts in natural product and bacterial profiles of the sponge Aplysina aerophoba?

The potential role of microbial consortia on sponge chemistry is well known. However, how environmental factors affect microbial and chemical profiles and how these shifts affect the sponge holobiont are far from being understood. This study experimentally investigated the effect of light on both the concentration of secondary metabolites and the bacterial assemblages of the sponge Aplysina aerophoba. We quantified major brominated alkaloids (BAs) using a high‐performance liquid chromatography system coupled with a UV‐detector and analysed sponge‐associated bacteria using denaturing gradient gel electrophoresis of 16S rRNA gene amplicons. We identified distinct chemical and bacterial profiles between the ectosome and the choanosome of A. aerophoba. The abundance of most secondary metabolites increased regardless of the illumination regime. We found that the probability of occurrence of three microbial phylotypes (operational taxonomic units 84, 86 and 87) was strongly associated with increasing concentrations of three brominated compounds (aerophobin‐1, aplysinamisin‐1 and isofistularin‐3). Although the role of these bacteria remains uninvestigated, these associations between natural products and specific microbial phylotypes outline further hypotheses that will improve our understanding of the organization and functioning of these complex host–symbiont interactions.     

HPLC analysis of algal pigments: a comparison exercise among laboratories and recommendations for improved analytical performance

Pure individual and mixed pigment standards were distributed among Joint Global Ocean Flux Study (JGOFS) pigment analysts to estimate the variability of their spectrophotometric and chromatographic systems. To monitor the integrity of the pigments during the comparison exercise, chlorophyll and carotenoid standards were archived and periodically analyzed by high-performance liquid chromatography (HPLC). Pigment standards stored in the dark under nitrogen at − 20 °C were found to be stable for periods of at least one year. Results from three separate intercalibration exercises document a better agreement for spectrophotometric analyses than for HPLC. For the spectrophotometric comparisons, 90% of the pigments analyzed by participant laboratories were within ± 6% of the mean “consensus” values. By contrast, 65 and 85% of the laboratories agreed to within ± 10 and ± 20%, respectively, when chromatographic analyses were compared. Chlorophyll absorption measurements obtained with a diode array-type spectrophotometer were 6–9% lower than those obtained with monochromator-type spectrophotometers. These underestimates probably result from chlorophyll fluorescence contamination associated with the optical configuration of the diode array spectrophotometer. It was also determined that HPLC methods which are not capable of separating monovinyl chlorophyll a from divinyl chlorophyll a can produce 15–25% overestimates of total chlorophyll a concentration in Prochlorococcus-dominated oceanic waters. A simple dichromatic approach is described for eliminating this variable source of error caused by co-elution of these structurally-related pigments. The use of internal standards and periodic calibration checks with external standards is highly recommended for improving analytical performance.