Sensitivity analysis of an ocean carbon cycle model in the North Atlantic: an investigation of parameters affecting the air-sea CO2 flux, primary production and export of detritus

The sensitivity of the biological parameters in a nutrient-phytoplankton-zooplankton-detritus (NPZD) model in the calculation of the air-sea CO2 flux, primary production and detrital export is analysed. We explore the effect on these outputs of variation in the values of the twenty parameters that control ocean ecosystem growth in a 1-D formulation of the UK Met Office HadOCC NPZD model used in GCMs. We use and compare the results from one-at-a-time and all-at-a-time perturbations performed at three sites in the EuroSITES European Ocean Observatory Network: the Central Irminger Sea (60° N 40° W), the Porcupine Abyssal Plain (49° N 16° W) and the European Station for Time series in the Ocean Canary Islands (29° N 15° W). Reasonable changes to the values of key parameters are shown to have a large effect on the calculation of the air-sea CO2 flux, primary production, and export of biological detritus to the deep ocean. Changes in the values of key parameters have a greater effect in more productive regions than in less productive areas. The most sensitive parameters are generally found to be those controlling well-established ocean ecosystem parameterisations widely used in many NPZD-type models. The air-sea CO2 flux is most influenced by variation in the parameters that control phytoplankton growth, detrital sinking and carbonate production by phytoplankton (the rain ratio). Primary production is most sensitive to the parameters that define the shape of the photosynthesis-irradiance curve. Export production is most sensitive to the parameters that control the rate of detrital sinking and the remineralisation of detritus.

The impact of diurnal variability in sea surface temperature on the atlantic air-sea CO2 flux

The effect of diurnal variations in sea surface temperature (SST) on the air-sea flux of CO2 over the central Atlantic ocean and Mediterranean Sea (60 S–60 N, 60 W–45 E) is evaluated for 2005–2006. We use high spatial resolution hourly satellite ocean skin temperature data to determine the diurnal warming (ΔSST). The CO2 flux is then computed using three different temperature fields – a foundation temperature (Tf, measured at a depth where there is no diurnal variation), Tf, plus the hourly ΔSST and Tf, plus the monthly average of the ΔSSTs. This is done in conjunction with a physically-based parameterisation for the gas transfer velocity (NOAA-COARE). The differences between the fluxes evaluated for these three different temperature fields quantify the effects of both diurnal warming and diurnal covariations. We find that including diurnal warming increases the CO2 flux out of this region of the Atlantic for 2005–2006 from 9.6 Tg C a−1 to 30.4 Tg C a−1 (hourly ΔSST) and 31.2 Tg C a−1 (monthly average of ΔSST measurements). Diurnal warming in this region, therefore, has a large impact on the annual net CO2 flux but diurnal covariations are negligible. However, in this region of the Atlantic the uptake and outgassing of CO2 is approximately balanced over the annual cycle, so although we find diurnal warming has a very large effect here, the Atlantic as a whole is a very strong carbon sink (e.g. −920 Tg C a−1 Takahashi et al., 2002) making this is a small contribution to the Atlantic carbon budget.

Mixed layer lateral eddy fluxes mediated by air-sea interaction

The modulation of air–sea heat fluxes by geostrophic eddies due to the stirring of temperature at the sea surface is discussed and quantified. It is argued that the damping of eddy temperature variance by such air–sea fluxes enhances the dissipation of surface temperature fields. Depending on the time scale of damping relative to that of the eddying motions, surface eddy diffusivities can be significantly enhanced over interior values. The issues are explored and quantified in a controlled setting by driving a tracer field, a proxy for sea surface temperature, with surface altimetric observations in the Antarctic Circumpolar Current (ACC) of the Southern Ocean. A new, tracer-based diagnostic of eddy diffusivity is introduced, which is related to the Nakamura effective diffusivity. Using this, the mixed layer lateral eddy diffusivities associated with (i) eddy stirring and small-scale mixing and (ii) surface damping by air–sea interaction is quantified. In the ACC, a diffusivity associated with surface damping of a comparable magnitude to that associated with eddy stirring (;500 m2 s21) is found. In frontal regions prevalent in the ACC, an augmentation of surface lateral eddy diffusivities of this magnitude is equivalent to an air–sea flux of 100 W m22 acting over a mixed layer depth of 100 m, a very significant effect. Finally, the implications for other tracer fields such as salinity, dissolved gases, and chlorophyll are discussed. Different tracers are found to have surface eddy diffusivities that differ significantly in magnitude.

An assessment of air-sea heat fluxes from ocean and coupled reanalyses

Sixteen monthly air–sea heat flux products from global ocean/coupled reanalyses are compared over 1993–2009 as part of the Ocean Reanalysis Intercomparison Project (ORA-IP). Objectives include assessing the global heat closure, the consistency of temporal variability, comparison with other flux products, and documenting errors against in situ flux measurements at a number of OceanSITES moorings. The ensemble of 16 ORA-IP flux estimates has a global positive bias over 1993–2009 of 4.2 ± 1.1 W m−2. Residual heat gain (i.e., surface flux + assimilation increments) is reduced to a small positive imbalance (typically, +1–2 W m−2). This compensation between surface fluxes and assimilation increments is concentrated in the upper 100 m. Implied steady meridional heat transports also improve by including assimilation sources, except near the equator. The ensemble spread in surface heat fluxes is dominated by turbulent fluxes (>40 W m−2 over the western boundary currents). The mean seasonal cycle is highly consistent, with variability between products mostly <10 W m−2. The interannual variability has consistent signal-to-noise ratio (~2) throughout the equatorial Pacific, reflecting ENSO variability. Comparisons at tropical buoy sites (10°S–15°N) over 2007–2009 showed too little ocean heat gain (i.e., flux into the ocean) in ORA-IP (up to 1/3 smaller than buoy measurements) primarily due to latent heat flux errors in ORA-IP. Comparisons with the Stratus buoy (20°S, 85°W) over a longer period, 2001–2009, also show the ORA-IP ensemble has 16 W m−2 smaller net heat gain, nearly all of which is due to too much latent cooling caused by differences in surface winds imposed in ORA-IP.

Air-sea interaction studies at the station occupied by R/V Sirius

The transfer formulas are used to compute hourly values of sensible and latent heat exchanges for the period July 29 to July 31, 1974. Also, values of eddy shearing stress at the sea surface, Bowen's ratio and dissipation of kinetic energy are computed. The data used cover part of the GATE period for station 20 occupied by the Brazilian Naval Ship R/V Sirius in the Equatorial Atlantic (0730 N 4000 W). The variations in the computed values are studied in relation to rainfall and the synoptic situation. © 1980 D. Reidel Publishing Co.

Air-Sea Interactions of Natural Long-Lived Greenhouse Gases (CO2, N2O, CH4)in a Changing Climate

[EN] Understanding and quantifying ocean-atmosphere exchanges of the long-lived greenhouse gases carbon dioxide (CO2), nitrous oxide (N2O) and methane (CH4) are important for understanding the global biogeochemical cycles of carbon and nitrogen in the context of ongoing global climate change. In this chapter we summarise our current state of knowledge regarding the oceanic distributions, formation and consumption pathways, and oceanic uptake and emissions of CO2, N2O and CH4, with a particular emphasis on the upper ocean. We specifically consider the role of the ocean in regulating the tropospheric content of these important radiative gases in a world in which their tropospheric content is rapidly increasing and estimate the impact of global change on their present and future oceanic uptake and/or emission. Finally, we evaluate the various uncertainties associated with the most commonly used methods for estimating uptake and emission and identify future research needs.

A study of air-sea interaction processes on water mass formation and upwelling in the Mediterranean sea

Air-sea interactions are a key process in the forcing of the ocean circulation and the climate. Water Mass Formation is a phenomenon related to extreme air-sea exchanges and heavy heat losses by the water column, being capable to transfer water properties from the surface to great depth and constituting a fundamental component of the thermohaline circulation of the ocean. Wind-driven Coastal Upwelling, on the other hand, is capable to induce intense heat gain in the water column, making this phenomenon important for climate change; further, it can have a noticeable influence on many biological pelagic ecosystems mechanisms. To study some of the fundamental characteristics of Water Mass Formation and Coastal Upwelling phenomena in the Mediterranean Sea, physical reanalysis obtained from the Mediterranean Forecating System model have been used for the period ranging from 1987 to 2012. The first chapter of this dissertation gives the basic description of the Mediterranean Sea circulation, the MFS model implementation, and the air-sea interaction physics. In the second chapter, the problem of Water Mass Formation in the Mediterranean Sea is approached, also performing ad-hoc numerical simulations to study heat balance components. The third chapter considers the study of Mediterranean Coastal Upwelling in some particular areas (Sicily, Gulf of Lion, Aegean Sea) of the Mediterranean Basin, together with the introduction of a new Upwelling Index to characterize and predict upwelling features using only surface estimates of air-sea fluxes. Our conclusions are that latent heat flux is the driving air-sea heat balance component in the Water Mass Formation phenomenon, while sensible heat exchanges are fundamental in Coastal Upwelling process. It is shown that our upwelling index is capable to reproduce the vertical velocity patterns in Coastal Upwelling areas. Nondimensional Marshall numbers evaluations for the open-ocean convection process in the Gulf of Lion show that it is a fully turbulent, three-dimensional phenomenon.

Dynamics of Bora wind over the Adriatic sea: atmospheric water balance and role of air-sea fluxes and orography

The Bora wind is a mesoscale phenomenon which typically affects the Adriatic Sea basin for several days each year, especially during winter. The Bora wind has been studied for its intense outbreak across the Dinaric Alps. The properties of the Bora wind are widely discussed in the literature and scientific papers usually focus on the eastern Adriatic coast where strong turbulence and severe gust intensity are more pronounced. However, the impact of the Bora wind can be significant also over Italy, not only in terms of wind speed instensity. Depending on the synoptic pressure pattern (cyclonic or anticyclonic Bora) and on the season, heavy snowfall, severe storms, storm surges and floods can occur along the Adriatic coast and on the windward flanks of the Apennines. In the present work five Bora cases that occurred in recent years have been selected and their evolution has been simulated with the BOLAM-MOLOCH model set, developed at ISAC-CNR in Bologna. Each case study has been addressed by a control run and by several sensitivity tests, performed with the purpose of better understanding the role played by air-sea latent and sensible heat fluxes. The tests show that the removal of the fluxes induces modifications in the wind approching the coast and a decrease of the total precipitation amount predicted over Italy. In order to assess the role of heat fluxes, further analysis has been carried out: column integrated water vapour fluxes have been computed along the Italian coastline and an atmospheric water balance has been evaluated inside a box volume over the Adriatic Sea. The balance computation shows that, although latent heat flux produces a significant impact on the precipitation field, its contribution to the balance is relatively minor. The most significant and lasting case study, that of February 2012, has been studied in more detail in order to explain the impressive drop in the total precipitation amount simulated in the sensitivity tests with removed heat fluxes with respect to the CNTRL run. In these experiments relative humidity and potential temperature distribution over different cross-sections have been examined. With respect to the CNTRL run a drier and more stable boundary layer, characterised by a more pronounced wind shear at the lower levels, has been observed to establish above the Adriatic Sea. Finally, in order to demonstrate that also the interaction of the Bora flow with the Apennines plays a crucial role, sensitivity tests varying the orography height have been considered. The results of such sensitivity tests indicate that the propagation of the Bora wind over the Adriatic Sea, and in turn its meteorological impact over Italy, is influenced by both the large air-sea heat fluxes and the interaction with the Apennines that decelerate the upstream flow.

Air-Sea Interactions During the Passage of a Winter Storm Over the Gulf Stream: A Three-Dimensional Coupled Atmosphere-Ocean Model Study

A three-dimensional, regional coupled atmosphere-ocean model with full physics is developed to study air-sea interactions during winter storms off the U. S. east coast. Because of the scarcity of open ocean observations, models such as this offer valuable opportunities to investigate how oceanic forcing drives atmospheric circulation and vice versa. The study presented here considers conditions of strong atmospheric forcing (high wind speeds) and strong oceanic forcing (significant sea surface temperature (SST) gradients). A simulated atmospheric cyclone evolves in a manner consistent with Eta reanalysis, and the simulated air-sea heat and momentum exchanges strongly affect the circulations in both the atmosphere and the ocean. For the simulated cyclone of 19-20 January 1998, maximum ocean-to-atmosphere heat fluxes first appear over the Gulf Stream in the South Atlantic Bight, and this results in rapid deepening of the cyclone off the Carolina coast. As the cyclone moves eastward, the heat flux maximum shifts into the region near Cape Hatteras and later northeast of Hatteras, where it enhances the wind locally. The oceanic response to the atmospheric forcing is closely related to the wind direction. Southerly and southwesterly winds tend to strengthen surface currents in the Gulf Stream, whereas northeasterly winds weaken the surface currents in the Gulf Stream and generate southwestward flows on the shelf. The oceanic feedback to the atmosphere moderates the cyclone strength. Compared with a simulation in which the oceanic model always passes the initial SST to the atmospheric model, the coupled simulation in which the oceanic model passes the evolving SST to the atmospheric model produces higher ocean-to-atmosphere heat flux near Gulf Stream meander troughs. This is due to wind-driven lateral shifts of the stream, which in turn enhance the local northeasterly winds. Away from the Gulf Stream the coupled simulation produces surface winds that are 5 similar to 10% weaker. Differences in the surface ocean currents between these two experiments are significant on the shelf and in the open ocean.

Chemical characterization of dissolved organic compounds from coastal sea surface microlayers (Baltic Sea, Germany)

The physicochemical properties of the sea surface microlayer (SML), i.e. the boundary layer between the air and the sea, and its impact on air-sea exchange processes have been investigated for decades. However, a detailed description about these processes remains incomplete. In order to obtain a better chemical characterization of the SML, in a case study three pairs of SML and corresponding bulk water samples were taken in the southern Baltic Sea. The samples were analyzed for dissolved organic carbon and dissolved total nitrogen, as well as for several organic nitrogen containing compounds and carbohydrates, namely aliphatic amines, dissolved free amino acids, dissolved free monosaccharides, sugar alcohols, and monosaccharide anhydrates. Therefore, reasonable analytical procedures with respect to desalting and enrichment were established. All aliphatic amines and the majority of the investigated amino acids (11 out of 18) were found in the samples with average concentrations between 53 ng/l and 1574 ng/l. The concentrations of carbohydrates were slightly higher, averaging 2900 ng/l. Calculation of the enrichment factor (EF) between the sea surface microlayer and the bulk water showed that dissolved total nitrogen was more enriched (EF: 1.1 and 1.2) in the SML than dissolved organic carbon (EF: 1.0 and 1.1). The nitrogen containing organic compounds were generally found to be enriched in the SML (EF: 1.9-9.2), whereas dissolved carbohydrates were not enriched or even depleted (EF: 0.7-1.2). Although the investigated compounds contributed on average only 0.3% to the dissolved organic carbon and 0.4% to the total dissolved nitrogen fraction, these results underline the importance of single compound analysis to determine SML structure, function, and its potential for a transfer of compounds into the atmosphere.

Sea-suface microlayer (SML) and underlying water (ULW) samples focusing on CDOM spectral characteristics during METEOR cruise M91

The coastal upwelling system off the coast of Peru is characterized by high biological activity and a pronounced subsurface oxygen minimum zone, as well as associated emissions of atmospheric trace gases such as N2O, CH4 and CO2. From 3 to 23 December 2012, R/V Meteor (M91) cruise took place in the Peruvian upwelling system between 4.59 and 15.4°S, and 82.0 to 77.5°W. During M91 we investigated the composition of the sea-surface microlayer (SML), the oceanic uppermost boundary directly subject to high solar radiation, often enriched in specific organic compounds of biological origin like chromophoric dissolved organic matter (CDOM) and marine gels. In the SML, the continuous photochemical and microbial recycling of organic matter may strongly influence gas exchange between marine systems and the atmosphere. We analyzed SML and underlying water (ULW) samples at 38 stations focusing on CDOM spectral characteristics as indicator of photochemical and microbial alteration processes. CDOM composition was characterized by spectral slope (S) values and excitation-emission matrix fluorescence (EEMs), which allow us to track changes in molecular weight (MW) of DOM, and to determine potential DOM sources and sinks. Spectral slope S varied between 0.012 to 0.043 1 nm-1 and was quite similar between SML and ULW, with no significant differences between the two compartments. Higher S values were observed in the ULW of the southern stations below 15°S. By EEMs, we identified five fluorescent components (F1-5) of the CDOM pool, of which two had excitation/emission characteristics of amino-acid-like fluorophores (F1, F4) and were highly enriched in the SML, with a median ratio SML : ULW of 1.5 for both fluorophores. In the study region, values for CDOM absorption ranged from 0.07 to 1.47 m-1. CDOM was generally highly concentrated in the SML, with a median enrichment with respect to the ULW of 1.2. CDOM composition and changes in spectral slope properties suggested a local microbial release of DOM directly in the SML as a response to light exposure in this extreme environment. In a conceptual model of the sources and modifications of optically active DOM in the SML and underlying seawater (ULW), we describe processes we think may take place (Fig. 1); the production of CDOM of higher MW by microbial release through growth, exudation and lysis in the euphotic zone, includes the identified fluorophores (F1, F2, F3, F4, F5). Specific amino-acid-like fluorophores (F1, F4) accumulate in the SML with respect to the ULW, as photochemistry may enhance microbial CDOM release by (a) photoprotection mechanisms and (b) cell-lysis processes. Microbial and photochemical degradation are potential sinks of the amino-acid-like fluorophores (F1, F4), and potential sources of reworked and more refractory humic-like components (F2, F3, F5). In the highly productive upwelling region along the Peruvian coast, the interplay of microbial and photochemical processes controls the enrichment of amino-acid-like CDOM in the SML. We discuss potential implications for air-sea gas exchange in this area.

Global air-sea flux of CO2: An estimate based on measurements of sea–air pCO2 difference

Approximately 250,000 measurements made for the pCO2 difference between surface water and the marine atmosphere, ΔpCO2, have been assembled for the global oceans. Observations made in the equatorial Pacific during El Nino events have been excluded from the data set. These observations are mapped on the global 4° × 5° grid for a single virtual calendar year (chosen arbitrarily to be 1990) representing a non-El Nino year. Monthly global distributions of ΔpCO2 have been constructed using an interpolation method based on a lateral advection–diffusion transport equation. The net flux of CO2 across the sea surface has been computed using ΔpCO2 distributions and CO2 gas transfer coefficients across sea surface. The annual net uptake flux of CO2 by the global oceans thus estimated ranges from 0.60 to 1.34 Gt-C⋅yr−1 depending on different formulations used for wind speed dependence on the gas transfer coefficient. These estimates are subject to an error of up to 75% resulting from the numerical interpolation method used to estimate the distribution of ΔpCO2 over the global oceans. Temperate and polar oceans of the both hemispheres are the major sinks for atmospheric CO2, whereas the equatorial oceans are the major sources for CO2. The Atlantic Ocean is the most important CO2 sink, providing about 60% of the global ocean uptake, while the Pacific Ocean is neutral because of its equatorial source flux being balanced by the sink flux of the temperate oceans. The Indian and Southern Oceans take up about 20% each.