Reconstruction of the Gulf Stream variability since 1940 using a variational inverse method and study of its interaction with the North Atlantic Oscillation

In this study, the Gulf Stream’s (GS) response to the North Atlantic oscillation (NAO) is investigated by generating an observation-based reconstruction of the GS path between 70° and 50°W since 1940. Using in situ data from WOD, SeaDataNet, ICES, Hydrobase3 and ARGO floats, a harmonized database of more than 40 million entries is created. A variational inverse method implemented in the software DIVA (Data-Interpolating Variational Analysis) allows the production of time series of monthly analyses of temperature and salinity over the North Atlantic (NA). These time series are used to derive two GS indices: the GS North Wall (GSNW) index for position and the GS Delta (GSD) index as a proxy of its transport. We find a significant correlation (0.37) between the GSNW and the NAO at a lag of 1 year (NAO preceding GS) since 1940 and significant correlations (0.50 and 0.43) between the GSD and the NAO at lags of 0 and 2 years between 1960–2014. We suggest this 2-year lag is due to Rossby waves, generated by NAO variability, that propagate westwards from the center of the NA. This is the first reconstruction of GS indices over a 75-year period based on an objective method using the largest in situ dataset so far. This enhanced tracking and quantification of the GS confirms and extends the temporal scope of this property: NAO+ phases lead to a stronger and more northward GS, and conversely for NAO− phases. The teleconnections between the NAO and the variability of the GS were extensively studied these last years, often exhibiting time delays between both phenomena. These time lags, usually ranging between 0–2 years, are sometimes explained by the hypothesis of baroclinic Rossby waves generated by the NAO in the central NA and travelling westward before interacting with the GS. In this study, we use a numerical hindcast at an eddy-resolving resolution (1/12°) from the DRAKKAR project to examine the occurrence and properties of such Rossby waves between 1970–2015, thus including a large pre-TOPEX/Poseidon period. Through the use of a two-dimensional Radon Transform (2D-RT) on Hovmöller diagrams of the Sea Surface Height (SSH), a methodology easily portable to other oceanic model outputs, we show evidence of baroclinic Rossby waves travelling at 39°N at a speed of 4.17 cm/s. This study extends the period over which Rossby waves have been found that far north to a much longer period, which reinforces the findings of previous works. These results are consistent with the time lags observed between the NAO and the GS transport while the GS latitudinal shifts might obey additional processes. The Barents Sea, located between the Norwegian Sea and the Arctic Ocean, is one of the main pathways of the Atlantic Meridional Overturning Circulation. Changes in the water mass transformations in the Barents Sea potentially affect the thermohaline circulation through the alteration of the dense water formation process. In order to investigate such changes, we present here a seasonal atlas of the Barents Sea including both temperature and salinity for the period 1965–2016. The atlas is built as a compilation of datasets from the World Ocean Database, the Polar Branch of Russian Federal Research Institute of Fisheries and Oceanography, and the Norwegian Polar Institute using the DIVA tool. DIVA allows for a minimization of the expected error variance with respect to the true field. The atlas is used to provide a volumetric analysis of water mass characteristics and an estimation of the ocean heat and freshwater contents. The results show a recent "Atlantification" of the Barents Sea, that is a general increase of both temperature and salinity, while its density remains stable. The atlas is made freely accessible as handy NetCDF files to encourage further research in the Barents Sea physics.SeaDataNet I

Web-based visualization of gridded dataset usings OceanBrowser

OceanBrowser is a web-based visualization tool for gridded oceanographic data sets. Those data sets are typically four-dimensional (longitude, latitude, depth and time). OceanBrowser allows one to visualize horizontal sections at a given depth and time to examine the horizontal distribution of a given variable. It also offers the possibility to display the results on an arbitrary vertical section. To study the evolution of the variable in time, the horizontal and vertical sections can also be animated. Vertical section can be generated by using a fixed distance from coast or fixed ocean depth. The user can customize the plot by changing the color-map, the range of the color-bar, the type of the plot (linearly interpolated color, simple contours, filled contours) and download the current view as a simple image or as Keyhole Markup Language (KML) file for visualization in applications such as Google Earth. The data products can also be accessed as NetCDF files and through OPeNDAP. Third-party layers from a web map service can also be integrated. OceanBrowser is used in the frame of the SeaDataNet project (http://gher-diva.phys.ulg.ac.be/web-vis/) and EMODNET Chemistry (http://oceanbrowser.net/emodnet/) to distribute gridded data sets interpolated from in situ observation using DIVA (Data-Interpolating Variational Analysis)

A volumetric census of the Barents Sea in a changing climate

The Barents Sea, located between the Norwegian Sea and the Arctic Ocean, is one of the main pathways of the Atlantic Meridional Overturning Circulation. Changes in the water mass transformations in the Barents Sea potentially affect the thermohaline circulation through the alteration of the dense water formation process. In order to investigate such changes, we present here a seasonal atlas of the Barents Sea including both temperature and salinity for the period 1965–2016. The atlas is built as a compilation of datasets from the World Ocean Database, the Polar Branch of the Russian Federal Research Institute of Fisheries and Oceanography and the Norwegian Polar Institute using the Data-Interpolating Variational Analysis (DIVA) tool. DIVA allows for a minimization of the expected error with respect to the true field. The atlas is used to provide a volumetric analysis of water mass characteristics and an estimation of the ocean heat and freshwater contents. The results show a recent “Atlantification” of the Barents Sea, that is a general increase in both temperature and salinity, while its density remains stable. The atlas is made freely accessible as user-friendly NetCDF files to encourage further research in the Barents Sea physics (https://doi.org/10.21335/NMDC-2058021735, Watelet et al., 2020).publishedVersio

Web-based workflows to produce ocean climatologies using DIVA (Data-Interpolating Variational Analysis) and Jupyter notebooks

The analysis tool DIVA (Data-Interpolating Variational Analysis) is designed to generate gridded fields or climatologies from in situ observations. The tool DIVA minimizes a cost function to ensure that the analysed field is relatively close to the observations and conforms at the same time to a set of dynamical constraints. In particular, DIVA naturally decouples water bodies which are not directly connected and it uses a (potentially spatial varying) correlation length to describe over which length-scale the analysed variable is correlated. In addition, DIVA can also take ocean currents into account to introduce a preferential direction for the correlation. The SeaDataCloud project aims to facilitate the access and use of ocean in situ data from 45 national oceanographic data centres and marine data centres from 35 countries riparian to all European seas. A central aspect is to provide web-based virtual research environment, where scientists can easily access and explore the data sets through the SeaDataCloud infrastructure. For users familiar with programming languages like Julia and Python, Jupyter (acronym for Julia, Python and R) notebooks provide an exciting way to analyse and to interact with ocean data. Jupyter notebooks are made up of cells that can be run individually and can contain text, formulas or code fragment. A complete notebook explains how to go from input data and parameters to a result, in this case a gridded field obtained executing DIVA. This presentation discusses this new web-based workflow for generating climatologies using DIVA. It explores its new possibilities in particular, in terms of improved ease of use and reproducibility of the results. The integration in the infrastructure of EUDAT is also addressed

A new global interior ocean mapped climatology: the 1° × 1° GLODAP version 2

We present a mapped climatology (GLODAPv2.2016b) of ocean biogeochemical variables based on the new GLODAP version 2 data product (Olsen et al., 2016; Key et al., 2015), which covers all ocean basins over the years 1972 to 2013. The quality-controlled and internally consistent GLODAPv2 was used to create global 1°  ×  1° mapped climatologies of salinity, temperature, oxygen, nitrate, phosphate, silicate, total dissolved inorganic carbon (TCO2), total alkalinity (TAlk), pH, and CaCO3 saturation states using the Data-Interpolating Variational Analysis (DIVA) mapping method. Improving on maps based on an earlier but similar dataset, GLODAPv1.1, this climatology also covers the Arctic Ocean. Climatologies were created for 33 standard depth surfaces. The conceivably confounding temporal trends in TCO2 and pH due to anthropogenic influence were removed prior to mapping by normalizing these data to the year 2002 using first-order calculations of anthropogenic carbon accumulation rates. We additionally provide maps of accumulated anthropogenic carbon in the year 2002 and of preindustrial TCO2. For all parameters, all data from the full 1972–2013 period were used, including data that did not receive full secondary quality control. The GLODAPv2.2016b global 1°  ×  1° mapped climatologies, including error fields and ancillary information, are available at the GLODAPv2 web page at the Carbon Dioxide Information Analysis Center (CDIAC; doi:10.3334/CDIAC/OTG.NDP093_GLODAPv2)

Reconstruction of the Gulf Stream from 1940 to the present and correlation with the North Atlantic Oscillation

AbstractIn this study, the Gulf Stream (GS)’s response to the North Atlantic Oscillation (NAO) is investigated by generating an observation-based reconstruction of the GS path between 70° and 50°W since 1940. Using in situ data from the World Ocean Database (WOD), SeaDataNet, International Council for the Exploration of the Sea (ICES), Hydrobase3, and Argo floats, a harmonized database of more than 40 million entries is created. A variational inverse method implemented in the software Data Interpolating Variational Analysis (DIVA) allows the production of time series of monthly analyses of temperature and salinity over the North Atlantic (NA). These time series are used to derive two GS indices: the GS north wall (GSNW) index for position and the GS delta (GSD) index as a proxy of its transport. This study finds a significant correlation (0.37) between the GSNW and the NAO at a lag of 1 year (NAO preceding GS) since 1940 and significant correlations (0.50 and 0.43) between the GSD and the NAO at lags of 0 and 2 years between 1960 and 2014. The authors suggest this 2-yr lag is due to Rossby waves, generated by NAO variability, that propagate westward from the center of the NA. This is the first reconstruction of GS indices over a 75-yr period based on an objective method using the largest in situ dataset so far.</jats:p

Software citation in ocean sciences: towards the reproducibility of results

ODI