The Baltic Algae Watch System – a remote sensing application for monitoring cyanobacterial blooms in the Baltic Sea
The aim of this paper is to describe our satellite based system for monitoring of cyanobacterial blooms and to present the method we used for evaluating and analyzing the dataset collected during the period 1997-2006. To compare different year and see trends normalized index for bloom duration, extent and intensity were developed.
Cyanobacterial blooms, usually dominated by the toxic Nodularia sp. and non toxic Aphanizomenon sp., are commonly found during summer in the Baltic Sea. Although blooms have occurred in the Baltic Sea for at least 7000 years B.P. [1,2], it has been debated if the areal extent and intensity of the blooms have increased because of anthropogenic eutrophication [3-5]. In recent years there has been an increasing interest from both the public and authorities since the peak of the bloom usually coincide with the summer holidays when people reside near the coast. The thick surface accumulation of cyanobacteria is a major nuisance, especially when the accumulations aggregate along the coast, releasing bad odour and prevent swimming. At worst, toxic substances are released from Nodularia that can kill pets and cattle drinking polluted water. Humans and especially children who get in contact with contaminated water usually suffer from skin irritation, stomach problems and symptoms similar to the flu [6,7].
Another environmental concern with the cyanobacterial blooms in the Baltic Sea is the increasing input of nitrogen to the water. Most cyanobacteria use dissolved molecular nitrogen (N2) as an additional nutrient source, which allows them to bloom during summer when growths of other phytoplankton are limited by nitrogen [8,9]. Estimates of nitrogen fixation by cyanobacteria may be as high as 400 000-700 000 ton/year . In 2000, rivers transported ~660 000 ton nitrogen to the Baltic Sea and ~300 000 ton/year due to atmospheric deposition . Hence, it is likely that nitrogen fixation is the single largest source of nitrogen to the Baltic Sea.
The signature of cyanobacterial blooms in the Baltic Sea is distinct even if the character can differ depending on the vertical position of blooms in the water column and its concentration. Cyanobacteria have the ability to move vertically to find favorable light or nutrient conditions, subsurface blooms often cause a discoloring. Towards the end of a bloom cyanobacteria can loose their buoyancy control and aggregate at the surface forming dense accumulations. The blooms are formed into meandering structures or patches (Langmuir circulation) by winds and currents (see Fig. 1).
Since cyanbacterial blooms have a high spatial and temporal variability it is difficult to provide accurate near real-time (NRT) information using only conventional shipborn sampling. Collecting data using remote sensing techniques enables monitoring of vast areas on a regular basis in NRT. The main disadvantage is that most sensors cannot ”see” through clouds, which limits the monitoring ability at high latitudes, which have a high presence of clouds . Remote sensing data is also limited to the top meters of the sea surface; blooms at greater depths are not visible .
Material and method
2.1 Detection of cyanobacteria blooms
Observations of cyanobacteria blooms forming heavy bloom or scum at or close to the sea surface in the Baltic Sea is associated with high optical scattering properties rather than discolouring by absorption properties. This experience emanates from experimental work using satellite data and airborne or ship based reconnaissance (visual inspection and photos). Early attempts to demonstrate this experience can be found in [3,14-18].
The scattering properties are however intrinsically difficult to measure with existing in situ instruments simultaneous with remote sensing data. Both spatial texture and near surface accumulations of heavy blooms are the main problem for in situ sampling, whereas atmospheric correction methods generally fail for satellite data . Hence, validation of in situ scattering and absorption characteristics of these heavy phytoplankton blooms are still unsolved, mainly due to lack of adequate methodology and technology.
We therefore conclude that heavy blooms of cyanobacteria close to or at the sea surface can be detected by their scattering characteristics in the red and near infrared part of the electromagnetic spectrum by satellite sensors such as NOAA-AVHRR, MODIS, MERIS and Hyperion (cf. ).
Detection of algae blooms using optical satellite data is restricted to cloud free areas, since direct sun radiation is absorbed and reflected by clouds, hindering radiation information from the upper water column to reach the satellite sensor. Time series of optical satellite data can therefore be biased and/or aliased, concerning spatial and temporal coverage of blooms. These effects remain to be proven but are subjects outside the scope of this investigation. Here we can provide information based on the ten years of time series data (1997 - 2006). We find that algae blooms are detected in 20 to 60 % of the cloud free satellite data. Typically cloud free days occur as often as cloud covered days during July to August, but varying from year to year. Furthermore, the temporal coverage of blooms varies between a few days up to 20 days, always considerably less than available number of cloud free days. We conclude that heavy blooms are not under sampled in cloud free areas and these areas are on average evenly distributed over the open Baltic Sea during summer months.
Our hypothesis is therefore that optical satellite data, available several times a day, in NRT, (e.g. NOAA-AVHHR) can be used to estimate the occurrence of these blooms in time and space.
The method to process NOAA-AVHRR data to detect cyanobacteria is basically a supervised classification algorithm that relies on multiple thresholding and difference in visible, near infrared and thermal channels to distinguish cyanobacteria accumulations from other common features as land, clear water, clouds, haze, sun glint and error pixels. The method includes manual interpretation and correction of each satellite scene .
A common problem when using NOAA/AVHRR data is the relatively coarse pixel size (~1 km2), which make bloom detection impossible near the coast. Hence, the method can only be applied for monitoring of offshore blooms. Hindcasting with MERIS or MODIS satellite data with a spatial resolution of 300m (Available at ESA, ENVISAT Web File Server) and 250-500m respectively (Available at NASA, Level 1 and Atmosphere Archive and Distribution System) are therefore used to visually validate the daily bloom maps.
To get an estimate of cloud-free observations in relation to bloom and cloud observations the cloud cover is screened for every analyzed image. A cloud mask produced by the NoWCasting Satellite Application Facility (NWC-SAF)  is loaded during the image processing and the cloud cover over the Baltic Sea is extracted. To further minimize the cloud problem hindcasting with composite maps of the daily bloom maps during the past 7 days are produced showing number of observations in each pixel. Hence, even if cloud cover makes it impossible to detect any blooms, the composite image can picture the algal situation .
2.2 The monitoring system
The web based monitoring system was constructed to provide a direct overview of the blooming situation at the default view but also offering more thoroughly information about the sea state. Maps showing the extent of cyanobacterial bloom are updated daily throughout the blooming season (June-September). During cloud free conditions one map per day is enough but if the cloud situation changes rapidly revealing new information maps can be updated up to three times per day. Enclosed to each map a short text is presented describing the current bloom situation and the 7 days composite. The latest sea surface temperature image from Ocean & SeaIce – Satellite Application Facility (O&SI-SAF) is also presented. The default view of the website is presented in Fig. 2.
The Swedish national marine monitoring programme performs monthly water sampling cruises in the Baltic Sea. In situ measurements of phytoplankton are analyzed onboard and an algae report is also produced shortly after the cruise giving a brief indication on the phytoplankton species composition, which can be used for sea-truthing of the satellite products. These reports together with cruise reports are available at the website.
The main sources of information assembled at the Baltic Algae Watch System can be found in Fig. 3.
To be able to compare the blooms between different years definitions of bloom normalized duration (T), extent (A) and intensity (I) have been developed. Based on the yearly summaries where the area (ai) is equal to the extent that is covered by surface accumulations of blooms during (i) number of days, the normalized duration, Eq. (1), and extent, Eq. (2), is given. Where (i) ranges from 1 to the maximum number of days with bloom observations during the current year. The intensity, Eq. (3) is given in “extent days” or km2days .
The highest normalized intensity was recorded during the bloom 2006, while both 2005 and 2006 had the longest duration. The largest extent was noted in 1998 but the bloom was interrupted by cold and windy weather and the duration of the bloom was short, see Table 1. Ten days of data is missing from the year 1999 due to antenna malfunction. It is likely that blooms were present during this time gap since blooms were detected both before and after the interruption which would have given both a larger extent and longer duration e.i. higher intensity. It has also been suggested  that the bloom during 1999 had the highest cumulative blooming area during the period 1979-1984, 1998-2006. Table 1. Presents a comparison of normalized extent, duration and intensity of cyanobacterial blooms between the period 1997-2006. The results are based on the yearly summaries. Year 2001 is missing due to satellite antenna problems. Year Extent (A) km2 Duration (T) days Intensity (I) km2 days 1997 8 009 6 49 821 1998 12 240 3 31 214 1999 9 621 7 68 342 2000 6 061 3 17 461 2002 5 448 5 26 474 2003 6 386 7 42 841 2004 5 384 3 17 726 2005 6 572 8 55 319 2006 9 149 8 73 452
Based on the annual summaries and the normalized indexes for extent, duration and intensity hindcasting products are produced after the blooming season. These products are part of the marine environmental Indicator Fact Sheets  that Helsinki Commission is publishing at http://www.helcom.fi to provide information about recent state of and trends in the Baltic marine environment.
Summary and conclusion
The Baltic Algal Watch System has shown that that the combination of satellite imagery, observations, forecasts, models and tools assembled at one website, is a powerful nowcasting tool for monitoring and prediction of cyanobacterial blooms. The system has been appreciated and is frequently used by the Marine Information Offices, authorities, environmental managers, media and public. The number of visitors of the website has increased rapidly since the system became operational in 2002. During 2004 there were about ~43 000 visitors, which increased during the strong bloom in 2005 to ~120 000 visitors. In 2006, which had the highest normalized intensity during the 10 years analyzed the number of visitor increased further too over 149 000. About 70-80% of the visits at the website were made during July and August when the blooms are most frequent and tourists are spending their holidays near the coasts of the Baltic Sea. Future work will focus on incorporating NRT satellite data from MERIS and MODIS with higher spatial and radiometric resolution than NOAA-AVHRR and obtaining homogenous data sets of the mapping procedure using different sensors. Algorithms for certain cyanobacterial pigments e.i. phycocyanin [27,28] yet not evaluated in the Baltic Sea could improve the mapping procedure. Further evaluations on the impact of clouds on the mapping procedure together with NRT in situ validations with ferrybox systems in the Baltic region  are required.
The main conclusions and experience from 5 years of operational monitoring of cyanobacterial blooms are:
• NOAA-AVHHR is suitable for nowcasting purposes of algal blooms because of high number (10-12) overpasses per day, NRT-access to data, availability of relatively long compiled time series were data continuously are being added. • The cyanobacterial blooms showed large inter-annual variability. Years with massive blooms can be followed by years with only minor blooms. • Weather conditions are an important factor for the formation of cyanobacterial blooms. Warm, sunny and calm weather in combination with high concentrations of phosphorus provides suitable conditions for cyanobacterial blooms. • Cloud cover, which is a major problem when using NOAA-AVHRR data, can be minimized using 7-days composite as a complement to the daily bloom maps. • The increase of the blooms both in area extent and intensity that have been debated is difficult to prove. Satellite data has been available for about 3 decades, and this period is too short to show trends or natural variability of the phenomena. Though, low N:P ratios observed in the Baltic Sea since 2002 could be one explanation to the high intensity bloom during recent years.
We would like to acknowledge MISTRA RESE-miljömål for financially supporting the development of the Baltic Algae Watch System and Ove Rud for providing satellite data from 1997-2000 and getting us started. We would also like to acknowledge the Swedish Coastguard for providing us with aerial photographs of cyanobacterial blooms.