The Use of Ocean-Colour Data to Estimate Chla Trends in European Seas

Ocean-colour remote-sensing products have been used to estimate Chl-a trends in European seas with the aim to develop a new indicator based on ocean-colour data for the European Environment Agency (EEA). The new indicator, called CSI023(+), derived from satellite ocean-colour products from the MyOcean Marine Core Service (www.myocean.eu) has been defined and calculated. In our analysis, we have used 3 MyOcean satellite products: 2 global satellite products (SeaWiFS and a merged product) and one regional (adjusted to specific regional Mediterranean conditions) ocean-colour product. We have evaluated the differences among the 3 different products in estimating Chl-a trends. CSI023(+) complements the EEA CSI023 indicator for eutrophication based on chlorophyll-a (Chl-a) in-situ observations. Analysis has revealed the potential of ocean colour as a CSI023(+) indicator to detect large-scale, and in some cases, even local-scale, changes and decreasing trends of Chl-a were observed throughout the Black Sea, the Eastern Mediterranean, the southern part of the Western Mediterranean, the English Channel and the north part of the North Sea. Large areas with increasing trends were observed in the Bay of Biscay, in the North-East Atlantic regions of Ireland and the UK, in the northern part of the North Sea, in the Kattegat and in the Baltic. Specific analysis has been performed in the Mediterranean coastal areas using regional products to investigate local scale results. Validation of ocean-colour products has been carried out through comparison with observations of the Eionet EEA database. The validation results highlight that regional products produced with regional algorithms are recommended for the future.


Introduction
A recent review [1] has presented an overview of eutrophication indicators for assessing environmental status within the European Marine Strategy Framework Directive (MSFD).Several studies presented the use Chl-a, used as a proxy for phytoplankton biomass, as an indicator of eutrophication [2][3][4][5][6][7][8][9][10][11][12][13].Chl-a alone is not enough to estimate eutrophication processes [1] and to provide a complete picture of eutrophic conditions should also be included (i.e.changes in community composition) but Chl-a is one of the relevant biological indicators [1] suitable to describe the status of eutrophication since Chl-a will increase due to increased nutrient availability.
Eutrophication in marine waters has been a management concern in Europe over recent decades.Legislative frameworks have been set up, including the Water Framework Directive (WFD) and the MSFD on transitional and marine waters.Moreover, several other EU Directives are aimed at reducing nutrient loads and impacts.These include the Nitrates Directive (91/676/EEC); the Urban Waste Water Treatment Directive (91/271/EEC) as well as the Integrated Pollution Prevention and Control Directive (96/61/EEC).Measures also arose from a number of other international initiatives and policies (i.e.UN Global Programme of Action for the Protection of the Marine Environment against Land-Based Activities; the Medi-terranean Action Plan (MAP); the Helsinki Convention (HELCOM) for the Baltic Sea; the OSPAR Convention for the North East Atlantic as well as the Bucharest Convention for the Black Sea.
In the European and international framework, EEA has set up an indicator based on in-situ Chl-a trends to monitor eutrophication in the European seas; this is referred to as CSI023 (Chlorophyll in transitional, coastal and marine waters) in the EEA system.For a complete overview of the indicator, please refer to the following web site: http://www.eea.europa.eu/data-and-maps/indicators/chlorophyll-in-transitional-coastal-and/chlorophyll-in-transitio nal-coastal-and-3.
The objective of the CSI023 indicator is to demonstrate the effects of policy measures taken to reduce loading of nitrogen and phosphates from rivers that affect primary production biomass in the coastal zones.CSI023 is calculated from Chl-a in-situ data.The last EEA assessment was performed in March 2013, and its results were presented in terms of concentration of Chl-a in the European seas and CSI023 from 1985 to 2010.
CSI023 is based on measurements collected using fluorometers are given at selected coastal stations, providing an accurate measure of Chl-a.In this regard, two main problems relating to in-situ Chl-a dataset were presented, which are: a low temporal and spatial resolution as well as the poor coverage, especially in Southern European seas.We have therefore identified the challenge of integrating the in-situ dataset with ocean-colour products to complement the former and provide EEA and Member States with valuable information for eutrophication assessment.
Measurements using satellite radiometers of waterleaving radiance in the visible range (ocean colour) can today be used to determine Chl-a concentration, which is an indicator of algal photosynthetic activity and thus related to phytoplankton biomass.Chl-a can now be estimated from ocean-colour data at daily frequencies and 250 m horizontal resolutions.
Ocean-colour satellite products are now available from SeaWiFS, MERIS-Envisat and MODIS-aQUA sensors.The future Sentinel-3 GMES satellite will also have an ocean-colour sensor ensuring continuous monitoring for the period 2015-2030.The first ocean-colour sensor, CZCS, was in operation from approximately 1980 to 1986, afterwards the MOS sensor was from approximately 1997 to 2004 and OCTS from October 1996 to June 1997, finally, SeaWiFS, only started in 1997.
Estimation of Chl-a from ocean colour is an integral value over the e-folding scale of light in water.The optical properties of oceanic waters can be classed into Case 1 or Case 2 waters [17][18][19].In our analysis we used global and regional (adjusted to specific regional Mediterranean conditions) ocean-colour products, both devel-oped for Case 1 waters since only these types of products were available.
A first attempt of evaluating Chl-a trend was performed by using CZCS products [20] and demonstrated the importance of ocean colour data in assessing Chl-a trends and studying the ocean carbon cycle, Recent studies underlined the importance of maintaining continuous, intercalibrated, climate-quality satellite data records for climate-change detection and attribution studies [21].Similar requirements are underlined also from our study that highlights the importance of intercalibrated and consistent merged products.
The analyses of temporal variability of derived apparent optical properties in European seas are conducted in [22] on the SeaWiFS dataset (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007) and the Chl-a trends presented in our study were consistent with the products presented in [22] such as reflectance band ratio trends.Moreover, inter-annual variations in the SeaWiFS global Chl-a concentration (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007) are presented in term of annual trends in [23] for the Global Ocean.A direct comparison with our results cannot be presented because we focused on the summer period trend analysis for consistency with existing EEA CSI023 and specifically on European Seas, while [23] considered the trends of the full year period.
Compared to previous work that only took into account remote sensing [22,23] or in situ [5,[13][14][15][16] the new approach presented in our study focuses on integrating in-situ and ocean colour data to assess the differences in estimating Chl-a concentration and trends and to better respond to the EEA need of providing an estimation of Chl-a trends in all European Seas.
This study aims to developing a methodology to complement the EEA CSI023 indicator with observations based on ocean colour from space and the indicator obtained which was named CSI023(+).
The purpose of the study is to:  Derive and analyse the ocean colour based trends in European seas;  Compare ocean-colour Chl-a data and trends with insitu ones;  Present the trends of Chl-a as CSI023(+) in the European seas for the period 1998-2009.
The study is organized as follows: Section 2 describes the satellite and in-situ data sources and defines the indicator and the Chl-a areas concept and methods; Section 3 is on the validation of datasets; Section 4 is the results; Section 5 concludes the study.
The first dataset "Global Ocean GSM-MyOcean" was developed within the framework of GlobColour (www.globcolour.info)and it is regularly produced and distributed by the Global Monitoring for Environment and Security (GMES) MyOcean project.This is a multi sensor merged product based on SeaWiFS, AQUA and MERIS data available from MyOcean at 4 km resolution.The Global Ocean GSM-MyOcean has been produced using the bio-optical model-based merging procedure observed in [24,25], which combines the normalized waterleaving radiances from different sensor data sets.Over each particular pixel of a geographical grid common to SeaWiFS, AQUA and MERIS, the spectral water-leaving radiance, NLw(λ) spectra from the available sensors at that pixel are selected and combined in a single, multisource, spectrum which is then used in the inversion of the GSM01 semi-analytical ocean colour model [24].The use of these three sensors contributed to reduce data gaps, increasing the coverage over ocean by a factor, which is nearly twice that of any single mission's observations [26].The Global Ocean GSM-MyOcean products have been validated by [26] and by comparing them to the data sets obtained from individual missions.This product has been available since September 1997 and the timeseries is constituted of daily products delivered by MyOcean.The time frame of the dataset chosen for this research work was between 1998-2009 based on consistent processing chains for MODIS and MERIS.Only later in 2010 and 2011 the processing chains were upgraded: the NASA MODIS L2 processing has been switched from R2009.1 to R2010.0 in June 2011 and the MERIS L2 processing switched from the 2nd to the 3rd MERIS reprocessing in autumn of 2012 (see MYO-OC-PUM manual available on line at www.myocean.eu).Moreover the SeaWiFS products are associated with the latest reproc-essing performed by NASA for SeaWiFS using SeaDAS 6.1 software.
The second dataset is the "Global Ocean SeaWiFS RAN-MyOcean" Chl-a dataset.The full SeaWiFStime series, based on SeaWiFS algorithm [27] reprocessed a consistent time series of L2 input data using OC4-V4 algorithm for chlorophyll retrieval.This MyOcean product is associated with the latest reprocessing performed by NASA for SeaWiFS using SeaDAS 6.1 software.Temporal characteristics of the dataset comprises of standard mapped image monthly mean global sea surface Chl-a maps at 9 km resolution (L3 product) and distributed by the GMES MyOcean project.This product has been evaluated in the framework of MyOcean.
The last satellite dataset is the "Med Regional SeaWiFS RAN-MyOcean"; it is the Mediterranean regional product based on SeaWiFS data using the Mediterranean regional ocean colour algorithm for chlorophyll retrieval (MedOC4, [28]).The MedOC4 algorithm has been validated with a large in-situ bio-optical dataset for the Mediterranean area, and its performance has been compared with global algorithms such as OC4v4 for SeaWiFS, and the results show that MedOC4 is the best algorithm for satellite chlorophyll estimates in the Mediterranean [28].The Med regional SeaWiFS RAN has been produced at once, reprocessing the entire SeaWiFS L1 time series with a single software configuration, using the latest version of calibration and ancillary data with the SeaWiFS Data Analysis System (SeaDAS) software package version 6.1 (available seadas.gsfc.nasa.gov).The description of these MyOcean Mediterranean products and the processing system and results on the products validation has been done in [19].
It is important to underline that the three MyOcean products are calibrated for open ocean waters and not specifically for coastal waters, and thus a lower performance in the coastal zone is therefore to be expected.The complexity of optical properties of the European coastal zone requires the use of specialized local algorithm and processing system.Presently, they are not available in the European scale coverage.Thus, we used the My-Ocean global and regional products to investigate trends in both open ocean and coastal areas of the European Seas.The availability of in-situ coastal data will be use at posterior to check potential use of MyOcean ocean colour products to estimate trends also in the coastal waters.The trends analysis using ocean colour products should take into account how these products have been created.For a series based on one sensor, it therefore means making sure that the series have been created with one processing chain, including a consistent calibration table and calibration history [30,31].This applies to the Med Regional SeaWiFS RAN-MyOcean and to Global Ocean SeaWiFS RAN-MyOcean.The possible problems connected with the GSM product are: the use of L2 data provided by the space agencies at the time of the operation update of the dataset that does not ensure the complete consistency of the time series and the method does not account for a specific inter-bias calibration effort in the production of the merged datasets.We are not facing the first problem because our dataset consider the period until 2009 and the upgrade of processing chains occurred later for MERIS and MODIS.Meanwhile, the second problem should be minimized because the methodology of production of the Global Ocean GSM-MyOcean merged product should take into account specific interbias calibration [26].
We have performed a specific analysis of comparing Global Ocean SeaWiFS RAN-MyOcean reprocessed dataset from 1998 to 2009 and Global Ocean GSM-MyOcean to evaluate the possible problems and the differences in the estimation of the Chl-a trends.This analysis is presented in terms of summer Chl-a trends for the European Seas calculated from the Global Ocean GSM-MyOcean and from the Global Ocean SeaWiFS RAN-MyOcean dataset.

In-Situ Data Sources
To validate the Global Ocean GSM-MyOcean and the Med Regional SeaWiFS RAN-MyOcean daily satellite products, we used the in-situ data collected in the framework of EEA-Eionet databank and distributed by ICES.The data consists of Chl-a profiles collected in European Seas.The in-situ dataset covers the period 1998-2009 and contains, after the removal of duplicates, 103170 profiles.For comparison with satellite dataset, we have selected the in-situ daily profiles that had a corresponding satellite data by regrinding the in-situ data on the two satellites grid Global Ocean GSM-MyOcean and the Med Regional SeaWiFS RAN-MyOcean.The total corresponding in-situ and Global Ocean GSM-MyOcean data are 8910.
In the case of the satellite regional product Med Regional SeaWiFS RAN-MyOcean in order to perform the comparison between point in-situ observations and oceancolour data we decide to search in a surrounding area of 4 grid points around each single in-situ observation.The MedOC4 algorithm was tested and validated up to value of 10 mg•m −3 of Chl-a and therefore we masked the Chl-a data higher than the value.
The total corresponding in-situ and Med Regional SeaWiFS RAN-MyOcean data are 13611.Table 2 presents the principal information on the in-situ dataset used in the comparison with the two satellite products.

Complementing the EEA CSI023 In-Situ Indicator
The EEA classical CSI023 indicator is calculated based on profiles estimated by fluorometer and averaged in the summer 3 and in the first 10 metres depth at each station of sampling.The CSI023 is defined as the significant increasing or decreasing temporal trend for each station.
The in-situ Chl-a estimates are provided to the EEA through the Eionet network (http://www.eionet.europa.eu/).The workflow chart illustrating the methodology developed to complement the EEA CSI023 indicator with observations based on ocean colour from space is presented in

CSI023(+) Indicator Definition
CSI023(+) is computed from MyOcean ocean-colour gridded data as a temporal trend at each grid point starting from 1998.CSI023(+) consists of significant Chl-a trends in Mediterranean Sea based on a regional ocean-colour dataset (Med Regional SeaWiFS RAN-MyOcean 7 product).

Statistical Analysis
This section describes the method used to calculate the trend and climatology (Figure 2 left panel) from oceancolour data.For each grid point a yearly time series of summer mean Chl-a concentration was calculated for the period of 1998-2009 (12 points maximum, some summer values can be missing at some grid points if not enough values where available to calculate the summer mean).
To calculate summer values the monthly mean was calculated and the summer values were computed.A threshold on the minimum number of days required to create a valid monthly value was fixed at the level of 10 days and a threshold on the minimum number of years required to create a valid time-series was fixed at the level of 9 years.
For each grid point Sen's method calculates the slope of all data value pairs: , for all i and j where i > j (1) where:  X i and X j are the summer mean Chl-a concentration for years i and j;  Q m is the linear slope between X i and X j summer values;  m is the index of slope estimates. 3Summer is defined as the months from June to September for stations north of 59 degrees in the Baltic Sea (Gulf of Bothnia and Gulf of Finland) and from May to September for all other stations. 4In each figure of the paper the Sea of Azov has been masked because it is very shallow (maximum 15 metres depth) and ocean-colour products would be strongly affected by bottom reflectance; the White Sea has been masked because it was out of our domain of interest. 5The Global Ocean GSM-MyOcean product is described in Section 2.1.1. 6Chl-a areas are defined in Section 2.5. 7The Med Regional SeaWiFS RAN-MyOcean dataset is described in Section 2.   If there are n summer values X i in the time series we 2) Histograms of percentages of positive, negative and not significant values of CSI023(+) in the Mediterranean Chl-a areas (Figures 6(a)-(f)).get as many as

Chl-a Areas Description
Sen's estimator of slope is the median of these M values of Q.For each grid point the M values of Q are ranked from the smallest to the largest and Sen's estimator is: In the Mediterranean Sea, where a regional ocean-colour product at high spatial resolution was available, a special application of CSI023(+) was performed for the period of 1998-2009.This application uses the concept of "Chl-a areas", defined on the basis of river basins and political borders, where Chl-a trends are calculated.This method utilizes the high spatial resolution of the colour images and thus produces a more robust trend estimate than the pan-European trend indicator.
Ocean colour Chl-a and in situ trends were calculated with this method (Figures 3 and 4).
Chl-a areas were defined in the Mediterranean (68 Chl-a areas), using information on the River Basin Districts (RBDs), and political borders when RBDs were not defined.Moreover, 18 open-ocean sub-basins 8 were identified in the Mediterranean and used to design the Chl-a areas.Within each Chl-a area two sub-areas are defined: a coastal one (IN), from the coast to a depth of 30 metres and an offshore one (OFF), from a depth of 30 to 200 metres.A name composed of 3 parts was associated with each single Chl-a area as following: 1) name of the RBD or name of the country; 2) name of the sub-basin; 3) "IN" if it is the inshore part of the Chl-a area or "OFF" if it is the offshore part of the Chl-a area.
A Mann-Kendal statistics check [16,17] is applied to each Chl-a grid point to identify the statistically significant trends values at a 95% confidence level.
Chl-a standard deviation (STD) values were calculated using a non-parametric approach (68% confidence interval is considered as 2σ).Non-parametric Sen's slope [15][16][17] is normalized by the non-parametric Chl-a STD values of each grid point, so that CSI023(+) is then calculated as: Once CSI023(+) values have been calculated at each ocean-colour data grid point they are presented in two different ways: For example, the offshore part of the Maltese Chl-a area in the Strait of Sicily will be named "Malta-SSI-OFF".In Annex 1 we present details of Chl-a areas with a table listing their names.1) In a map over the European seas for values above the 95% confidence level (Figure 5(a)).

Validation against in-Situ Data
In the comparison between Global Ocean GSM-My-Ocean and in-situ data we proceeded as follows: for each in-situ data of a specific day we check if there is a corresponding satellite daily data for that day.Then we build a datasets of daily match-up in-situ and satellite data.Finally, summer and annual mean values are calculated by averaging all the corresponding (same day) in-situ and satellite data for that summer and year respectively.Summer and yearly mean satellite values are then compared with the corresponding in-situ values.
The validation of the ocean-colour products used to calculate CSI023(+) consists of the comparison of the oceancolour products with in-situ values from the EEA-Eionet databank.
In-situ Eionet observations are re-gridded on the satellite products spatial grids (1/24˚ when using the Global Ocean GSM-My-Ocean and 1 km × 1 km when using Med Regional SeaWiFS RAN-MyOcean).
In the case of the Med Regional SeaWiFS RAN-MyOcean and in-situ data comparison we enlarged the area for searching the corresponding satellite and in-situ data due to the lower coverage of the regional product than the global merged product and the spatial resolution is higher decreasing the chance of finding corresponding in-situ and satellite data.In order to perform the comparison we searched in a surrounding area of 4 grid points around each single in-situ observation.Then we build a datasets of daily match-up in-situ and satellite data.Finally summer and annual mean values were calculated by averaging all the corresponding daily in-situ and satellite data for that summer and year respectively.Summer and yearly mean satellite values were then compared with the corresponding in-situ values.
The MedOC4 algorithm was tested and validated up to value of 10 mg•m −3 of Chl-a and no higher value was present in the satellite dataset.Therefore we removed from the comparison the in-situ Chl-a data higher than this value.
The number of corresponding in-situ and ocean-colour daily values for the different satellite products is reported in Table 3.
The number of corresponding summer and annual mean values used in the comparison were reported in Table 3 together with the validation results.
We did not compare Global Ocean SeaWiFS RAN-MyOcean, because only monthly ocean-colour values were available.
Global Ocean GSM-MyOcean product validation results were presented in Figures 7 to 10, while Med Regional SeaWiFS RAN-MyOcean product validation was presented in Figure 10.The correlations were performed taking the logarithm on a decimal basis of the Chl-a values.The validation results were presented in Table 3.
The correlation between in-situ data and Global Ocean GSM-MyOcean ocean colour Chl-a concentration was relatively high (r 2 = 0.53) when analysis was carried out using all the data covering the entire European seas domain (Figure 7), although bias is equal to 1.10 (mg•m −3 ) and RMSD is 4.46 (mg•m −3 ).At the basin scales correlation values are lower, for the Mediterranean (Figure 8) r 2 is equal to 0.31 in the summer analysis (0.34 in the annual analysis) with a low bias (0.36 (mg•m −3 )) and RMSD is equal to 1.95 (mg•m −3 ).In the annual analysis the Mediterranean bias and RMSD increase to 0.42 (mg•m −3 ) and 2.15 (mg•m −3 ) respectively.In the Black Sea (Figure 8) we register the minimum number of data availability (46 for the full year and 23 for the summer period) and results showed that r 2 = 0.39 in the annual analysis and r 2 = 0.38 in the summer analysis.Bias is equal to −0.62 (mg•m −3 ) in the summer analysis indicating an overestimation of the in-situ data by the satellite ones, while RMSD is equal to 1.00 (mg•m −3 ).In the annual analysis of Black Sea data Bias is positive and equal to 0.41 (mg•m −3 ), while RMSD is equal to 1.86 (mg•m −3 ).The North-east Atlantic (Figure 9) shows relative high r 2 = 0.45 in summer and r 2 = 0.46 in the yearly analysis.Consequently, the Bias is low in the annual analysis (0.02 (mg•m −3 )) and in the summer period 0.23 (mg•m −3 ).Thus, RMSD for the North-east Atlantic is 1.43 in the annual analysis and 1.44 in the summer analysis.While in the North Sea (Figure 9) we found r 2 values similar to the other basins (0.40 in the annual analysis and 0.33 in the summer period analysis) showing that the higher Bias RMSD equal to 2.12 (mg•m −3 ) in the annual analysis (2.03 (mg•m −3 ) in the summer analysis) and RMSD is     The comparison of Med Regional SeaWiFS RAN-MyOcean with in-situ data in the Mediterranean Sea shows a correlation r 2 = 0.41 in summer analysis (Table 3 and Figure 10) and 0.36 in the yearly analysis, which are higher than that observed with the Global Ocean GSM-MyOcean data (Figure 8 and Table 3; r 2 = 0.31 for the summer analysis and 0.34 for the yearly analysis).Bias of the Mediterranean regional products (−0.04 (mg•m −3 ) in the summer analysis and 0.15 (mg•m −3 ) in the yearly analysis) and RMSD (1.44 (mg•m −3 ) in the summer analysis and 1.36 (mg•m −3 ) in the yearly analysis) are lower than that of the Global product.
In addition to the comparison of in Chl-a concentrations presented above we also propose a comparison of in-situ and satellite trends.To perform this comparison we identified locations at which there were at least 9 summer Chl-a mean values of corresponding in-situ and satellite data.The identified time series are 48 in the Global Ocean GSM-MyOcean and in-situ trends comparison.The results of this comparison were presented in a scatter plot in Figure 11 and in a map in Figure 12.
Results show that there are 16 locations in which both satellite and in-situ products detect positive trends, 17 locations in which both satellite and in-situ products detected negative trends, 6 locations in which in-situ products detected positive trends while satellite products

Results and Discussion
Ocean colour Chl-a concentrations and trends were compared with in-situ ones.Paragraph 4.1 presents the results of this comparison.The mean summer (May-September) Chl-a concentration (mg•m −3 ) in coastal European seas for the period of 1998-2009 was presented as calculated from the Global Ocean GSM-MyOcean ocean-colour product (Figure 2(a)) and in-situ data (Figure 2(b)).CSI023(+) was shown in the two formats: 1) CSI023(+) pan-European trend indicator and 2) CSI023(+) Chl-a areas trend indicator in the Mediterranean Sea.Paragraphs 4.2 to 4.4 present the abovementioned results.

Validation against In-Situ Data
In some locations the ocean-colour products differ greatly from the in-situ observations; one reason could be that in some cases in-situ data quality could be low (i.e.problem with calibration of instruments is expected due to the large number of data providers; some geographical coordinates of the single profiles were provided without precise rounding at the minutes value, and therefore with possible uncertainties on the geographical location of the sampling position) and this information is not provided with the Eionet datasets (i.e.data are not flagged for different quality levels).Moreover as explained in paragraph 2.1, ocean colour data used in the study was based on open ocean algorithms and therefore was not as accurate as needed in the very coastal zone where part of the in-situ data have been collected.For these reasons we expect that the validation will be worse in some coastal regions.
The comparison of Med Regional SeaWiFS RAN-MyOcean with in-situ data in the Mediterranean Sea shows a correlation which are higher than that observed with the Global Ocean GSM-MyOcean data (Figure 8 and Table 3).Bias of the Mediterranean regional products presented in Table 3 is lower than that of the Global product.The regional product which is also refers to as Med Regional SeaWiFS RAN-MyOcean thus seems to be preferable with respect to Global Ocean GSM-MyOcean in comparison with in-situ data.
The comparison of in-situ and satellite Chl-a trends shows that satellite products are able to capture the sign of the in-situ trend in the majority of the stations (68%) but the intensity of trends are not well captured by the satellite products.

Spatial Distribution of Chl-a Concentration in Coastal European Seas (Chl-a Summer Climatology)
The European seas summer Chl-a mean (mg•m −3 ) over the period of 1998-2009 (Figure 2 left panel) indicated the highest values in the Baltic Sea, in the southern North Sea and the Western Black Sea.The lowest values were reached in the Mediterranean Sea open-ocean areas.The Irish Sea, Bay of Biscay and Portuguese Atlantic areas also show high Chl-a values, although these are lower than the Baltic and North Sea.In the Mediterranean the low Chl-a values show a negative West-to-east gradient: this is very well known for this area [35,36].The Black Sea shows highest values in the North-western part of the basin (Danube area).The climatology based on ocean colour products (Global Ocean GSM-MyOcean) presented in this paper shows similar patterns of Chl-a compared to the ones based on SeaWIFS data [22] but lower concentrations of Chl-a are detected.
The in-situ climatology (Figure 2 right panel) reveals the sampling problem of the in-situ dataset: large areas with no data are present in the Black Sea, Mediterranean Sea, Atlantic Ocean, North Sea and Norwegian Sea.The poor coverage of the in-situ observations is due to the lack of monitoring programs and to the possible problems in the reporting systems.
The in-situ and ocean colour climatologies show similar values in most of the Baltic Sea with the exceptions of the Gulf of Finland and Gulf of Riga where ocean-colour products underestimate in-situ values.The southern part of the North Sea in-situ and ocean colour show similar mean values, although there is still an evident underestimation of the values by the satellite products.In the northern part of the Adriatic Sea, the in-situ and oceancolour products show similar values, while in the southern areas of the Adriatic ocean colour once again underestimates the in-situ values.In like manner, in the Ligurian Sea, the ocean-colour products underestimate in-situ observations, as they also do in the Tyrrhenian Sea.The fact that the in situ and ocean colour climatologies show similar values is encouraging because it will allow to complement the in-situ climatology with the ocean colour one thus providing a larger coverage and a full Pan-European picture of Chl-a mean values in the European Seas.

Comparison of Global Ocean GSM-MyOcean, Global Ocean SeaWiFS RAN-MyOcean and Med Regional SeaWiFS RAN-MyOcean Trends
To control the results of CSI023(+) obtained with Global Ocean GSM-MyOcean we decide to perform a comparison between the Chl-a trends obtained with Global Ocean GSM-MyOcean and the ones obtained with Global Ocean SeaWiFS RAN-MyOcean.
In Figure 3 we show both significant and non-significant trends to allow a better comparison of the two trends.

Global Ocean GSM-MyOcean and Global Ocean
SeaWiFS RAN-MyOcean trends show similar patterns but the intensity of trends detected using Global Ocean SeaWiFS RAN-MyOcean time series are higher than the ones obtained using Global Ocean GSM-MyOcean.Both products show positive trends in the northern part of the North Sea, in the North East Atlantic, in the central and southern part of the Baltic Sea, as well as in the Bay of Biscay, in the Alboran Sea, in the Ligurian Sea, in the southern part of the Gulf of Lion, in the southern Adriatic Sea and in the North-East Ionian Sea.Moreover negative trends were detected by the two datasets in the Black Sea, in most of the Mediterranean Sea, in part of the North-East Atlantic, in the English Channel and in the central part of the North Sea.Although, the main differences are found in the Gulf of Botnia and Gulf of Finland where Global Ocean SeaWiFS RAN-MyOcean shows negative trends while Global Ocean GSM-MyOcean shows positive trends.Besides, in the Skagerrak, Norwegian coasts and, in the English Channel, in the Po River mouth and in the western part of the Black Sea SeaWiFS showed similar negative sign but stronger negative trends in respect to GMS product.The positive trend in the central part of the Baltic Sea is stronger in Global Ocean SeaWiFS RAN-MyOcean than in the Global Ocean GSM-MyOcean dataset.
The comparison of the Global Ocean SeaWiFS RAN-MyOcean and the Global Ocean GSM-MyOcean reveals a general similarity among the two datasets and therefore the usage of the merged product (Global Ocean GSM-MyOcean), which offers longer time series, could be taken into account for continuing the update of the indicator also after the end of Sea-WIFS.On the contrary the differences highlighted in the Chl-a trends should be taken into account to identify areas of larger uncertainties where the satellite products might be less appropriate to estimate Chl-a trends (i.e.Baltic Sea).
The Chl-a trend analysis was also carried out in the Mediterranean Sea using the regional products Med Regional SeaWiFS RAN-MyOcean ocean-colour dataset for the period of 1998-2009 and thus was presented in Figure 4. Results obtained with the regional product appear similar to the ones obtained using the global product Global Ocean GSM-MyOcean presented in Figure 3(b).In the western Mediterranean, global and regional satellite products show similar trends.Some differences were detected in the Tyrrhenian Sea along the eastern coasts of Sardinia and Corsica where the global product presents positive trends not shown in the regional product.In the southern part of the Tyrrhenian basin, the regional product shows positive trends areas that appear larger than the ones of the global product.The Adriatic Sea trends appear very similar in the two datasets.In the Ionian Sea the two products show similar trends with the exception of a positive trend area south of Sicily that appears larger in the regional product, the Gulf of Gabes in which the regional product shows a negative trend and the Gulf of Sirte where the regional product highlights a positive trend that was not clearly visible in the global product.In the central part of the Levantine Basin, in southern part of the Aegean Sea and in front of the Nile River mouth the regional product show stronger positive trends respect to the global product.The differences among the two products Med Regional SeaWiFS RAN-MyOcean and Global Ocean GSM-MyOcean might be due to the fact that the global product consists of a merged product and uses different satellites (MERIS, MODIS and SEAWIFS) while the regional product is based on the SeaWiFS satellite product.Moreover the SeaWiFS regional product is based on a regional algorithm specifically calibrated for the Mediterranean Sea.
The results of the regional product Med Regional SeaWiFS RAN-MyOcean were also compared with the ones obtained in the Mediterranean with the Global Ocean SeaWiFS RAN-MyOcean.Results appear very similar in the entire basin; although small differences were detected in front of the Po River mouth where the regional product shows weaker negative trend values, while along the Italian Adriatic coast the negative trend appears stronger in the regional product.In the Gulf of Gabes the regional product shows a negative trend that which was detected as positive trend in the global Sea-WiFS product.In front of the Nile River mouth the regional product detects a stronger positive trend with respect to the one shown in the global product.
The main difference among the global (Global Ocean SeaWiFS RAN-MyOcean) and regional (Med Regional SeaWiFS RAN-MyOcean) is that the first is produced using the global ocean algorithm while the second is produced using the regional alrogithm (MedOC4, [28]).This difference might be the main reason for the differences detected in the trends.The use of regional product (Med Regional SeaWiFS RAN-MyOcean) is recommended because the algorithm (MedOC4, [28]) used to calculate it is calibrated for the specific sea (i.e.Mediterranean Sea in our study).

CSI023(+)
In this section we present CSI023(+) analysis in terms of pan-European trends map (based on the Global Ocean GSM-MyOcean product) and Chl-a area trends in the Mediterranean Sea (based on Med Regional SeaWiFS RAN-MyOcean product).CSI023(+) shows areas with decreasing Chl-a concentrations in the Black Sea, the Eastern Mediterranean, the southern part of the Western Mediterranean and the English Channel.Areas with increasing trends were observed in the northern Ionian Sea, the off-shore area of the Bay of Biscay, the North-East Atlantic, the North Sea, the Kattegat and the Baltic, trend values only above the 95% confidence level were presented, and therefore white areas in the map correspond to grid points with no significant trend.
The CSI023(+) positive trends in the central part of the Baltic Sea are consistent to the results presented in [13] but CSI023(+) seems to overestimate CHl-a trends in the Gulf of Bothnia.In interpreting the results in the Baltic Sea we need to consider that it is known that the global ocean products used in this study are not calibrated for the Case2 waters of the Baltic Sea and there fore large error are expected in this area.
In the Mediterranean Sea CSI023(+) did not show significant Chl-a trends in most of the basin.Negative values of CSI023(+) were detected in several off-shore regions (South of Sardinia, in the Ionian Sea, and in the Southern-Eastern Levantine Sea) and in the coastal reagion of the Eastern Mediterranean.The negative trend of Chla might be related to the general warming of Mediterranean waters [37,38] that leads to an intensification of water column stratification that consequently limits the uplift of nutrients from the deeper water column below the nutricline.
In the analysis performed with the regional product (Med Regional SeaWiFS RAN-MyOcean) the results of the slope of the trends presented in Figure 4 (statistical significance was not calculated, both significant and non significant trends are shown in this figure) some of the off-shore positive trends of Chl-a (i.e. in the Northern Ionian Sea, in the Western and Central Levantine) are possibly correlated with the intensification of cyclonic activities identified in the same areas in a similar period 1997-2006 [39] that could contribute to upwelling processes.In the Most of the off-shore negative trends are detected in areas (i.e.Algerian Basin, Central and Southern Ionian Sea) that are characterized by anti-cyclonic activities [39].
The negative trend in the Black Sea has been previously described in previous studies [22].The variations of the ecosystem of the Black Sea in the last decades have been reported [40][41][42] and they are related to complex set of pressures, including impacts from large scale climatic variability [43], variations in the nutrient and pollutant inputs [44,45].Also the appearance and evolu-tion of invasive species [46,47] was very relevant for the Black Sea ecosystem changes.
The Chl-a standard deviation over the period of 1998-2009 presented in Figure 5(b) shows that natural variability is high in the Baltic, the North Sea, the North-East Atlantic shelf areas and the northern regions of the Mediterranean and the Black Sea.Chl-a natural variability appears low in the Mediterranean, with a minimum in the Eastern Mediterranean basin.

CSI023(+) Chl-a Area Trends
In the Mediterranean, a special application of CSI023(+) was performed for the period of 1998-2009 when a regional ocean-colour product at high spatial resolution was available.CSI023(+) is calculated for each Chl-a area and results for this special application of CSI023(+) are presented with the classic CSI023 indicator mapping showing the percentage of decreasing, increasing and no-trend stations (Figures 6).The CSI023(+) Chl-a areas indicator for the Mediterranean Sea shows that about 52% of Chl-a areas do not show significant trends (i.e., significant trends found for less than 10% of the grid points in each area).Increasing Chl-a trends (i.e., found for equal or more than 10% of the grid points in each area) were detected in 9% ( 12) of the Chl-a areas (in Egypt, Greece, Tunisia, Malta and Turkey coast lines) see Figures 6(e) and (f).The increasing Chl-a trend in the coastal areas of Egypt, higher in the in-shore Chl-a area than in the off-shore Chl-a area (Figure 6(f)), could be related to the discharge of nutrients from the Nile.The increasing trend in Greece is detected in the Epirus areas (North West Greece) and might be related to nutrient river inputs such as the ones from Kalamas River in which eutrophication signs were obvious during the warm period [48].
Decreasing significant Chl-a trends (i.e., found for equal or more than 10% of the grid points in each area) were detected in the 39% (53) Ch-a areas: Italian coasts in the Adriatic, Slovenian coasts (Figures 6(d The decreasing Chl-a trends in the coastal areas in the Mediterranean Sea are possibly related to a combination of factors such as the reduction of river flow due to agriculture and industrial use [49], impact of climate change (i.e.reduction of rivers flow and increase of ocean temperature) and the implementation of EU policies for the reduction of nutrients loads in the water bodies (i.e.WFD).In this sense the decreasing trend of Chl-a in the Ebro coastal area (see Chl-a area Ebro-IBS IN and OFF) could be related to the decrease in nutrients concentration and flow [49].Also in the Northern Adriatic Sea our results confirm the decrease in Chl-a concentration as discussed in [50] that might be related to the reduced freshwater discharges, together with phosphorus being banned by Italian law in the mid-1980s and the general improvement in sewage treatment.

Conclusions
The CSI023(+) indicator based on the ocean colour products has been developed to contribute to the state-ofthe-environment assessment and, in particular, to monitor eutrophication trends.We performed the comparison of Chl-a concentrations and trends estimated through remote sensing with in-situ ones.The comparison presented indicated that the global ocean-colour algorithm seems to underestimate Chl-a concentration (Figure 2).To evaluate the performance of satellite products to estimate Chla trends at Pan-European Scale, we compared different satellite products (i.e.Global Ocean GSM-MyOcean and Global Ocean SeaWiFS RAN-MyOcean), and found that the trends intensity detected using SeaWiFS time series were higher with the GSM once, nevertheless, GSM and SeaWiFS trends show similar patterns.The differences highlighted in the northern Baltic Sea point out the lower confidence of the satellite products in detecting Chl-a trends in this region.The comparison of the global products (i.e.Global Ocean GSM-MyOcean) with the regional product available in the Mediterranean Sea Med Regional SeaWiFS RAN-MyOcean show similarities in the general patterns while specific differences have been highlighted.The comparison of Global Ocean SeaWiFS RAN-MyOcean and the corresponding regional product Med Regional SeaWiFS RAN-MyOcean show strong similarities.Thus, some local differences were discussed.
The CSI023(+) indicator based on the Global Ocean GSM-MyOcean shows the capability of detecting significant negative and positive terms for Chl-a trends in European seas, thereby allowing us to complement the in-situ CSI023 covering only part of the European coastal areas.The CSI023(+) indicator also shows a large area with decreasing Chl-a concentrations in the Black Sea, the Mediterranean, the English Channel and the northern part of the North Sea; whereas a large area with increasing trends were observed in the Bay of Biscay and the Baltic.Trends estimated by ocean colour products are compared with trends estimated by in-situ data (Figures 11 and 12).Hence, they show a clear differences as well as the capability of ocean colour products to capture the sign of the in-situ trend in the majority of the comparison locations.Moreover, the Med Regional SeaWiFS RAN-MyOcean and CNR regional daily dataset used in the Mediterranean Sea seems to show a good comparison with in-situ data (Table 3 and Figure 10).Therefore, further investigations are needed to be carried out to compare the trends at the level of Chl-a areas when more coastal in-situ Chl-a time series will be available and when new satellite products built upon coastal algorithm will be available.Validation should continue in the future, perhaps with a dedicated data-collection exercise.Chl-a area analysis offers the possibility of synthesizing CSI023 (+) information thus focuses attention on coastal areas.The analysis also revealed the need for regional oceancolour products availability to develop support of the EEA indicator.In addition there is potential in a long-term trend analysis based on ocean colour, as large-scale, and in some cases even regional-scale, changes appear to be captured by the satellite images.It is clear, however, that in order to build confidence in this analysis it needs to be based on the best possible regional products.As not all the used MyOcean products consist of complete reprocessed data, it is planned that as soon as new regional datasets are available they will be used to calculate CSI023(+).
Although Global Ocean GSM-MyOcean and Med Regional SeaWiFS RAN-MyOcean and CNR were available as daily products, not all the information from ocean colour has been used yet, because, CSI023(+) using summer mean values in analogy with classical in-situ CSI023.However, for the future we plan to use full daily satellite Chl-a estimate resolution to evaluate changes in the statistical significance of the results.

Figure 1 .
Eleven steps were identified starting from the ocean colour and in-situ data download (Step 1) to the presentation of CSI023(+) indicator (Step 11).The different aspects of the methodology illustrated in Figure 1 are presented in this Section 2 of the paper.

Figure 1 .
Figure 1.The workflow chart presenting the methodology developed to complement the EEA CSI023 indicator with observations based on ocean colour from space.Eleven steps are defined starting from the data download (Step 1) to the delivery of CSI023(+) (Step 11).The column on the left presents the activities performed on the ocean colour data while the column on the right presents the activities related to the in-situ data.The comparison of in-situ and ocean colour Chl-a concentration (Step 2) and Chl-a trends (Step 7) are identified as common activities.

Figure 2 .
Figure 2. Summer (May-Sept.)mean Chl-a concentration in European seas for the period of 1998-2009 from the MyOcean ocean-colour global dataset (a) and the Eionet in-situ dataset (b).

Figure 4 .
Figure 4. Mediterranean trend values (mg•m −3 •y −1 ) (Sen's slope) for the period of 1998-2009 calculated using the regional product Med Regional SeaWiFS RAN-MyOcean.Both significant and non-significant trends are displayed.

Figure 7 .
Figure 7. Summer Chl-a in-situ concentration and the ocean-colour Global Ocean GSM-MyOcean Chl-a concentration for all selected stations in the European seas (Log10 mg•m −3 ).The black line is the best linear approximation between ocean colour (y) and in-situ values (x).

Figure 8 .
Figure 8. Summer Chl-a in-situ concentration and the ocean-colour Global Ocean GSM-MyOcean Chl-a concentration for all selected stations in the Mediterranean Sea (Log10 mg•m −3 ).The black line is the best linear approximation between ocean colour (y) and in-situ values (x) in the Mediterranean Sea and the green line is the one for the Black Sea.

Figure 9 .
Figure 9. Summer Chl-a in-situ concentration and the ocean-colour Global Ocean GSM-MyOcean Chl-a concentration (Log10 mg•m −3 ) for all selected stations in the Baltic Sea (blue diamond), North-east Atlantic (orange circle) and North Sea (green cross).The best linear approximation between ocean colour (y) and in-situ values (x) is in blue for the Baltic Sea, black for the North-east Atlantic and yellow for the North Sea.

Figure 10 .
Figure 10.Summer log10 Chl-a in-situ concentration (x axis) compared with the ocean-colour Med Regional SeaWiFS RAN-MyOcean log10 Chl-a concentration (y axis).Log10 Chl-a values were expressed in mg•m −3 .The black line is the best linear approximation between ocean colour (y) and in-situ values (x).equals to 6.80 (mg•m −3 ) in the annual analysis (5.99 (mg•m −3 ) in the summer analysis).The Baltic Sea (Figure9) the correlation noticeably drops (r 2 equal to 0.02 in the annual analysis and 0.15 in the summer analysis) and Bias is equals to 1.55 (mg•m −3 ) in the annual analysis and 1.26 (mg•m −3 ) in the summer analysis, while RMSD is equal to 5.10 (mg•m −3 ) in the annual analysis and 4.74 (mg•m −3 ) in the summer analysis.The comparison of Med Regional SeaWiFS RAN-MyOcean with in-situ data in the Mediterranean Sea shows a correlation r 2 = 0.41 in summer analysis (Table3and Figure10) and 0.36 in the yearly analysis, which are higher than that observed with the Global Ocean GSM-MyOcean data (Figure8and Table3; r 2 = 0.31 for the summer analysis and 0.34 for the yearly analysis).Bias of the Mediterranean regional products (−0.04 (mg•m −3 ) in the summer analysis and 0.15 (mg•m −3 ) in the yearly analysis) and RMSD (1.44 (mg•m −3 ) in the summer analysis and 1.36 (mg•m −3 ) in the yearly analysis) are lower than that of the Global product.In addition to the comparison of in Chl-a concentrations presented above we also propose a comparison of in-situ and satellite trends.To perform this comparison we identified locations at which there were at least 9 summer Chl-a mean values of corresponding in-situ and satellite data.The identified time series are 48 in the Global Ocean GSM-MyOcean and in-situ trends comparison.The results of this comparison were presented in a scatter plot in Figure11and in a map in Figure12.Results show that there are 16 locations in which both satellite and in-situ products detect positive trends, 17 locations in which both satellite and in-situ products detected negative trends, 6 locations in which in-situ products detected positive trends while satellite products

Figure 11 .
Figure 11.Summer Chl-a in-situ trends and the ocean-colour Global Ocean GSM-MyOcean Chl-a trends for all selected stations in the European Seas (mg•m −3 •y −1 ).The black line is the best linear approximation between ocean colour (y) trends and in-situ trend values (x).The significance of the trends is was presented with different symbols: In-situ (IS) and Ocean Colour (OC) not significant-blue diamond; IS significant (95%) and OC not significant-Black square; IS and OC both significant (95%) green triangle; IS not significant and OC significant (95%) blue triangle.

Figure 12 .
Figure 12.Comparison of Chl-a trends calculated estimated from ocean colour data Global Ocean GSM-MyOcean in the European seas (a) for the period of 1998-2009 (mg•m −3 •y −1 ) and from the in-situ dataset (b).Large circles correspond to significant (95% confidence level) trends while small circles correspond to non-significant trends.The palette indicates different trend's range (mg•m −3 ) with different colours.detectnegative trends and 9 locations in which in-situ products detected negative trends while satellite products detected positive trends.Position of the locations where trends were identified was shown in Figure12 together

4. 4 . 1 .
CSI023(+) Pan-European Trends CSI023(+) over the period of 1998-2009 was presented in Figure 5 (a).It is worthy to note that the STD of Chl-a (mg•m −3 ) is naturally high in regions of high Chl-a.With regard to relative changes, CSI023(+), provides a Pan-European picture that highlights the spatial variability of the trends and minimizes the natural Chl-a concentration differences among the European Seas.

Table 3 . Validation of the ocean-colour data products used in the calculation of CSI023(+) 9 .
9Correlation coefficient (r 2 ), slope and intercept of the regression lines correspond to the one presented in Figures7-10.Numbers of summer/annual mean values compared are reported in the table for each of the seas and European Sea and for each of the aggregation period (summer or year).Root Mean Square Distance (RMSD) and BIAS are calculated for Chl-a absolute values (mg•m 3 ), not logarithmic.The Mediterranean Sea results are reported for the Global Ocean GSM-MyOcean product (GSM) and Med Regional SeaWiFS RAN-MyOcean product (Med Reg).