Interactions spatiales et auto-organisation des végétations semi-arides

Les recherches récapitulées dans cette thèse de doctorat ont porté sur les causes de l’organisation spatiale des végétations périodiques. Ces structures paysagères aux motifs réguliers, tachetés, tigrés ou labyrinthiques, d’échelle décamétrique à hectométrique, couvrant des étendues considérables sur au moins trois continents, constituent un cas d’école dans l’étude des processus endogènes présidant à l’hétérogénéité du couvert végétal. Ces structures prennent place sur un substrat homogène, mis à part la rétroaction du couvert lui-même, et sont marquées par des écotones abrupts et la persistance d’une proportion considérable de sol nu. Plusieurs modèles ont mis en avant l’existence possible d’un phénomène d’auto-organisation du couvert, qui verrait une structure d’ensemble émerger des interactions locales entre individus. Ces modèles se basent sur le jeu simultané de la consommation de la ressource (compétition) et de l’amélioration de l’un ou l’autre des éléments du bilan de la même ressource par le couvert (facilitation). La condition à l’existence d’une structure d’ensemble spatialement périodique et stable réside dans une différence entre la portée de la compétition (plus grande) et celle de la facilitation. L’apparition de ces structures est modulée par le taux de croissance biologique, qui est le reflet des contraintes extérieures telles que l’aridité, le pâturage ou la coupe de bois. Le modus operandi des interactions spatiales supposées entre individus reste largement à préciser.

Nos recherches ont été menées au sud-ouest de la République du Niger, à l’intérieur et dans les environs du parc Régional du W. Trois axes ont été explorés :(i) Une étude de la dépendance spatiale entre la structure de la végétation (biovolumes cartographiés) et les paramètres du milieu abiotique (relief, sol), sur base d’analyses spectrales et cross-spectrales par transformée de Fourier (1D et 2D). (ii) Une étude diachronique (1956, 1975 et 1996) à large échelle (3000 km²) de l’influence de l’aridité et des pressions d’origine anthropique sur l’auto-organisation des végétations périodiques, basée sur la caractérisation de la structure spatiale des paysages sur photos aériennes via la transformée de Fourier en 2D. (iii) Trois études portant sur les interactions spatiales entre individus :En premier lieu, via l’excavation des systèmes racinaires (air pulsé) ;Ensuite, par un suivi spatio-temporel du bilan hydrique du sol (blocs de gypse) ;Enfin, via le marquage de la ressource par du deutérium.

Nous avons ainsi pu établir que les végétations périodiques constituent bien un mode d’auto-organisation pouvant survenir sur substrat homogène et modulé par les contraintes climatiques et anthropiques. Un ajustement rapide entre l’organisation des végétations périodiques et le climat a pu être montrée en zone protégée. La superficie et l’organisation des végétations périodiques y ont tour à tour progressé et régressé en fonction d’épisodes secs ou humides. Par contre, en dehors de l’aire protégée, la possibilité d’une restauration du couvert semble fortement liée au taux d’exploitation des ressources végétales. Ces résultats ont d’importantes implications quant à la compréhension des interactions entre climat et écosystèmes et à l’évaluation de leurs capacités de charge. La caractérisation de la structure spatiale des végétations arides, notamment par la transformée de Fourier d’images HR, devrait être généralisée comme outil de monitoring de l’état de ces écosystèmes. Nos études portant sur les modes d’interactions spatiales ont permis de confirmer l’existence d’une facilitation à courte portée du couvert végétal sur la ressource. Cependant, cette facilitation ne semble pas s’exercer sur le terme du bilan hydrique traditionnellement avancé, à savoir l’infiltration, mais plutôt sur le taux d’évaporation (deux fois moindre à l’ombre des canopées). Ce mécanisme exclut l’existence de transferts diffusifs souterrains entre sols nu et fourrés. Des transferts inverses semblent d’ailleurs montrés par le marquage isotopique. L’étude du bilan hydrique et la cartographie du micro-relief, ainsi que la profondeur fortement réduite de la zone d’exploitation racinaire, jettent de sérieux doutes quant au rôle communément admis des transferts d’eau par ruissellement/diffusion de surface en tant que processus clé dans la compétition à distance entre les plantes. L’alternative réside dans l’existence d’une compétition racinaire de portée supérieure aux canopées. Cette hypothèse trouve une confirmation tant par les rhizosphères excavées, superficielles et étendues, que dans le marquage isotopique, montrant des contaminations d’arbustes situés à plus de 15 m de la zone d’apport. De même, l’étude du bilan hydrique met en évidence les influences simultanées et contradictoires (facilitation/compétition) des ligneux sur l’évapotranspiration.

/

This PhD thesis gathers results of a research dealing with the causes of the spatial organisation of periodic vegetations. These landscape structures, featuring regular spotted, labyrinthine or banded patterns of decametric to hectometric scale, and extending over considerable areas on at least three continents, constitute a perfect study case to approach endogenous processes leading to vegetation heterogeneities. These patterns occur over homogeneous substratum, except for vegetation’s own feedbacks, and are marked by sharp ecotones and the persistence of a considerable amount of bare soil. A number of models suggested a possible case of self-organized patterning, in which the general structure would emerge from local interactions between individuals. Those models rest on the interplay of competitive and facilitative effects, relating to soil water consumption and to soil water budget enhancement by vegetation. A general necessary condition for pattern formation to occur is that negative interactions (competition) have a larger range than positive interactions (facilitation). Moreover, all models agree with the idea that patterning occurs when vegetation growth decreases, for instance as a result of reduced water availability, domestic grazing or wood cutting, therefore viewing patterns as a self-organised response to environmental constraints. However the modus operandi of the spatial interactions between individual plants remains largely to be specified.

We carried out a field research in South-West Niger, within and around the W Regional Park. Three research lines were explored: (i) The study of the spatial dependency between the vegetation pattern (mapped biovolumes) and the factors of the abiotic environment (soil, relief), on the basis of spectral and cross-spectral analyses with Fourier transform (1D and 2D). (ii) A broad scale diachronic study (1956, 1975, 1996) of the influence of aridity and human induced pressures on the vegetation self-patterning, based on the characterisation of patterns on high resolution remote sensing data via 2D Fourier transform. (iii) Three different approaches of the spatial interactions between individuals: via root systems excavation with pulsed air; via the monitoring in space and time of the soil water budget (gypsum blocks method); and via water resource labelling with deuterated water.

We could establish that periodic vegetations are indeed the result of a self-organisation process, occurring in homogeneous substratum conditions and modulated by climate and human constraints. A rapid adjustment between vegetation patterning and climate could be observed in protected zones. The area and patterning of the periodic vegetations successively progressed and regressed, following drier or wetter climate conditions. On the other hand, outside protected areas, the restoration ability of vegetation appeared to depend on the degree of vegetation resource exploitation. These results have important implications regarding the study of vegetation-climate interactions and the evaluation of ecosystems’ carrying capacities. Spatial pattern characterisation in arid vegetations using Fourier transform of HR remote sensing data should be generalised for the monitoring of those ecosystems. Our studies dealing with spatial interaction mechanisms confirmed the existence of a short range facilitation of the cover on water resource. However, this facilitation does not seem to act through the commonly accepted infiltration component, but rather on the evaporative rate (twice less within thickets). This mechanism excludes underground diffusive transfers between bare ground and vegetation. Inverse transfers were even shown by deuterium labelling. Water budget study and micro-elevation mapping, along with consistent soil shallowness, together cast serious doubts on the traditional mechanism of run-off/diffusion of surface water as a key process of the long range competition between plants. An alternative explanation lies in long range root competition. This hypothesis find support as well in the excavated root systems, shallow and wide, as in isotopic labelling, showing contaminations of shrubs located up to 15 m of the irrigated area. Water budget study also evidenced simultaneous contradictory effects (facilitation/competition) of shrubs on evapotranspiration.

Case history and persistence of the non-indigenous diatom Coscinodiscus wailesii in the north-east Atlantic

The introduction of non-indigenous marine plankton species can have a considerable ecological and economic effect on regional systems. Their presence, however, can go unnoticed until they reach nuisance status and as a consequence few case histories exist containing information on their initial appearance and their spatio-temporal patterns. Here we report on the occurrence of the non-indigenous diatom Coscinodiscus wailesii in 1977 in the English Channel, its subsequent geographical spread into European shelf seas, and its persistence as a significant member of the diatom community in the north-east Atlantic from 1977-1995.

Relating sardine recruitment in the Northern Benguela to satellite-derived sea surface height using a neural network pattern recognition approach

Processes of enrichment, concentration and retention are thought to be important for the successful recruitment of small pelagic fish in upwelling areas, but are difficult to measure. In this study, a novel approach is used to examine the role of spatio-temporal oceanographic variability on recruitment success of the Northern Benguela sardine Sardinops sagax. This approach applies a neural network pattern recognition technique, called a self-organising map (SOM), to a seven-year time series of satellite-derived sea level data. The Northern Benguela is characterised by quasi-perennial upwelling of cold, nutrient-rich water and is influenced by intrusions of warm, nutrient-poor Angola Current water from the north. In this paper, these processes are categorised in terms of their influence on recruitment success through the key ocean triad mechanisms of enrichment, concentration and retention. Moderate upwelling is seen as favourable for recruitment, whereas strong upwelling, weak upwelling and Angola Current intrusion appear detrimental to recruitment success. The SOM was used to identify characteristic patterns from sea level difference data and these were interpreted with the aid of sea surface temperature data. We found that the major oceanographic processes of upwelling and Angola Current intrusion dominated these patterns, allowing them to be partitioned into those representing recruitment favourable conditions and those representing adverse conditions for recruitment. A marginally significant relationship was found between the index of sardine recruitment and the frequency of recruitment favourable conditions (r super(2) = 0.61, p = 0.068, n = 6). Because larvae are vulnerable to environmental influences for a period of at least 50 days after spawning, the SOM was then used to identify windows of persistent favourable conditions lasting longer than 50 days, termed recruitment favourable periods (RFPs). The occurrence of RFPs was compared with back-calculated spawning dates for each cohort. Finally, a comparison of RFPs with the time of spawning and the index of recruitment showed that in years where there were 50 or more days of favourable conditions following spawning, good recruitment followed (Mann-Whitney U-test: p = 0.064, n = 6). These results show the value of the SOM technique for describing spatio-temporal variability in oceanographic processes. Variability in these processes appears to be an important factor influencing recruitment in the Northern Benguela sardine, although the available data time series is currently too short to be conclusive. Nonetheless, the analysis of satellite data, using a neural network pattern-recognition approach, provides a useful framework for investigating fisheries recruitment problems.

The CPR Survey (1948-1997): a gridded database browser of plankton abundance in the North Sea.

The Continuous Plankton Recorder (CPR) survey provides a unique multi- decadal dataset on the abundance of plankton in the North Sea and North Atlantic and is one of only a few monitoring programmes operating at a large spatio- temporal scale. The results of all samples analysed from the survey since 1946 are stored on an Access Database at the Sir Alister Hardy Foundation for Ocean Science (SAHFOS) in Plymouth. The database is large, containing more than two million records (~80 million data points, if zero results are added) for more than 450 taxonomic entities. An open data policy is operated by SAHFOS. However, the data are not on-line and so access by scientists and others wishing to use the results is not interactive. Requests for data are dealt with by the Database Manager. To facilitate access to the data from the North Sea, which is an area of high research interest, a selected set of data for key phytoplankton and zooplankton species has been processed in a form that makes them readily available on CD for research and other applications. A set of MATLAB tools has been developed to provide an interpolated spatio-temporal description of plankton sampled by the CPR in the North Sea, as well as easy and fast access to users in the form of a browser. Using geostatistical techniques, plankton abundance values have been interpolated on a regular grid covering the North Sea. The grid is established on centres of 1 degree longitude x 0.5 degree latitude (~32 x 30 nautical miles). Based on a monthly temporal resolution over a fifty-year period (1948-1997), 600 distribution maps have been produced for 54 zooplankton species, and 480 distribution maps for 57 phytoplankton species over the shorter period 1958-1997. The gridded database has been developed in a user-friendly form and incorporates, as a package on a CD, a set of options for visualisation and interpretation, including the facility to plot maps for selected species by month, year, groups of months or years, long-term means or as time series and contour plots. This study constitutes the first application of an easily accessed and interactive gridded database of plankton abundance in the North Sea. As a further development the MATLAB browser is being converted to a user- friendly Windows-compatible format (WinCPR) for release on CD and via the Web in 2003.

Summarising spatial and temporal information in CPR data

The Continuous Plankton Recorder survey provides pan-oceanic data on geographic distribution, species composition, seasonal cycles of abundance, and long-term change during the last 70 years. In this paper we compare and contrast some of the historic data-analytic protocols of the survey, focusing primarily on the various means by which spatio-temporal information in CPR data has been exposed. Relative strengths and limitations are assessed, followed by suggestions for future approaches to the visualisation and summarising of CPR data.

The meso-scale response of subarctic North Pacific seabird community structure to lower trophic level abundance and diversity

Understanding the mechanisms that structure communities and influence biodiversity are fundamental goals of ecology. To test the hypothesis that the abundance and diversity of upper-trophic level predators (seabirds) is related to the underlying abundance and diversity of their prey (zooplankton) and ecosystem-wide energy availability (primary production), we initiated a monitoring program in 2002 that jointly and repeatedly surveys seabird and zooplankton populations across a 7,500 km British Columbia-Bering Sea-Japan transect. Seabird distributions were recorded by a single observer (MH) using a strip-width technique, mesozooplankton samples were collected with a Continuous Plankton Recorder, and primary production levels were derived using the appropriate satellite parameters and the Vertically Generalized Production Model (Behrenfeld and Falkowski 1997). Each trophic level showed clear spatio-temporal patterns over the course of the study. The strongest relationship between seabird abundance and diversity and the lower trophic levels was observed in March/April ('spring') and significant relationships were also found through June/July ('summer'). No discernable relationships were observed during the September/October ('fall') months. Overall, mesozooplankton abundance and biomass explained the dominant portion of seabird abundance and diversity indices (richness, Simpson's Index, and evenness), while primary production was only related to seabird richness. These findings underscore the notion that perturbations of ocean productivity and lower trophic level ecosystem constituents influenced by climate change, such as shifts in timing (phenology) and synchronicity (match-mismatch), could impart far-reaching consequences throughout the marine food web.

Marine plankton phenology and life history in a changing climate: current research and future directions

Increasing availability and extent of biological ocean time series (from both in situ and satellite data) have helped reveal significant phenological variability of marine plankton. The extent to which the range of this variability is modified as a result of climate change is of obvious importance. Here we summarize recent research results on phenology of both phytoplankton and zooplankton. We suggest directions to better quantify and monitor future plankton phenology shifts, including (i) examining the main mode of expected future changes (ecological shifts in timing and spatial distribution to accommodate fixed environmental niches vs. evolutionary adaptation of timing controls to maintain fixed biogeography and seasonality), (ii) broader understanding of phenology at the species and community level (e.g. for zooplankton beyond Calanus and for phytoplankton beyond chlorophyll), (iii) improving and diversifying statistical metrics for indexing timing and trophic synchrony and (iv) improved consideration of spatio-temporal scales and the Lagrangian nature of plankton assemblages to separate time from space changes.

Effects of environmental conditions on the seasonal distribution of phytoplankton biomass in the North Sea

We study the spatial and seasonal variability of phytoplankton biomass (as phytoplankton color) in relation to the environmental conditions in the North Sea using data from the Continuous Plankton Recorder survey. By using only environmental fields and location as predictor variables we developed a nonparametric model (generalized additive model) to empirically explore how key environmental factors modulate the spatio-temporal patterns of the seasonal cycle of algal biomass as well as how these relate to the ,1988 North Sea regime shift. Solar radiation, as manifest through changes of sea surface temperature (SST), was a key factor not only in the seasonal cycle but also as a driver of the shift. The pronounced increase in SST and in wind speed after the 1980s resulted in an extension of the season favorable for phytoplankton growth. Nutrients appeared to be unimportant as explanatory variables for the observed spatio-temporal pattern, implying that they were not generally limiting factors. Under the new climatic regime the carrying capacity of the whole system has been increased and the southern North Sea, where the environmental changes have been more pronounced, reached a new maximum.

Seasonal plankton dynamics along a cross-shelf gradient off northern Spain

The development of population models able to reproduce the dynamics of zooplankton is a central issue when trying to understand how a changing environment would affect zooplankton in the future. Using 10 years of monthly data on phytoplankton and zooplankton abundance in the Bay of Biscay from the IEO's RADIALES time-series programme, we built non-parametric Generalized Additive Models (GAMs) able to reproduce the dynamics of plankton on the basis of environmental factors (nutrients, temperature, upwelling and photoperiod). We found that the interaction between these two plankton components is approximately linear, whereas the effects of environmental factors are non-linear. With the inclusion of the environmental variability, the main seasonal and inter-annual dynamic patterns observed within the studied plankton assemblage indicate the prevalence of bottom-up regulatory control. The statistically deduced models were used to simulate the dynamics of the phytoplankton and zooplankton. A good agreement between observations and simulations was obtained, especially for zooplankton. We are presently developing spatio-temporal GAM models for the North Sea based on the Continuous Plankton Recorder database.

Non-linearities, regime shifts and recovery: The recent influence of climate on Black Sea chlorophyll

The Black Sea ecosystem experienced severe eutrophication-related degradation during the 1970s and 1980s. However, in recent years the Black Sea has shown some signs of recovery which are often attributed to a reduction in nutrient loading. Here, SeaWiFS chlorophyll a (chl a), a proxy for phytoplankton biomass, is used to investigate spatio-temporal patterns in Black Sea phytoplankton dynamics and to explore the potential role of climate in the Black Sea's recovery. Maps of chl a anomalies, calculated relative to the 8 year mean, emphasize spatial and temporal variability of phytoplankton biomass in the Black Sea, particularly between the riverine-influenced Northwest Shelf and the open Black Sea. Evolution of phytoplankton biomass has shown significant spatial variability of persistence of optimal bloom conditions between three major regions of the Black Sea. With the exception of 2001, chl a has generally decreased during our 8 year time-series. However, the winter of 2000–2001 was anomalously warm with low wind stress, resulting in reduced vertical mixing of the water column and retention of nutrients in the photic zone. These conditions were associated with anomalously high levels of chl a throughout much of the open Black Sea during the following spring and summer. The unusual climatic conditions occurring in 2001 may have triggered a shift in the Black Sea's chl a regime. The long-term significance of this recent shift is still uncertain but illustrates a non-linear response to climate forcing that makes future ecosystem changes in the pelagic Black Sea ecosystem difficult to predict.

Bottom-up effects of climate on fish populations: data from the Continuous Plankton Recorder

The Continuous Plankton Recorder (CPR) dataset on fish larvae has an extensive spatio-temporal coverage that allows the responses of fish populations to past changes in climate variability, including abrupt changes such as regime shifts, to be investigated. The newly available dataset offers a unique opportunity to investigate long-term changes over decadal scales in the abundance and distribution of fish larvae in relation to physical and biological factors. A principal component analysis (PCA) using 7 biotic and abiotic parameters is applied to investigate the impact of environmental changes in the North Sea on 5 selected taxa of fish larvae during the period 1960 to 2004. The analysis revealed 4 periods of time (1960–1976; 1977–1982; 1983–1996; 1997–2004) reflecting 3 different ecosystem states. The larvae of clupeids, sandeels, dab and gadoids seemed to be affected mainly by changes in the plankton ecosystem, while the larvae of migratory species such as Atlantic mackerel responded more to hydrographic changes. Climate variability seems more likely to influence fish populations through bottom-up control via a cascading effect from changes in the North Atlantic Oscillation (NAO) impacting on the hydro dynamic features of the North Sea, in turn impacting on the plankton available as prey for fish larvae. The responses and adaptability of fish larvae to changing environmental conditions, parti cularly to changes in prey availability, are complex and species-specific. This complexity is enhanced with fishing effects interacting with climate effects and this study supports furthering our under - standing of such interactions before attempting to predict how fish populations respond to climate variability

Identifying four phytoplankton functional types from space: An ecological approach

Deriving maps of phytoplankton taxa based on remote sensing data using bio-optical properties of phytoplankton alone is challenging. A more holistic approach was developed using artificial neural networks, incorporating ecological and geographical knowledge together with ocean color, bio-optical characteristics, and remotely sensed physical parameters. Results show that the combined remote sensing approach could discriminate four major phytoplankton functional types (diatoms, dinoflagellates, coccolithophores, and silicoflagellates) with an accuracy of more than 70%. Models indicate that the most important information for phytoplankton functional type discrimination is spatio-temporal information and sea surface temperature. This approach can supply data for large-scale maps of predicted phytoplankton functional types, and an example is shown.

Scale-dependent foraging ecology of a marine top predator modelled using passive acoustic data

1.Understanding which environmental factors drive foraging preferences is critical for the development of effective management measures, but resource use patterns may emerge from processes that occur at different spatial and temporal scales. Direct observations of foraging are also especially challenging in marine predators, but passive acoustic techniques provide opportunities to study the behaviour of echolocating species over a range of scales. 2.We used an extensive passive acoustic data set to investigate the distribution and temporal dynamics of foraging in bottlenose dolphins using the Moray Firth (Scotland, UK). Echolocation buzzes were identified with a mixture model of detected echolocation inter-click intervals and used as a proxy of foraging activity. A robust modelling approach accounting for autocorrelation in the data was then used to evaluate which environmental factors were associated with the observed dynamics at two different spatial and temporal scales. 3.At a broad scale, foraging varied seasonally and was also affected by seabed slope and shelf-sea fronts. At a finer scale, we identified variation in seasonal use and local interactions with tidal processes. Foraging was best predicted at a daily scale, accounting for site specificity in the shape of the estimated relationships. 4.This study demonstrates how passive acoustic data can be used to understand foraging ecology in echolocating species and provides a robust analytical procedure for describing spatio-temporal patterns. Associations between foraging and environmental characteristics varied according to spatial and temporal scale, highlighting the need for a multi-scale approach. Our results indicate that dolphins respond to coarser scale temporal dynamics, but have a detailed understanding of finer-scale spatial distribution of resources.

Basin-scale phenology and effects of climate variability on global timing of initial seaward migration of Atlantic salmon (Salmo salar)

Migrations between different habitats are key events in the lives of many organisms. Such movements involve annually recurring travel over long distances usually triggered by seasonal changes in the environment. Often, the migration is associated with travel to or from reproduction areas to regions of growth. Young anadromous Atlantic salmon (Salmo salar) emigrate from freshwater nursery areas during spring and early summer to feed and grow in the North Atlantic Ocean. The transition from the freshwater (parr') stage to the migratory stage where they descend streams and enter salt water (smolt') is characterized by morphological, physiological and behavioural changes where the timing of this parr-smolt transition is cued by photoperiod and water temperature. Environmental conditions in the freshwater habitat control the downstream migration and contribute to within- and among-river variation in migratory timing. Moreover, the timing of the freshwater emigration has likely evolved to meet environmental conditions in the ocean as these affect growth and survival of the post-smolts. Using generalized additive mixed-effects modelling, we analysed spatio-temporal variations in the dates of downstream smolt migration in 67 rivers throughout the North Atlantic during the last five decades and found that migrations were earlier in populations in the east than the west. After accounting for this spatial effect, the initiation of the downstream migration among rivers was positively associated with freshwater temperatures, up to about 10 degrees C and levelling off at higher values, and with sea-surface temperatures. Earlier migration occurred when river discharge levels were low but increasing. On average, the initiation of the smolt seaward migration has occurred 2.5days earlier per decade throughout the basin of the North Atlantic. This shift in phenology matches changes in air, river, and ocean temperatures, suggesting that Atlantic salmon emigration is responding to the current global climate changes.

Comparison of new and primary production models using SeaWiFS data in contrasting hydrographic zones of the northern North Atlantic.

The accuracy of two satellite models of marine primary (PP) and new production (NP) were assessed against 14C and 15N uptake measurements taken during six research cruises in the northern North Atlantic. The wavelength resolving model (WRM) was more accurate than the Vertical General Production Model (VGPM) for computation of both PP and NP. Mean monthly satellite maps of PP and NP for both models were generated from 1997 to 2010 using SeaWiFS data for the Irminger basin and North Atlantic. Intra- and inter-annual variability of the two models was compared in six hydrographic zones. Both models exhibited similar spatio-temporal patterns: PP and NP increased from April to June and decreased by August. Higher values were associated with the East Greenland Current (EGC), Iceland Basin (ICB) and the Reykjanes Ridge (RKR) and lower values occurred in the Central Irminger Current (CIC), North Irminger Current (NIC) and Southern Irminger Current (SIC). The annual PP and NP over the SeaWiFS record was 258 and 82 gC m-2 yr-1 respectively for the VGPM and 190 and 41 gC m-2 yr-1 for the WRM. Average annual cumulative sum in the anomalies of NP for the VGPM were positively correlated with the North Atlantic Oscillation (NAO) in the EGC, CIC and SIC and negatively correlated with the multivariate ENSO index (MEI) in the ICB. By contrast, cumulative sum of the anomalies of NP for the WRM were significantly correlated with NAO only in the EGC and CIC. NP from both VGPM and WRM exhibited significant negative correlations with Arctic Oscillation (AO) in all hydrographic zones. The differences in estimates of PP and NP in these hydrographic zones arise principally from the parameterisation of the euphotic depth and the SST dependence of photo-physiological term in the VGPM, which has a greater sensitivity to variations in temperature than the WRM. In waters of 0 to 5C PP using the VGPM was 43% higher than WRM, from 5 to 10C the VGPM was 29% higher and from 10 to 15C the VGPM was 27% higher.