Biogeochemistry of greenhouse gases in coastal upwelling systems: Processes and sensitivity to global change

Major coastal upwelling systems are among the most productive marine ecosystems in the world. They contribute disproportionately to the cycling of carbon and nutrients in the ocean and influence marine biogeochemistry beyond their productive regions. Characterized by intense microbial respiration (both aerobic and anaerobic), major coastal upwelling systems are also hotspots for the production and outgassing of potent greenhouse gases (GHG) such as CO2, N2O, and CH4. Quantifying and understanding these roles in the context of a changing climate is therefore a subject of great interest. Here we provide a short synthesis of the current knowledge of the contributions of major coastal upwelling systems to the cycling of GHG. Despite variations within and among different systems, low-latitude coastal upwelling systems typically act as a net carbon source to the atmosphere, while those at higher latitudes function as weak sinks or remain neutral regarding atmospheric CO2. These systems also significantly contribute to oceanic N2O and CH4 emissions, although the extent of their contribution to the latter remains poorly constrained. We also overview recent and future changes to upwelling systems in the context of a warmer climate and discuss uncertainties and implications for GHG production. Although rapid coastal warming is anticipated in all major coastal upwelling systems, the future changes in upwelling-favorable winds and their implications within the context of increased stratification are uncertain. Finally, we examine the major challenges that impede our ability to accurately predict how major coastal upwelling systems will respond to future climate change, and present recommendations for future research to better capture ongoing changes and disentangle natural and forced variability

A comparative study of biological production in eastern boundary upwelling systems using an artificial neural network

Eastern Boundary Upwelling Systems (EBUS) are highly productive ocean regions. Yet, substantial differences in net primary production (NPP) exist within and between these systems for reasons that are still not fully understood. Here, we explore the leading physical processes and environmental factors controlling NPP in EBUS through a comparative study of the California, Canary, Benguela, and Humboldt Current systems. The NPP drivers are identified with the aid of an artificial neural network analysis based on self-organizing-maps (SOM). Our results suggest that in addition to the expected NPP enhancing effect of stronger equatorward alongshore wind, three factors have an inhibiting effect: (1) strong eddy activity, (2) narrow continental shelf, and (3) deep mixed layer. The co-variability of these 4 drivers defines in the context of the SOM a continuum of 100 patterns of NPP regimes in EBUS. These are grouped into 4 distinct classes using a Hierarchical Agglomerative Clustering (HAC) method. Our objective classification of EBUS reveals important variations of NPP regimes within each of the four EBUS, particularly in the Canary and Benguela Current systems. Our results show that the Atlantic EBUS are generally more productive and more sensitive to upwelling favorable winds because of weaker factors inhibiting NPP. Perturbations of alongshore winds associated with climate change may therefore lead to contrasting biological responses in the Atlantic and the Pacific EBUS.ISSN:1726-4170ISSN:1726-417

Rôle des tourbillons de méso-échelle océaniques dans la distribution et les flux air-mer de CO2 anthropique à l'échelle globale

Les tourbillons de méso-échelle jouent un rôle essentiel dans la dynamique océanique et particulièrement dans l'océan Austral. Dans cette thèse nous explorons le rôle des tourbillons de méso-échelle dans l'estimation du puits océanique de CO2 anthropique grâce à une approche multi-traceurs et multi-résolutions. Nos analyses montrent que les tourbillons résolus réduisent sensiblement les flux air-mer du CO2 anthropique permettant ainsi un meilleur accord avec les observations. Les mécanismes de l action des tourbillons sur les flux air-mer des traceurs sont identifiés. L'augmentation de la résolution horizontale affecte également la pénétration verticale des traceurs dans l'océan Austral, et particulièrement la ventilation des couches intermédiaires. Les tourbillons induisent également une modification substantielle de la distribution spatiale des régions source des Eaux Antarctiques Intermédiaires. Enfin, l'efficacité de la paramétrisation de Gent & McWilliams (1990) à représenter les effets des tourbillons dans les simulations de traceurs est évaluée.PARIS-BIUSJ-Thèses (751052125) / SudocPARIS-BIUSJ-Physique recherche (751052113) / SudocSudocFranceF

What controls biological production in coastal upwelling systems? Insights from a comparative modeling study

The magnitude of net primary production (NPP) in Eastern Boundary Upwelling Systems (EBUS) is traditionally viewed as directly reflecting the wind-driven upwelling intensity. Yet, different EBUS show different sensitivities of NPP to upwelling-favorable winds (Carr and Kearns, 2003). Here, using a comparative modeling study of the California Current System (California CS) and Canary Current System (Canary CS), we show how physical and environmental factors, such as light, temperature and cross-shore circulation modulate the response of NPP to upwelling strength. To this end, we made a series of eddy-resolving simulations of the two upwelling systems using the Regional Oceanic Modeling System (ROMS), coupled to a nitrogen-based Nutrient-Phytoplankton-Zooplankton-Detritus (NPZD) ecosystem model. Using identical ecological/biogeochemical parameters, our coupled model simulates a level of NPP in the California CS that is 50 % smaller than that in the Canary CS, in agreement with observationally based estimates. We find this much lower NPP in the California CS despite phytoplankton in this system having nearly 20 % higher nutrient concentrations available to fuel their growth. This conundrum can be explained by: (1) phytoplankton having a faster nutrient-replete growth in the Canary CS relative to the California CS; a consequence of more favorable light and temperature conditions in the Canary CS, and (2) the longer nearshore water residence times in the Canary CS, which permit a larger buildup of biomass in the upwelling zone, thereby enhancing NPP. The longer residence times in the Canary CS appear to be a result of the wider continental shelves and the lower mesoscale activity characterizing this upwelling system. This results in a weaker offshore export of nutrients and organic matter, thereby increasing local nutrient recycling and reducing the spatial decoupling between new and export production in the Canary CS. Our results suggest that climate change-induced perturbations such as upwelling favorable wind intensification might lead to contrasting biological responses in the California CS and the Canary CS, with major implications for the biogeochemical cycles and fisheries in these two ecosystems.ISSN:1726-4170ISSN:1726-417

Spatiotemporal variability of pCO2 and airsea CO2 uxes in the California Current System: An eddy-resolving modeling study

We quantify the CO2 source/sink nature of the California Current System (CalCS) and determine the drivers and processes behind the mean and spatiotemporal variability of the partial pressure of CO2 (pCO2) in the surface ocean. To this end, we analyze eddy-resolving, climatological simulations of a coupled physical-ecosystem-biogeochemical5 ocean model on the basis of the Regional Oceanic Modeling System (ROMS). The model-simulated pCO2 agrees very well with in situ observations over the entire domain with virtually no bias, but the model overestimates pCO2 in the nearshore 100km, and underestimates the observed temporal variability. In the annual mean, the entire CalCS within 800 km of the coast and from ∼ 33◦ N to 1046 ◦ N is essentially neutral with regard to atmospheric CO2. The model simulates an integrated uptake flux of − 0.9 TgCyr −1, corresponding to a very small average flux density of −0.05 molCm−2 yr−1, with an uncertainty of the order of ± 0.20 molCm −2 yr −1.This near zero flux is a consequence of an almost complete regional compensation between the strong outgassing in the nearshore region (first 100 km), with flux densities of more than 3molCm−2 yr−1 and a weaker, but more widespread uptake flux in the offshore region with an average flux density of −0.17 molCm−2yr−1. This pattern is primarily a result of the interaction between upwelling in the nearshore that brings waters with high concentrations of dissolved inorganic carbon (DIC) to the surface, and an intense biological drawdown of this DIC, driven by the nutrients that are upwelled20 together with the DIC. The biological drawdown occurs too slowly to prevent the escape of a substantial amount of CO2 into the atmosphere, but this is compensated by the biological generation of undersaturated conditions offshore of 100 km, permitting the CalCS to take up most of the escaped CO2. Thus, the biological pump over the entire CalCS is essentially 100 % efficient, making the preformed DIC and nutrient25 concentrations of the upwelled waters a primary determinant of the overall source/sink nature of the CalCS. The comparison of the standard simulation with one for preindustrial conditions show that the CalCS is taking up anthropogenic CO2 at a rate of about −1 molCm −2yr−1, implying that the region was a small source of CO2 to the atmosphere in preindustrial times. The air–sea CO2 fluxes vary substantially in time, both on seasonal and subseasonal timescales, largely driven by variations in surface ocean pCO2. There are important differences among the subregions. Notably, the total variance of the fluxes5 in the central nearshore CalCS is roughly 4–5 times larger than elsewhere. Most of the variability in pCO2 is associated with the seasonal cycle, except in the nearshore, where sub-seasonal variations driven by mesoscale processes dominate. In the regions off shore of 100km, changes in surface temperature are the main driver, while in the nearshore region, changes in surface temperature, as well as anomalies in DIC and 10alkalinity (Alk) owing to changes in circulation, biological productivity and air–sea CO2 fluxes dominate. The dominance of eddy-driven variability in the nearshore 100km leads to a complex spatiotemporal mosaic of surface ocean pCO2 and air–sea CO2 fluxes that require a substantial observational effort to determine the source/sink nature of this region reliably.ISSN:1810-6277ISSN:1810-628

Shrinking of the Arabian Sea oxygen minimum zone with climate change projected with a downscaled model

International audienceIn Arabian Sea (AS), landlocked northern boundary and strong seasonal productivity lead to the formation of one of the most intense open ocean Oxygen Minimum Zones (OMZs). Presence of this perennial OMZ has significant consequences on adjacent coastal fisheries and ecosystem. Simulations from CMIP5 suggest significant weakening of both monsoonal winds and productivity under high emission scenario. But the fate of AS OMZ in this scenario-whether it will expand or shrink-still remains elusive, mainly due to poor representation of extent and strength of AS OMZ in CMIP5 presentday simulations. To address this, we analyze the distribution of O 2 in AS from a subset of three contrasted CMIP5 simulations, and complemented with a set of regional downscaled model experiments which we forced at surface and open boundaries using information from those three CMIP5 models. We tested two regional downscaling approaches-with and without correction of CMIP5 biases with respect to observations. Using a set of sensitivity experiments, we disentangle the contributions of local (atmospheric) forcing vs. remote (at the lateral boundaries) forcing in driving the future projected O 2 changes. While CMIP5 projects either shrinking or expansion of the AS OMZ depending on the model, our downscaling experiments consistently project a shrinking of AS OMZ. We show that projected O 2 changes in OMZ layer are affected by both local and remote processes. In the southern AS, the main response to climate change is oxygenation that originates from the boundaries, and hence downscalled and CMIP5 model responses are similar. In contrast, in northern AS, downscaling yields a substantial reduction in O 2 projection discrepancies because of a minimal influence of remote forcing there leading to a stronger sensitivity to improved local physics and improved model representation of present-day conditions. We find that when corrected for present-day biases, projected deoxygenation in the northern AS is shallower. Our findings indicate the importance of downscaling of global models in regions where local forcing is dominant, and the need for correcting global model biases with respect to observations to reduce uncertainties in future O 2 projections

Spatiotemporal variability and drivers of pCO2 and air–sea CO2 fluxes in the California Current System: an eddy-resolving modeling study

We quantify the CO2 source/sink nature of the California Current System (CalCS) and determine the drivers and processes behind the mean and spatiotemporal variability of the partial pressure of CO2 (pCO2) in the surface ocean. To this end, we analyze eddy-resolving, climatological simulations of a coupled physical–biogeochemical oceanic model on the basis of the Regional Oceanic Modeling System (ROMS). In the annual mean, the entire CalCS within 800 km of the coast and from ∼33° N to 46° N is essentially neutral with regard to atmospheric CO2: the model simulates an integrated uptake flux of −0.9 ± 3.6 Tg C yr−1, corresponding to an average flux density of −0.05 ± 0.20 mol C m−2 yr−1. This near zero flux is a consequence of an almost complete regional compensation between (i) strong outgassing in the nearshore region (first 100 km) that brings waters with high concentrations of dissolved inorganic carbon (DIC) to the surface and (ii) and a weaker, but more widespread uptake flux in the offshore region due to an intense biological reduction of this DIC, driven by the nutrients that are upwelled together with the DIC. The air–sea CO2 fluxes vary substantially in time, both on seasonal and sub-seasonal timescales, largely driven by variations in surface ocean pCO2. Most of the variability in pCO2 is associated with the seasonal cycle, with the exception of the nearshore region, where sub-seasonal variations driven by mesoscale processes dominate. In the regions offshore of 100 km, changes in surface temperature are the main driver, while in the nearshore region, changes in surface temperature, as well as anomalies in DIC and alkalinity (Alk) owing to changes in circulation, biological productivity and air–sea CO2 fluxes dominate. The prevalence of eddy-driven variability in the nearshore 100 km leads to a complex spatiotemporal mosaic of surface ocean pCO2 and air–sea CO2 fluxes that require a substantial observational effort to determine the source/sink nature of this region reliably.ISSN:1726-4170ISSN:1726-417

Strong Intensification of the Arabian Sea Oxygen Minimum Zone in Response to Arabian Gulf Warming

International audienceThe highly saline, oxygen-saturated waters of the Arabian Gulf (hereafter the Gulf) sink to intermediate depths (200-300 m) when they enter the Arabian Sea, ventilating the World's thickest oxygen minimum zone (OMZ). Here, we investigate the impacts of a warming of the Gulf consistent with climate change projections on the intensity of this OMZ. Using a series of eddy-resolving model simulations, we show that the warming of the Gulf waters increases their buoyancy and hence limits their contribution to the ventilation of intermediate depths. This leads to an intensification of the OMZ and an increase in denitrification that depletes subsurface nitrate and limits deoxygenation at depth. The projected future concomitant increase of Gulf salinity only partially reduces the OMZ intensification. Our findings highlight the importance of the Arabian marginal seas for the biogeochemistry of the North Indian Ocean and stress the need for improving their representation in global climate models

Intensification and deepening of the Arabian Sea Oxygen Minimum Zone in response to increase in Indian monsoon wind intensity

International audienceThe decline in oxygen supply to the ocean associated with global warming of sea-surface temperatures is expected to expand the oxygen minimum zones (OMZs). This global trend can be attenuated or amplified by regional processes. In the Arabian Sea, the World's thickest OMZ is highly vulnerable to changes in the Indian monsoon wind. Evidence from paleo records and future climate projections indicate strong variations of the Indian monsoon wind intensity over climatic timescales. Yet, the response of the OMZ to these wind changes remains poorly understood and its amplitude and timescale unexplored. Here, we investigate the impacts of perturbations in Indian monsoon wind intensity (from -50% to +50%) on the size and intensity of the Arabian Sea OMZ, and examine the biogeochemical and ecological implications of these changes. To this end, we conducted a series of eddy-resolving simulations of the Arabian Sea using the Regional Oceanic Modeling System (ROMS) coupled to a nitrogen based Nutrient-Phytoplankton-Zooplankton-Detritus (NPZD) ecosystem model that includes a representation of the O2 cycle. We show that the Arabian Sea productivity increases and its OMZ expands and deepens in response to monsoon wind intensification. These responses are dominated by the perturbation of the summer monsoon wind, whereas the changes in the winter monsoon wind play a secondary role. While the productivity responds quickly and nearly linearly to wind increase (i.e., on a timescale of years), the OMZ response is much slower (i.e., a timescale of decades). Our analysis reveals that the OMZ expansion at depth is driven by increased oxygen biological consumption, whereas its surface weakening is induced by increased lateral ventilation. The enhanced lateral ventilation favors episodic intrusions of oxic waters in the lower epipelagic zone (100-200m) of the western and central Arabian Sea, leading to intermittent expansions of habitats and a more frequent alternation of hypoxic and oxic conditions there. The increased productivity and deepening of the OMZ also lead to a strong intensification of denitrification at depth, resulting in a substantial amplification of fixed nitrogen depletion in the Arabian Sea. We conclude that changes in the Indian monsoon can affect, on longer timescales, the large-scale biogeochemical cycles of nitrogen and carbon, with a positive feedback on climate change in the case of stronger winds