Sedimentary environments and climate change: a case study (late Miocene, central Crete)

In this study, conditions of deposition and stratigraphical architecture of Neogene (Tortonian, 11-6,7Ma) sediments of southern central Crete were analysed. In order to improve resolution of paleoclimatic data, new methods were applied to quantify environmental parameters and to increase the chronostratigraphic resolution in shallow water sediments. A relationship between paleoenvironmental change observed on Crete and global processes was established and a depositional model was developed. Based on a detailed analysis of the distribution of non geniculate coralline red algae, index values for water temperature and water depth were established and tested with the distribution patterns of benthic foraminifera and symbiont-bearing corals. Calcite shelled bivalves were sampled from the Algarve coast (southern Portugal) and central Crete and then 87Sr/86Sr was measured. A high resolution chronostratigraphy was developed based on the correlation between fluctuations in Sr ratios in the measured sections and in a late Miocene global seawater Sr isotope reference curve. Applying this method, a time frame was established to compare paleoenvironmental data from southern central Crete with global information on climate change reflected in oxygen isotope data. The comparison between paleotemperature data based on red algae and global oxygen isotope data showed that the employed index values reflect global change in temperature. Data indicate a warm interval during earliest Tortonian, a second short warm interval between 10 and 9,5Ma, a longer climatic optimum between 9 and 8Ma and an interval of increasing temperatures in the latest Tortonian. The distribution of coral reefs and carpets shows that during the warm intervals, the depositional environment became tropical while temperate climates prevailed during the cold interval. Since relative tectonic movements after initial half-graben formation in the early Tortonian were low in southern central Crete, sedimentary successions strongly respond to global sea-level fluctuation. A characteristic sedimentary succession formed during a 3rd order sea-level cycle: It comprises mixed siliciclastic-limestone deposited during sea-level fall and lowstand, homogenous red algal deposits formed during sea-level rise and coral carpets formed during late rise and highstand. Individual beds in the succession reflect glacioeustatic fluctuations that are most prominent in the mixed siliciclastic-limestone interval. These results confirm the fact that sedimentary successions deposited at the critical threshold between temperate and tropical environments develop characteristic changes in depositional systems and biotic associations that can be used to assemble paleoclimatic datasets.

External drivers of biogeochemical cycles in a tropical montane forest in Ecuador

Successful conservation of tropical montane forest, one of the most threatened ecosystems on earth, requires detailed knowledge of its biogeochemistry. Of particular interest is the response of the biogeochemical element cycles to external influences such as element deposition or climate change. Therefore the overall objective of my study was to contribute to improved understanding of role and functioning of the Andean tropical montane forest. In detail, my objectives were to determine (1) the role of long-range transported aerosols and their transport mechanisms, and (2) the role of short-term extreme climatic events for the element budget of Andean tropical forest. In a whole-catchment approach including three 8-13 ha microcatchments under tropical montane forest on the east-exposed slope of the eastern cordillera in the south Ecuadorian Andes at 1850-2200 m above sea level I monitored at least in weekly resolution the concentrations and fluxes of Ca, Mg, Na, K, NO3-N, NH4-N, DON, P, S, TOC, Mn, and Al in bulk deposition, throughfall, litter leachate, soil solution at the 0.15 and 0.3 m depths, and runoff between May 1998 and April 2003. I also used meteorological data from my study area collected by cooperating researchers and the Brazilian meteorological service (INPE), as well as remote sensing products of the North American and European space agencies NASA and ESA. My results show that (1) there was a strong interannual variation in deposition of Ca [4.4-29 kg ha-1 a-1], Mg [1.6-12], and K [9.8-30]) between 1998 and 2003. High deposition changed the Ca and Mg budgets of the catchments from loss to retention, suggesting that the additionally available Ca and Mg was used by the ecosystem. Increased base metal deposition was related to dust outbursts of the Sahara and an Amazonian precipitation pattern with trans-regional dry spells allowing for dust transport to the Andes. The increased base metal deposition coincided with a strong La Niña event in 1999/2000. There were also significantly elevated H+, N, and Mn depositions during the annual biomass burning period in the Amazon basin. Elevated H+ deposition during the biomass burning period caused elevated base metal loss from the canopy and the organic horizon and deteriorated already low base metal supply of the vegetation. Nitrogen was only retained during biomass burning but not during non-fire conditions when deposition was much smaller. Therefore biomass burning-related aerosol emissions in Amazonia seem large enough to substantially increase element deposition at the western rim of Amazonia. Particularly the related increase of acid deposition impoverishes already base-metal scarce ecosystems. As biomass burning is most intense during El Niño situations, a shortened ENSO cycle because of global warming likely enhances the acid deposition at my study forest. (2) Storm events causing near-surface water flow through C- and nutrient-rich topsoil during rainstorms were the major export pathway for C, N, Al, and Mn (contributing >50% to the total export of these elements). Near-surface flow also accounted for one third of total base metal export. This demonstrates that storm-event related near-surface flow markedly affects the cycling of many nutrients in steep tropical montane forests. Changes in the rainfall regime possibly associated with global climate change will therefore also change element export from the study forest. Element budgets of Andean tropical montane rain forest proved to be markedly affected by long-range transport of Saharan dust, biomass burning-related aerosols, or strong rainfalls during storm events. Thus, increased acid and nutrient deposition and the global climate change probably drive the tropical montane forest to another state with unknown consequences for its functions and biological diversity.

Biogenic volatile organic compounds in tropical, temperate and boreal forest ecosystems

Biogene flüchtige organische Verbindungen (BFOV) werden in großen Mengen aus terrestrischenrnÖkosystemen, insbesondere aus Wäldern und Wiesen, in die untere Troposphäre emittiert. Austausch-rnFlüsse von BFOVs sind in der troposphärischen Chemie wichtig, weil sie eine bedeutende Rolle in derrnOzon- und Aerosolbildung haben. Trotzdem bleiben die zeitliche und räumliche Änderung der BFOVrnEmissionen und ihre Rolle in Bildung und Wachstum von Aerosolen ungewiss.rnDer Fokus dieser Arbeit liegt auf der in-situ Anwendung der Protonen Transfer ReaktionsrnMassenspektrometrie (PTR-MS) und der Messung von biogenen flüchtigen organischen Verbindungen inrnnordländischen, gemäßigten und tropischen Waldökosystemen während drei unterschiedlicherrnFeldmesskampagnen. Der Hauptvorteil der PTR-MS-Technik liegt in der hohen Messungsfrequenz,rnwodurch eine eventuelle Änderung in der Atmosphäre durch Transport, Vermischung und Chemiernonline beobachtet werden kann. Die PTR-MS-Messungen wurden zweimal am Boden aus und einmalrnvon einem Forschungsflugzug durchgeführt.rnIn Kapitel 3 werden die PTR-MS-Daten, gesammelt während der Flugmesskampagne über demrntropischen Regenwald, vorgelegt. Diese Studie zeigt den Belang der Grenzschichtdynamik für diernVerteilung von Spurengasen mittels eines eindimensionalen Säule - Chemie und KlimaModells (SCM).rnDer Tagesablauf von Isopren zeigte zwischen 14:00 und 16:15 Uhr lokaler Zeit einen Mittelwert vonrn5.4 ppbv auf der Höhe der Baumspitzen und von 3.3 ppbv über 300 m Höhe. Dies deutet darauf hin, dassrnsowohl der turbulente Austausch als auch die hohe Reaktionsfähigkeit von Isopren mit den OxidantienrnOH und Ozon eine wichtige Rolle spielen. Nach dem Ausgleich von chemischen Verlusten undrnEntrainment (Ein- und Ausmischung von Luft an der Grenzschicht), wurde ein Fluss vonrn8.4 mg Isopren m-2h-1 unter teilweise bewölkten Bedingungen für den tropischen Regenwald in derrnGuyanregion abgeschätzt. Dies entspricht einem täglichen Emissionsfluss von 28 mg Isopren prornQuadratmeter.rnIm Kapitel 4 werden die Messungen, welche auf einer Hügellage in gemäßigter Breite inrnsüddeutschland stattgefunden haben, diskutiert. Bei diesem Standort ist die Grenzschicht nachts unter diernStandorthöhe abgefallen, was den Einsatzort von Emissionen abgesondert hatte. Während diernGrenzschicht morgens wieder über die Höhe des Einsatzortes anstieg, konnten die eingeschlossenenrnnächtlichen Emissionen innerhalb der bodennahen Schicht beobachtet werden. Außerdem wurde einrndeutlicher Anstieg von flüchtigen organischen Verbindungen gemessen, wenn die Luftmassen überrnMünchen geführt wurden oder wenn verschmutzte Luftmassen aus dem Po-Tal über die Alpen nachrnDeutschland transportiert wurden. Daten von dieser Kampagne wurden genutzt, um die Änderungen inrndem Mischungsverhältnis der flüchtigen organischen Verbindungen, verbunden mit dem Durchfluss vonrnwarmen und kalten Wetterfronten sowie bei Regen zu untersuchen.rnIm Kapitel 5 werden PTR-MS-Messungen aus dem nördlichen Nadelwaldgürtel beschrieben. Starkernnächtliche Inversionen mit einer niedrigen Windgeschwindigkeit fingen die Emissionen vonrnnahegelegenen Kiefernwäldern und andere BFOV-Quellen ab, was zu hohen nächtlichen BFOVMischverhältnissenrnführte. Partikelereignisse wurden für Tag und Nacht detailliert analysiert. Diernnächtlichen Partikelereignisse erfolgten synchron mit starken extremen von Monoterpenen, obwohl dasrnzweite Ereignis Kernbildung einschloss und nicht mit Schwefelsäure korrelierte. Die MonoterpenrnMischungsverhältnisse von über 16 ppbv waren unerwartet hoch für diese Jahreszeit. NiedrigernWindgeschwindigkeiten und die Auswertung von Rückwärtstrajektorien deuten auf eine konzentrierternQuelle in der Nähe von Hyytiälä hin. Die optische Stereoisomerie von Monoterpenen hat bestätigt, dassrndie Quelle unnatürlich ist, da das Verhältnis von [(+)-α-pinen]/[(−)-α-pinen] viel höher ist als dasrnnatürliches Verhältnis der beiden Enantiomeren.

Lateglacial and early Holocene glaciation in the tropical Andes caused by La Niña-like conditions

The response of the tropics to North Atlantic cold events, such as Heinrich Event I (H-I, ∼ 17–15 ka) and the Younger Dryas (YD, 12.7–11.5 ka), is still one of the most tantalizing, yet unresolved issues in paleoclimatology. The advent of surface exposure dating has therefore instigated the establishment of glacial chronologies in the tropical Andes to investigate potential climate teleconnections. Here, we present new exposure ages from the Cordillera Cochabamba (17°17′S), Bolivia, that reveal glacial advances during H-I and YD, as well as during the Early Holocene. Our chronology correlates well with cold sea surface temperatures in the eastern tropical Pacific, which indicates that La Niña-like conditions, i.e. forcings intrinsic to the tropics, played a key role for moisture advection and glaciation in the tropical Andes.

Sensitivity of atmospheric CO2 and climate to explosive volcanic eruptions

Impacts of low-latitude, explosive volcanic eruptions on climate and the carbon cycle are quantified by forcing a comprehensive, fully coupled carbon cycle-climate model with pulse-like stratospheric aerosol optical depth changes. The model represents the radiative and dynamical response of the climate system to volcanic eruptions and simulates a decrease of global and regional atmospheric surface temperature, regionally distinct changes in precipitation, a positive phase of the North Atlantic Oscillation, and a decrease in atmospheric CO2 after volcanic eruptions. The volcanic-induced cooling reduces overturning rates in tropical soils, which dominates over reduced litter input due to soil moisture decrease, resulting in higher land carbon inventories for several decades. The perturbation in the ocean carbon inventory changes sign from an initial weak carbon sink to a carbon source. Positive carbon and negative temperature anomalies in subsurface waters last up to several decades. The multi-decadal decrease in atmospheric CO2 yields a small additional radiative forcing that amplifies the cooling and perturbs the Earth System on longer time scales than the atmospheric residence time of volcanic aerosols. In addition, century-scale global warming simulations with and without volcanic eruptions over the historical period show that the ocean integrates volcanic radiative cooling and responds for different physical and biogeochemical parameters such as steric sea level or dissolved oxygen. Results from a suite of sensitivity simulations with different magnitudes of stratospheric aerosol optical depth changes and from global warming simulations show that the carbon cycle-climate sensitivity γ, expressed as change in atmospheric CO2 per unit change in global mean surface temperature, depends on the magnitude and temporal evolution of the perturbation, and time scale of interest. On decadal time scales, modeled γ is several times larger for a Pinatubo-like eruption than for the industrial period and for a high emission, 21st century scenario.

Atlantic hurricanes and associated insurance loss potentials in future climate scenarios: limitations of high-resolution AGCM simulations

Potential future changes in tropical cyclone (TC) characteristics are among the more serious regional threats of global climate change. Therefore, a better understanding of how anthropogenic climate change may affect TCs and how these changes translate in socio-economic impacts is required. Here, we apply a TC detection and tracking method that was developed for ERA-40 data to time-slice experiments of two atmospheric general circulation models, namely the fifth version of the European Centre model of Hamburg model (MPI, Hamburg, Germany, T213) and the Japan Meteorological Agency/ Meteorological research Institute model (MRI, Tsukuba city, Japan, TL959). For each model, two climate simulations are available: a control simulation for present-day conditions to evaluate the model against observations, and a scenario simulation to assess future changes. The evaluation of the control simulations shows that the number of intense storms is underestimated due to the model resolution. To overcome this deficiency, simulated cyclone intensities are scaled to the best track data leading to a better representation of the TC intensities. Both models project an increased number of major hurricanes and modified trajectories in their scenario simulations. These changes have an effect on the projected loss potentials. However, these state-of-the-art models still yield contradicting results, and therefore they are not yet suitable to provide robust estimates of losses due to uncertainties in simulated hurricane intensity, location and frequency.

Climate-induced interannual variability of marine primary and export production in three global coupled climate carbon cycle models

Fully coupled climate carbon cycle models are sophisticated tools that are used to predict future climate change and its impact on the land and ocean carbon cycles. These models should be able to adequately represent natural variability, requiring model validation by observations. The present study focuses on the ocean carbon cycle component, in particular the spatial and temporal variability in net primary productivity (PP) and export production (EP) of particulate organic carbon (POC). Results from three coupled climate carbon cycle models (IPSL, MPIM, NCAR) are compared with observation-based estimates derived from satellite measurements of ocean colour and results from inverse modelling (data assimilation). Satellite observations of ocean colour have shown that temporal variability of PP on the global scale is largely dominated by the permanently stratified, low-latitude ocean (Behrenfeld et al., 2006) with stronger stratification (higher sea surface temperature; SST) being associated with negative PP anomalies. Results from all three coupled models confirm the role of the low-latitude, permanently stratified ocean for anomalies in globally integrated PP, but only one model (IPSL) also reproduces the inverse relationship between stratification (SST) and PP. An adequate representation of iron and macronutrient co-limitation of phytoplankton growth in the tropical ocean has shown to be the crucial mechanism determining the capability of the models to reproduce observed interactions between climate and PP.

Spectral biases in tree-ring climate proxies

External forcing and internal dynamics result in climate system variability ranging from sub-daily weather to multi-centennial trends and beyond1, 2. State-of-the-art palaeoclimatic methods routinely use hydroclimatic proxies to reconstruct temperature (for example, refs 3, 4), possibly blurring differences in the variability continuum of temperature and precipitation before the instrumental period. Here, we assess the spectral characteristics of temperature and precipitation fluctuations in observations, model simulations and proxy records across the globe. We find that whereas an ensemble of different general circulation models represents patterns captured in instrumental measurements, such as land–ocean contrasts and enhanced low-frequency tropical variability, the tree-ring-dominated proxy collection does not. The observed dominance of inter-annual precipitation fluctuations is not reflected in the annually resolved hydroclimatic proxy records. Likewise, temperature-sensitive proxies overestimate, on average, the ratio of low- to high-frequency variability. These spectral biases in the proxy records seem to propagate into multi-proxy climate reconstructions for which we observe an overestimation of low-frequency signals. Thus, a proper representation of the high- to low-frequency spectrum in proxy records is needed to reduce uncertainties in climate reconstruction efforts.

LakeMIP Kivu: evaluating the representation of a large, deep tropical lake by a set of one-dimensional lake models

The African great lakes are of utmost importance for the local economy (fishing), as well as being essential to the survival of the local people. During the past decades, these lakes experienced fast changes in ecosystem structure and functioning, and their future evolution is a major concern. In this study, for the first time a set of one-dimensional lake models are evaluated for Lake Kivu (2.28°S; 28.98°E), East Africa. The unique limnology of this meromictic lake, with the importance of salinity and subsurface springs in a tropical high-altitude climate, presents a worthy challenge to the seven models involved in the Lake Model Intercomparison Project (LakeMIP). Meteorological observations from two automatic weather stations are used to drive the models, whereas a unique dataset, containing over 150 temperature profiles recorded since 2002, is used to assess the model’s performance. Simulations are performed over the freshwater layer only (60 m) and over the average lake depth (240 m), since salinity increases with depth below 60 m in Lake Kivu and some lake models do not account for the influence of salinity upon lake stratification. All models are able to reproduce the mixing seasonality in Lake Kivu, as well as the magnitude and seasonal cycle of the lake enthalpy change. Differences between the models can be ascribed to variations in the treatment of the radiative forcing and the computation of the turbulent heat fluxes. Fluctuations in wind velocity and solar radiation explain inter-annual variability of observed water column temperatures. The good agreement between the deep simulations and the observed meromictic stratification also shows that a subset of models is able to account for the salinity- and geothermal-induced effects upon deep-water stratification. Finally, based on the strengths and weaknesses discerned in this study, an informed choice of a one-dimensional lake model for a given research purpose becomes possible.

Ecosystem services, biodiversity and environmental change in a tropical mountain ecosystem of south Ecuador

An interdisciplinary research unit consisting of 30 teams in the natural, economic and social sciences analyzed biodiversity and ecosystem services of a mountain rainforest ecosystem in the hotspot of the tropical Andes, with special reference to past, current and future environmental changes. The group assessed ecosystem services using data from ecological field and scenario-driven model experiments, and with the help of comparative field surveys of the natural forest and its anthropogenic replacement system for agriculture. The book offers insights into the impacts of environmental change on various service categories mentioned in the Millennium Ecosystem Assessment (2005): cultural, regulating, supporting and provisioning ecosystem services. Examples focus on biodiversity of plants and animals including trophic networks, and abiotic/biotic parameters such as soils, regional climate, water, nutrient and sediment cycles. The types of threats considered include land use and climate changes, as well as atmospheric fertilization. In terms of regulating and provisioning services, the emphasis is primarily on water regulation and supply as well as climate regulation and carbon sequestration. With regard to provisioning services, the synthesis of the book provides science-based recommendations for a sustainable land use portfolio including several options such as forestry, pasture management and the practices of indigenous peoples. In closing, the authors show how they integrated the local society by pursuing capacity building in compliance with the CBD-ABS (Convention on Biological Diversity - Access and Benefit Sharing), in the form of education and knowledge transfer for application.

Are Greenhouse Gas Signals of Northern Hemisphere winter extra-tropical cyclone activity dependent on the identification and tracking algorithm?

For Northern Hemisphere extra-tropical cyclone activity, the dependency of a potential anthropogenic climate change signal on the identification method applied is analysed. This study investigates the impact of the used algorithm on the changing signal, not the robustness of the climate change signal itself. Using one single transient AOGCM simulation as standard input for eleven state-of-the-art identification methods, the patterns of model simulated present day climatologies are found to be close to those computed from re-analysis, independent of the method applied. Although differences in the total number of cyclones identified exist, the climate change signals (IPCC SRES A1B) in the model run considered are largely similar between methods for all cyclones. Taking into account all tracks, decreasing numbers are found in the Mediterranean, the Arctic in the Barents and Greenland Seas, the mid-latitude Pacific and North America. Changing patterns are even more similar, if only the most severe systems are considered: the methods reveal a coherent statistically significant increase in frequency over the eastern North Atlantic and North Pacific. We found that the differences between the methods considered are largely due to the different role of weaker systems in the specific methods.

Protecting Tropical Forests through Public Procurement of Emission Reductions

The rate of destruction of tropical forests continues to accelerate at an alarming rate contributing to an important fraction of overall greenhouse gas emissions. In recent years, much hope has been vested in the emerging REDD+ framework under the UN Framework Convention on Climate Change (UNFCCC), which aims at creating an international incentive system to reduce emissions from deforestation and forest degradation. This paper argues that in the absence of an international consensus on the design of results-based payments, “bottom-up” initiatives should take the lead and explore new avenues. It suggests that a call for tender for REDD+ credits might both assist in leveraging private investments and spending scarce public funds in a cost-efficient manner. The paper discusses the pros and cons of results-based approaches, provides an overview of the goals and principles that govern public procurement and discusses their relevance for the purchase of REDD+ credits, in particular within the ambit of the European Union.

Stalagmite water content as a proxy for drip water supply in tropical and subtropical areas

In this pilot study water was extracted from samples of two Holocene stalagmites from Socotra Island, Yemen, and one Eemian stalagmite from southern continental Yemen. The amount of water extracted per unit mass of stalagmite rock, termed "water yield" hereafter, serves as a measure of its total water content. Based on direct correlation plots of water yields and δ18Ocalcite and on regime shift analyses, we demonstrate that for the studied stalagmites the water yield records vary systematically with the corresponding oxygen isotopic compositions of the calcite (δ18Ocalcite). Within each stalagmite lower δ18Ocalcite values are accompanied by lower water yields and vice versa. The δ18Ocalcite records of the studied stalagmites have previously been interpreted to predominantly reflect the amount of rainfall in the area; thus, water yields can be linked to drip water supply. Higher, and therefore more continuous drip water supply caused by higher rainfall rates, supports homogeneous deposition of calcite with low porosity and therefore a small fraction of water-filled inclusions, resulting in low water yields of the respective samples. A reduction of drip water supply fosters irregular growth of calcite with higher porosity, leading to an increase of the fraction of water-filled inclusions and thus higher water yields. The results are consistent with the literature on stalagmite growth and supported by optical inspection of thin sections of our samples. We propose that for a stalagmite from a dry tropical or subtropical area, its water yield record represents a novel paleo-climate proxy recording changes in drip water supply, which can in turn be interpreted in terms of associated rainfall rates.

Provenance versus weathering control on the composition of tropical river mud (southern Africa)

This study presents an integrated mineralogical-geochemical data base on fine-grained sediments transported by all major rivers of southern Africa, including the Zambezi, Okavango, Limpopo, Olifants, Orange and Kunene. Clay mineralogy, bulk geochemistry, Sr and Nd isotopic signatures of river mud, considered as proxy of suspended load, are used to investigate the influence of source-rock lithology and weathering intensity on the composition of clay and silt produced in subequatorial to subtropical latitudes. Depletion in mobile alkali and alkaline-earth metals, minor in arid Namibia, is strong in the Okavango, Kwando and Upper Zambezi catchments, where recycling is also extensive. Element removal is most significant for Na, and to a lesser extent for Sr. Depletion in K, Ca and other elements, negligible in Namibia, is moderate elsewhere. The most widespread clay minerals are smectite, dominant in muds derived from Karoo or Etendeka flood basalts, or illite and chlorite, dominant in muds derived from metasedimentary rocks of the Damara Orogen or Zimbabwe Craton. Kaolinite represents 30-40% of clay minerals only in Okavango and Upper Zambezi sediments sourced in humid subequatorial Angola and Zambia. After subtracting the effects of recycling and of local accumulation of authigenic carbonates in soils, the regional distribution of clay minerals and chemical indices consistently reflect weathering intensity primarily controlled by climate. Bulk geochemistry identifies most clearly volcaniclastic sediments and mafic sources in general, but cannot discriminate the other sources of detritus in detail. Instead, Sr and Nd isotopic fingerprints are insensitive to weathering, and thus mirror faithfully the tectonic structure of the southern African continent. Isotopic tools thus represent a much firmer basis than bulk geochemistry or clay mineralogy in the provenance study of mudrocks.

Multiyear predictability of tropical marine productivity

With the emergence of decadal predictability simulations, research toward forecasting variations of the climate system now covers a large range of timescales. However, assessment of the capacity to predict natural variations of relevant biogeochemical variables like carbon fluxes, pH, or marine primary productivity remains unexplored. Among these, the net primary productivity (NPP) is of particular relevance in a forecasting perspective. Indeed, in regions like the tropical Pacific (30°N–30°S), NPP exhibits natural fluctuations at interannual to decadal timescales that have large impacts on marine ecosystems and fisheries. Here, we investigate predictions of NPP variations over the last decades (i.e., from 1997 to 2011) with an Earth system model within the tropical Pacific. Results suggest a predictive skill for NPP of 3 y, which is higher than that of sea surface temperature (1 y). We attribute the higher predictability of NPP to the poleward advection of nutrient anomalies (nitrate and iron), which sustain fluctuations in phytoplankton productivity over several years. These results open previously unidentified perspectives to the development of science-based management approaches to marine resources relying on integrated physical-biogeochemical forecasting systems.