Dec 14, 2006 - short-term temperature fluctuations (Bernard 1963,. Haq 1972). .... dishing DH (1978) Biological effects of climate change. Rapp. Pv Reun Cons ...
Vol. 129: 269-281.1995
MARINE ECOLOGY PROGRESS SERIES Mar Ecol Prog Ser
Published December 14
Hydroclimatic relationships with planktonic time series from 1975 to 1992 in the North Sea off Gravelines, France Genevieve Le Fevre-~ehoerff Frederic Ibanez2, Paola Poniz2. Jean-Marc ~ r o m e n t i n ~ 'l*,
'IFREMER Centre de Brest, BP 70, F-29280 Plouzane, France 'Observatoire Oceanologique, Station Zoologique, BP 28, F-06230 Villefranche-sur-Mer, France
ABSTRACT: Planktonic time series collected since 1975 off Gravelines, France, were analysed in the context of climatic changes. Principal component analysis and the cumulative sums method were applied to 30 data sets which included meteorological and hydrological parameters and planktonic descriptors. Eight groups of biological descriptors were linked to external parameters. The study of winds highlighted the importance of cyclonic and anticyclonic influences. Seven year cycles appeared clearly in air temperature and in species abundance with different trophic levels for species such as Ternora longicornis, Oikopleura dioica and Sagitta setosa, whereas an 11 yr cycle appeared for Pseudocalanus minutus. A 3 to 4 yr cycle was detected for atmospheric pressure, precipitation and some species. This cycle may or may not be superimposed on another (for example, on an 11 yr cycle).
KEY WORDS: Time series . Zooplankton. Meteorology Hydrology. Southern North Sea
INTRODUCTION
At the turn of the century, British researchers had already set up a system for climatic and biological monitoring on a regular basis for the North Atlantic, Channel and North Sea regions, such as the E l station located in Plymouth, England, followed by a series of CPR (Continuous Plankton Recorder) measurements. These enabled long-term fluctuations to be monitored (Russell 1935) and then summarised to highlight the links between climatic changes (Lamb 1969, 1985) and biological changes in terms of fish stocks, plankton and the intertidal benthos (Cushing 1961, 1978, Southward 1963, Russell et al. 1971, Cushing & Dickson 1976).The periodic phenomena described as the 'Russell Cycle' are still referred to today. Since plankton recorders have remained in operation, various simple or complex mathematical processing methods have enabled assessment of the results and comparisons to be drawn
O Inter-Research 1995
Resale of full article not permitted
(Colebrook 1981, 1982, Colebrook & Taylor 1984, Ibanez 1984, Ibanez & Etienne 1992, Fromentin et al. 1993, Ibanez et al. 1993). The plankton record from Gravelines, France, presented here covers part of the chronological data gathered from 1975 onwards. The study was begun following the decision to build a high-power (5500 MW) coastal nuclear power plant. The objective was first to assess site characteristics (e.g. hydrology, microbiology, plankton, fish, benthos) before the plant was put into operation (1980) and then to monitor its impact. The plankton data were used to study seasonal species variations (Toularastel et al. 1977, Le Fevre-Lehoerff & Quintin 1981, Woehrling & Le Fevre-Lehoerff 1993), for the impact study (Le Fevre-Lehoerff & Woehrling 1991, Le Fevre-Lehoerff et al. 1993b) and for spatial comparison between Gravelines and other Channel sites (Le Fevre-Lehoerff et al. 1983, 1993a). The data used within the context of regional climatic changes (Fromentin & Ibanez 1994) concern only the area outside of the power plant's impact area.
27 0
Mar Ecol Prog Ser 129: 269-281, 1995
MATERIALS AND METHODS Sampling area. Gravelines is located between Calais and Dunkerque (Fig 1) in the southern part of the North Sea on a low, sandy coastline. Semi-diurnal alternating tidal currents run parallel to the coastline and the general residual current flows toward the North Sea. Coastal counter-currents are exceptional, and can divert water from large northern European rivers (e.g. the Escaut, Belgium) southwards. Irregular flows of 350000 m3 d-' issue from a small canalised river in France called the Aa. Sampling frequency. The sampling approach adopted was determined by the geographical situation and research objectives. The first aspect was a multidisciplinary approach during sampling surveys. The second involved an assessment of the environment's natural variability by replicate samples and samples taken during tidal cycles. Man-made constructions implemented from 1973 to 1980 modified the shoreline, and consequently, our research strategy. From 1980 on, impact studies based on spatial comparisons rather than on time-based variations (Chardy & Menesguen 1984) were mainly used in this research. Surveys were conducted between September 1974 a n d October 1986 at a rate of 10 to 12 times yr-' Their number then dropped to 3 yr-' Since measurements are less accurate offshore, weekly measurements were taken at the entrance to the plant. Few measurement points were used here. The non-impacted zone succes-
Fig. 1. Geographical zone of study
sively shifted on a geographical scale. This does not appear to be of great importance, considering the movement of water masses with the tide. All measurements are included in 1 file list: Points 9 and 17 from September 1974 to December 1975, Point 2 from August 1976 to May 1978 and from April 1982 to September 1993, and Point 5 from August 1976 to February 1982 (Fig. l ) . Measurement techniques. Conventional oceanographic techniques (Aminot & Chaussepled 1983) were applied in field studies and laboratory analyses. The more specific techniques employed are described in a collective publication (Belsher et al. 1986). Samples were collected using Niskin bottles and a WP2 plankton net (200 pm). Hydrologically speaking, the shallow water (20 to 25 m) only required samples taken at a few levels (3) and for vertical tows of zooplankton. The filtered volume was evaluated using a TSK flowmeter attached to the net opening (Le Fevre-Lehoerff 1985). To estimate the suspended matter (SM) in the sea water, samples were filtered through a precombusted 0.45 pm Whatman filter. The samples of SM, nutrients, chlorophyll and zooplankton biomass were frozen as soon as possible after being collected. The zooplankton samples which were to be used for counting and determination of species were put in formalin (8% final concentration). Nutrients were analysed using a Technlcon Autoanalyser and expressed in mm01 1-' for each nutrient, NO3-, NO2-, NH,+, PO,-, SiOH-. The chlorophyll biomass was assessed using a fluorometer and expressed in mg m-3. The zooplankton biomass was dried for 48 h at 60°C, weighed (mg m-3), then homogenised and analysed with a Carlo Erba elemental analyser; results were given as carbon and nitrogen percentages of dry weight. Densities (no. m-3) were estimated for the species determined. Data. Meteorological parameters: The following parameters were available: air temperature, atmospheric pressure, precipitation, number of days with east and west winds. Readings of monthly values were recorded in Dunkerque from January 1974 to December 1992, except for wind, which was measured from 1976 onward. Hydrological parameters: Conservative hydrological parameters such as water temperature ("C) and salinity (PSU) were available for the period 1976 to 1992. Other non-conservative parameters such as nitrogen salts NO3-, NO2-, NH,' (from 1976 to 1988) as well a s PO, and SiOH- (from 1976 to 1986) were expressed in mm01 I-'. Overall parameters were estimated for SM (from 1976 to 1986) in mg 1 - l , chlorophyll (mg 1 - l ) from 1976 to 1992, phaeopigments (mg 1-l) (1978, 1986), and zooplankton biomass in dry weight (mg m-3) from 1976 to 1992.
Le Fevre-Lehoerff et al.: Hydroclimatic relationships with planktonic time series
27 1
Biological descriptors: The total zooplankton Table 1. Mean, standard deviation, maximum (all in ind. per 10 m3) and percentage of zero ( X )in counted categories species present at Gravelines comprised approximately 150 taxa (Le Fevre-Lehoerff et al. 1976). This inventory was drawn up at the start of the Category Mean SD Max. study; subsequently only about 40 taxa were reguSagitta selosa ( S A G I ) 65 277 3687 larly counted. This selection also supplied the 14 Lanice conchilega (LANI) 256 1079 10625 biological descriptors (Table l ) used in the analysis. Pleurobrachia pileus (PLEU) 7 27 229 These species were well represented in terms of Paracalanus parvus (PARA) 292 899 7476 Pseudocalanus minutus (PSEU) 424 1347 11565 frequency and number of individuals. Other abunCentropages hama tus ( C E N T ) 1890 94 17 129090 dant species were rejected because their presence Temora longicornis (TEMO) 6324 19433 163158 was too erratic. Many observations were lacking Cyphonaut larvae bryozoa (CYPH) 184 763 8853 for the time series (110 monthly sampled values out Cirriped nauplii (CINA) 172 564 6382 of 220). This was especially true for later years, Cirriped cypris (CICY) 140 776 11625 Acartia clausi (ACAR) 1544 4817 46010 when sampling dropped to 4 times in 1986, then to Euterpina acutifrons (EUTE) 686 6504 112960 2 in 1989. Abundance data were log-transformed in Oncaea sp. (ONCA) 122 625 9674 order to stabilise the variance. To render the data Oikopleura dioica (OIKO) 897 2728 22116 more regular, median observation values within each month were used (sampling unit: 1 mo). Numerical analyses. Cumulative sums: Chronologidescriptors. The 11 hydrological parameters were cal series in oceanography often contain missing valentered as supplementary variables. Correlations with ues. In spite of this it is important to detect the general the first 3 components graphically displayed the proximities of the 25 variables (Dauvin & Ibanez 1986). trend in order to compare biological and environmental conditions. Toward this aim the simple method of the cumulated function (Ibanez et al. 1993),which does RESULTS not require special conditions, can be applied. For a chronological series with data X, sampled for each t ( t varying between 1 and n), a reference value k is choYear-to-year changes in meteorological and hydrological parameters sen (which could be e.g. the mean of the series, or an indicative value). After subtracting k from each data Meteorological parameters point, the residuals are added successively:
For successive values equal to k, the curve will be horizontal, and for successive values lower than k, the slope will be negative, and vice versa. If an erratic, high value exists (large peak or dip), this will also show up as a marked anomaly on the curve of cumulative sums. However, small fluctuations are smoothed. For information purposes, Bravais-Pearson correlations were calculated using series of cumulative sums. The r significance value for 100 degrees of freedom was 0.19 at the 5 % threshold and 0.25 at the 1 % threshold. However, this conventional test could not be used since the terms of cumulative sums are autocorrelated. An extremely scrupulous approach would have included a Monte Carlo simulation producing an r value distribution. However, this would have been excessive since the cited correlations taken for our research were all considerably higher than the previous statistical thresholds. Principal component analysis (PCA): An analysis was conducted by principal components, using the series of cumulative sums from the 14 biological
Air temperature at Dunkerque (Fig. 2a) showed various phases of trends over several years: stable oscillation from 1974 to 1980, followed by a slight rise from 1980 to 1984. The recent drop in temperature, characterised by cold winters from 1985 to 1987 and the rise in temperatures after 1987, are meteorological events which correspond to periodicities highlighted in studies of data from the last 50 yr (1949 to 1992), and particularly to 7-8 yr cycles (Fromentin & Ibanez 1994). Atmospheric pressure (Fig. 2b) showed a very smooth profile from 1974 to 1986 and an overall decreasing trend, except for the anomaly of the 1976 anticyclonic period. Following a period of relative stability from 1986 to 1989, an increasing trend was seen up to 1992. Study of the last 50 yr revealed a significant 32 mo cycle. Precipitation (Fig. 2c) also presented a major cycle lasting 64 mo. A 32 mo harmonic cycle is detected if all observations over the last half century are taken into account. Anomalies can be seen for the period 1974 to 1992, with high winter rainfall (1977, 1982, 1988) and summer drought (1976 and at the end of the series), to which a long-term, increasing trend observed between 1978 and 1984 should be added.
Mar Ecol Prog Ser 129: 269-281, 1995
272
74
76
78
80
82
84
86
88
90
92
Year ..... .........-....................... ....... . ...
Year
74
76
78
80
82
Year
84
86
88
90
92
Fig. 2. Chronological series of meteorological parameters: (a) air temperature; (b) atmospheric pressure, (c)precipitations; (d)number of days with strong western wind. Dotted line: raw data; solid line data after cumulative sums treatment
Le Fevre-Lehoerff et al.. Hydroclimatic relationsh~pswith planktonic time series
East winds were estimated as no. of days mo-' with winds over 7 m ss1 (following recommendations set out by Lagadeuc 1992). Two periods of instability appeared in 1979 and 1986, when a decreasing trend was observed. A characteristic anticyclonic condition period was seen in 1988-1992. West winds (Fig. 2d) from the sea are generally stronger than east winds; they tended to increase until 1986, followed by a decrease. For information purposes, correlations between parameters were calculated from the cumulative sums. These correlations were all highly significant (with the reservation that observations were not independently performed). Overall, the results were largely coherent. High temperatures (T,,,, in "C) and anticyclonic conditions coincided with high atmospheric pressure (atm P) (Tai,/atm P, r = 0.911"), high precipitation and low atmospheric pressure are correlated (r = -0.595") and correlated with west winds (atm P/west winds, r = -0.845").
Conservative hydrological parameters Water temperature (T,,,,,,) (Fig. 3a) during the yearly cycle fluctuated between 5OC in January and 18.8"C in August (mean monthly temperatures calculated for the period 1975 to 1991). The general annual mean temperature was 11.27"C (Le Fevre-Lehoerff & Woehrling 1991). The cumulative sums of water temperature trace a U-shaped curve; this movement reflects the influence of climatic processes in Dunkerque, especially those of air temperature and atmospheric pressure. Thus calculated correlations for , (r = 0.941 ") and atm P/T , (r = 0.876") T,,,/T were noted. From 1976 onward, the series showed a decreasing trend, reversing to a general warming trend from 1986 onward. Fig. 3b presents the salinity (S) variations. Mean S ranged from 33.63 PSU yearly to 32.96 PSU for March and 34.08 PSU for September. Although seasonal variations occasionally show correlations between temperature and sal.inity, long-term series are more difficult to understand, as the correlation coefficient is calculated on the cumulative sums (T/S = 0.4"). High and low S anomalies and general variations result from the combined action of general circulation phenomena, evaporation, precipitation, thermal flux, wind force and direction (Dickson 1971, Taylor et al. 1981, Maddock & Pingree 1982, Taylor & Stevens 1983). The curve of cun~ulativesums reflected high values at the beginning and end of the series. This tallies with results from other authors who noted high S values in the Channel and the North Sea in recent years (Becker et al. 1992). The causal factors in this research, as in the studies
,,,,
273
previously mentioned, were not yet demonstrated, nor was a relation with wind direction. The only suggested relation is that of low S In water coming from all the coastal and northern European rivers along with north winds which could push these waters against the coastline (Salomon & Breton 1993). The only coherent results observed were the coincidence, in the middle of the period studied, of an increasing trend for precipitation and a decreasing trend for salinity.
Non-conservative chemical parameters At Gravelines, very high values were seen for nutrient salt concentrations; over 20 mm01 1-' of NO3- and over 10 mm01 1-' of SiOH- for a winter monthly mean calculated for the period 1975 to 1991 (Le FevreLehoerff et al. 1993a). On cumulative sum curves, yearly cycles can be recognised, conventionally with maxima in winter. Over the entire period studied, there were 2 notable peaks in 1979 and 1983, 4 yr apart. This was the case for all nutrients (Fig. 3d, e , f). This phenomenon may be related to circulation in the southern North Sea, especially the coastal counter-currents coming from the large northern European rivers, since large declines in salinity and low winter temperatures were also observed during these periods. Correlations were calculated between nutrients and water temperature: N03/T,,,,,,, r = -0.819 " . Correlations between different nutrients were strong: N03/SiOH, r = 0.552".
Overall parameters Chlorophyll values were high, reaching 10 to 40 mg 1-' in spring, depending on the year. In Fig. 3c, seasonal variations are easily recognised. For the period 1976 to 1989, the cumulative sum curve followed, overall, the changes in water temperature. Phaeopigments developed quite similarly. Correlations were calculated on the cumulative values chl a/T,,,,,,,, r = 0.806" and chl a/phaeopigments, r = 0.798". Although the yearly cycle showed increased values for spring, when the water was still cool, long-term changes appeared to indicate another phenomenon. The positive combined action on chlorophyll in the warmest years, anticyclonic periods, and good periods of sunshine were certainly favourable factors (at the end of the period, after 1985). SM followed the same trend as seen for nitrates, with 2 marked reversals, one in 1980 and the other in 1983, giving maximum values on the cumulative sum curves. Over the entire period studied, the increases mentioned were most likely linked to the same phenome-
78
80
82
Year
84
86
FICJ 3 C h r o n o l o g ~ c a lserlcs of h y d r o b ~ o l o g ~ c d l parameters: ( a ) water t e r n p e ~ a l u r e (; b ) s a l ~ n ~ t y ; ( C ) chlorophyll a ; (cl) n~tr-atesNO;, ( e ) phos: phates P O , ; ( f ) s l l ~ c a t e sS I O H Dotted I ~ n erdw d a t a ; s o l ~ dI ~ n e ,d a t a alter cumulat.~vesunls treatment
76
76
88
78
90
80
92
Yeas
82
I
84
76
'
l I
l
78
86
.
.
I
I
76
.
I
.
I
80
'
1
I
'
1
78
.
I
1
.
76
'
82
1
~
#
'
I
84
'
1
1
80
78
'
'
Year
.
I
1
I
.
'
I
I
'
86
I
'
80
'
Year
82
.
I
I
I
I
I
84
'
'
'
I
82
88
I
Year
'
'
'
I
I
'
~
'
I
I
86
l
.
1
84
.
'
90
l
l
'
'
.
I
1
88
l
.
'
1
L
L
86
.
92
~
Le Fevre-Lehoerff et al.. Hydrocllmatic relationships with planktonic time series
non of southern North Sea water circulation already mentioned for nutrient salt variations.
Biological descriptors
.,;
,+,...........
Results of PCA PCA was done using 14 biological descriptors, with hydroclimatic parameters introduced as complementary variables (Fig. 4). All variables were studied as cumulative sums. By plotting ordinates on the 3 mean axes, it was possible to determine 8 groups of biological descriptors, by classification. In order to interpret changes over time, only 1 key category from each group is described here. These were chosen on the basis of the strongest link with a hydrological parameter. All changes over time for the other descriptors of the group were very similar to those of the key category. The following 8 groups were used: H 1: Oikopleura dioica, Temora longicornis; H2: Pseudocalanus minutus, Oncaea sp.; H3: cirriped cypris larvae, Centropages
.
A RA
i,W
",
NO3 O
'"-'~u'
WESTERN WIND
CHLOO 8
PHAE 0
/ ./"".'.. '
' SAG1 0
.:'
,/ j i'.O........,. lK0 ;
;
.......... !,,,./."'..'-*--..
,
AIR To
i
PLEU
. 'j
{ACAR '.,,,
-....
AXIS l
hama tus, Lanice conchilega; H4: Euterpina acutifrons; H5: cirriped nauplii larvae; H6: Bryozoa cyphonaut larvae, Pleuro-
brachia pileus, Acartia clausi; H?: Paracalanus parvus; H8: Sagitta setosa. Axis 1 very clearly indicates the opposition between atmospheric pressure, air and water temperatures a n d chlorophyll (positive correlations) and precipitation, along with west winds, dominant winds and nitrates (negative correlations). Therefore, it represents the opposition between periods of stability and instability. The greatest variability appears to be d u e to meteorological conditions. The organisms which stability most favoured were Temora longjcornis and Oikopleura dioica ( H I ) whereas those most favoured by instability were Pseudocalanus minutus, Oncaea sp. (H2) and Paracalan us parvus (H?). Axis 2 demonstrates the influence of typical hydrological characteristics on population changes over time. High values of salinity correlate positively on Axis 2, as opposed to precipitation. The cirriped
f /'B[2/
AIR To
1
D.WEIGHT,
EA TERN 4
i ?
AXIS 2
Fig. 4 . Principal component analysis results conducted with meteorological, hydrological parameters a n d 14 biological descriptors. Codes as in Table 1
1 c. Cirriped cyuris larvae
I
74
'
~
.
76
I
.
I
.
78
I
.
I
.
80
I
.
I
I
82
.
I
.
Year
.
1 f. Bryozoa cyplionaut larvae
l
I
84
.
I
.
I
86
.
I
.
I
88
'
I
'
I
'
90
I
.
I
.
92
L
Le Fevre-Lehoerff e t al.. Hydroclimatic relationships with plankton~ctime series
277
cypris larvae, Lanice conchilega a n d Centropages hamatus (H3) were strongly linked to episodic high temperatures and salinity, with east winds, whereas the most abundant species in low salinity water situations were those from groups H2 a n d H6: Pseudocalanus mjnutus. Oncaea sp., and Bryozoa cyphonaut larvae, Pleurobrachia pileus, and Acartia clausi. Axis 3 shows how great an effect biological variables have and highlights the opposition between various groups. For instance, the opposition between cirriped nauplii (H5) and chaetognaths Sagitta setosa (H8) is clearly seen.
Changes in biological descriptors over time
Cumulat~vedeviations ,
2
m
.
s
-
,
n
o
Group H I : Oikopleura dioica (Fig. 5a) is one of the coastal appendicularians, a food source for flat fish on sandy bottoms (Wyatt 1973). They can reach sizes which represent 20% of zooplankton dry weight. Thanks to its 'house' system, this animal filters phytoplankton and fine organic particles. At Gravelines, a high correlation was seen (correlations calculated on a yearly cycle) with high temperatures: r = 0.800" and large amounts of phytoplankton (Toularastel et al. 1977). The seasonal variation can be easily identified at the start of the series. O n the curve of cumulative sums 3 clearly recognizable periods can be distinguished. An overall view of the 7 yr cycle is also provided. There are very sharp positive slopes at the end of 1979, in early 1983 and in spring 1986. These peaks follow nitrate values, with a 3 to 5 mo lag. O n the cumulative sum curves, a highly significant correlation is seen with chlorophyll ( r = 0.871 " ) with temperature ( r = 0.763") and a negative correlation with respect to nitrate (r = -0.646"). The copepod Temora longicornis can reach a population level of 104 ind. and make up 80% of zooplankton abundance in spring. The seasonal variation is easily seen, as for Oikopleura diojca. The correlations calculated on cumulative sums a r e strongly linked to food (phytoplankton, represented schematically by chlorophyll). Therefore, the H1 group represented by Oikopleura dioica a n d Temora longicornis was distinguished by species favoured by chlorophyll-rich waters, according to the herbivore food index. The 7 yr cycle was easily recognised for this group a n d correlations between Temora a n d Ojkopleura were strong, r = 0.638". Group H2: For Pseudocalanus minutus (Fig. 5b) the cumulative sums showed very clear trends: from 1974 to 1977, a decreasing trend; between 1977 and 1988, marked population growth, especially in 1979 a n d 1986. From 1988 on, a decreasing trend was observed. The parallel with temperature changes is obvious. The
278
Mar Ecol Prog Ser
trend showed a continuous growth gradient over 11 yr from 1977 to 1988. By presenting long time senes of P. minutus and other species from 1948 to 1980, Colebrook highlighted the possibility of detecting several combined cycles at the same time: 3 to 4 yr, 5 to 6 yr and 11 yr cycles (Colebrook 1986) and the persistence of a species stock from year to year (Colebrook 1981, Colebrook & Taylor 1984).This recurring phenomenon can be seen in planktonic populations as well as in benthic ones. For the Oncaea species in the same group, the relation to temperature was not as strong; however, it was very strong with low salinity. Thus the H2 group, represented by both P. minutus and Oncaea sp. copepods, featured a n affinity with low salinity, low temperature and could easily withstand a low chlorophyll content in the water. An 11 yr cycle with 3 to 4 yr fluctuations was evident. The 2 copepods which represent this group are considered to be omnivorous. Their capability of adapting to various food forms, ranging from large phytoplankton cells to small flagellates (Smetacek 1985) or even living or non-living matter in variable proportions depending on the season (Poulet 1973, 1974, 1976), enhances their dominance over other species for living in winter periods in cold water, sometimes with low pnmary production in the North Sea and the North A t l a ~ t i c . Group H3: The cirriped cypris (Fig. 5c) are not naturally abundant on sandy coasts, but their larvae can colonise man-made structures such as seawalls and industrial or harbour installations. At the Gravelines site, chlorine is injected into the power plant's pumped sea water cooling circults to limit proliferation of cirriped populations (and other fouling groups such as the bryozoa, molluscs a n d hydrozoa). Therefore, changes in larval abundance at the nauplii and cypris stages were closely monitored. Cypris larvae (Fig. 5c) showed very high seasonal variability owing to their life span; however, the year-to-year changes were far from random. The evolution precisely followed that of salinity, with a very high correlation (r = 0.751"). Salinity may well be a limiting factor for this larval stage, or at the least, its influence may be greater than that of the environment's nutrient resources. The correlations which were calculated between taxa in the same group were r = 0.622" for cypris/Centropages and r = 0.676" for cypris/Lanice. Lanice are sand-dwelling annelids generally abundant in low salinity waters, and large numbers of them are carned to Gravelines by the north European rivers (J.M. Dewarumez pers. comm.).Thus, in the first analysis, it may seem curious that the curve of cumulative sums revealed a positive relation with salinity. The long-term phenomenon must take longterm changes in salinity into account over a wider scale as well as the effect on benthic populations from either salinity (Dickson et al. 1988) or temperature anomalies,
particularly in winter, which can have a recurrent effect on benthic populations (Beukema 1989). For larvae groups, fluctuations also depend on all conditions encountered during their long pelagic life (Marcaco 1990) and buoyancy characteristics of aulophore larva (Bhaud & Cazaux 1990). Summer periods are favourable for keeping larvae on the surface, both for cypris (Le Fevre & Bourget 1991) and for Lanice, and the combined effect of hlgh temperature and high salinity cannot be distinguished (r = 0.541 "). The Centropages hamatus copepod was favoured by increases in temperature, sallnity and chlorophyll (CentropageslT ,,, r = 0.509 ). Group H4 (Fig. 5d):The harpacticoid copepod Euterpina acutifrons produces many generations per year in its seasonal cycle and its development time is linked to short-term temperature fluctuations (Bernard 1963, Haq 1972). At Gravelines, maximal abundance in the seasonal cycle is seen in late summer. The geographical distribution of E. acutifronsis vast, ranging from the North Sea to the Mediterranean and the Black Sea. It is an omnivorous, even carnivorous feeder. On the cumulative sum curves (Fig. 5d) 3 rather distinct periods can be seen, with a reverse in trends in 1983. The calculated correlations give positive values with nitrates and salinity and negative values with chlorophyll (r = -0.545 ' ) . The ecological characteristics mentioned above agree with these results. It can be observed that this carnivore's reaction was opposite to that of the herbivore O~kopleuradioica, and slope changes appeared in 1979, 1983 and 1986, as for nitrates. ( E ~ f e r p i n a l N O ~ , r = 0.652"). It could be interesting to compare the H4 group represented by Euterpina acutifrons in a subsequent study. Seeing its wide geographical distribution, this could be one of the species enabling long-term changes in general climatic factors to be compared over a large geographical area. Group H5: Cirripeda nauplii. Considering the extreme variety in the category (Fig. 5e) which is comprised of many very small larvae in early cirriped stages and distributed in swarms, the general population curve is quite similar to noise. Only the seasonal alternation due to the disappearance of this young stage can be recognised. The same was not true for the cumulative sum curve where a 7 yr cycle could clearly be distinguished. This periodicity was already revealed in the year-to-year evolution of meteorological parameters represented by air temperature. The relation with T,,, was negative (r = 0 . 3 0 6 " ) . Relations to salinity (r = 0.468") and suspended matter (r = 0.591 " ) were calculated. Group H6: Bryozoa cyphonauts (Fig. 5f). For this group, which also lncludes the ctenarian Pleurobrachia and the copepod Acartia clausi, neither a clear cycle nor a period could be determined on the senes. High "
Le Fevre-Lehoerff et al.: Hydroclimatic relationships with planktonic time serles
growth periods were noted in 1983 and 1987, along with a relationship with low salinity, although correlations were rather low. Group H?: Paracalanus parvus (Fig. 5g) is a copepod with a wide geographical distribution, abundant in the nel-itic plankton of the Ivory Coast (Binet 1977). Although this species was never dominant at Gravelines, in autumn it is regularly abundant, from September until December. In Fig. 5g the marked seasonal variation is evident, but neither a regular year-to-year trend nor segments of homogeneous periods can be seen. The species is omnivorous and even detritivorous. It seems logical to highlight a negative relation between the cumulative deviations of abundance and chlorophyll (r = -0.566") which this species has more or less in common with Euterpina. There was a significant link with nitrates (r = 0.238 "). Group H8: Sagitta setosa (Fig. 5h) showed a clearly recognizable yearly cycle. The geographical distribution of this chaetognath is more coastal (Oresland 1983) than that of Sagitta elegans, which is also considered to be a good descriptor of long-term climatic fluctuations, from the research conducted by Russell (1935, 1936, 1939) and Southward (1962, 1980, 1984). This coastal species is widely distributed: it is found in the North Sea and the Channel, as well as in the Mediterranean Sea (Dallot 1978) and the Black Sea (Furnestin 1958). As for Euterpina acutifrons, mentioned previously, the advantage of choosing this species can thus be emphasized for long-term monitoring to compare sites on a large geographical scale. Sagitta setosa is carnivorous, and its predation on the copepod stock, especially Paracalanusparvus, was considerable at Gravelines. The year-to-year evolution (Fig. 5h) showed a 7 yr cycle. This cycle reflected that highlighted for meteorological parameters such as atmospheric temperature (Fromentin & Ibanez 1994). The correlation calculated on cumulative sums reached r = 0.551 and the Sagitta/T,",, = 0.425". For the period 1974 to 1982 in the cycle under study, a reverse trend appeared with salinity, but this was no longer a constant relation at the end of the period. On a year-toyear scale, opposition between S. setosa and cirriped nauplii could be noted. On the cumulative sums, the correlation cirriped nauplii/Sagitta was r = -0.609". "
DISCUSSION AND CONCLUSION
Numerical processing of cumulative sums applied to Gravelines plankton time records and of hydroclimatic parameters in the region has enabled cycles to be detected. These cycles are either common to several parameters or biological descriptors or are more specific. The methods used allow data gathered irregu-
larly over time to be processed. They provide a useful description of changes which have occurred and reveal the phenomena which often do not stand out on seasonal variation curves. They also report the intervals where col-relations between biological and physical processes are very marked. These methods render a plotting approach, such as that for PCA, more heuristic by reinforcing the variance percentage explained by the first principal components (high correlations at the start). Seven year cycles appeared clearly on the series studied from 1975 to 1992, both in meteorological parameters such as air temperature and in biological species whose life spans and trophic levels can be very different. These include cirriped larvae, the copepod Temora longicornis, the appendicularian Oikopleura dioica and the chaetognath Sagitta setosa. Therefore, not only did the 7 yr cycle appear for both groups, but a phase opposition between cirriped nauplii and chaetognaths also appeared. This strong relation is of interest; however, it should be noted that the cirriped nauplii group is not homogeneous, but rather made up of larval stages of different species. However, the spring species Balanus balanoides did predominate. Inter-site comparisons remain absolutely indispensable, for this group as for others. In fact most research carried out in the Western Channel, mainly from the Plymouth laboratory, shows that 10 to 11 yr cycles are very important in the ecosystem and are linked to climatic fluctuations (sunspot variation link). However, shorter harmonic cycles also exist (Southward et al. 1975, Colebrook 1982) and this makes it interesting to use the same species over a longer number of years than the period available at Gravelines to compare long-term variations in order to distinguish local phenomena from overall phenomena on a larger spatial scale. Atmospheric pressure and precipitation appeared in 32 mo cycles over the past 50 yr (Fromentin & Ibanez 1994). Studies of strong winds showed the distinction between anticyclonic periods (characterised by high air temperatures) and cyclic conditions combining west winds and precipitation. The 11 yr cycle shown by the Plymouth data was only seen at Gravelines in the Pseudocalanus rnlnutus group and Oncaea sp. The Gravelines time series was too short to conclude an 11 yr cycle periodicity. For some years the abrupt changes in parameter values were noted, especially 1979, 1983 and 1986. The instability during 1979 and 1986 was great. During those years, changes occurred for 1 or several parameters. In 1986, an important change took place, and since 1987 a long period of anticyclonic conditions and a rise in temperatures were observed. Conservative and nonconservative parameters and biological descriptors do not change at the same speed.
Mar Ecol Prog Ser 129: 269-281, 1995
Several cycles can be superimposed (generally 11 yr a n d 3 to 4 yr). Other authors (Colebrook & Taylor 1984, Aebischer et al. 1990) have already observed this in t h e North Atlantic a n d the North Sea. Their research, based on a 30 yr time series, would indicate that several cycles lasting 10 to 12, 5 to 6 and 3 to 4 yr can be superimposed. Acknowledgements. The authors acknowledge the assistance received from all involved in sampling surveys, laboratory studies and data analysis since 1975 This work was supported by a n 'Electricite d e France' contract for the monitoring of French coastal power plants. This work is also part of the PNOC (Programme National d'oceanographie C6tiere) research programme's time series topic. We thank Dr C , Parrish and Mrs J . Heard Carnot for comments and help in presenting the manuscript.
LITERATURE CITED Aebischer NJ, Coulson J C , Colebrook J M (1990) Parallel long-term trends across four marine trophic levels and weather Nature 347-753-755 Aminot A, Chaussepied M (1983) Manuel des analyses chimiques e n milieu marin. Centre National pour 1'Exploitation des Oceans. CNEXO, BNDO DOC, Brest Becker G, Corten AD, Dooley H (1992) Recent high salinity in the English Channel and southern North Sea. Comm Meet int Coun Explor Sea (Hydrography Committee) CM-ICES 1992/C6 Belsher T, Dulesmont R, Degros N, Dewarumez JM, Grossel H, Le Fevre-Lehoerff G (1986) Recueil des techniques utilisees pour I'etude ecologique des sites de centrales nucleates, sur les cotes francaises de Manche el Atlantique. DERO-8624-EL, Rapp IFREMER, Brest Bernard M (1963) Le cycle vital e n laboratoire d'un copepode pelagique d e Mediterranee: Euterpina acutifrons Claus Pelagos. Bull Inst Oceanogr Alger 1 35-48 Beukema J J (1989) Long-term changes in macrozoobenthic abundance on the tidal flats of the western part of the Dutch Wadden Sea. Helgolander Meeresunters 43: 405-415 Bhaud MR, Cazaux CP (1990) Buoyancy characteristics of Lanice conchilega (Pallas) larvae (Terebellidae). Implications for settlement. J exp mar Biol Ecol 141:31-45 Binet D (977) Contributions a la connaissance d u zooplancton neritique ivoirien. Ecologie descnptive et dynamique These Doct es Sciences Naturelles, Univ P et M Curie, Paris VI Chardy P, Menesguen A (1984) Ecological monitoring and assessment in a coastal power plant impact study: alternative strategy. Thalassia jugosl 20.1 15-126 Colebrook JM (1981.) Continuous plankton records: persistence in time-series of annual means of abundance of zooplankton. Mar Biol 61:143-149 Colebrook J M (1982) Continuous plankton records: seasonal variations in the distribution and abundance of plankton in the North Atlantic Ocean and the North Sea. J Plankton Res 4:435-462 Colebrook J M (1986) Environmental influences on long-term variability in marine plankton. Hydrobiologia 142: 309-325
Colebrook JM, Taylor AH (1984) Significant time scales of long term variability in the plankton and the environment Rapp Pv Reun Cons int Explor Mer 183:20-26 Cushing DH (1961) O n the failure of the Plymouth herring fishery. J mar biol Ass UK 41:799-816 d i s h i n g DH (1978) Biological effects of climate change. Rapp Pv Reun Cons int Explor Mer 173:107-116 Cushing DH, Dickson RR (1976) The biological response in the sea to climatic changes. Adv mar Biol 14:l-122 Dallot S (1978) Sur la presence du chaetognathe planctonique Sagitta setosa Miiller 1897, dans les eaux neritiques d e Castellon, Inv Pesq 42:33-52 Dauvin JC, Ibanez F (19861 Variations a long terme (1977-1985) du peuplement des sables fins d e la Pierre Noire (Baie d e Morlaix, Manche occidentale): analyse statistique d e l'evolution structurale. Hydrobiologia 142: 171-186 Dickson RR (1971) A recurrent and persistent pressure-anomaly pattern a s the principal cause of intermediate-scale hydrographic variation in the European shelf seas. Dtsch hydrogr Z 24:97-119 Dickson R R , Meincke J , Malmberg SVA, Lee AJ (1988) The great salinity anomaly in the northern North Atlantic 1968-1982. Prog Oceanogr 20:103-151 Fromentin JM, Ibanez F (1994) Year to year changes in meteorological features of the French coast during the last half century. Examples of two biological responses. Oceanol Acta 17:285-296 Fromentin JM, Ibanez F, Legendre P (1993) A phytosociological method for interpreting plankton data. Mar Ecol Prog Ser 93-285-306 Furnestin ML (1958) Chaetognathes recoltes en Mediterranee orientals et en Mer Noire par la 'Calypso' (campagne 1955).Rapp Pv Reun Comm int Explor scient Mer mediterr 14:202-209 Haq SM (1972) Breeding of Euterpina aciitifrons, a harpacticold copepod with special reference to dimorphic males Mar Biol 15:221-235 Ibanez F (1984) Sur la segmentation des series chronologiques planctoniques multivariables. Oceanol Acta 7(4): 481-491 Ibanez F, Etienne M (1992) Le filtrage des series chronologiques par l'analyse en composantes principales d e processus (ACPP).J Rech Oceanogr 16:66-72 Ibanez F, Fromentin JM, Castel J (1993) Application d e la methode des sommes cumulees a l'analyse des series chronologiques en oceanographie Cr Acad Sci Paris (Sci Vie) 316:745-748 Lagadeuc Y (1992) Transport larvaire e n Manche. Exemple de Pectinaria korenj (Malmgren) annehde polychete en bale de Seine. Oceanol Acta 15(4):383-395 Lamb HH (1969) The new look of climatology Nature 223: 1209-1215 Lamb HH (1985) Climatic history and the future. Princeton University Press, Princeton, NJ Le Fevre J , Bourget E (1991) Neustonic niche for cirnpede larvae as a possible adaptation to long-range dispersal. Mar Ecol Prog Ser 74:185-194 Le Fevre-Lehoerff G (1985) Techniques for zooplankton studies at nuclear power station sites on the English Channel and Atlantic coasts of France DERO-8712-EL, IFREMER, Brest, p 1-24 Le Fevre-Lehoerff G, Erard-Le Denn E, Arzul G (1993a) Planktonic ecosystems in the Channel. Trophic relations. Oceanol Acta 1993 16661-670 Le Fevre-Lehoerff G , Grossel H,Derrien A (1983) Evolution des populations planctoniques animales en Manche et au
Le Fevre-Lehoerff et al.: Hydroclimatic 1-elationships with planktonic tlme serles
281
-
sud d e la Mer du Nord. Proc 17th Eur Mar Biol Symp, Brest, Sept-Oct 1982. Oceanol Actcl v01 spec:131-135 Le Fevre-Lehoerff G, Grossel H, \iVoehrling D , Dewarumez JM, Delesmont R (1993b) La centrale nucleaire littorale de Gravelines: ecologie, environnement et impact. J Rech Oceanogr 18.33-38 Le Fevre-Lehoerff G. Quintin JY 11981)Etucle c o m ~ a r a t i v ede la sensibiliti. d e differentes ispeces de copipodes aux variations de la t e m ~ e r a t u r ee n Manche. Relations entre la taille des individus et les facteurs du milieu. Influence des rejets thermiques sur le m ~ l i e uvivant e n mer et e n estuaire. In- 2emes JournCes de la Thermoecologie, Nantes, 14-15 nov 1979. EDF Equlpement. Actes du Colloque 1.71-86 Le Fevre-Lehoerff G, Toularastel F, Arnal O (1976) Zooplancton. In: Etude ecologique, site de Gravelmes. Rapport final, Vol 1. Rapport interne des CNEXO pour EDF. p 105-194 Le Fevre-Lehoerff G, Woehrling D (1991) Amenagement industriel du littoral et surveillance de l'environnement: la centrale nucleaire d e Gravelines (1975 a 1989). Actes du Colloque International sur I'environnement des mers epicont~nentales,Lille, 20-22 mars 1990 Oceanol Acta 11 299-311 Maddock L, Pingree RD (1982) Mean heat and salts budgets for the eastern English Channel and southern bight of the North Sea. J mar biol Ass UK 62:559-575 Marcaco G (1990) Study of competence and recruitment in Lanice conchilega larvae (Polychaeta Terebellidae). Oceanis 16:191-206 Bresland V (1983) Abundance breeding and temporal size distrlbutlon of the chaetognath Sagitta selosa in the Kattegat. J Plankton Res 5:425-443 Poulet S (1973) Grazing of Pseudocalanus minutus on naturally occurring particulate matter Llmnol Oceanogr 18. 564-573 Poulet S (1974) Seasonal grazing of Pseudocalanus rninutus on particles. Mar Biol 25:109-123 Poulet S (1976) Feeding of Pseudocalanus rninutus on living and non living particles. Mar Biol 34.1 17-125 Russell FS (1935) On the value of certain plankton animals as indicators of water movements In the English Channel and North Sea. J mar biol Ass U K 20:309-332 Russell FS (1936) The importance of certain plankton anlmals as ind~catorsof water movements in the western end of the English Channel. Rapp Pv Cons Perm int Explor Mer 100: 7-10
Russell FS (1939) Hydrographical and biological conditions in the North Sea a s indicated by plankton organisms. J Cons int Explor Mer 14.171-192 Russcll FS, Southward AJ, Boalch GT, Butler E1 (1971) Changes in b~ologicalconditions in the Engl~shChannel off Plymouth durlng the last half century. Nature 234: 468-470 Salomon JC, Breton M (1993)An atlas of long term currents in the Channel. Oceanol Acta 1993 16:439-448 Smetacek V (1985) Thc annual cycle of Kiel Bight plankton: a long term analysis. lstuaries 8:145-157 Southward AJ (1962) The distribution of some plankton animals in the English Channel and approaches 11. Surveys with the Gulf. 111 High speed sample 1958-1960. J mar biol Ass UK 42:275-375 Southward AJ (1963) The distribution of some plankton animals in the English channel and approaches. 111. Theories about long term biological changes including fish. J mar biol Ass UK 43: 1-29 Southward AJ (1980) The western English Channel - a n inconstant ecosystem? Nature 285:361-366 Southward AJ (1984) Fluctuations in the 'indicator' chaetognaths Sagjtta elegans and Sagitta setosa in the western channel. Oceanol Acta 7:229-239 Southward AJ, Butler El, Pennycuick L (1975) Recent cyclic changes in climate and in abundance of marine life. Nature 253:714-717 Taylor AH, Reid PC, Marsh TJ, Jonas TD, Stephens JA (1981) Year to year changes in the salinity of the eastern English channel, 1948-1973: a budget. J mar biol Ass UK 61: 489-507 Taylor AH, Stephens JA (1983) Seasonal and year to year changes in the temperatures of the English Channel and the southern North Sea, 1961-1976- a budget. Oceanol Acta 653-72 Toularastel F, Le Fevre-Lehoerff G, Arnal O (1977) L'interet du zooplancton dans l'elaboration d'un etat d e reference. Influence des rejets thermiques sur le milieu vivant e n mer et e n estuaire. In: leres Journees d e la Thermoecologie, Brest. EDF Equipement. Actes du Colloque 1. 461-480 Woehrling D, Le F&vre-Lehoerff G (1993) Fluctuations pluriannuelles de la ponte de la sole Solea solea (L.) sur le littoral du sud de la Mer du Nord. J Rech Oceanogr 18:74-79 Wyatt T (1973) The biology of Oikopleura dioica and Fritil/aria borealis in the southern bight. Mar B101 22:137-158
This article was submitted to the editor
Manuscript first received: December 19, 1994 Revised version accepted: J u n e 26, 1995