Trends in Ocean Colour and Chlorophyll Concentration from 1889 to 2000, Worldwide

Marine primary productivity is an important agent in the global cycling of carbon dioxide, a major ‘greenhouse gas’, and variations in the concentration of the ocean's phytoplankton biomass can therefore explain trends in the global carbon budget. Since the launch of satellite-mounted sensors globe-wide monitoring of chlorophyll, a phytoplankton biomass proxy, became feasible. Just as satellites, the Forel-Ule (FU) scale record (a hardly explored database of ocean colour) has covered all seas and oceans – but already since 1889. We provide evidence that changes of ocean surface chlorophyll can be reconstructed with confidence from this record. The EcoLight radiative transfer numerical model indicates that the FU index is closely related to chlorophyll concentrations in open ocean regions. The most complete FU record is that of the North Atlantic in terms of coverage over space and in time; this dataset has been used to test the validity of colour changes that can be translated to chlorophyll. The FU and FU-derived chlorophyll data were analysed for monotonously increasing or decreasing trends with the non-parametric Mann-Kendall test, a method to establish the presence of a consistent trend. Our analysis has not revealed a globe-wide trend of increase or decrease in chlorophyll concentration during the past century; ocean regions have apparently responded differentially to changes in meteorological, hydrological and biological conditions at the surface, including potential long-term trends related to global warming. Since 1889, chlorophyll concentrations have decreased in the Indian Ocean and in the Pacific; increased in the Atlantic Ocean, the Mediterranean, the Chinese Sea, and in the seas west and north-west of Japan. This suggests that explanations of chlorophyll changes over long periods should focus on hydrographical and biological characteristics typical of single ocean regions, not on those of ‘the’ ocean.


Introduction
Oceanographic data archaeology (salvaging historical ocean data) is, according to Parker [1], a critical requirement for climate change research when the data (e.g., temperature and salinity) cover long time periods. Changing climate will affect marine life, and this, in turn, plays a crucial role in climate control [2]. As a matter of fact, phytoplankton accounts for nearly half of Earth's total primary productivity [3]. Phytoplankton biomass is nowadays readily derived from satellite data [4]; indeed, chlorophyll provides the best index of phytoplankton biomass for primary productivity studies, as was recently confirmed by Huot et al. [5]. Long-term changes in ocean primary production may well have important consequences for the ocean food chain, as well as for the global carbon cycle [6].
One of the so-called 'Essential Climate Variables' listed by the World Meteorological Organization (WMO) to detect biological activity in the ocean's surface layer is ocean colour [7]. More recently an important aspect of climate research has been the detection of ocean colour from space with derived products such as water transparency that is influenced by coloured dissolved organic matter absorption and chlorophyll concentration. Water transparency can be used to predict the depth of the upper ocean mixed layer [8]. This layer plays a critical role in the flux of energy between atmosphere and ocean.
Changes in water colour are caused by a change in the composition of optically active substances [9]: suspended particulate matter, pigments (chlorophyll-a, b and c and carotenoids) in algae [10], and dissolved organic matter [11]. The blue colour of oligotrophic oceans is caused by domination of selective absorption and scattering of the water molecules. Chlorophyll colours the water green (or even red by the carotenoids peridinin in dinoflagellates), and the presence of coloured dissolved organic matter (CDOM) will, in high concentrations, account for the absorption of most of the blue part of the incoming sunlight, resulting in brownish-coloured water.
The colour of seas and oceans is, next to water transparency, salinity and temperature, one of the few oceanographic parameters that have been recorded for over a century [12]. One outstanding long-term record (since 1931) has been produced by the Continuous Plankton Recorder Surveys in the North Sea and the North Atlantic [13], [14]; other records are the Hawaii-and Bermuda Atlantic Ocean Time-Series (HOTS and BATS [15], [16]) and the California Cooperative Oceanic Fisheries Investigations (CalCOFI) [17]. All these series have from the launch of the Coastal Zone Color Scanner been supplemented with satellite-derived oceanographic products for temperature change and, as to biological characteristics, for ocean colour and its derivative, chlorophyll concentration. From the late 1990s more sophisticated imaging ocean colour satellite sensors such as the Sea viewing Wide angle Field of view Sensor SeaWiFS, the MEdium Resolution Imaging Spectrometer MERIS and the Moderate Resolution Imaging Spectroradiometer MODIS sample the ocean's colour. Since the launch of SeaWiFS, ocean colour from Long-Term Trends in Ocean Colour and Chlorophyll PLOS ONE | www.plosone.org space with rigid measurement protocols [18] has been archived as a standard oceanographic parameter. With revisit times of 1 to 3 days, a swath width of around 1,000 to 2,000 kilometres and a spatial resolution of 250 to 1,000 meters, the ocean has never before been sampled optically on such a scale [19].
The record of chlorophyll changes over time is limited; only a few authors have presented analyses of trends in chlorophyll over limited numbers of years. Venrick et al. [20] found, by means of trend analysis of standing-stock chlorophyll, a doubling of phytoplankton chlorophyll in the central North Pacific gyre between 1964 and 1985. We must bear in mind that the analysis of chlorophyll has become more accurate with the introduction of HPLC in the late seventies of the past century [21]. Falkowski and Wilson [22] found areas (10 by 10 degrees boxes) with significant positive and negative average annual changes in Secchi depth since 1900, but overall their estimated annual change averaged over the whole North-Pacific Ocean's chlorophyll concentration is small and not significant. A similar transformation was recently used by Boyce et al. [23] in order to derive world-wide changes in ocean phytoplankton biomass since 1899.
Since 1889 a colour comparator method has been used to establish the colour of the sea, through Forel-Ule scale observations [24], [25]. Besides Forel and Ule, the inventors of the scale, Krümmel (1893), Luksch (1901) and Von Drygalski (1897) [26], [27], [28] were among the first to use this method in the open sea. Krümmel made his observations during the plankton expedition of the Humboldt Society in 1889. At that time only the blue to bluegreen part of the scale (scale numbers 1 to 11) was used to classify the colour of open water. An example of Krümmel's contoured sea colour map is shown in Figure 1. This map shows surprisingly well the oligotrophic Sargasso Sea gyre ('deep cobalt blue, FU-1'), as well as the outflow of the Congo River ('greenish-blue, FU2-3'), and the Canary and Greenland Currents ('dark green, FU4-7'). The data were collected in the Atlantic during the Plankton Expedition of 1889 and the colour of the sea was expressed as a percentage of a yellow potassium chromate solution added to a blue copper-sulphate solution. Luksch collected observations over eight years during the 'Pola' expedition (1890-1898) in the Mediterranean, Aegean-and Red Sea and Von Drygalski performed observations during the Greenland Expedition of the Geographical Society of 1891-1893. Since 1890 the Forel-Ule scale became the most commonly used and most simple scale to determine, through comparison, the colour of seas, lakes and rivers.
Over the years only a limited number of geographical maps, based on interpolation of sets of Forel-Ule data, became available. The U.S. Navy Hydrographic Office published in the 1950s three atlases on marine geography, including Forel-Ule contours of the Seas of Japan, Korea and Indochina [29], [30], [31]. A more comprehensive effort was made by Frederick [32], who presented contoured Forel-Ule maps. She had around 24,000 Forel-Ule observations at her disposal to create contour maps of averaged sea colours, but she did not analyse trends. Wernand and Van der Woerd [33] have recently found, in the North Pacific Ocean, a greening trend from 1930 until 1955, and since then a bluing, based on an analysis of 17,171 FU-observations that are available for that ocean.
The aim of the work presented here was to visualise changes in ocean colour, notably its chlorophyll component, over the longest possible period of time by an analysis of the 'forgotten' globe-wide ocean colour obtained by the Forel-Ule method. The temporal variation of ocean colour has been analysed from 1889 till the year 2000, based on subsets in a total of 221,110 FU observations collected globe-wide. First the data selection and quality control procedures are described and a model is introduced that provides simple relations for the conversion from FU-scale number to Chlorophyll concentration. Subsequently, the data are grouped in 28 seas and oceans. All this should eventually allow speculation on mechanisms behind basin-scale external forcing of the plankton abundance near the ocean surface.

Materials and Methods
With the Forel-Ule method the colour of natural waters is compared to a scale, that consists of twenty-one tubes filled with mixtures of coloured chemical solutions. The method has recently been described by Wernand and van der Woerd [34]. Most FU observations were retrieved from oceanographic and meteorological databases archived by NOAA-NODC [35]. This dataset contains, besides date, geographical position and Secchi depth (SD), the FU index (coded 1 to 21; see introduction). Two new attributes were added for each entry, indicating the meteorological season and the sea area. The geographical naming follows the conventions of the U.S. Geological Survey (USGS) that is based upon the latitude and longitude at which the observation took place. This dataset contains 220,440 observations collected from 1907 to 1999, inclusive.
From The re-formatted dataset is referred to as dataset-0 and contains, after the first quality check, 221,110 FU globally collected observations of open water, lakes and rivers. Figure 2 shows a bar chart of the total number of FU observations collected within a decade and per season. The period from 1889 to 1899 period covers eleven years in order to gain as much available information of the earliest period of ocean colour observation. Data filtering Subsequently, dataset-0 was differentiated into three geographical areas by applying Basin Masks (BM). The Basin Masks are used to filter out observations within land boundaries (lakes and rivers), and observations close (100 km or 500 km) to the coast. The Basin Masks were created and executed in ArcGIS 9.3 software [36] by means of the ArcGIS clip function.
First, BM1 was created to extract data from dataset-0 free of observations within the land boundaries just mentioned. This dataset is referred to as dataset-1 and contains 219,839 sea observations. Table 1 shows the distribution of the FU observations per sea area with details on the total and seasonal numbers of FU observations, also expressed as a percentage of the total. The definition of seasons follows the meteorological convention. We highlight some results of the major oceans and seas. The selected oceans and seas are the Barents Sea, Bering Sea, East China Sea, Indian Ocean, Mediterranean, North-and South Atlantic, Northand South Pacific, Norwegian Sea, Pacific Coast, Philippine Sea, South China Sea, Sea of Japan, Sea of Okhotsk and the Yellow Sea.
Second, BM2 was chosen pragmatically to extract open sea observations from dataset-0, at a distance of more than 100 km off-coast, to avoid effects on water colouration of anthropogenic pressure in the coastal zones such as locally increased nutrient loading (eutrophication) that tends to enhance phytoplankton biomass or high sediment loading caused by changes in land use or erosion. These phenomena result in short-term, drastic colour changes, no doubt influencing the open sea FU values. This dataset is referred to as dataset-2 and contains 61.434 open sea observations. The last mask, BM3, was again chosen pragmatically to extract observations from dataset-0 at a distance of over 500 km from the coast to include the oceans but at the same time avoid the influences of mixing with the differently coloured water of nearby seas. This dataset will be referred to as dataset-3 and contains 21.971 open ocean observations. An example of the shapes of the three data extraction masks is given in Figure 3. One must bear in mind that by clipping the data to the dimensions of the mentioned masks the number of observations per sea area diminish by roughly a factor of 4 in case of BM2 and a factor 10 in case of BM3. However, to establish a realistic colour of the bulk water of a sea or ocean, not influenced by its nearby coastal or surrounding seawater (with much higher FU numbers) this action is unavoidable. The data extraction masks, datasets name and the number of observations included in the dataset are summarized in Table 2. From dataset-1, a global map of the FU-colour distribution was established through an Inverse Distance Weighted interpolation (IDW) [37], [38]. IDW interpolation explicitly implements the assumption that observations that are close to one another are more alike than those that are further apart. The IDW determines cell values by calculating a linearly weighted combination of nearby FU-observations. The weight is calculated by taking the inverse of the distance between the cell and FU-observation, raised to a power (usually equal to two). This power option controls the significance of known points on the interpolated values, based on their distance from the output point. Datasets-2 and -3 are statistically analysed per sea area and per year. Data were binned per year to eliminate any influence of a seasonal cycle in the colour of the sea. The arithmetic mean FU values FU with the confidence interval of the mean (95%) and of the observations (95%) [39] of each sea were calculated. Linear trends over time were calculated by a, to the number of observations, weighted-linear regression of the (FU) values. The  weighted linear regression takes into account the number of observations as weights [40].
To check the presence of a monotonously increasing or decreasing trend in the observations the nonparametric Mann-Kendall trend test per sea and ocean has been applied as proposed by Mann [41] and further studied by Kendall [42] and Sen [43]. The test was found to be a good tool for trend detection and has been utilized by many other investigators in the analysis of various types of environmental data [44], [45], [46]. For the statistical analyses including the Mann-Kendall test we applied the Excel add-in XLSTAT-software (http://www.xlstat.com). In the trend analysis all results are presented under the assumption that the FUmeasurements have been made independently. The potential influence of temporal and spatial autocorrelation in the dataset on the simple linear regression analysis has not been taken into account, mainly because this is a very complex and certainly not well defined exercise [47]. Therefore the provided temporal trends should be treated with caution; we are aware of the fact that they provide a simplified description of reality.

Bio-optical Modelling
An FU colour change from a lower to a higher scale index indicates a greening of seawater caused by the presence of lightabsorbing substances, notably chlorophyll and CDOM. In order to simulate the relation of chlorophyll or CDOM and the FU index the HydroLight-EcoLight [48] modelling software was adopted. This model solves the radiation-transfer equations of light in natural waters and computes radiance or irradiance distributions and derived quantities, like remote sensing reflectance R RS and FU [49]. By varying the model's input parameters, like chlorophyll, CDOM or mineral concentrations, a relation can be established between concentrations and R RS or the FU index. HydroLight-EcoLight employs mathematically sophisticated invariant imbedding techniques to solve the radiative transfer equation. Details of this solution method can be found in Mobley's Specific absorption -a p *(z, l) a CDOM *(z, l) a p *(l) for calcareous sand   book on Light and Water [50] and in his 'Comparison of numerical models for the computation of underwater light fields' [51]. Both models, HydroLight and EcoLight, generate the same output; notice the difference, namely that HydroLight computes the full 3D-radiance distribution, whereas EcoLight computes only irradiances. Two simulation runs were made to establish the influence of either chlorophyll or CDOM on the colour of seawater i.e. on the remote sensing reflectance R RS and on the FU number. For modelling the oceanic waters, the new IOP model ''NEW'' case 1 was used. This model is based on a combination of i) the results of modelled particle absorption as given by Bricaud et al. [52] and ii) the results of modelled particle absorption described in more recent publications by Morrison and Nelson [53] and Vasilkov et al. [54]. The ''NEW'' case 1 two-component IOP model consists of 1) pure water [55] and 2) chlorophyll-bearing particles with covarying CDOM and detritus. One of the model's output parameters is the remote sensing reflectance R RS (l) which was automatically convoluted with CIE1931 curves into a chromaticity coordinate set x, y and accordingly into a FU scale number (see [34] for more details). In this way a relation between chlorophyll and FU could be established for oceanic waters alone, and is used in the trend analysis to convert FU back into a chlorophyll concentration.
To simulate the chlorophyll and mineral concentration in open sea water the values of these properties were set respectively at 0.1 mg m 23 and 0 or 0.2 g m 23 . The CDOM absorption (a 440 ) was varied between 0.01 and 1.0 m 21 . Again, the calculated remote sensing reflectance R RS (l) was convoluted with CIE1931 curves into a chromaticity coordinate set x, y and into a FU scale number [34]. In Table 3 the EcoLight model parameters, including the scattering and absorption properties of the constituents are tabulated for both models. Figure 4 shows the geographical positions of the observations extracted under respectively BM3 (.500 km, black triangles) and    analysis was also done for each separate season, the number of seasonally binned observations is given too. In Figure 6 these numbers are provided for all seas that were treated under BM2 (.100 km off-coast).

Results
In order to understand the oceanic trend analysis of chlorophyll, we first present the Ecolight modelling results of the conversion between FU and chlorophyll concentration (chapter 3.1).
Hereafter, a global view of averaged ocean colour of the oceans and seas utilizing all observations from 1889 to 2000 is presented in chapter 3.2. The chapter is followed by a regional trend analysis, first per ocean and then per sea.

FU Modelling
Varying chlorophyll concentrations were used as input for the two-component case 1 water model, for the four-component case 2 water model CDOM as input was varied (a 440 between 0 and 1 m 21 ) with fixed values for components three and four mimicking open sea bio-optical values. For each run the EcoLight modelling software output parameters are the spectral remote sensing reflectance R RS (l), as shown in Figure 7, and the FU number shown in Figure 8. Figure 7a shows R RS spectra generated by the case 1 model with varying chlorophyll, Figure 7b shows R RS spectra generated by the case 2 model with varying CDOM absorption and chlorophyll and mineral concentrations set to 0. Figures 7c and 7d show R RS spectra generated by the case 2 model with varying CDOM and chlorophyll set to 0.1 mg m 23 and minerals set to 0 and to 0.2 g m 23 respectively.
The results presented in Figure 8A show a change in chlorophyll concentration from 0.1 to 1 mg m 23 resulting in a shift in FU index from 1 to 4. An exponential fit to the model results provides the following relation:  This relation takes into account the fact that in the open ocean chlorophyll bearing particles are co-varying with CDOM and detritus, whereas detritus contributes only 5% to 20% to the nonwater absorption coefficient [59]. Equation 1 is used to transform the acquired FU into an average chlorophyll concentration per year for oceanic regions. After having consulted a chlorophyll composite of an entire MODIS mission, with maximum chlorophyll concentration of a few mg m 23 , we decided to set the upper limit of observed FU to 7 (<6 mg m 23 chlorophyll) to avoid extreme outliers which could greatly influence the results; remember the exponential relation between FU and chlorophyll (Equation 1). In this way, observations taken above extremely dense phytoplankton blooms (e.g., coccolithophorids) could be avoided. As an example (BM2), for the Indian Ocean no data were found .FU7, for the North Atlantic 25 observations were omitted and for the North Pacific (total of around 23,000 obs.) we omitted 240 observations for our chlorophyll calculations.
For the seas, results will be presented as FU only, because they generally have more optically active constituents and no single relation between FU and chlorophyll can be established. However, to demonstrate the influence of CDOM, on the FU scale index we present case 2 model results for 3 modelled parameter settings.   Table 4). The Barents Sea is the most greenish sea (bottom) and the Equatorial Pacific is the most bluish ocean (top). doi:10.1371/journal.pone.0063766.g011 Long-Term Trends in Ocean Colour and Chlorophyll PLOS ONE | www.plosone.org The parameter settings are identical to the settings used to generate the results in Figures 7b, 7c and 7d. In Figure 8B we see that a change in CDOM absorption from 0.01 to 0.1 m 21 results in a shift over the first 4 FU scale colours in the case that CDOM, next to water itself, is the only absorber of light. Through an exponential fit CDOM 440 absorption can be calculated with a coefficient of determination R 2 = 0.999 according to Setting both the background chlorophyll concentration to 0.1 mg m 23 and the mineral concentration to zero the FU numbers change from FU = 3 to FU<8 (open circles in Figure 8B), which implies that the low chlorophyll concentration (set at 0.1 mg m 23 ) already colours the water to a FU value of 3. This outcome is somewhat different from the case 1 modelling results, shown in Figure 8A, where the same chlorophyll concentration corresponds to FU = 1. We must bear in mind that in the case 1 model, different from the case 2 model, CDOM is co-varying with chlorophyll. Through a power law fit CDOM 440 absorption can be calculated with a coefficient of determination R 2 = 0.96 according to In the last set of results of the case 2 Ecolight model, the background chlorophyll concentration was set to 0.1 mg m 23 and the mineral concentration was set to 0.2 g m 23 . For FU numbers The CDOM-FU relations above are presented to help with the interpretation and comparison of a derived FU in seas for which more background knowledge of the area-specific inherent optical properties are available.
In the following we show the significance i.e. accuracy and applicability of the FU scale observations presented here by visualizing i) an overall (1889-1999) monthly FU and by ii) visualizing an overall monthly chlorophyll concentration applying From the findings presented above we conclude that FU is a good proxy for chlorophyll in the open ocean. We have already argued that besides water itself chlorophyll (mostly inside the cells of phytoplankton) gives the ocean its colour [9]; but we must bear in mind that next to chlorophyll the usually non-covarying CDOM can, through absorption, have a similar influence on the colour of the sea and therefore on FU.

Century averaged ocean colour
To establish a global view of the FU data collected between 1889 and 1999 all observations contained within dataset-1 were interpolated according to the IDW technique as described in material and methods. The inverse distance weight function was set to a power law of index 2. The search radius, that limits the number of FU observations used for calculating each interpolated FU value, was set to respectively 2.5 and 4 degrees with an output grid size of 0.5 degrees. The difference between results of interpolation using a search radius of 2.5 and 4 degrees turned out to be negligible, although for under-sampled seas the applied interpolation technique results in FU values for areas where no observation were performed. A map of IDW interpolated FU values, with a search radius of 4 degrees, is presented in Figure 10. The southern oceans, highly under-sampled, are shown in white.
A first analysis of Figure 10 shows that the lowest FU numbers are found in the Equatorial Pacific and the Indian Ocean (FU,2) and in the Equatorial Atlantic (FU,3). Note the extremely high FU values (16 to 21) east of the Ob and Yenisei estuaries in the Kara Sea (75uN, 70uE), most likely caused by very high CDOM values [61]. Also in the Baltic, where CDOM is the major light absorber [62], FU numbers of around 14 are encountered. The Atlantic shows increasing FU values (up to FU = 5) towards higher  latitudes, most likely caused by the influence of more extensive phytoplankton blooms. Of all data extracted under either BM3 or BM2, the mean FU 1889-1999 per sea area was calculated, with 95% confidence intervals for the mean and for the observations. The results of the statistical analysis show that the Equatorial Pacific and Indian Ocean are the bluest oceans of our globe with a mean FU of 1.6 and 1.7, respectively. With the results of our statistical analysis one can classify and rank the oceans and greater seas in terms of their mean FU colour, presented in Figure 11 (similar to the IDW interpolation of Figure 10, however, less visible in this coloured graph). This figure shows that in the 20 th century the North Atlantic Ocean has been the greenest ocean of our planet while the Barents Sea was the greenest open sea of all.
In order to support these classifications with recently collected global data (after 1999), we show in Figure 12 an entire MODIS mission composite (4 July 2002 to 30 Jun 2010) for the CDOM index [63] and for the chlorophyll-a concentrations. Note that in case 1 waters a mean relationship exists between the CDOM content and the chlorophyll concentration, anomalies in this relation, for both case 1 and case 2 waters are given by this socalled CDOM index W [63]. In the hyper-oligotrophic Equatorial Pacific gyre [64], located within the greater white box, we see a W,1 and very low chlorophyll concentrations, that are expected in blue oligotrophic oceans (around 0.02 mg m 23 ). Within the smaller white box, indicating the Barents Sea, we see a CDOM index around 3.5 and chlorophyll values around 5 mg m 23 colouring the sea green. Areas with an increasing CDOM index and or increasing chlorophyll concentration ( Figure 12) are generally compatible with an increasing FUindex (Figure 11).

Yearly averaged regional chlorophyll and ocean colour
To establish the influence of the choice of Basin Mask on the trend analysis for open ocean water, we compared the results of dataset-2 (BM2) and dataset-3 (BM3). Because the conversion from FU index to chlorophyll concentration depends on the CDOM absorption, only oceanic waters with a CDOM index W,2.5 were converted with Equation 1; see also the global CDOM composite in Figure 12. For the remaining seas, as indicated in bold in Table 1, we applied the trend analyses on BM2 extracted data. For all data presented here, per sea or per ocean, the overall trend covers the whole period of data collection. The least-squares regression lines, indicated by a full blue or green line, indicate either a bluing or greening of the ocean/sea under investigation. For the oceans bluing/greening means decrease/ increase in chlorophyll, for the rest of the seas this can mean a decrease/increase in FU which can either mean a decrease/ increase in chlorophyll or CDOM or a combination of both.
In the next paragraphs the results of our statistical analysis are presented per year per sea area for illustrational purposes only. For the Mann-Kendall test all single observations were used (no binning) In some cases, the BM2 extracted data is enclosed within the boundary of a neighbouring sea or ocean, hence a second data selection was needed, done with the ArcGIS data selection tool, to exclude overlapping data. The total number of extracted FU observations per sea or ocean can therefore be slightly different from the number of observations given in Table 2 and Figures 5  and 6. Data extracted under the BM3 mask concerns only the oceans, data extracted under the BM2 mask concerns both oceans and world seas.
Oceans. In Figure 13 the results of the statistical analysis are presented for the Atlantic and Pacific Ocean and in Figure 14 for the Indian Ocean, per year between 1889 and 1999. The FU has been converted to chlorophyll according to Equation 1. For all oceans analysed, the yearly FU values, including the number of observations, regression coefficients and trend lines with the 95% confidence interval curves for the mean and for the observations, are presented in the graphs. At the left-hand side the graphs are based on data extracted under the BM2, graphs on the right are based on data extracted under BM3. The number of observations is indicated near the yearly means. The trend line with regression coefficients, where y is the chlorophyll concentration in mg m 23 and x is the year is indicated in each graph. For the oceans Table 4 shows the regression coefficient a with standard deviation delta a, Kendall's tau coefficient (represents the degree of concordance between two columns of ranked data) and the Mann-Kendall The results for the Indian Ocean, presented in Figure 14, show a bluing ocean for both BM2 and BM3 data with an average decay of 20.0021 mg m 23 per year. Between 1898 and 1993 the Indian Ocean showed significant decline in chlorophyll, from 0.37 to 0.18 mg m 23 , at the end of the observational period.
Seas. For all seas analysed, the yearly FU values, including the number of observations, regression coefficients and trend lines with the 95% confidence interval curves for the mean and for the observations, are presented in Figure 15 and Figure 16. Table 5 shows the general trend of a sea, bluing or greening, the weighted (no. of observations) linear regression coefficient a with standard deviation delta a, Kendall's tau and the two-tailed p-value for the analysed seas. The last two columns of the table shows the period over which observations were collected and the start and end value of the FU number. BM2 was used as data extraction mask. The FU values are not converted into chlorophyll concentrations in waters with a CDOM index W.2.5 (Figure 12), which occurs in most of the seas. In these seas a change in CDOM as well as chlorophyll can cause shifts in sea colour. Figure 15 shows for both the northern Barents and Bering Sea a bluing trend. The Barents Sea, the greenest open sea of our planet, shows a dramatic bluing (p = 0.036) with a factor of 20.059 FU per year from < 9.0 (1925) to <5. 1 (1992). Between 1935 and 1998 the Bering Sea is bluing (p,0.0001) with a factor of 20.0211 FU per year from <5.2 to <4.6. The Norwegian Sea, geographically situated below the Barents Sea, has an average FU value of 5.05 and shows no significant colour change (p = 0.302) over the period .
The Mediterranean was geographically split into a western and eastern part, based on the known existence of an east-west oligotrophy gradient at 15u East [65], [66], [67]. In total 824 FU observations (BM1) between 1889 and 1999 were collected. From these only 237 observations (BM2) were analysed. Our analyses shows that the western part is slightly greener compared to the eastern part; FU MED:West = FU MED:East +0.8FU. BM1 and BM2 related FU both show a significant greening between 1890 and 1991. The results for the BM2 dataset are plotted in Figure 15: the whole Mediterranean has been greening significantly over the past century, with an increase of 0.016FU per year since 1890. (p,0.0001). Interestingly, this greening rate is almost similar to the results found for the North Atlantic (BM3) with an increase of 0.02FU per year. As an indication of the potentially related changes in biomass, the FU data of the Mediterranean were converted into chlorophyll, assuming its water being Case1 (see Table 6). With a FU = 1. For the Sea of Japan an overall greening is found between 1933 and 1999 (p,0.0001). The Sea of Okhotsk, a bluish-green sea, connected with the Sea of Japan, shows an overall greening between 1935 and 1995 with an undulating pattern in FU not seen so pronounced in other seas (minor trend, p = 0.4).
We emphasise that the trend analysis presented here is only a simplified description of the data. It was already shown in a previous paper on FU measurements made in the North Pacific [33] that oceans can show a much more complex temporal behaviour.

Discussion and Conclusions
We conclude that the ocean colour dataset of Forel-Ule scale observations that has been established over the period 1889 to 1999 offers a unique opportunity to investigate worldwide temporal changes of colour at the sea and ocean surface in the 20 th century, from 1889 to present; long-term colour changes that we show can safely be related to concentrations of chlorophyll, the most used phytoplankton biomass proxy. Bio-optical modelling, based on the state-of-the-art radiation transfer code called 'EcoLight' that is based on a two-component and four-component model has demonstrated how the concentration of chlorophyll and CDOM relates to ocean colour (FU scale index). However, we must keep in mind that the modelled chlorophyll behaves exponentially. This means that in case a FU observation was done over an oceanic plankton bloom (FU5-FU7, personally observed), the model generates for a FU5 observation 1.7 mg m 23 of chlorophyll and for a FU7 observation 6.5 mg m 23 of chlorophyll. Indeed, for the clearest oceanic waters (FU = 1 to 4), where algal pigments are the dominant factor of light absorption, a highly significant exponential relation was found between the FU index and chlorophyll concentration ( Figure 8A). Thus, significant changes in the average chlorophyll concentration per year could be reconstructed (Figures 13 and 14). Although, for water closer to land (,200 km off-coast) the FU-Chlorophyll relation can be improved by a better estimate of the CDOM index and specific inherent absorption by phytoplankton in each ocean, the observed variation in colour is real and indicates a change in phytoplankton biomass at the surface.
The Forel-Ule classification of natural waters, that was started long ago, in the late nineteenth century, can, as we have shown, be regarded as an oceanographic record as reliable for interpretation as the long-term records of salinity, temperature and water transparency that are often used to relate ocean water properties to long-term climate changes that receive so much attention since the IPCC reporting on this phenomenon. The ocean colour record of the FU-scale observations can hardly be replaced by other data if only because it goes back all the way to 1889. We have shown here that observations made with this colour scale are not simply 'subjective', as earlier has been stated when Secchi disc visibility records were evaluated [68].
The Forel-Ule classification accuracy is highlighted by our analysis of the North Atlantic dataset; the annual cycle of Forel-Ule ocean colour and derived chlorophyll results (Figure 9b) are fully in line with the annual cycles established by satellite ocean colour sensors [60]. Also for the Northern Pacific the FU-based seasonal changes of calculated chlorophyll concentration compare well with the values of chlorophyll changes reported for this region [33] (Wernand and van der Woerd, 2010b).
We are well aware of the fact that in some seas the number of observations available and used for our analyses is rather low, but they are at the same time the only data related to phytoplankton biomass change in these ocean regions. The Mann-Kendall analysis was performed to evaluate the ocean colour or chlorophyll trend of an ocean or sea. The outcome of the analysis indicates that the North Atlantic and Equatorial Atlantic have been greening over the whole past century. The increase in North Atlantic chlorophyll after the 1950's is consistent with the trend of increase in the Phytoplankton Colour Index (PCI) obtained during the Continuous Plankton Recorder (CPR) survey over the past 50 years [14] [69], [70], [71], [72]. Further evidence of a greening Atlantic is given by Edwards et al. [73] in a 'SAHFOS ecological status report' of 2007/2008 ('there has been a large increase in phytoplankton since the late 1980s in most regional area of the North Atlantic'). For the five oceans under investigation a significant trend was found (p,0.0001, significance level of 0.05). Out of eleven investigated seas seven show trends with a p,0.0001.
Our findings do not indicate a global trend in the ocean's colour as established on the basis of the FU record, and therewith in its chlorophyll contents. We cannot confirm the recent report of an ocean-wide phytoplankton decrease [23]; Boyce et al. found a global decline of phytoplankton over the past century that they related to increasing sea surface temperature. Differences in data filtering techniques may have been the cause of the considerable discrepancy between the results of Boyce et al. [23] and the results we present here. Boyce et al. based their conclusions on datasets from which only data obtained in waters ,25 m deep or ,1 km from the coast were excluded. This means that their results close to coasts may have been biased by chlorophyll-rich eutrophic coastal zones. This will have an impact on trends in oligotrophic regions, where chlorophyll varies much less. In contrast to Boyce et al. [23], we have omitted data ,100 km (seas and oceans) or ,500 km (oceans) from the coast.
Henson et al. [74] recently concluded that detection of climate change-driven trends in the satellite data record is confounded by the relatively short time series and large inter-annual and decadal variability in productivity. According to Henson et al., recently observed changes in chlorophyll, primary production and the size of the oligotrophic gyres cannot be explicitly ascribed to the impact of global climate change. They suggest that time series of at least 40 years are needed to distinguish trends from natural variability.
Gregg and Conkright [60] re-analysed the CZCS global ocean chlorophyll product, using SeaWiFS compatible atmospheric correction methods, and made a first quantitative comparison of the decadal trends in global ocean chlorophyll over the periods 1979-1986 (CZCS) and 1997-2000 (SeaWiFS). They reported a decrease in global chlorophyll concentrations from the CZCS records to present of about 6%; they found larger reductions in the northern high latitudes and an increase in chlorophyll in the low latitudes. Their conclusions are not in line with ours but we must bear in mind that results are based upon the results of two relatively short periods -much shorter than the shortest period presented by us. Antoine et al. [19] also used CZCS and SeaWiFS data; they reported an overall increase of the world ocean average chlorophyll concentration by about 22%.

Recommendations
The dataset presented here contains 220,000 Forel-Ule observations; 160,000 remain unexploited for the time being. It would be of great interest to investigate this latter dataset, collected in coastal and shelf seas (,100 km off-shore) globe-wide, to establish the possible influence of coloured near shore waters on open-ocean colour; moreover, with a refinement of the model results presented here it should be possible to classify water quality not only in the open ocean but also in seas where coloured dissolved organic matter (CDOM) is not co-varying with chlorophyll, usually sea regions where human influence plays a role, notably eutrophication by nutrient input of rivers, and by sediment and CDOM load changes. However, we must bear in mind that such a refinement is not an easy task as chlorophyll and CDOM absorption can have a similar influence on sea and ocean colour and therefore on FU values. Future work should concentrate on explanations of driving forces behind the decadal trends in open-ocean colour presented here, and anomalies in the ocean colour record over long terms of up to a century. Hindcasting is the only way to understand and predict changes in the ocean's ecosystem in its relation to climate change.
Long-term series are limited; a literature search indicates that 'long-term' generally means: between 5 and at most 30 years. It is therefore not only challenging but also necessary to combine archived long-term data series of basic water quality parameters such as Secchi disc visibility depth, chlorophyll, or phytoplankton abundance (e.g., the colour index recorded since 1948 during the Continuous Plankton Recorder Survey of SAHFOS) to establish inter-annual or inter-decadal means from which results can be compared to our results. An attempt to combine historic water quality datasets has been presented recently by Boyce et al. [75] (2012). Data of other long-term time series, namely the Hawaiiand Bermuda Atlantic Ocean Time-Series (HOTS and BATS, both since 1988) and the California Cooperative Oceanic Fisheries Investigations (CalCOFI since 1949) should help to achieve this goal.
Finally, we like to make a plea, substantiated in the present paper, for the reintroduction of the FU-scale to classify the colour of seas and oceans. Ocean colour classification by means of a single numerical value instead of a hyper-spectral classification will facilitate the interpretation of long-term ocean colour data series and at the same time facilitate a connection between the present and the past. MERIS satellite Forel-Ule mapping can play a role here, through validation and coupling of historic observations to the present era of satellite-based monitoring.