Radon 222 and Radium 226 measured on water bottle samples during POLARSTERN cruise ARK-XXVI/3 (TransArc)

Air-sea gas exchange plays a key role in the cycling of greenhouse and other biogeochemically important gases. Although air-sea gas transfer is expected to change as a consequence of the rapid decline in summer Arctic sea ice cover, little is known about the effect of sea ice cover on gas exchange fluxes, especially in the marginal ice zone. During the Polarstern expedition ARK-XXVI/3 (TransArc, August/September 2011) to the central Arctic Ocean, we compared 222Rn/226Ra ratios in the upper 50 m of 14 ice-covered and 4 ice-free stations. At three of the ice-free stations, we find 222Rn-based gas transfer coefficients in good agreement with expectation based on published relationships between gas transfer and wind speed over open water when accounting for wind history from wind reanalysis data. We hypothesize that the low gas transfer rate at the fourth station results from reduced fetch due to the proximity of the ice edge, or lateral exchange across the front at the ice edge by restratification. No significant radon deficit could be observed at the ice-covered stations. At these stations, the average gas transfer velocity was less than 0.1 m/d (97.5% confidence), compared to 0.5-2.2 m/d expected for open water. Our results show that air-sea gas exchange in an ice-covered ocean is reduced by at least an order of magnitude compared to open water. In contrast to previous studies, we show that in partially ice-covered regions, gas exchange is lower than expected based on a linear scaling to percent ice cover.

Seawater carbonate chemistry and calcification during tropical reef studies at Moorea (French Polynesia), 1996

Community metabolism was investigated using a Lagrangian flow respirometry technique on 2 reef flats at Moorea (French Polynesia) during austral winter and Yonge Reef (Great Barrier Reef) during austral summer. The data were used to estimate related air-sea CO2 disequilibrium. A sine function did not satisfactorily model the diel light curves and overestimated the metabolic parameters. The ranges of community gross primary production and respiration (Pg and R; 9 to 15 g C m-2 d-1) were within the range previously reported for reef flats, and community net calcification (G; 19 to 25 g CaCO3 m-2 d-1) was higher than the 'standard' range. The molar ratio of organic to inorganic carbon uptake was 6:1 for both sites. The reef flat at Moorea displayed a higher rate of organic production and a lower rate of calcification compared to previous measurements carried out during austral summer. The approximate uncertainty of the daily metabolic parameters was estimated using a procedure based on a Monte Carlo simulation. The standard errors of Pg,R and Pg/R expressed as a percentage of the mean are lower than 3% but are comparatively larger for E, the excess production (6 to 78%). The daily air-sea CO2 flux (FCO2) was positive throughout the field experiments, indicating that the reef flats at Moorea and Yonge Reef released CO2 to the atmosphere at the time of measurement. FCO2 decreased as a function of increasing daily irradiance.

Seawater carbonate chemistry and calcification during tropical reef studies at Younge reef (Great Barrier Reef), 1996

Community metabolism was investigated using a Lagrangian flow respirometry technique on 2 reef flats at Moorea (French Polynesia) during austral winter and Yonge Reef (Great Barrier Reef) during austral summer. The data were used to estimate related air-sea CO2 disequilibrium. A sine function did not satisfactorily model the diel light curves and overestimated the metabolic parameters. The ranges of community gross primary production and respiration (Pg and R; 9 to 15 g C m-2 d-1) were within the range previously reported for reef flats, and community net calcification (G; 19 to 25 g CaCO3 m-2 d-1) was higher than the 'standard' range. The molar ratio of organic to inorganic carbon uptake was 6:1 for both sites. The reef flat at Moorea displayed a higher rate of organic production and a lower rate of calcification compared to previous measurements carried out during austral summer. The approximate uncertainty of the daily metabolic parameters was estimated using a procedure based on a Monte Carlo simulation. The standard errors of Pg,R and Pg/R expressed as a percentage of the mean are lower than 3% but are comparatively larger for E, the excess production (6 to 78%). The daily air-sea CO2 flux (FCO2) was positive throughout the field experiments, indicating that the reef flats at Moorea and Yonge Reef released CO2 to the atmosphere at the time of measurement. FCO2 decreased as a function of increasing daily irradiance.

The South Atlantic carbon isotope record of planktic foraminifera

We reviewed the paleoceanographic application of the carbon isotope composition of planktic foraminifera. Major controls on the distribution of d13C of dissolved CO2 (d13CSCO2) in the modern ocean are photosynthesis-respiration cycle, isotopic fractionation during air-sea exchange, and circulation. The carbon isotope composition of surface waters is not recorded without perturbations by planktic foraminifera. Besides d13CSCO2 of the surrounding seawater, the d13C composition of planktic foraminifera is affected by vital effects, the water depth of calcification and postdepositional dissolution. We compared several high-resolution (>10cm/ka) carbon isotope records from the Southern Ocean, the Benguela upwelling system, and the tropical Atlantic. In the Southern Ocean, carbon isotope values are about 1.2 per mil lower during the LGM and up to 1.7 per mil lower during the last deglaciation, when compared to the Holocene. These depletions might be explained with a combination of a subsurface nutrient enrichment and reduced air-sea exchange due to an increased stratification of surface waters. In the Benguela Upwelling system, waters originating in the south are upwelled. While the deglacial minimum is transferred and recorded in its full extent in the d13C record of Globigerina bulloides, glacial values show only little changes. This might suggest, that the lower glacial d13C values of high-latitude surface waters are not upwelled off Namibia, or that G. bulloides records post-upwelling conditions, when increased seasonal production has already increased surface-water d13C. Synchronous to the d13C depletions in high latitudes, low d13C values were recorded in Globigerinoides sacculifer during the LGM and during the last deglaciation in the nutrient-depleted western equatorial Atlantic. Hence, part of the glacial-interglacial variability presumably transferred from high to low latitudes seems to be related to changes in thermodynamic fractionation. The variability in d13C is lowest in the northernmost core M35003-4 from the eastern Caribbean, implying that the Antarctic Intermediate Water might have acted as a conduit to transfer the deglacial minimum to tropical surface waters.

Underway measurements of dissolved gases along POLARSTERN cruise track ANT-XXVII/2

We present measurements of pCO2, O2 concentration, biological oxygen saturation (Delta O2/Ar) and N2 saturation (Delta N2) in Southern Ocean surface waters during austral summer, 2010-2011. Phytoplankton biomass varied strongly across distinct hydrographic zones, with high chlorophyll a (Chla) concentrations in regions of frontal mixing and sea-ice melt. pCO2 and Delta O2 /Ar exhibited large spatial gradients (range 90 to 450 µatm and -10 to 60%, respectively) and co-varied strongly with Chla. However, the ratio of biological O2 accumulation to dissolved inorganic carbon (DIC) drawdown was significantly lower than expected from photosynthetic stoichiometry, reflecting the differential time-scales of O2 and CO2 air-sea equilibration. We measured significant oceanic CO2 uptake, with a mean air-sea flux (~ -20 mmol m-2 d-1) that significantly exceeded regional climatological values. N2 was mostly supersaturated in surface waters (mean Delta N2 of +2.5 %), while physical processes resulted in both supersaturation and undersaturation of mixed layer O2 (mean Delta O2phys = 2.1 %). Box model calculations were able to reproduce much of the spatial variability of Delta N2 and Delta O2phys along the cruise track, demonstrating significant effects of air-sea exchange processes (e.g. atmospheric pressure changes and bubble injection) and mixed layer entrainment on surface gas disequilibria. Net community production (NCP) derived from entrainment-corrected surface Delta O2 /Ar data, ranged from ~ -40 to > 300 mmol O2 m-2 d-1 and showed good coherence with independent NCP estimates based on seasonal mixed layer DIC deficits. Elevated NCP was observed in hydrographic frontal zones and regions of sea-ice melt with shallow mixed layer depths, reflecting the importance of mixing in controlling surface water light and nutrient availability.

Breaking Waves on the Ocean Surface

Thesis (Ph.D.)--University of Washington, 2016-08

FluxEngine: A flexible processing system for calculating atmosphere-ocean carbon dioxide gas fluxes and climatologies

The air-sea flux of greenhouse gases (e.g. carbon dioxide, CO2) is a critical part of the climate system and a major factor in the biogeochemical development of the oceans. More accurate and higher resolution calculations of these gas fluxes are required if we are to fully understand and predict our future climate. Satellite Earth observation is able to provide large spatial scale datasets that can be used to study gas fluxes. However, the large storage requirements needed to host such data can restrict its use by the scientific community. Fortunately, the development of cloud-computing can provide a solution. Here we describe an open source air-sea CO2 flux processing toolbox called the ‘FluxEngine’, designed for use on a cloud-computing infrastructure. The toolbox allows users to easily generate global and regional air-sea CO2 flux data from model, in situ and Earth observation data, and its air-sea gas flux calculation is user configurable. Its current installation on the Nephalae cloud allows users to easily exploit more than 8 terabytes of climate-quality Earth observation data for the derivation of gas fluxes. The resultant NetCDF data output files contain >20 data layers containing the various stages of the flux calculation along with process indicator layers to aid interpretation of the data. This paper describes the toolbox design, the verification of the air-sea CO2 flux calculations, demonstrates the use of the tools for studying global and shelf-sea air-sea fluxes and describes future developments.

Dissolved inorganic carbon budgets in the eastern subpolar North Atlantic in the 2000s from in situ data

The subpolar North Atlantic (SPNA) is important in the global carbon cycle because of the deep water ventilation processes that lead to both high uptake of atmospheric CO2 and large inventories of anthropogenic CO2 (C-ant). Thus, it is crucial to understand its response to increasing anthropogenic pressures. In this work, the budgets of dissolved inorganic carbon (DIC), C-ant and natural DIC (DICnat) in the eastern SPNA in the 2000s, are jointly analyzed using in situ data. The DICnat budget is found to be in steady state, confirming a long-standing hypothesis from in situ data for the first time. The biological activity is driving the uptake of natural CO2 from the atmosphere. The C-ant increase in the ocean is solely responsible of the DIC storage rate which is explained by advection of C-ant from the subtropics (65%) and C-ant air-sea flux (35%). These results demonstrate that the C-ant is accumulating in the SPNA without affecting the natural carbon cycle.

Tracking the attenuation and nonbreaking dissipation of swells using altimeters

A method for systematically tracking swells across oceanic basins is developed by taking advantage of high-quality data from space-borne altimeters and wave model output. The evolution of swells is observed over large distances based on 202 swell events with periods ranging from 12 to 18 s. An empirical attenuation rate of swell energy of about 4 × 10−7 m−1 is estimated using these observations, and the nonbreaking energy dissipation rates of swells far away from their generating areas are also estimated using a point source model. The resulting acceptance range of nonbreaking dissipation rates is −2.5 to 5.0 × 10−7 m−1, which corresponds to a dissipation e-folding scales of at least 2000 km for steep swells, to almost infinite for small-amplitude swells. These resulting rates are consistent with previous studies using in-situ and synthetic aperture radar (SAR) observations. The frequency dispersion and angular spreading effects during swell propagation are discussed by comparing the results with other studies, demonstrating that they are the two dominant processes for swell height attenuation, especially in the near field. The resulting dissipation rates from these observations can be used as a reference for ocean engineering and wave modeling, and for related studies such as air-sea and wind-wave-turbulence interactions.

Swell dissipation from 10 years of Envisat ASAR in wave mode

Swells are found in all oceans and strongly influence the wave climate and air-sea processes. The poorly known swell dissipation is the largest source of error in wave forecasts and hindcasts. We use synthetic aperture radar data to identify swell sources and trajectories, allowing a statistically significant estimation of swell dissipation. We mined the entire Envisat mission 2003–2012 to find suitable storms with swells (13 < T < 18 s) that are observed several times along their propagation. This database of swell events provides a comprehensive view of swell extending previous efforts. The analysis reveals that swell dissipation weakly correlates with the wave steepness, wind speed, orbital wave velocity, and the relative direction of wind and waves. Although several negative dissipation rates are found, there are uncertainties in the synthetic aperture radar-derived swell heights and dissipation rates. An acceptable range of the swell dissipation rate is −0.1 to 6 × 10−7 m−1 with a median of 1 × 10−7 m−1.

Oceanic control of multidecadal variability in an idealized coupled GCM

Idealized ocean models are known to develop intrinsic multidecadal oscillations of the meridional overturning circulation (MOC). Here we explore the role of ocean–atmosphere interactions on this low-frequency variability. We use a coupled ocean–atmosphere model set up in a flat-bottom aquaplanet geometry with two meridional boundaries. The model is run at three different horizontal resolutions (4°, 2° and 1°) in both the ocean and atmosphere. At all resolutions, the MOC exhibits spontaneous variability on multidecadal timescales in the range 30–40 years, associated with the propagation of large-scale baroclinic Rossby waves across the Atlantic-like basin. The unstable region of growth of these waves through the long wave limit of baroclinic instability shifts from the eastern boundary at coarse resolution to the western boundary at higher resolution. Increasing the horizontal resolution enhances both intrinsic atmospheric variability and ocean–atmosphere interactions. In particular, the simulated atmospheric annular mode becomes significantly correlated to the MOC variability at 1° resolution. An ocean-only simulation conducted for this specific case underscores the disruptive but not essential influence of air–sea interactions on the low-frequency variability. This study demonstrates that an atmospheric annular mode leading MOC changes by about 2 years (as found at 1° resolution) does not imply that the low-frequency variability originates from air–sea interactions.

The EMSO-ERIC Pan-European Consortium: Data Benefits and Lessons Learned as the Legal Entity Forms

The European Multidisciplinary Seafloor and water-column Observatory (EMSO) European Research Infrastructure Consortium (ERIC) provides power, communications, sensors, and data infrastructure for continuous, high-resolution, (near-)real-time, interactive ocean observations across a multidisciplinary and interdisciplinary range of research areas including biology, geology, chemistry, physics, engineering, and computer science, from polar to subtropical environments, through the water column down to the abyss. Eleven deep-sea and four shallow nodes span from the Arctic through the Atlantic and Mediterranean, to the Black Sea. Coordination among the consortium nodes is being strengthened through the EMSOdev project (H2020), which will produce the EMSO Generic Instrument Module (EGIM). Early installations are now being upgraded, for example, at the Ligurian, Ionian, Azores, and Porcupine Abyssal Plain (PAP) nodes. Significant findings have been flowing in over the years; for example, high-frequency surface and subsurface water-column measurements of the PAP node show an increase in seawater pCO2 (from 339 μatm in 2003 to 353 μatm in 2011) with little variability in the mean air-sea CO2 flux. In the Central Eastern Atlantic, the Oceanic Platform of the Canary Islands open-ocean canary node (aka ESTOC station) has a long-standing time series on water column physical, biogeochemical, and acidification processes that have contributed to the assessment efforts of the Intergovernmental Panel on Climate Change (IPCC). EMSO not only brings together countries and disciplines but also allows the pooling of resources and coordination to assemble harmonized data into a comprehensive regional ocean picture, which will then be made available to researchers and stakeholders worldwide on an open and interoperable access basis.

First day of an oil spill on the open sea: early mass transfers of hydrocarbons to air and water.

During the first hours after release of petroleum at sea, crude oil hydrocarbons partition rapidly into air and water. However, limited information is available about very early evaporation and dissolution processes. We report on the composition of the oil slick during the first day after a permitted, unrestrained 4.3 m(3) oil release conducted on the North Sea. Rapid mass transfers of volatile and soluble hydrocarbons were observed, with >50% of ≤C17 hydrocarbons disappearing within 25 h from this oil slick of <10 km(2) area and <10 μm thickness. For oil sheen, >50% losses of ≤C16 hydrocarbons were observed after 1 h. We developed a mass transfer model to describe the evolution of oil slick chemical composition and water column hydrocarbon concentrations. The model was parametrized based on environmental conditions and hydrocarbon partitioning properties estimated from comprehensive two-dimensional gas chromatography (GC×GC) retention data. The model correctly predicted the observed fractionation of petroleum hydrocarbons in the oil slick resulting from evaporation and dissolution. This is the first report on the broad-spectrum compositional changes in oil during the first day of a spill at the sea surface. Expected outcomes under other environmental conditions are discussed, as well as comparisons to other models.