British Oceanographic Data Centre

Abstract. Many scientists have made use of the wavelet method in analyzing time series, often using popular free software. However, at present there are no similar easy to use wavelet packages for analyzing two time series together. We discuss the cross wavelet transform and wavelet coherence for examining relationships in time frequency space between two time series. We demonstrate how phase angle statistics can be used to gain confidence in causal relationships and test mechanistic models of physical relationships between the time series. As an example of typical data where such analyses have proven useful, we apply the methods to the Arctic Oscillation index and the Baltic maximum sea ice extent record. Monte Carlo methods are used to assess the statistical significance against red noise backgrounds. A software package has been developed that allows users to perform the cross wavelet transform and wavelet coherence (www.pol.ac.uk/home/research/waveletcoherence/).

A new 30-arc second resolution global topography/bathymetry grid (SRTM30_PLUS) has been developed from a wide variety of data sources. Land and ice topography comes from the SRTM30 and ICESat topography, respectively. Ocean bathymetry is based on a new satellite-gravity model where the gravity-to-topography ratio is calibrated using 298 million edited soundings. The main contribution of this study is the compilation and editing of the raw soundings, which come from NOAA, individual scientists, SIO, NGA, JAMSTEC, IFREMER, GEBCO, and NAVOCEANO. The gridded bathymetry is available for ftp download in the same format as the 33 tiles of SRTM30 topography. There are 33 matching tiles of source identification number to convey the provenance of every grid cell. The raw sounding data, converted to a simple common format, are also available for ftp download.

The International Bathymetric Chart of the Arctic Ocean (IBCAO) released its first gridded bathymetric compilation in 1999. The IBCAO bathymetric portrayals have since supported a wide range of Arctic science activities, for example, by providing constraint for ocean circulation models and the means to define and formulate hypotheses about the geologic origin of Arctic undersea features. IBCAO Version 3.0 represents the largest improvement since 1999 taking advantage of new data sets collected by the circum‐Arctic nations, opportunistic data collected from fishing vessels, data acquired from US Navy submarines and from research ships of various nations. Built using an improved gridding algorithm, this new grid is on a 500 meter spacing, revealing much greater details of the Arctic seafloor than IBCAO Version 1.0 (2.5 km) and Version 2.0 (2.0 km). The area covered by multibeam surveys has increased from ∼6% in Version 2.0 to ∼11% in Version 3.0.

Abstract General Bathymetric Chart of the Oceans (GEBCO) has released the GEBCO_2014 grid, a new digital bathymetric model of the world ocean floor merged with land topography from publicly available digital elevation models. GEBCO_2014 has a grid spacing of 30 arc sec and updates the 2010 release (GEBCO_08) by incorporating new versions of regional bathymetric compilations from the International Bathymetric Chart of the Arctic Ocean, the International Bathymetric Chart of the Southern Ocean, the Baltic Sea Bathymetry Database, and data from the European Marine Observation and Data network bathymetry portal, among other data sources. Approximately 33% of ocean grid cells (not area) have been updated in GEBCO_2014 from the previous version, including both new interpolated depth values and added soundings. These updates include large amounts of multibeam data collected using modern equipment and navigation techniques, improving portrayed details of the world ocean floor. Of all nonland grid cells in GEBCO_2014, approximately 18% are based on bathymetric control data, i.e., primarily multibeam and single‐beam soundings or preprepared grids which may contain some interpolated values. The GEBCO_2014 grid has a mean and median depth of 3897 m and 3441 m, respectively. Hypsometric analysis reveals that 50% of the Earth's surface is composed of seafloor located 3200 m below mean sea level and that ~900 ship years of surveying would be needed to obtain complete multibeam coverage of the world's oceans.

Holgate, S.J.; Matthews, A.; Woodworth, P.L.; Rickards, L.J.; Tamisiea, M.E.; Bradshaw, E.; Foden, P.R.; Gordon, K.M.; Jevrejeva, S., and Pugh, J., 2013. New data systems and products at the Permanent Service for Mean Sea Level. Sea-level rise remains one of the most pressing societal concerns relating to climate change. A significant proportion of the global population, including many of the world's large cities, are located close to the coast in potentially vulnerable regions such as river deltas. The Permanent Service for Mean Sea Level (PSMSL) continues to evolve and provide global coastal sea-level information and products that help to develop our understanding of sea-level and land motion processes. Its work aids a range of scientific research, not only in long-term change, but also in the measurement and understanding of higher frequency variability such as storm surges and tsunamis. The PSMSL has changed considerably over the past 10 years, and the aim of this paper is to update the community about these changes as well as provide an overview of our continuing work.

A total of 954 continuous GPS position time series from 414 individual sites in nine different GPS solutions were analyzed for noise content using maximum likelihood estimation (MLE). The lengths of the series varied from around 16 months to over 10 years. MLE was used to analyze the data in two ways. In the first analysis the noise was assumed to be white noise only, a combination of white noise plus flicker noise, or a combination of white noise plus random walk noise. For the second analysis the spectral index and amplitude of the power law noise were estimated simultaneously with the white noise. In solutions where the sites were globally distributed, the noise can be best described by a combination of white noise plus flicker noise. Both noise components show latitude dependence in their amplitudes (higher at equatorial sites) together with a bias to larger values in the Southern Hemisphere. In the regional solutions, where a spatially correlated (common mode) signal has been removed, the noise is significantly lower. The spectral index of the power law in regional solutions is more varied than in the global solutions and probably reflects a mixture of local effects. A significant reduction in noise can be seen since the first continuous GPS networks began recording in the early 1990s. A comparison of the noise amplitudes to the different monument types in the Southern California Integrated GPS Network suggests that the deep drill braced monument is preferred for maximum stability.

A quasi‐global sea level data set from tide gauges has been used to investigate extreme sea level events and their spatial and temporal variabilities. Modern methods based on a nonstationary extreme value analysis have been applied to the maxima of the total elevations and surges for the period of 1970 and onward. A subset of the data was used to study changes over the 20th century. The analyses demonstrate the magnitude and timing of the seasonal cycle of extreme sea level occurrence, the magnitude of long‐term trends in extreme sea levels, the evidence for perigean and nodal astronomical tidal components in the extremes, and the relationship of the interannual variability in high water levels to other ocean and atmosphere variations as represented by climate indices. The subtraction from the extreme sea levels of the corresponding annual median sea level results in a reduction in the magnitude of trends at most stations, leading to the conclusion that much of the change in the extremes is due to change in the mean values.

Despite many of years of mapping effort, only a small fraction of the world ocean’s seafloor has been sampled for depth, greatly limiting our ability to explore and understand critical ocean and seafloor processes. Recognizing this poor state of our knowledge of ocean depths and the critical role such knowledge plays in understanding and maintaining our planet, GEBCO and the Nippon Foundation have joined forces to establish the Nippon Foundation GEBCO Seabed 2030 Project, an international effort with the objective of facilitating the complete mapping of the world ocean by 2030. The Seabed 2030 Project will establish globally distributed regional data assembly and coordination centers (RDACCs) that will identify existing data from their assigned regions that are not currently in publicly available databases and seek to make these data available. They will develop protocols for data collection (including resolution goals) and common software and other tools to assemble and attribute appropriate metadata as they assimilate regional grids using standardized techniques. A Global Data Assembly and Coordination Center (GDACC) will integrate the regional grids into a global grid and distribute to users world-wide. The GDACC will also act as the central focal point for the coordination of common data standards and processing tools as well as the outreach coordinator for Seabed 2030 efforts. The GDACC and RDACCs will collaborate with existing data centers and bathymetric compilation efforts. Finally, the Nippon Foundation GEBCO Seabed 2030 Project will encourage and help coordinate and track new survey efforts and facilitate the development of new and innovative technologies that can increase the efficiency of seafloor mapping and thus make the ambitious goals of Seabed 2030 more likely to be achieved.

Storm surges are the sea level response to meteorological conditions . Scientists and engineers need to understand the interaction of surges with the tide in order to provide better estimates of extreme sea level for use in coastal defense. Using data from five tide gauges, spaced equally along the North Sea coastline around the UK, we show that the mode of peak residual occurrence is everywhere 3 to 5 hours before the nearest high water. We reveal a previously unobserved mode that falls 1 to 2 hours prior to high water, although this cluster is not associated with the highest residuals. A simple mathematical explanation for surge clustering on the rising tide is presented. The phase shift of the tidal signal is combined with the modulation of surge production due to water depth in a model that provides a good description of the residual data set. The results contain several features of interest for flood risk management. We show that large, locally generated surges are precluded close to high water. For physically realistic arrival times of any travelling surge component, the residual peak will avoid high water for any finite tidal phase shift. Furthermore, increasing the tidal range reduces the risk of residual peaks near high water. We draw attention to the existence of critical time and space scales for surge development and decay. For reliable operational forecasts of sea level, coastal numerical models need to reproduce both tides and surges with improved accuracy.

We present a reconstruction of global sea level (GSL) since 1700 calculated from tide gauge records and analyse the evolution of global sea level acceleration during the past 300 years. We provide observational evidence that sea level acceleration up to the present has been about 0.01 mm/yr 2 and appears to have started at the end of the 18th century. Sea level rose by 6 cm during the 19th century and 19 cm in the 20th century. Superimposed on the long‐term acceleration are quasi‐periodic fluctuations with a period of about 60 years. If the conditions that established the acceleration continue, then sea level will rise 34 cm over the 21st century. Long time constants in oceanic heat content and increased ice sheet melting imply that the latest Intergovernmental Panel on Climate Change (IPCC) estimates of sea level are probably too low.

Regional sea level changes can deviate substantially from those of the global mean, can vary on a broad range of timescales, and in some regions can even lead to a reversal of long-term global mean sea level trends. The underlying causes are associated with dynamic variations in the ocean circulation as part of climate modes of variability and with an isostatic adjustment of Earth's crust to past and ongoing changes in polar ice masses and continental water storage. Relative to the coastline, sea level is also affected by processes such as earthquakes and anthropogenically induced subsidence. Present-day regional sea level changes appear to be caused primarily by natural climate variability. However, the imprint of anthropogenic effects on regional sea level-whether due to changes in the atmospheric forcing or to mass variations in the system-will grow with time as climate change progresses, and toward the end of the twenty-first century, regional sea level patterns will be a superposition of climate variability modes and natural and anthropogenically induced static sea level patterns. Attribution and predictions of ongoing and future sea level changes require an expanded and sustained climate observing system.

In the past two decades, the Argo Program has collected, processed and distributed over two million vertical profiles of temperature and salinity from the upper two kilometers of the global ocean. A similar number of subsurface velocity observations near 1000 dbar have also been collected. This paper recounts the history of the global Argo Program, from its aspiration arising out of the World Ocean Circulation Experiment, to the development and implementation of its instrumentation and telecommunication systems, and the various technical problems encountered. We describe the Argo data system and its quality control procedures, and the gradual changes in the vertical resolution and spatial coverage of Argo data from 1999 to 2019. The accuracies of the float data have been assessed by comparison with high-quality shipboard measurements, and are concluded to be 0.002°C for temperature, 2.4 dbar for pressure, and 0.01 PSS-78 for salinity, after delayed-mode adjustments. Finally, the challenges faced by the vision of an expanding Argo Program beyond 2020 are discussed.

Variability in time series of ice conditions in the Baltic Sea is examined within the context of atmospheric circulation represented by the North Atlantic Oscillation (NAO) and Arctic Oscillation (AO) winter indices using the wavelet approach. We develop methods of assessing statistical significance and confidence intervals of cross‐wavelet phase and wavelet coherence. Cross‐wavelet power for the time series indicates that the times of largest variance in ice conditions are in excellent agreement with significant power in the AO at 2.2–3.5, 5.7–7.8, and 12–20 year periods, similar patterns are also seen with the Southern Oscillation Index (SOI) and Niño3 sea surface temperature (Niño3) series. Wavelet coherence shows in‐phase linkages between the 2.2–7.8 and 12–20 year period signals in both tropical and Arctic atmospheric circulation and also with ice conditions in the Baltic Sea. These results are consistent with GCM simulations showing dynamical connections between high‐latitude surface conditions, tropical sea surface temperatures mediated by tropical wave propagation, the wintertime polar vortex, and the AO and with models of sea ice and oceanic feedbacks at decadal periods.

The GEOTRACES Intermediate Data Product 2017 (IDP2017) is the second publicly available data product of the international GEOTRACES programme, and contains data measured and quality controlled before the end of 2016. The IDP2017 includes data from the Atlantic, Pacific, Arctic, Southern and Indian oceans, with about twice the data volume of the previous IDP2014. For the first time, the IDP2017 contains data for a large suite of biogeochemical parameters as well asaerosol and rain data characterising atmospherictrace element and isotope (TEI) sources. The TEI data in the IDP2017 are quality controlled by careful assessment of intercalibration results and multi-laboratory data comparisons at crossover stations. The IDP2017 consists of two parts: (1) a compilation of digital data for more than 450 TEIs as well as standard hydrographic parameters, and (2) the eGEOTRACES Electronic Atlas providing an on-line atlas that includes more than 590 section plots and 130 animated 3D scenes. The digital data are provided in several formats, including ASCII, Excel spreadsheet, netCDF, and Ocean Data View collection. Users can download the full data packages or make their own custom selections with a new on-line data extraction service. In addition to the actual data values, the IDP2017 also contains data quality flags and 1-σ data error values where available. Quality flags and error values are useful for data filtering and for statistical analysis. Metadata about data originators, analytical methods and original publications related to the data are linked in an easily accessible way. The eGEOTRACES Electronic Atlas is the visual representation of the IDP2017 as section plots and rotating 3D scenes. The basin-wide 3D scenes combine data from many cruises and provide quick overviews of large-scale tracer distributions. These 3D scenes provide geographical and bathymetric context that is crucial for the interpretation and assessment of tracer plumes near ocean margins or along ridges. The IDP2017 is the result of a truly international effort involving 326 researchers from 25 countries. This publication provides the critical reference for unpublished data, as well as for studies that make use of a large cross-section of data from the IDP2017.

Abstract. The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The ocean component is CLIO3, which consists of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is of 3° by 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the ocean carbon cycle is represented by LOCH, a comprehensive model that takes into account both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the ice-atmosphere and ice-ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, and an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.

Introduction of active chlorophyll a fluorescence<br/>protocols, in particular fast repetition rate<br/>(FRR) fluorometry, to oceanography and limnology<br/>15 yr ago has enabled rapid assessment of photosynthetic<br/>physiology in situ. The FRR protocol generates<br/>simultaneous measurements of Photosystem II (PSII)<br/>effective absorption cross sections (termed ?PSII) and<br/>photochemical efficiency (termed Fv/Fm). Both Fv/Fm<br/>and ?PSII measurements have been utilised to examine<br/>the effects of physiological stress on the photosynthetic<br/>apparatus of phytoplankton in an ever growing number<br/>of fluorescence-based studies. However, it is now<br/>becoming clearer that in situ values of Fv/Fm and ?PSII<br/>also contain taxonomic information. Here, we present<br/>a synthesis of previously unpublished and published<br/>data, which show that Fv/Fm and ?PSII vary principally<br/>with broad-scale changes in community structure.<br/>These patterns observed in situ conform to trends<br/>observed in laboratory-grown cultures of a range of<br/>phytoplankton taxa. The magnitudes of variability in<br/>Fv/Fm and ?PSII driven by changes in phytoplankton<br/>community structure often exceed that induced by<br/>nutrient limitation (as determined from controlled<br/>nutrient addition experiments). An exception to this<br/>general trend occurs in high-nutrient, low-chlorophyll<br/>a (HNLC) regions, where strong phenotypic changes<br/>in Fv/Fm and ?PSII have been repeatedly demonstrated<br/>on relief of iron limitation. Overall, FRR fluorescence<br/>measurements of both Fv/Fm and ?PSII in natural populations<br/>represent a combination of the taxonomic ‘signature’<br/>(values of Fv/Fm and ?PSII determined by the<br/>taxa present) within the phytoplankton community<br/>that is further modified according to the (photo-) physiological<br/>status. As such, fluorescence-based investigations<br/>of mixed populations must account for potential<br/>variations in phytoplankton community structure before<br/>interpretations of physiological status are made.

We analyze the Permanent Service for Mean Sea Level (PSMSL) database of sea level time series using a method based on Monte Carlo Singular Spectrum Analysis (MC‐SSA). We remove 2–30 year quasi‐periodic oscillations and determine the nonlinear long‐term trends for 12 large ocean regions. Our global sea level trend estimate of 2.4 ± 1.0 mm/yr for the period from 1993 to 2000 is comparable with the 2.6 ± 0.7 mm/yr sea level rise calculated from TOPEX/Poseidon altimeter measurements. However, we show that over the last 100 years the rate of 2.5 ± 1.0 mm/yr occurred between 1920 and 1945, is likely to be as large as the 1990s, and resulted in a mean sea level rise of 48 mm. We evaluate errors in sea level using two independent approaches, the robust bi‐weight mean and variance, and a novel “virtual station” approach that utilizes geographic locations of stations. Results suggest that a region cannot be adequately represented by a simple mean curve with standard error, assuming all stations are independent, as multiyear cycles within regions are very significant. Additionally, much of the between‐region mismatch errors are due to multiyear cycles in the global sea level that limit the ability of simple means to capture sea level accurately. We demonstrate that variability in sea level records over periods 2–30 years has increased during the past 50 years in most ocean basins.

Sea level rise over the last 55 years is estimated to have been 1.7 ± 0.2 mm yr −1 , based upon 177 tide gauges divided into 13 regions with near global coverage and using a Glacial Isostatic Adjustment (GIA) model to correct for land movements. We present evidence from altimeter data that the rate of sea level rise around the global coastline was significantly in excess of the global average over the period 1993–2002. We also show that the globally‐averaged rate of coastal sea level rise for the decade centered on 1955 was significantly larger than any other decade during the past 55 years. In some models of sea level rise, enhanced coastal rise is a pre‐cursor of global average rise. It remains to be seen whether the models are correct and whether global‐average rates in the future reflect the high rates of coastal rise observed during the 1990s.

Over 20 global ocean tide models have been developed since 1994, primarily as a consequence of analysis of the precise altimetric measurements from TOPEX/POSEIDON and as a result of parallel developments in numerical tidal modeling and data assimilation. This paper provides an accuracy assessment of 10 such tide models and discusses their benefits in many fields including geodesy, oceanography, and geophysics. A variety of tests indicate that all these tide models agree within 2–3 cm in the deep ocean, and they represent a significant improvement over the classical Schwiderski 1980 model by approximately 5 cm rms. As a result, two tide models were selected for the reprocessing of TOPEX/POSEIDON Geophysical Data Records in late 1995. Current ocean tide models allow an improved observation of deep ocean surface dynamic topography using satellite altimetry. Other significant contributions include their applications in an improved orbit computation for TOPEX/POSEIDON and other geodetic satellites, to yield accurate predictions of Earth rotation excitations and improved estimates of ocean loading corrections for geodetic observatories, and to allow better separation of astronomical tides from phenomena with meteorological and geophysical origins. The largest differences between these tide models occur in shallow waters, indicating that the current models are still problematic in these areas. Future improvement of global tide models is anticipated with additional high‐quality altimeter data and with advances in numerical techniques to assimilate data into high‐resolution hydrodynamic models.

Bathymetry (seafloor depth), is a critical parameter providing the geospatial context for a multitude of marine scientific studies. Since 1997, the International Bathymetric Chart of the Arctic Ocean (IBCAO) has been the authoritative source of bathymetry for the Arctic Ocean. IBCAO has merged its efforts with the Nippon Foundation-GEBCO-Seabed 2030 Project, with the goal of mapping all of the oceans by 2030. Here we present the latest version (IBCAO Ver. 4.0), with more than twice the resolution (200 × 200 m versus 500 × 500 m) and with individual depth soundings constraining three times more area of the Arctic Ocean (∼19.8% versus 6.7%), than the previous IBCAO Ver. 3.0 released in 2012. Modern multibeam bathymetry comprises ∼14.3% in Ver. 4.0 compared to ∼5.4% in Ver. 3.0. Thus, the new IBCAO Ver. 4.0 has substantially more seafloor morphological information that offers new insights into a range of submarine features and processes; for example, the improved portrayal of Greenland fjords better serves predictive modelling of the fate of the Greenland Ice Sheet.