Stirring across the Antarctic Circumpolar Current's southern boundary at the prime meridian, Weddell Sea - Ria Oelerich; Karen J. Heywood; Gillian M. Damerell; Marcel du Plessis; Louise C. Biddle; Sebastiaan Swart (Ocean Science, 2023)
At the southern boundary of the Antarctic Circumpolar Current (ACC), relatively warm ACC waters encounter the colder waters surrounding Antarctica. Strong density gradients across the southern boundary indicate the presence of a frontal jet and are thought to modulate the southward heat transport across the front. In this study, the southern boundary in the Weddell Sea sector at the prime meridian is surveyed for the first time in high resolution over 2 months during an austral summer with underwater gliders occupying a transect across the front on five occasions. The five transects show that the frontal structure (i.e. hydrography, velocities and lateral density gradients) varies temporally. The results demonstrate significant, transient (a few weeks) variability of the southern boundary and its frontal jet in location, strength and width. A mesoscale cold-core eddy is identified to disrupt the southern boundary’s frontal structure and strengthen lateral density gradients across the front. The front's barrier properties are assessed using mixing length scales and potential vorticity to establish the cross-frontal exchange of properties between the ACC and the Weddell Gyre. The results show that stronger lateral density gradients caused by the mesoscale eddy strengthen the barrier-like properties of the front through reduced mixing length scales and pronounced gradients of potential vorticity. In contrast, the barrier-like properties of the southern boundary are reduced when no mesoscale eddy is influencing the density gradients across the front. Using satellite altimetry, we further demonstrate that the barrier properties over the past decade have strengthened as a result of increased meridional gradients of absolute dynamic topography and increased frontal jet speeds in comparison to previous decades. Our results emphasise that locally and rapidly changing barrier properties of the southern boundary are important to quantify the cross-frontal exchange, which is particularly relevant in regions where the southern boundary is located near the Antarctic shelf break (e.g. in the West Antarctic sector).

Technical note: Determining Arctic Ocean halocline and cold halostad depths based on vertical stability - Enrico P. Metzner; Marc Salzmann (Ocean Science, 2023)
The Arctic Ocean halocline separates the cold surface mixed layer from the underlying warm Atlantic Water (AW), and thus provides a precondition for sea ice formation. Here, we introduce a new method in which the halocline base depth is determined from vertical stability and compare it to two existing methods. We also propose a novel method for detecting the cold halostad, a layer characterized by a small vertical salinity gradient, which is formed by the Pacific Winter Water in the Canada Basin or by meltwater off the eastern coast of Greenland and off Svalbard. Our main motivation for determining the halocline base depth depending on vertical stability was that vertical stability is closely related to vertical mixing and heat exchange. Vertical stability is a crucial parameter for determining whether the halocline can prevent vertical heat exchange and protect sea ice from warm AW. When applied to measurements from ice-tethered profilers, ships, and moorings, the new method for estimating the halocline base depth provides robust results with few artifacts. Analyzing a case in which water previously homogenized by winter convection was capped by fresh water at the surface suggests that the new method captured the beginning of new halocline formation in the Eurasian Basin. Comparatively large differences between the methods for detecting the halocline base depth were found in warm AW inflow regions for which climate models predict halocline thinning and increased net surface energy fluxes from the ocean to the atmosphere.

Intraseasonal and interannual variability of sea temperature in the Arabian Sea Warm Pool - Na Li; Xueming Zhu; Hui Wang; Shouwen Zhang; Xidong Wang (Ocean Science, 2023)
The Arabian Sea Warm Pool (ASWP) is a part of the Indian Ocean Warm Pool, formed in the Arabian Sea before the onset of the Indian summer monsoon. The ASWP has a significant impact on climate change in the Indian Peninsula and globally. In this study, we examined the intraseasonal and interannual variability of sea temperature in the ASWP using the latest Simple Ocean Data Assimilation (SODA) reanalysis dataset. We quantified the contributions of sea surface heat flux forcing, horizontal advection, and vertical entrainment to the sea temperature using the mixed-layer heat budget analysis method. We also used a lead–lag correlation method to examine the relationship between the interannual variability of the ASWP and various large-scale modes in the Indo-Pacific Ocean. We found that the ASWP formed in April and decayed in June; its formation and decay processes were asymmetrical, with the decay rate being twice as fast as the formation rate. During the ASWP development phase, the sea surface heat flux forcing had the largest impact on the mixed-layer temperature with a contribution of up to 85 %. Its impact was divided into the net surface heat flux (0.41–0.50 ∘C per 5 d) and the shortwave radiation loss penetrating the mixed layer (from −0.08 ∘C per 5 d to −0.17 ∘C per 5 d). During the decay phase, the cooling effect of the vertical entrainment on the temperature variation increased (from −0.05 ∘C per 5 d to −0.18 ∘C per 5 d) and dominated the temperature variation jointly with the sea surface heat flux forcing. We also found that the ASWP has strong interannual variability related to the basin warming of the Indian Ocean. The lead–lag correlation indicated that the ASWP had a good synchronous correlation with the Indian Ocean Dipole. The ASWP had the largest correlation coefficient at a lag of 5–7 months of the Niño3.4 index, showing the characteristics of modulation by the El Niño–Southern Oscillation (ENSO). When the El Niño (La Niña) event peaked in the winter of the previous year, the ASWP that occurred before the summer monsoon was more significant (insignificant) in the following year.

Ocean 2D eddy energy fluxes from small mesoscale processes with SWOT - Elisa Carli; Rosemary Morrow; Oscar Vergara; Robin Chevrier; Lionel Renault (Ocean Science, 2023)
We investigate ocean dynamics at different scales in the Agulhas Current system, a region of important interocean exchange of heat and energy. While ocean observations and some of the most advanced climate models capture the larger mesoscale dynamics (> 100 km), the smaller-scale fronts and eddies are underrepresented. The recently launched NASA–CNES Surface Water and Ocean Topography (SWOT) wide-swath altimeter mission observes the smaller ocean geostrophic scales down to 15 km in wavelength globally. Here we will analyse different eddy diagnostics in the Agulhas Current region and quantify the contributions from the larger mesoscales observable today and the smaller scales to be observed with SWOT. Surface geostrophic diagnostics of eddy kinetic energy, strain, and energy cascades are estimated from modelled sea surface height (SSH) fields of the Massachusetts Institute of Technology general circulation model (MITgcm) latitude–longitude polar cap (LLC4320) simulation subsampled at 1/10∘. In this region, the smaller scales (

Evaluating altimetry-derived surface currents on the south Greenland shelf with surface drifters - Arthur Coquereau; Nicholas P. Foukal (Ocean Science, 2023)
The pathways and fate of freshwater in the East Greenland Coastal Current (EGCC) are crucial to the climate system. The EGCC transports large amounts of freshwater in close proximity to sites of deep open-ocean convection in the Labrador and Irminger seas. Many studies have attempted to analyze this system from models and various observational platforms, but the modeling results largely disagree with one another, and observations are limited due to the harsh conditions typical of the region. Altimetry-derived surface currents, constructed from remote-sensing observations and applying geostrophic equations, provide a continuous observational data set beginning in 1993. However, these products have historically encountered difficulties in coastal regions, and thus their validity must be checked. In this work, we use a comprehensive methodology to compare these Eulerian data to a Lagrangian data set of 34 surface drifter trajectories and demonstrate that the altimetry-derived surface currents are surprisingly capable of recovering the spatial structure of the flow field on the south Greenland shelf and can mimic the Lagrangian nature of the flow as observed from surface drifters.

Bivariate sea-ice assimilation for global-ocean analysis–reanalysis - Andrea Cipollone; Deep Sankar Banerjee; Doroteaciro Iovino; Ali Aydogdu; Simona Masina (Ocean Science, 2023)
In the last decade, various satellite missions have been monitoring the status of the cryosphere and its evolution. Besides sea-ice concentration data, available since the 1980s, sea-ice thickness retrievals are now ready to be used in global operational prediction and global reanalysis systems. Nevertheless, while univariate algorithms are commonly used to constrain sea-ice area or volume, multivariate approaches have not yet been employed due to the highly non-Gaussian distribution of sea-ice variables together with the low accuracy of thickness observations. This study extends a 3DVar system, called OceanVar, which is routinely employed in the production of global/regional operational/reanalysis products, to process sea-ice variables. The tangent/adjoint versions of an anamorphosis operator are used to locally transform the sea-ice anomalies into Gaussian control variables and back, minimizing in the latter space. The benefit achieved by such a transformation is described. Several sensitivity experiments are carried out using a suite of diverse datasets. The sole assimilation of the CryoSat-2 provides a good spatial representation of thickness distribution but still overestimates the total volume that requires the inclusion of Soil Moisture and Ocean Salinity (SMOS) mission data to converge towards the observation estimates. The intermittent availability of thickness data can lead to potential jumps in the evolution of the volume and requires a dedicated tuning. The use of the merged L4 product CS2SMOS shows the best skill score when validated against independent measurements during the melting season when satellite data are not available. This new sea-ice module is meant to simplify the future coupling with ocean variables.

Spatial and temporal variability in mode-1 and mode-2 internal solitary waves from MODIS-Terra sun glint off the Amazon shelf - Carina Regina de Macedo; Ariane Koch-Larrouy; José Carlos Bastos da Silva; Jorge Manuel Magalhães; Carlos Alessandre Domingos Lentini; Trung Kien Tran; Marcelo Caetano Barreto Rosa; Vincent Vantrepotte (Ocean Science, 2023)
The Amazon shelf is a key region for intense internal tides (ITs) and nonlinear internal solitary wave (ISWs) generation associated with them. The region shows well-marked seasonal variability (from March to July, MAMJJ, and from August to December, ASOND) of the circulation and stratification, which can both induce changes in the ISW physical characteristics. The description of the seasonal and neap–spring tidal variability in the ISWs off the Amazon shelf is performed for the first time using a meaningful data set composed of 140 MODIS-Terra imagery from 2005 to 2021, where about 500 ISW signatures were identified in the sun glint region. Previous studies have documented the existence of mode-1 ISWs, but the region appears as a newly described hotspot for mode-2 ISWs. ISW packets separated by typical mode-1 (95–170 km; 2.1–3.8 m s−1) and mode-2 (46–85 km; 1.0–1.9 m s−1) IT wavelengths have been identified and mapped coming from different IT generation sites. For each ISW, a group of waves (3 to 10) is generally follows the largest crest. The intra-packet distance between each wave in the group is about 10 to 20 km. Regions of higher occurrence of ISWs are spaced by a IT mode-1 wavelength. We make the assumption that it might correspond to the IT reflection beams at the surface, which may generate newer ISWs. The mean mode-1 and mode-2 inter-packet distances do not show significant differences according to their IT generation sites. The ISW activity is higher (more than 60 % of signatures) during spring tides than neap tides. In the region under the influence of the North Equatorial Counter Current (NECC), ISWs are separated by a mean mode-1 IT wavelength which is 14.3 % higher during ASOND than during MAMJJ due to a deeper thermocline and the reinforcement of the NECC. These ISWs are also characterized by a wider inter-packet distance distribution (higher standard deviation) that may be related to the stronger eddy kinetic energy (EKE) during ASOND compared to MAMJJ. The mean inter-packet distance of mode-2 ISWs remains almost unchanged during the two seasons, but the inter-packet distance distribution is wider in ASOND than in MAMJJ as for mode 1. Note that these results need to be treated with caution, as only few occurrences of mode-2 waves were found during MAMJJ. In the region of the NECC, the direction of propagation for all modes is very similar in MAMJJ (about 30∘ clockwise from the north), whereas, for ASOND, the ISWs propagate in a wider pathway (from 0 to 60∘ clockwise from the north), due to a much larger eddy activity. During ASOND, as the background flux goes further east, the inter-packet distances become larger (4 % for mode 1 and 7.8 % for mode 2). These results show that the reinforcement of the NECC in ASOND appears to play a role in diverting the waves towards the east, increasing their phase velocities and their eastern traveling direction component when compared to MAMJJ. Calculations of the IT velocities using the Taylor–Goldstein equation supported our results regarding the presence of ISWs associated with mode-2 ITs and additionally the IT seasonal variability.

Coastal and regional marine heatwaves and cold spells in the northeastern Atlantic - Amélie Simon; Coline Poppeschi; Sandra Plecha; Guillaume Charria; Ana Russo (Ocean Science, 2023)
The latest Intergovernmental Panel on Climate Change (IPCC) report describes an increase in the number and intensity of marine heatwaves (MHWs) and a decrease in marine cold spells (MCSs) in the global ocean. However, these reported changes are not uniform on a regional to local basis, and it remains unknown if coastal areas follow the open-ocean trends. Surface ocean temperature measurements collected by satellites (from 1982–2022) and 13 coastal buoys (from 1990–2022) are analyzed in the northeastern Atlantic and three subregions: the English Channel, Bay of Brest and Bay of Biscay. The activity metric, combining the number of events, intensity, duration and spatial extent, is used to evaluate the magnitude of these extreme events. The results from in situ and satellite datasets for each of the studied regions are quite in agreement, although the satellite dataset underestimates the amplitude of activity for both MHWs and MCSs. This supports the applicability of the method to both in situ and satellite data, albeit with caution on the amplitude of these events. Also, this localized study in European coastal northeastern Atlantic water highlights that similar changes are being seen in coastal and open oceans regarding extreme events of temperature, with MHWs being more frequent and longer and extending over larger areas, while the opposite is seen for MCSs. These trends can be explained by changes in both the mean of and variance in sea-surface temperature. In addition, the pace of evolution and dynamics of marine extreme events differ among the subregions. Among the three studied subregions, the English Channel is the region experiencing the strongest increase in summer MHW activity over the last 4 decades. Summer MHWs were very active in the English Channel in 2022 due to long events, in the Bay of Biscay in 2018 due to intense events and in the Bay of Brest in 2017 due to a high occurrence of events. Winter MCSs were the largest in 1987 and 1986 due to long and intense events in the English Channel. Finally, our findings suggest that at an interannual timescale, the positive North Atlantic Oscillation favors the generation of strong summer MHWs in the northeastern Atlantic, while low-pressure conditions over northern Europe and a high off the Iberian Peninsula in winter dominate for MCSs. A preliminary analysis of air–sea heat fluxes suggests that, in this region, reduced cloud coverage is a key parameter for the generation of summer MHWs, while strong winds and increased cloud coverage are important for the generation of winter MCSs.

Regional modeling of internal-tide dynamics around New Caledonia – Part 1: Coherent internal-tide characteristics and sea surface height signature - Arne Bendinger; Sophie Cravatte; Lionel Gourdeau; Laurent Brodeau; Aurélie Albert; Michel Tchilibou; Florent Lyard; Clément Vic (Ocean Science, 2023)
The southwestern tropical Pacific exhibits a complex bathymetry and represents a hot spot of internal-tide generation. Based on a tailored high-resolution regional model, we investigate for the first time the internal-tide field around the New Caledonia islands through energy budgets that quantify the coherent internal-tide generation, propagation, and dissipation. A total of 15.27 GW is converted from the barotropic to the baroclinic M2 tide with the main conversion sites associated with the most prominent bathymetric structures such as continental slopes and narrow passages in the north (2.17 GW) and ridges and seamounts south of New Caledonia (3.92 GW). The bulk of baroclinic energy is generated in shallow waters around 500 m depth and on critical to supercritical slopes, highlighting the limitations of linear semi-analytical models in those areas. Despite the strongly dominant mode-1 generation, more than 50 % of the locally generated energy either dissipates in the near field close to the generation sites or loses coherence. The remaining baroclinic energy propagates within well-defined tidal beams with baroclinic energy fluxes of up to 30 kW m−1 toward the open ocean. The New Caledonia site represents a challenge for SWOT (Surface Water and Ocean Topography) observability of balanced motion in the presence of internal tides with sea surface height (SSH) signatures >6 cm at similar wavelengths. We show for our study region that a correction of SSH for the coherent internal tide potentially increases the observability of balanced motion from wavelengths >160 km to well below 100 km.

Nouméa: a new multi-mission calibration and validation site for past and future altimetry missions? - Clémence Chupin; Valérie Ballu; Laurent Testut; Yann-Treden Tranchant; Jérôme Aucan (Ocean Science, 2023)
Today, monitoring the evolution of sea level in coastal areas is of importance, since almost 11 % of the world's population lives in low-lying areas. Reducing uncertainties in sea level estimates requires a better understanding of both altimetry measurements and local sea level dynamics. In New Caledonia, the Nouméa lagoon is an example of this challenge, as altimetry, coastal tide gauge, and vertical land motions from global navigation satellite systems (GNSSs) do not provide consistent information. The GEOCEAN-NC 2019 field campaign addresses this issue with deployments of in situ instruments in the lagoon (GNSS buoy, pressure gauge, etc.), with a particular focus on the crossover of one Jason-series track and two Sentinel-3A missions tracks. In this study, we propose a method to virtually transfer the Nouméa tide gauge at the altimetry crossover point, using in situ data from the field campaign. Following the philosophy of calibration and validation (Cal/Val) studies, we derive absolute altimeter bias time series over the entire Jason and Sentinel-3A periods. Overall, our estimated altimeter mean biases are slightly larger by 1–2 cm compared to Corsica and Bass Strait results, with inter-mission biases in line with those of Bass Strait site. Uncertainties still remain regarding the determination of our vertical datum, only constrained by the three days of the GNSS buoy deployment. With our method, we are able to re-analyse about 20 years of altimetry observations and derive a linear trend of −0.2 ± 0.1 mm yr−1 over the bias time series. Compared to previous studies, we do not find any significant uplift in the area, which is more consistent with the observations of inland permanent GNSS stations. These results support the idea of developing Cal/Val activities in the lagoon, which is already the subject of several experiments for the scientific calibration phase of the SWOT wide-swath altimetry mission.

Potential artifacts in conservation laws and invariants inferred from sequential state estimation - Carl Wunsch; Sarah Williamson; Patrick Heimbach (Ocean Science, 2023)
In sequential estimation methods often used in oceanic and general climate calculations of the state and of forecasts, observations act mathematically and statistically as source or sink terms in conservation equations for heat, salt, mass, and momentum. These artificial terms obscure the inference of the system's variability or secular changes. Furthermore, for the purposes of calculating changes in important functions of state variables such as total mass and energy or volumetric current transports, results of both filter and smoother-based estimates are sensitive to misrepresentation of a large variety of parameters, including initial conditions, prior uncertainty covariances, and systematic and random errors in observations. Here, toy models of a coupled mass–spring oscillator system and of a barotropic Rossby wave system are used to demonstrate many of the issues that arise from such misrepresentations. Results from Kalman filter estimates and those from finite interval smoothing are analyzed. In the filter (and prediction) problem, entry of data leads to violation of conservation and other invariant rules. A finite interval smoothing method restores the conservation rules, but uncertainties in all such estimation results remain. Convincing trend and other time-dependent determinations in “reanalysis-like” estimates require a full understanding of models, observations, and underlying error structures. Application of smoother-type methods that are designed for optimal reconstruction purposes alleviate some of the issues.

The Iceland–Faroe warm-water flow towards the Arctic estimated from satellite altimetry and in situ observations - Bogi Hansen; Karin M. H. Larsen; Hjálmar Hátún; Steffen M. Olsen; Andrea M. U. Gierisch; Svein Østerhus; Sólveig R. Ólafsdóttir (Ocean Science, 2023)
The inflow of warm and saline Atlantic water to the Arctic Mediterranean (Nordic Seas and Arctic Ocean) between Iceland and the Faroes (IF inflow) is the strongest Atlantic inflow branch in terms of volume transport and is associated with a large transport of heat towards the Arctic. The IF inflow is monitored in a section east of the Iceland–Faroe Ridge (IFR) by use of sea level anomaly (SLA) data from satellite altimetry, a method that has been calibrated by in situ observations gathered over 2 decades. Monthly averaged surface velocity anomalies calculated from SLA data were strongly correlated with anomalies measured by moored acoustic Doppler current profilers (ADCPs) with consistently higher correlations when using the reprocessed SLA data released in December 2021 rather than the earlier version. In contrast to the earlier version, the reprocessed data also had the correct conversion factor between sea level slope and surface velocity required by geostrophy. Our results show that the IF inflow crosses the IFR in two separate branches. The Icelandic branch is a jet over the Icelandic slope with average surface speed exceeding 20 cm s−1, but it is narrow and shallow with an average volume transport of less than 1 Sv (106 m3 s−1). Most of the Atlantic water crosses the IFR close to its southernmost end in the Faroese branch. Between these two branches, water from the Icelandic branch turns back onto the ridge in a retroflection with a recirculation over the northernmost bank on the IFR. Combining multi-sensor in situ observations with satellite SLA data, monthly mean volume transport of the IF inflow has been determined from January 1993 to December 2021. The IF inflow is part of the Atlantic Meridional Overturning Circulation (AMOC), which is expected to weaken under continued global warming. Our results show no weakening of the IF inflow. Annually averaged volume transport of Atlantic water through the monitoring section had a statistically significant (95 % confidence level) increasing trend of (0.12±0.10) Sv per decade. Combined with increasing temperature, this caused an increase of 13 % in the heat transport, relative to 0 ∘C, towards the Arctic of the IF inflow over the 29 years of monitoring. The near-bottom layer over most of the IFR is dominated by cold water of Arctic origin that may contribute to the overflow across the ridge. Our observations confirm a dynamic link between the overflow and the Atlantic water flow above. The results also provide support for a previously posed hypothesis that this link may explain the difficulties in reproducing observed transport variations in the IF inflow in numerical ocean models, with consequences for its predictability under climate change.

Assessment of Indonesian Throughflow transports from ocean reanalyses with mooring-based observations - Magdalena Fritz; Michael Mayer; Leopold Haimberger; Susanna Winkelbauer (Ocean Science, 2023)
The transport of heat and freshwater from the Pacific Ocean to the Indian Ocean via the Indonesian seas is commonly referred to as the Indonesian Throughflow (ITF). The interaction between the ITF and large-scale phenomena occurring from intraseasonal to decadal timescales reflects its connection to the global ocean and the climate system, indicating the need for monitoring the ITF region. In situ observations in this region are highly valuable, but they are temporally and spatially insufficient for near-real-time monitoring. Ocean reanalyses have the potential to serve as near-real-time monitoring tools and to extend time series backward in time, yet a comprehensive quality assessment of their realism in this region with challenging bathymetry has been lacking so far. We focus on oceanic transports diagnosed from the Copernicus Marine Service (CMEMS) Global Reanalysis Ensemble Product (GREP) and the higher-resolution product GLORYS12V1, totaling six reanalysis products. They are validated against in situ observations taken from two different monitoring programs, namely International Nusantara Stratification and Transport (INSTANT 2004–2006) and Monitoring the Indonesian Throughflow (MITF 2006–2011 and 2013–2017), resulting in a total time series of about 11.5 years in the major inflow passage of the Makassar Strait and shorter sampled time series in the Lombok Strait, the Ombai Strait, and the Timor Passage. Results show that there is reasonable agreement between reanalysis-based transports and observations in terms of means, seasonal cycles, and variability, although some shortcomings stand out. The lower-resolution products do not represent the spatial structure of the flow accurately. They also tend to underestimate the integrated net flow in the narrower straits of Lombok and Ombai, an aspect that is improved in GLORYS12V1. Reanalyses tend to underestimate the effect of seasonal Kelvin waves on the transports, which leads to errors in the mean seasonal cycle. Interannual variations of reanalyzed transports agree well with observations, but uncertainties are much larger on sub-annual variability. Finally, as an application of physically consistent and observationally constrained fields provided by ocean reanalyses, we study the impact of the vertically varying pressure gradient on the vertical structure of the ITF to better understand an apparent two-layer regime of the flow.

Sensitivity of gyre-scale marine connectivity estimates to fine-scale circulation - Saeed Hariri; Sabrina Speich; Bruno Blanke; Marina Lévy (Ocean Science, 2023)
We investigated the connectivity properties of an idealized western boundary current system separating two ocean gyres, where the flow is characterized by a well-defined mean circulation as well as energetic fine-scale features (i.e., mesoscale and submesoscale currents). We used a time-evolving 3D flow field from a high-resolution (HR-3D) ocean model of this system. In order to evaluate the role of the fine scales in connectivity estimates, we computed Lagrangian trajectories in three different ways: using the HR-3D flow, using the same flow but filtered on a coarse-resolution grid (CR-3D), and using the surface layer flow only (HR-SL). We examined connectivity between the two gyres along the western boundary current and across it by using and comparing different metrics, such as minimum and averaged values of transit time between 16 key sites, arrival depths, and probability density functions of transit times. We find that when the fine-scale flow is resolved, the numerical particles connect pairs of sites faster (between 100 to 300 d) than when it is absent. This is particularly true for sites that are along and near the jets separating the two gyres. Moreover, the connectivity is facilitated when 3D instead of surface currents are resolved. Finally, our results emphasize that ocean connectivity is 3D and not 2D and that assessing connectivity properties using climatologies or low-resolution velocity fields yields strongly biased estimates.

Geostrophic adjustment on the midlatitude β plane - Itamar Yacoby; Nathan Paldor; Hezi Gildor (Ocean Science, 2023)
Analytical and numerical solutions of the linearized rotating shallow water equations are combined to study the geostrophic adjustment on the midlatitude β plane. The adjustment is examined in zonal periodic channels of width Ly=4Rd (narrow channel, where Rd is the radius of deformation) and Ly=60Rd (wide channel) for the particular initial conditions of a resting fluid with a step-like height distribution, η0. In the one-dimensional case, where η0=η0(y), we find that (i) β affects the geostrophic state (determined from the conservation of the meridional vorticity gradient) only when b=cot(ϕ0)RdR≥0.5 (where ϕ0 is the channel's central latitude, and R is Earth's radius); (ii) the energy conversion ratio varies by less than 10 % when b increases from 0 to 1; (iii) in wide channels, β affects the waves significantly, even for small b (e.g., b=0.005); and (iv) for b=0.005, harmonic waves approximate the waves in narrow channels, and trapped waves approximate the waves in wide channels. In the two-dimensional case, where η0=η0(x), we find that (i) at short times the spatial structure of the steady solution is similar to that on the f plane, while at long times the steady state drifts westward at the speed of Rossby waves (harmonic Rossby waves in narrow channels and trapped Rossby waves in wide channels); (ii) in wide channels, trapped-wave dispersion causes the equatorward segment of the wavefront to move faster than the northern segment; (iii) the energy of Rossby waves on the β plane approaches that of the steady state on the f plane; and (iv) the results outlined in (iii) and (iv) of the one-dimensional case also hold in the two-dimensional case.

Observed multi-decadal trends in subsurface temperature adjacent to the East Australian Current - Michael P. Hemming; Moninya Roughan; Neil Malan; Amandine Schaeffer (Ocean Science, 2023)
Sea surface temperature observations have shown that western boundary currents, such as the East Australian Current (EAC), are warming faster than the global average. However, we know little about coastal temperature trends inshore of these rapidly warming regions, particularly below the surface. In addition to this, warming rates are typically estimated linearly, making it difficult to know how these rates have changed over time. Here we use long-term in situ temperature observations through the water column at five coastal sites between approximately 27.3–42.6∘ S to estimate warming trends between the ocean surface and the bottom. Using an advanced trend detection method, we find accelerating warming trends at multiple depths in the EAC extension region at 34.1 and 42.6∘ S. We see accelerating trends at the surface and bottom at 34.1∘ S but similar trends in the top 20 m at 42.6∘ S. We compare several methods, estimate uncertainty, and place our results in the context of previously reported trends, highlighting that magnitudes are depth-dependent, vary across latitude, and are sensitive to the data time period chosen. The spatial and temporal variability in the long-term temperature trends highlight the important role of regional dynamics against a background of broad-scale ocean warming. Moreover, considering that recent studies of ocean warming typically focus on surface data only, our results show the necessity of subsurface data for the improved understanding of regional climate change impacts.

Impact of sea level changes on future wave conditions along the coasts of western Europe - Alisée A. Chaigneau; Stéphane Law-Chune; Angélique Melet; Aurore Voldoire; Guillaume Reffray; Lotfi Aouf (Ocean Science, 2023)
Wind waves and swells are major drivers of coastal environment changes and coastal hazards such as coastal flooding and erosion. Wave characteristics are sensitive to changes in water depth in shallow and intermediate waters. However, wave models used for historical simulations and projections typically do not account for sea level changes whether from tides, storm surges, or long-term sea level rise. In this study, the sensitivity of projected changes in wave characteristics to the sea level changes is investigated along the Atlantic European coastline. For this purpose, a global wave model is dynamically downscaled over the northeastern Atlantic for the 1970–2100 period under the SSP5–8.5 climate change scenario. Twin experiments are performed with or without the inclusion of hourly sea level variations from regional 3D ocean simulations in the regional wave model. The largest impact of sea level changes on waves is located on the wide continental shelf where shallow-water dynamics prevail, especially in macro-tidal areas. For instance, in the Bay of Mont-Saint-Michel in France, due to an average tidal range of 10 m, extreme historical wave heights were found to be up to 1 m higher (+30 %) when sea level variations are included. At the end of the 21st century, extreme significant wave heights are larger by up to +40 % (+60 cm), mainly due to the effect of tides and mean sea level rise. The estimates provided in this study only partially represent the processes responsible for the sea-level–wave non-linear interactions due to model limitations in terms of resolution and the processes included.

Variability in coastal downwelling circulation in response to high-resolution regional atmospheric forcing off the Pearl River estuary - Wenfeng Lai; Jianping Gan (Ocean Science, 2023)
We investigated the variabilities in coastal circulation and dynamics in response to spatiotemporally variable high-resolution atmospheric forcing off the Pearl River estuary during the downwelling wind. Our investigation focused on the dynamics of coastal downwelling circulation in response to variable atmospheric forcing of (1) single-station observation, (2) global reanalysis data, and (3) a high-resolution regional atmospheric model. We found that the high-resolution atmospheric model significantly improved the representations of the near-surface wind and air temperature, and the ocean model driven by the high-resolution and spatially variable atmospheric forcing improved the circulation and associated hydrographic properties in the coastal ocean. Momentum and vorticity analyses further revealed that the cross-isobath water exchange was primarily governed by the along-isobath pressure gradient force (PGF), which was influenced by different components of the atmospheric forcing. The spatial–temporal variability in high-resolution wind forcing determined the strength and structure of coastal circulation and improved estimates of cross-isobath transport and the associated PGF by refining the net stress curl and nonlinear advection of relative vorticity in the simulation. The high-resolution heat forcing can greatly improve the sea surface temperature simulation and adjust the nonlinear advection of relative vorticity, resulting in changes in cross-isobath transport.

The Weddell Gyre heat budget associated with the Warm Deep Water circulation derived from Argo floats - Krissy Anne Reeve; Torsten Kanzow; Olaf Boebel; Myriel Vredenborg; Volker Strass; Rüdiger Gerdes (Ocean Science, 2023)
The Weddell Gyre plays an important role in the global climate system by supplying heat to underneath the ice shelves and in the formation of deep and bottom water masses, which have been subject to widespread warming over past decades. In this study, we investigate the re-distribution of heat throughout the Weddell Gyre by diagnosing the terms of the heat conservation equation for a 1000 m thick layer of water encompassing the core of Warm Deep Water. The spatial distributions of the different advective and diffusive terms in terms of heat tendencies are estimated using gridded climatologies of temperature and velocity, obtained from Argo floats in the Weddell Gyre from 2002 to 2016. While the results are somewhat noisy on the grid scale and the representation of the effects of eddy mixing is highly uncertain due to the need to parameterise them by means of turbulent diffusion, the heat budget (i.e. the sum of all terms) closes (within the uncertainty range) when integrated over the open inflow region in the southern limb, whereas the interior circulation cell remains unbalanced. There is an overall balance in the southern limb between the mean horizontal advection and horizontal turbulent diffusion of heat, whereas the vertical terms contribute comparatively little to the heat budget throughout the Weddell Gyre. Heat convergence due to mean horizontal advection balances with divergence due to horizontal turbulent diffusion in the open southern limb of the Weddell Gyre. In contrast, heat divergence due to mean horizontal advection is much weaker than convergence due to horizontal turbulent diffusion in the interior circulation cell of the Weddell Gyre, due to large values in the latter along the northern boundary due to large meridional temperature gradients. Heat is advected into the Weddell Gyre along the southern limb, some of which is turbulently diffused northwards into the interior circulation cell, while some is likely turbulently diffused southwards towards the shelf seas. This suggests that horizontal turbulent diffusion plays a role in transporting heat both towards the gyre interior where upwelling occurs and towards the ice shelves. Horizontal turbulent diffusion is also a mechanism by which heat can be transported into the Weddell Gyre across the open northern boundary. Temporal deviations from the mean terms are not included due to study limitations. In order to appreciate the role of transient eddying processes, a continued effort to increase the spatial and temporal coverage of observations in the eastern Weddell Sea is required.

Mode-1 N 2 internal tides observed by satellite altimetry - Zhongxiang Zhao (Ocean Science, 2023)
Satellite altimetry provides a unique technique for observing the sea surface height (SSH) signature of internal tides from space. Previous studies have constructed empirical internal tide models for the four largest constituents M2, S2, K1, and O1 by satellite altimetry. Yet no empirical models have been constructed for minor tidal constituents. In this study, we observe mode-1 N2 internal tides (the fifth largest constituent) using about 100 satellite years of SSH data from 1993 to 2019. We employ a recently developed mapping procedure that includes two rounds of plane wave analysis and a two-dimensional bandpass filter in between. The results show that mode-1 N2 internal tides have millimeter-scale SSH amplitudes. Model errors are estimated from background internal tides that are mapped using the same altimetry data but with a tidal period of 12.6074 h (N2 minus 3 min). The global mean error variance is about 25 % that of N2, suggesting that the mode-1 N2 internal tides can overcome model errors in some regions. We find that the N2 and M2 internal tides have similar spatial patterns and that the N2 amplitudes are about 20 % of the M2 amplitudes. Both features are determined by the N2 and M2 barotropic tides. The mode-1 N2 internal tides are observed to propagate hundreds to thousands of kilometers in the open ocean. The globally integrated N2 and M2 internal tide energies are 1.8 and 30.9 PJ, respectively. Their ratio of 5.8 % is larger than the theoretical value of 4 % because the N2 internal tides contain relatively larger model errors. Our mode-1 N2 internal tide model is evaluated using independent satellite altimetry data in 2020 and 2021. The results suggest that the model can make internal tide correction in regions where the model variance is greater than twice the error variance. This work demonstrates that minor internal tidal constituents can be observed using multiyear multi-satellite altimetry data and dedicated mapping techniques.

Coupled hydrological and hydrodynamic modelling application for climate change impact assessment in the Nemunas river watershed–Curonian Lagoon–southeastern Baltic Sea continuum - Rasa Idzelytė; Natalja Čerkasova; Jovita Mėžinė; Toma Dabulevičienė; Artūras Razinkovas-Baziukas; Ali Ertürk; Georg Umgiesser (Ocean Science, 2023)
We analyse the cumulative impacts of climate change in a complex basin–lagoon–sea system continuum, which covers the Nemunas river basin, Curonian Lagoon, and the southeastern part of the Baltic Sea. A unique, state-of-the-art coupled modelling system was developed using hydrological and hydrodynamic models. The results of four regional downscaled models from the Rossby Centre high-resolution regional atmospheric climate model have been bias-corrected using in situ measurements and were used as forcing to assess the changes that the continuum will undergo until the end of this century. Results show that the Curonian Lagoon will be subjected to higher river discharges that in turn increase the outgoing fluxes into the Baltic Sea. Through these higher fluxes, both the water residence time and saltwater intrusion into the lagoon event frequency will decrease. Most of these changes will be more pronounced in the northern part of the lagoon, which is more likely to be influenced by the variations in the Nemunas river discharge. Its delta area may be susceptible to flooding as a result of the elevated discharge during winter. The southern part of the lagoon will experience lesser changes. While water temperatures in the entire lagoon and the southeastern Baltic Sea will steadily increase and salinity will decrease, the foreseen changes in the physical characteristics will not cause significant shifts in the ecosystem functioning but may affect the nutrient retention capacity. However, some ecosystem services such as ice fishing are expected to vanish completely due to the loss of ice cover.

Consistent picture of the horizontal circulation of the Atlantic Ocean over 3 decades - Verónica Caínzos; M. Dolores Pérez-Hernández; Daniel Santana-Toscano; Cristina Arumí-Planas; Alonso Hernández-Guerra (Ocean Science, 2023)
The circulation in the Atlantic Ocean is marked by the complex system of pathways of the Atlantic Meridional Overturning Circulation (AMOC). These currents change meridionally due to the interaction with nearby water masses. Hydrographic data provide the opportunity to characterize these currents for the whole water column with high-resolution data over the last 30 years. Moreover, inverse methods enable the quantification of absolute zonal transports across these sections, determining the strength of each current at a certain latitude in terms of mass, heat, and freshwater, as well as their transport-weighted temperature and salinity. Generally, no changes can be found among decades for each of the currents in terms of transport or their properties. In the South Atlantic, the circulation describes the subtropical gyre affected by several recirculations. There are nearly 61 Sv entering from the Southern and Indian oceans at 45∘ S. The South Atlantic subtropical gyre exports 17.0 ± 1.2 Sv and around 1 PW northward via the North Brazil Current, as well as −55 Sv southward at 45∘ S into the Antarctic Circumpolar Current. In the North Atlantic, most of the transport is advected northward via the western boundary currents, which reduce their strength as they take part in convection processes in the subpolar North Atlantic, also reflected in the northward progress of mass and heat transport. Deep layers carry waters southward along the western boundary, maintaining similar values of mass and heat transport until the separation into an eastern branch crossing the mid-Atlantic Ridge in the South Atlantic. Abyssal waters originating in the Southern Ocean are distributed along the South Atlantic mainly through its western subbasin, flowing northward up to 24.5∘ N, subjected to an increasing trend in their temperature with time.

The acceleration of sea-level rise along the coast of the Netherlands started in the 1960s - Iris Keizer; Dewi Le Bars; Cees de Valk; André Jüling; Roderik van de Wal; Sybren Drijfhout (Ocean Science, 2023)
The global acceleration of sea-level rise (SLR) during the 20th century is now established. On the local scale, this is harder to establish as several drivers of SLR play a role, which can mask the acceleration. Here, we study the rate of SLR along the coast of the Netherlands from the average of six tide gauge records covering the period 1890–2021. To isolate the effects of the wind field variations and the nodal tide from the local sea-level trend, we use four generalised additive models (GAMs) which include different predictive variables. From the sea-level trend estimates, we obtain the continuous evolution of the rate of SLR and its uncertainty over the observational period. The standard error in the estimation of the rate of SLR is reduced when we account for nodal-tide effects and is reduced further when we also account for the wind effects, meaning these provide better estimates of the rate of SLR. A part of the long-term SLR is due to wind forcing related to a strengthening and northward shift of the jet stream, but this SLR contribution decelerated over the observational period. Additionally, we detect wind-forced sea-level variability on multidecadal timescales with an amplitude of around 1 cm. Using a coherence analysis, we identify both the North Atlantic Oscillation and the Atlantic Multidecadal Variability as its drivers. Crucially, accounting for the nodal-tide and wind effects changes the estimated rate of SLR, unmasking an SLR acceleration that started in the 1960s. Our best-fitting GAM, which accounts for nodal and wind effects, yields a rate of SLR of about 1.72.21.3 mm yr−1 in 1900–1919 and 1.51.91.2 mm yr−1 in 1940–1959 compared to 2.93.52.4 mm yr−1 in 2000–2019 (where the lower and upper bounds denote the 5th and 95th percentiles). If we discount the nodal tide, wind and fluctuation effects and assume a constant rate of SLR, then the probability (p value) of finding a rate difference between 1940–1959 and 2000–2019 of at least our estimate is smaller than 1 %. Consistent with global observations and the expectations based on the physics of global warming, our results show unequivocally that SLR along the Dutch coast has accelerated since the 1960s.

Global variability of high-nutrient low-chlorophyll regions using neural networks and wavelet coherence analysis - Gotzon Basterretxea; Joan S. Font-Muñoz; Ismael Hernández-Carrasco; Sergio A. Sañudo-Wilhelmy (Ocean Science, 2023)
We examine 20 years of monthly global ocean color data and modeling outputs of nutrients using self-organizing map (SOM) analysis to identify characteristic spatial and temporal patterns of high-nutrient low-chlorophyll (HNLC) regions and their association with different climate modes. The global nitrate-to-chlorophyll ratio threshold of NO3 : Chl > 17 (mmol NO3 mg Chl−1) is estimated to be a good indicator of the distribution limit of this unproductive biome that, on average, covers 92 × 106 km2 (∼ 25 % of the ocean). The trends in satellite-derived surface chlorophyll (0.6 ± 0.4 % yr−1 to 2 ± 0.4 % yr−1) suggest that HNLC regions in polar and subpolar areas have experienced an increase in phytoplankton biomass over the last decades, but much of this variation, particularly in the Southern Ocean, is produced by a climate-driven transition in 2009–2010. Indeed, since 2010, the extent of the HNLC zones has decreased at the poles (up to 8 %) and slightly increased at the Equator (

Intense anticyclones at the global Argentine Basin array of the Ocean Observatory Initiative - Camila Artana; Christine Provost (Ocean Science, 2023)
We analyzed physical oceanic parameters gathered by a mooring array at mesoscale spatial sampling deployed in the Argentine Basin within the Ocean Observatory Initiative, a National Science Foundation major research facility. The array was maintained at 42∘ S and 42∘ W, a historically sparsely sampled region with small ocean variability, over 34 months from March 2015 to January 2018. The data documented four anticyclonic extreme-structure events in 2016. The four anticyclonic structures had different characteristics (size, vertical extension, origin, lifetime and Rossby number). They all featured near-inertial waves (NIWs) trapped at depth and low Richardson values well below the mixed layer. Low Richardson values suggest favorable conditions for mixing. The anticyclonic features likely act as mixing structures at the pycnocline, bringing heat and salt from the South Atlantic Central Water to the Antarctic Intermediate Waters. The intense structures were unique in the 29-year-long satellite altimetry record at the mooring site. The Argentine Basin is populated with many anticyclones, and mixing associated with trapped NIWs probably plays an important role in setting up the upper-water-mass characteristics in the basin.

Dense water formation in the eastern Mediterranean under a global warming scenario - Iván M. Parras-Berrocal; Rubén Vázquez; William Cabos; Dimitry V. Sein; Oscar Álvarez; Miguel Bruno; Alfredo Izquierdo (Ocean Science, 2023)
Dense water formation in the eastern Mediterranean (EMed) is essential in sustaining the Mediterranean overturning circulation. Changes in the sources of dense water in the EMed point to changes in the circulation and water properties of the Mediterranean Sea. Here we examine with a regional climate system model the changes in the dense water formation in the EMed through the 21st century under the RCP8.5 emission scenario. Our results show a shift in the dominant source of Eastern Mediterranean Deep Water (EMDW) from the Adriatic Sea to the Aegean Sea in the first half of the 21st century. The projected dense water formation is reduced by 75 % for the Adriatic Sea, 84 % for the Aegean Sea, and 83 % for the Levantine Sea by the end of the century. The reduction in the intensity of deep water formation is related to hydrographic changes in surface and intermediate water that strengthen the vertical stratification, hampering vertical mixing and thus convection. Those changes have an impact on the water that flows through the Strait of Sicily to the western Mediterranean and therefore on the whole Mediterranean system.

Factors influencing the meridional width of the equatorial deep jets - Swantje Bastin; Martin Claus; Richard J. Greatbatch; Peter Brandt (Ocean Science, 2023)
Equatorial deep jets (EDJs) are vertically alternating, stacked zonal currents that flow along the Equator in all three ocean basins at intermediate depth. Their structure can be described quite well by the sum of high-baroclinic-mode equatorial Kelvin and Rossby waves. However, the EDJ meridional width is larger by a factor of 1.5 than inviscid theory predicts for such waves. Here, we use a set of idealised model configurations representing the Atlantic Ocean to investigate the contributions of different processes to the enhanced EDJ width. Corroborated by the analysis of shipboard velocity sections, we show that widening of the EDJs by momentum loss due to irreversible mixing or other processes contributes more to their enhanced time mean width than averaging over meandering of the jets. Most of the widening due to meandering can be attributed to the strength of intraseasonal variability in the jets' depth range, suggesting that the jets are meridionally advected by intraseasonal waves. A slightly weaker connection to intraseasonal variability is found for the EDJ widening by momentum loss. These results enhance our understanding of the dynamics of the EDJs and, more generally, of equatorial waves in the deep ocean.

Assessing the capability of three different altimetry satellite missions to observe the Northern Current by using a high-resolution model - Alice Carret; Florence Birol; Claude Estournel; Bruno Zakardjian (Ocean Science, 2023)
Over the last 3 decades, satellite altimetry has observed sea surface height variations, providing a regular monitoring of the surface ocean circulation. Altimetry measurements have an intrinsic signal-to-noise ratio that limits the spatial scales of the currents that can be captured. However, the recent progress made on both altimetry sensors and data processing allows us to observe smaller geophysical signals, offering new perspectives in coastal areas where these structures are important. In this methodological study, we assess the ability of three altimeter missions with three different technologies to capture the Northern Current (northwestern Mediterranean Sea) and its variability, namely Jason-2 (Ku-band low-resolution-mode altimeter, launched in 2008), SARAL/AltiKa (Ka-band low-resolution-mode altimeter, launched in 2013) and Sentinel-3A (synthetic aperture radar altimeter, launched in 2016). Therefore, we use a high-resolution regional model as a reference. We focus along the French coast of Provence, where we first show that the model is very close to the observations of high-frequency radars and gliders in terms of surface current estimates. In the model, the Northern Current is observed 15–20 km from the coast on average, with a mean core velocity of 0.39 m s−1. Its signature in terms of sea level consists of a drop whose mean value at 6.14∘ E is 6.9 cm, extending over 20 km. These variations show a clear seasonal pattern, but high-frequency signals are also present most of the time. In comparison, in 1 Hz altimetry data, the mean sea level drop associated with the Northern Current is overestimated by 3.0 cm for Jason-2, but this overestimation is significantly less with SARAL/AltiKa and Sentinel-3A (0.3 and 1.4 cm respectively). In terms of corresponding sea level variability, Jason-2 and SARAL altimetry estimates are larger than the model reference (+1.3 and +1 cm respectively), whereas Sentinel-3A shows closer values (−0.4 cm). When we derive geostrophic surface currents from the satellite sea level variations without any data filtering, in comparison to the model, the standard deviations of the velocity values are also very different from one mission to the other (3.7 times too large for Jason-2 but 2.4 and 2.9 times too large for SARAL and Sentinel-3A respectively). When low-pass filtering altimetry sea level data with different cutoff wavelengths, the best agreement between the model and the altimetry distributions of velocity values are obtained with a 60, 30 and 40–50 km cutoff wavelength for Jason-2, SARAL and Sentinel-3A data respectively. This study shows that using a high-resolution model as a reference for altimetry data allows us not only to illustrate how the advances in the performances of altimeters and in the data processing improve the observation of coastal currents but also to quantify the corresponding gain.

Improving the thermocline calculation over the global ocean - Emmanuel Romero; Leonardo Tenorio-Fernandez; Esther Portela; Jorge Montes-Aréchiga; Laura Sánchez-Velasco (Ocean Science, 2023)
According to the typical thermal structure of the ocean, the water column can be divided into three layers: the mixed layer, the thermocline and the deep layer. In this study, we provide a new methodology, based on a function adjustment to the temperature profile, to locate the minimum and maximum depths of the strongest thermocline. We first validated our methodology by comparing the mixed layer depth obtained with the method proposed here with three other methods from previous studies. Since we found a very good agreement between the four methods we used the function adjustment to compute the monthly climatologies of the maximum thermocline depth and the thermocline thickness and strength in the global ocean. We also provide an assessment of the regions of the ocean where our adjustment is valid, i.e., where the thermal structure of the ocean follows the three-layer structure. However, there are ocean regions where the water column cannot be separated into three layers due to the dynamic processes that alter it. This assessment highlights the limitations of the existing methods to accurately determine the mixed layer depth and the thermocline depth in oceanic regions that are particularly turbulent such as the Southern Ocean and the northern North Atlantic, among others. The method proposed here has shown to be robust and easy to apply.

Unsupervised classification identifies coherent thermohaline structures in the Weddell Gyre region - Dani C. Jones; Maike Sonnewald; Shenjie Zhou; Ute Hausmann; Andrew J. S. Meijers; Isabella Rosso; Lars Boehme; Michael P. Meredith; Alberto C. Naveira Garabato (Ocean Science, 2023)
The Weddell Gyre is a major feature of the Southern Ocean and an important component of the planetary climate system; it regulates air–sea exchanges, controls the formation of deep and bottom waters, and hosts upwelling of relatively warm subsurface waters. It is characterised by low sea surface temperatures, ubiquitous sea ice formation, and widespread salt stratification that stabilises the water column. Observing the Weddell Gyre is challenging, as it is extremely remote and largely covered with sea ice. At present, it is one of the most poorly sampled regions of the global ocean, highlighting the need to extract as much value as possible from existing observations. Here, we apply a profile classification model (PCM), which is an unsupervised classification technique, to a Weddell Gyre profile dataset to identify coherent regimes in temperature and salinity. We find that, despite not being given any positional information, the PCM identifies four spatially coherent thermohaline domains that can be described as follows: (1) a circumpolar class, (2) a transition region between the circumpolar waters and the Weddell Gyre, (3) a gyre edge class with northern and southern branches, and (4) a gyre core class. PCM highlights, in an objective and interpretable way, both expected and underappreciated structures in the Weddell Gyre dataset. For instance, PCM identifies the inflow of Circumpolar Deep Water (CDW) across the eastern boundary, the presence of the Weddell–Scotia Confluence waters, and structured spatial variability in mixing between Winter Water and CDW. PCM offers a useful complement to existing expertise-driven approaches for characterising the physical configuration and variability of oceanographic regions, helping to identify coherent thermohaline structures and the boundaries between them.

High interannual surface p CO 2 variability in the southern Canadian Arctic Archipelago's Kitikmeot Sea - Richard P. Sims; Mohamed M. M. Ahmed; Brian J. Butterworth; Patrick J. Duke; Stephen F. Gonski; Samantha F. Jones; Kristina A. Brown; Christopher J. Mundy; William J. Williams; Brent G. T. Else (Ocean Science, 2023)
Warming of the Arctic due to climate change means the Arctic Ocean is now free from ice for longer, as sea ice melts earlier and refreezes later. Yet, it remains unclear how this extended ice-free period will impact carbon dioxide (CO2) fluxes due to scarcity of surface ocean CO2 measurements. Baseline measurements are urgently needed to understand spatial and temporal air–sea CO2 flux variability in the changing Arctic Ocean. There is also uncertainty as to whether the previous basin-wide surveys are representative of the many smaller bays and inlets that make up the Canadian Arctic Archipelago (CAA). By using a research vessel that is based in the remote Inuit community of Ikaluqtuutiak (Cambridge Bay, Nunavut), we have been able to reliably survey pCO2 shortly after ice melt and access previously unsampled bays and inlets in the nearby region. Here we present 4 years of consecutive summertime pCO2 measurements collected in the Kitikmeot Sea in the southern CAA. Overall, we found that this region is a sink for atmospheric CO2 in August (average of all calculated fluxes over the four cruises was −4.64 mmol m−2 d−1), but the magnitude of this sink varies substantially between years and locations (average calculated fluxes of +3.58, −2.96, −16.79 and −0.57 mmol m−2 d−1 during the 2016, 2017, 2018 and 2019 cruises, respectively). Surface ocean pCO2 varied by up to 156 µatm between years, highlighting the importance of repeat observations in this region, as this high interannual variability would not have been captured by sparse and infrequent measurements. We find that the surface ocean pCO2 value at the time of ice melt is extremely important in constraining the magnitude of the air–sea CO2 flux throughout the ice-free season. However, further constraining the air–sea CO2 flux in the Kitikmeot Sea will require a better understanding of how pCO2 changes outside of the summer season. Surface ocean pCO2 measurements made in small bays and inlets of the Kitikmeot Sea were ∼ 20–40 µatm lower than in the main channels. Surface ocean pCO2 measurements made close in time to ice breakup (i.e. within 2 weeks) were ∼ 50 µatm lower than measurements made > 4 weeks after breakup. As previous basin-wide surveys of the CAA have focused on the deep shipping channels and rarely measure close to the ice breakup date, we hypothesize that there may be an observational bias in previous studies, leading to an underestimate of the CO2 sink in the CAA. These high-resolution measurements constitute an important new baseline for gaining a better understanding of the role this region plays in the uptake of atmospheric CO2.

Validating the spatial variability in the semidiurnal internal tide in a realistic global ocean simulation with Argo and mooring data - Gaspard Geoffroy; Jonas Nycander; Maarten C. Buijsman; Jay F. Shriver; Brian K. Arbic (Ocean Science, 2023)
The autocovariance of the semidiurnal internal tide (IT) is examined in a 32 d segment of a global run of the HYbrid Coordinate Ocean Model (HYCOM). This numerical simulation, with 41 vertical layers and 1/25∘ horizontal resolution, includes tidal and atmospheric forcing, allowing for the generation and propagation of ITs to take place within a realistic eddying general circulation. The HYCOM data are in turn compared with global observations of the IT around 1000 dbar, from Argo float park-phase data and mooring records. HYCOM is found to be globally biased low in terms of the IT variance and decay of the IT autocovariance over timescales shorter than 32 d. Except in the Southern Ocean, where limitations in the model cause the discrepancy with in situ measurements to grow poleward, the spatial correlation between the Argo and HYCOM tidal variance suggests that the generation of low-mode semidiurnal ITs is globally well captured by the model.

DUACS DT2021 reprocessed altimetry improves sea level retrieval in the coastal band of the European seas - Antonio Sánchez-Román; M. Isabelle Pujol; Yannice Faugère; Ananda Pascual (Ocean Science, 2023)
More than 29 years of altimeter data have been recently reprocessed by the multi-satellite Data Unification and Altimeter Combination System (DUACS) and made available under the name of DT2021 through the Copernicus Marine Service (CMEMS) and the Copernicus Climate Change Service (C3S). New standards have been applied and various geophysical correction parameters have been updated compared to the previous release in order to improve the product quality. This paper describes the assessment of this new release through the comparison of both the all satellites and the two satellites product with external in situ tide gauge measurements in the coastal areas of the European seas for a time period from 1 January 1993 to 31 May 2020. The aim is to quantify the improvements on the previous DT2018 processing version for the retrieval of sea level in the coastal zone. The results confirmed that the CMEMS product in the new DT2021 processing version better solves the signal in the coastal band. The all satellites dataset showed a reduction of 3 % in errors when compared with tide gauges and of 5 % in the variance of the differences between the datasets compared to DT2018 reprocessing. Moreover, the all satellites dataset provided more accurate sea level measurements when making a comparison with tide gauges with respect to the climatic two satellites dataset due to the better performance of the former for the assessment of higher than climatic frequency signals. By contrast, the two satellite dataset is the most suitable product for the assessment of long-term sea level sea surface height (SSH) trends in the coastal zone due to its larger stability to the detriment of the all satellites dataset.

Seasonal overturning variability in the eastern North Atlantic subpolar gyre: a Lagrangian perspective - Oliver John Tooth; Helen Louise Johnson; Chris Wilson; Dafydd Gwyn Evans (Ocean Science, 2023)
Both observations and ocean reanalyses show a pronounced seasonality in the strength of the Atlantic meridional overturning circulation (MOC) within the eastern North Atlantic subpolar gyre (eSPG). However, attributing this overturning seasonality to seasonal dense water formation remains challenging owing to the wide distribution of recirculation timescales within the Iceland and Irminger basins. Here, we investigate the nature of seasonal overturning variability using Lagrangian water parcel trajectories initialised across the Overturning in the Subpolar North Atlantic Program (OSNAP) East section within an eddy-permitting ocean sea ice hindcast. By adopting a Lagrangian perspective, we show that the seasonal minimum of the Eulerian overturning at OSNAP East in autumn results from a combination of enhanced stratification and increased southward transport within the upper East Greenland Current. This convergence of southward transport within the MOC upper limb is explained by decreasing water parcel recirculation times in the upper Irminger Sea, consistent with a gyre-scale response to seasonal wind forcing. To account for the diversity of recirculation times within the eSPG, we also quantify the Lagrangian overturning (LMOC) as the total dense water formation along water parcel trajectories. The majority of water parcels, sourced from the central and southern branches of the North Atlantic Current, fail to return to OSNAP East prior to experiencing wintertime diapycnal transformation into the lower limb, and thus they determine the mean strength of the LMOC within the eSPG (8.9 ± 2.2 Sv). The strong seasonality of the LMOC is explained by a small collection of upper-limb water parcels, circulating rapidly (≤ 8.5 months) in the upper Irminger and central Iceland basins, whose along-stream transformation is determined by their month of arrival at OSNAP East.

Mixing and air–sea buoyancy fluxes set the time-mean overturning circulation in the subpolar North Atlantic and Nordic Seas - Dafydd Gwyn Evans; N. Penny Holliday; Sheldon Bacon; Isabela Le Bras (Ocean Science, 2023)
The overturning streamfunction as measured at the OSNAP (Overturning in the Subpolar North Atlantic Program) mooring array represents the transformation of warm, salty Atlantic Water into cold, fresh North Atlantic Deep Water (NADW). The magnitude of the overturning at the OSNAP array can therefore be linked to the transformation by air–sea buoyancy fluxes and mixing in the region north of the OSNAP array. Here, we estimate these water mass transformations using observational-based, reanalysis-based and model-based datasets. Our results highlight that air–sea fluxes alone cannot account for the time-mean magnitude of the overturning at OSNAP, and therefore a residual mixing-driven transformation is required to explain the difference. A cooling by air–sea heat fluxes and a mixing-driven freshening in the Nordic Seas, Iceland Basin and Irminger Sea precondition the warm, salty Atlantic Water, forming subpolar mode water classes in the subpolar North Atlantic. Mixing in the interior of the Nordic Seas, over the Greenland–Scotland Ridge and along the boundaries of the Irminger Sea and Iceland Basin drive a water mass transformation that leads to the convergence of volume in the water mass classes associated with NADW. Air–sea buoyancy fluxes and mixing therefore play key and complementary roles in setting the magnitude of the overturning within the subpolar North Atlantic and Nordic Seas. This study highlights that, for ocean and climate models to realistically simulate the overturning circulation in the North Atlantic, the small-scale processes that lead to the mixing-driven formation of NADW must be adequately represented within the model's parameterisation scheme.

Improving sea surface temperature in a regional ocean model through refined sea surface temperature assimilation - Silje Christine Iversen; Ann Kristin Sperrevik; Olivier Goux (Ocean Science, 2023)
Infrared (IR) and passive microwave (PMW) satellite sea surface temperature (SST) retrievals are valuable to assimilate into high-resolution regional ocean forecast models. Still, there are issues related to these SSTs that need to be addressed to achieve improved ocean forecasts. Firstly, satellite SST products tend to be biased. Assimilating SSTs from different providers can thus cause the ocean model to receive inconsistent information. Secondly, while PMW SSTs are valuable for constraining models during cloudy conditions, the spatial resolution of these retrievals is rather coarse. Assimilating PMW SSTs into high-resolution ocean models will spatially smooth the modeled SST and consequently remove finer SST structures. In this study, we implement a bias correction scheme that corrects satellite SSTs before assimilation. We also introduce a special observation operator, called the supermod operator, into the Regional Ocean Modeling System (ROMS) four-dimensional variational data assimilation algorithm. This supermod operator handles the resolution mismatch between the coarse observations and the finer model. We test the bias correction scheme and the supermod operator using a setup of ROMS covering the shelf seas and shelf break off Norway. The results show that the validation statistics in the modeled SST improve if we apply the bias correction scheme. We also find improvements in the validation statistics when we assimilate PMW SSTs in conjunction with the IR SSTs. However, our supermod operator must be activated to avoid smoothing the modeled SST structures on spatial scales smaller than twice the PMW SST footprint. Both the bias correction scheme and the supermod operator are easy to apply, and the supermod operator can easily be adapted for other observation variables.

Ocean color algorithm for the retrieval of the particle size distribution and carbon-based phytoplankton size classes using a two-component coated-sphere backscattering model - Tihomir S. Kostadinov; Lisl Robertson Lain; Christina Eunjin Kong; Xiaodong Zhang; Stéphane Maritorena; Stewart Bernard; Hubert Loisel; Daniel S. F. Jorge; Ekaterina Kochetkova; Shovonlal Roy; Bror Jonsson; Victor Martinez-Vicente; et al. (Ocean Science, 2023)
The particle size distribution (PSD) of suspended particles in near-surface seawater is a key property linking biogeochemical and ecosystem characteristics with optical properties that affect ocean color remote sensing. Phytoplankton size affects their physiological characteristics and ecosystem and biogeochemical roles, e.g., in the biological carbon pump, which has an important role in the global carbon cycle and thus climate. It is thus important to develop capabilities for measurement and predictive understanding of the structure and function of oceanic ecosystems, including the PSD, phytoplankton size classes (PSCs), and phytoplankton functional types (PFTs). Here, we present an ocean color satellite algorithm for the retrieval of the parameters of an assumed power-law PSD. The forward optical model considers two distinct particle populations: phytoplankton and non-algal particles (NAPs). Phytoplankton are modeled as coated spheres following the Equivalent Algal Populations (EAP) framework, and NAPs are modeled as homogeneous spheres. The forward model uses Mie and Aden–Kerker scattering computations, for homogeneous and coated spheres, respectively, to model the total particulate spectral backscattering coefficient as the sum of phytoplankton and NAP backscattering. The PSD retrieval is achieved via spectral angle mapping (SAM), which uses backscattering end-members created by the forward model. The PSD is used to retrieve size-partitioned absolute and fractional phytoplankton carbon concentrations (i.e., carbon-based PSCs), as well as particulate organic carbon (POC), using allometric coefficients. This model formulation also allows the estimation of chlorophyll a concentration via the retrieved PSD, as well as percent of backscattering due to NAPs vs. phytoplankton. The PSD algorithm is operationally applied to the merged Ocean Colour Climate Change Initiative (OC-CCI) v5.0 ocean color data set. Results of an initial validation effort are also presented using PSD, POC, and picophytoplankton carbon in situ measurements. Validation results indicate the need for an empirical tuning for the absolute phytoplankton carbon concentrations; however these results and comparison with other phytoplankton carbon algorithms are ambiguous as to the need for the tuning. The latter finding illustrates the continued need for high-quality, consistent, large global data sets of PSD, phytoplankton carbon, and related variables to facilitate future algorithm improvements.

Seasonal and El Niño–Southern Oscillation-related ocean variability in the Panama Bight - Rafael R. Torres; Estefanía Giraldo; Cristian Muñoz; Ana Caicedo; Ismael Hernández-Carrasco; Alejandro Orfila (Ocean Science, 2023)
In the Panama Bight, two different seasonal surface circulation patterns coincide with a strong mean sea level variation, as observed from 27 years of absolute dynamic topography (ADT) and the use of self-organizing maps. From January to April, a cyclonic gyre with a strong southwestward Panama Jet Surface Current (PJSC) dominates the basin circulation, forced by the Panama surface wind jet that also produces upwelling, reducing sea surface temperature (SST), increasing sea surface salinity (SSS) and causing an ADT decrease. From June to December, the Choco surface wind jet enhances SST, precipitation and river runoff, which reduces SSS, causing an ADT rise, which in turn forces a weak circulation in the bight, vanishing the PJSC. Interannual variability in the region is strongly affected by El Niño–Southern Oscillation (ENSO); however this climatic variability does not modify the seasonal circulation patterns in the Panama Bight. In contrast, the positive (negative) ENSO phase increases (decreases) SST and ADT in the Panama Bight, with a mean annual difference of 0.9 ∘C and 9.6 cm, respectively, between the two conditions, while its effect on SSS is small. However, as the strong seasonal SST, SSS and ADT ranges are up to 2.2 ∘C, 2.59 g kg−1 and 28.3 cm, the seasonal signal dominates over interannual variations in the bight.

Sudden, local temperature increase above the continental slope in the southern Weddell Sea, Antarctica - Elin Darelius; Vår Dundas; Markus Janout; Sandra Tippenhauer (Ocean Science, 2023)
Around most of Antarctica, the Circumpolar Deep Water (CDW) shows a warming trend. At the same time, the thermocline is shoaling, thereby increasing the potential for CDW to enter the shallow continental shelves and ultimately increase basal melt in the ice shelf cavities that line the coast. Similar trends, on the order of 0.05 ∘C and 3 m per decade, have been observed in the Warm Deep Water (WDW), the slightly cooled CDW derivative found at depth in the Weddell Sea. Here, we report on a sudden, local increase in the temperature maximum of the WDW above the continental slope north of the Filchner Trough (74∘ S, 25–40∘ W), a region identified as a hotspot for both Antarctic Bottom Water formation (AABW) and potential changes in the flow of WDW towards the large Filchner–Ronne Ice Shelf. New conductivity–temperature–depth profiles, obtained in summer 2021, and recent (2017–2021) mooring records show that the temperature of the warm-water core increased by about 0.1 ∘C over the upper part of the slope (700–2750 m depth) compared with historical (1973–2018) measurements. The temperature increase occurred relatively suddenly in late 2019 and was accompanied by an unprecedented (in observations) freshening of the overlying winter water. The AABW descending down the continental slope from Filchner Trough is sourced by dense ice shelf water and consists to a large degree (60 %) of entrained WDW. The observed temperature increase can hence be expected to imprint directly on deep-water properties, increasing the temperature of newly produced bottom water (by up to 0.06 ∘C) and reducing its density.

Multi-model analysis of the Adriatic dense-water dynamics - Petra Pranić; Cléa Denamiel; Ivica Janeković; Ivica Vilibić (Ocean Science, 2023)
This study aims to enhance our understanding of the bora-driven dense-water dynamics in the Adriatic Sea using different state-of-the-art modelling approaches during the 2014–2015 period. Practically, we analyse and compare the results of the following four different simulations: the latest reanalysis product for the Mediterranean Sea, a recently evaluated fine-resolution atmosphere–ocean Adriatic Sea climate model, and a long-time-running Adriatic Sea atmosphere–ocean forecast model used in both hindcast and data assimilation (with 4 d cycles) modes. As a first step, we evaluate the resolved physics in each simulation by focusing on the performance of the models. Then, we derive the general conditions in the ocean and the atmosphere during the investigated period. Finally, we analyse in detail the numerical reproduction of the dense-water dynamics as seen by the four simulations. The likely prerequisites for proper modelling of the ocean circulation in the Adriatic basin, including a kilometre-scale atmosphere–ocean approach, non-hydrostatic atmospheric models, fine vertical resolutions in both atmosphere and ocean, and the location and forcing of the open boundary conditions, are thus discussed in the context of the different simulations. In conclusion, a 31-year-long run of the fine-resolution Adriatic Sea climate model is found to be able to outperform most aspects of the reanalysis product, the short-term hindcast, and the data-assimilated simulation in reproducing the dense-water dynamics in the Adriatic Sea.

Response of the sea surface temperature to heatwaves during the France 2022 meteorological summer - Thibault Guinaldo; Aurore Voldoire; Robin Waldman; Stéphane Saux Picart; Hervé Roquet (Ocean Science, 2023)
The summer of 2022 was memorable and record-breaking, ranking as the second hottest summer in France since 1900, with a seasonal surface air temperature average of 22.7 ∘C. In particular, France experienced multiple record-breaking heatwaves during the meteorological summer. As the main heat reservoir of the Earth system, the oceans are at the forefront of events of this magnitude which enhance oceanic disturbances such as marine heatwaves (MHWs). In this study, we investigate the sea surface temperature (SST) of French maritime basins using remotely sensed measurements to track the response of surface waters to the atmospheric heatwaves and determine the intensity of such feedback. Beyond the direct relationship between SSTs and surface air temperatures, we explore the leading atmospheric parameters affecting the upper-layer ocean heat budget. Despite some gaps in data availability, the SSTs measured during the meteorological summer of 2022 were record-breaking, the mean SST was between 1.3 and 2.6 ∘C above the long-term average (1982–2011), and the studied areas experienced between 4 and 22 d where the basin-averaged SSTs exceeded the maximum recorded basin-averaged SSTs from 1982 to 2011. We found a significant SST response during heatwave periods with maximum temperatures measured locally at 30.8 ∘C in the north-western Mediterranean Sea. Our results show that in August 2022 (31 July to 13 August), France experienced above-average surface solar radiation correlated with below-average total cloud cover and negative wind speed anomalies. Our attribution analysis based on a simplified mixed-layer heat budget highlights the critical role of ocean–atmosphere fluxes in initiating abnormally warm SSTs, while ocean mixing plays a crucial role in the cessation of such periods. We find that the 2 m temperatures and specific humidity that are consistently linked to the advection of warm and moist air masses are key variables across all the studied regions. Our results reveal that the influence of wind on heatwaves is variable and of secondary importance. Moreover, we observe that the incident solar radiation has a significant effect only on the Bay of Biscay (BB) and the English Channel (EC) areas. Our study findings are consistent with previous research and demonstrate the vulnerability of the Mediterranean Sea to the increasing frequency of extreme weather events resulting from climate change. Furthermore, our investigation reveals that the recurring heatwave episodes during the summer of 2022 had an undeniable impact on all the surveyed maritime areas in France. Our study therefore provides valuable insights into the complex mechanisms underlying the ocean–atmosphere interaction and demonstrates the need for an efficient and sustainable operational system combining polar-orbiting and geostationary satellites to monitor the alterations that threaten the oceans in the context of climate change.

Origins of mesoscale mixed-layer depth variability in the Southern Ocean - Yu Gao; Igor Kamenkovich; Natalie Perlin (Ocean Science, 2023)
Mixed-layer depth (MLD) exhibits significant variability, which is important for atmosphere–ocean exchanges of heat and atmospheric gases. The origins of the mesoscale MLD variability in the Southern Ocean are studied here in an idealised regional ocean–atmosphere model (ROAM). The main conclusion from the analysis of the upper-ocean buoyancy budget is that, while the atmospheric forcing and oceanic vertical mixing, on average, induce the mesoscale variability of MLD, the three-dimensional oceanic advection of buoyancy counteracts and partially balances these atmosphere-induced vertical processes. The relative importance of advection changes with both season and average MLD. From January to May, when the mixed layer is shallow, the atmospheric forcing and oceanic mixing are the most important processes, with the advection playing a secondary role. From June to December, when the mixed layer is deep, both atmospheric forcing and oceanic advection are equally important in driving the MLD variability. Importantly, buoyancy advection by mesoscale ocean current anomalies can lead to both local shoaling and deepening of the mixed layer. The role of the atmospheric forcing is then directly addressed by two sensitivity experiments in which the mesoscale variability is removed from the atmosphere–ocean heat and momentum fluxes. The findings confirm that mesoscale atmospheric forcing predominantly controls MLD variability in summer and that intrinsic oceanic variability and surface forcing are equally important in winter. As a result, MLD variance increases when mesoscale anomalies in atmospheric fluxes are removed in winter, and oceanic advection becomes a dominant player in the buoyancy budget. This study highlights the importance of oceanic advection and intrinsic ocean dynamics in driving mesoscale MLD variability and underscores the importance of MLD in modulating the effects of advection on upper-ocean dynamics.

Extension of the general unit hydrograph theory for the spread of salinity in estuaries - Huayang Cai; Bo Li; Junhao Gu; Tongtiegang Zhao; Erwan Garel (Ocean Science, 2023)
From both practical and theoretical perspectives, it is essential to be able to express observed salinity distributions in terms of simplified theoretical models, which enable qualitative assessments to be made in many problems concerning water resource utilization (such as intake of fresh water) in estuaries. In this study, we propose a general and analytical salt intrusion model inspired by Guo's general unit hydrograph theory for flood hydrograph prediction in a watershed. To derive a simple, general and analytical model of salinity distribution, we first make four hypotheses on the longitudinal salinity gradient based on empirical observations; we then derive a general unit hydrograph for the salinity distribution along a partially mixed or well-mixed estuary. The newly developed model can be well calibrated using a minimum of three salinity measurements along the estuary axis and does converge towards zero when the along-estuary distance approaches infinity asymptotically. The theory has been successfully applied to reproduce the salt intrusion in 21 estuaries worldwide, which suggests that the proposed method can be a useful tool for quickly assessing the spread of salinity under a wide range of riverine and tidal conditions and for quantifying the potential impacts of human-induced and natural changes.

Physical processes and biological productivity in the upwelling regions of the tropical Atlantic - Peter Brandt; Gaël Alory; Founi Mesmin Awo; Marcus Dengler; Sandrine Djakouré; Rodrigue Anicet Imbol Koungue; Julien Jouanno; Mareike Körner; Marisa Roch; Mathieu Rouault (Ocean Science, 2023)
In this paper, we review observational and modelling results on the upwelling in the tropical Atlantic between 10∘ N and 20∘ S. We focus on the physical processes that drive the seasonal variability of surface cooling and the upward nutrient flux required to explain the seasonality of biological productivity. We separately consider the equatorial upwelling system, the coastal upwelling system of the Gulf of Guinea and the tropical Angolan upwelling system. All three tropical Atlantic upwelling systems have in common a strong seasonal cycle, with peak biological productivity during boreal summer. However, the physical processes driving the upwelling vary between the three systems. For the equatorial regime, we discuss the wind forcing of upwelling velocity and turbulent mixing, as well as the underlying dynamics responsible for thermocline movements and current structure. The coastal upwelling system in the Gulf of Guinea is located along its northern boundary and is driven by both local and remote forcing. Particular emphasis is placed on the Guinea Current, its separation from the coast and the shape of the coastline. For the tropical Angolan upwelling, we show that this system is not driven by local winds but instead results from the combined effect of coastally trapped waves, surface heat and freshwater fluxes, and turbulent mixing. Finally, we review recent changes in the upwelling systems associated with climate variability and global warming and address possible responses of upwelling systems in future scenarios.

Modelling the barotropic sea level in the Mediterranean Sea using data assimilation - Marco Bajo; Christian Ferrarin; Georg Umgiesser; Andrea Bonometto; Elisa Coraci (Ocean Science, 2023)
This paper analyses the variability of the sea level barotropic components in the Mediterranean Sea and their reproduction using a hydrodynamic model with and without data assimilation. The impact of data assimilation is considered both in reanalysis and short-forecast simulations. We used a two-dimensional finite element model paired with an ensemble Kalman filter, which assimilated hourly sea level data from 50 stations in the Mediterranean basin. The results brought about a significant improvement given by data assimilation in the reanalysis of the astronomical tide, the surge, and the barotropic total sea level, even in coastal areas and far from the assimilated stations (e.g. the southeastern Mediterranean Sea). As with the reanalysis simulations, the forecast simulations, which start from analysis states, improve, especially on the first day (37 % average error reduction) and when seiche oscillations are triggered. Since seiches are free barotropic oscillations that depend only on the initial state, their reproduction improves very effectively with data assimilation. Finally, we estimate the periods and the energy of these oscillations by means of spectral analysis, both in the Adriatic Sea, where they have been extensively studied, and in the Mediterranean Sea, where the present documentation is scarce. While the periods are well reproduced by the model even without data assimilation, their energy shows a good improvement when using it.

Joint observation–model mixed-layer heat and salt budgets in the eastern tropical Atlantic - Roy Dorgeless Ngakala; Gaël Alory; Casimir Yélognissè Da-Allada; Olivia Estelle Kom; Julien Jouanno; Willi Rath; Ezinvi Baloïtcha (Ocean Science, 2023)
In this study, we use a joint observation–model approach to investigate the mixed-layer heat and salt annual mean as well as seasonal budgets in the eastern tropical Atlantic. The regional PREFCLIM (PREFACE Climatology) observational climatology provides the budget terms with a relatively low spatial and temporal resolution compared to the online NEMO (Nucleus for European Modeling of the Ocean; Madec, G., 2014) model, and this is later resampled as in PREFCLIM climatology. In addition, advection terms are recomputed offline from the model as PREFCLIM gridded advection computation. In the Senegal, Angola, and Benguela regions, the seasonal cycle of mixed-layer temperature is mainly governed by surface heat fluxes; however, it is essentially driven by vertical heat diffusion in the equatorial region. The seasonal cycle of mixed-layer salinity is largely controlled by freshwater flux in the Senegal and Benguela regions; however, it follows the variability of zonal and meridional salt advection in the equatorial and Angola regions, respectively. Our results show that the time-averaged spatial distribution of NEMO offline heat and salt advection terms compares much better to PREFCLIM horizontal advection terms than the online heat and salt advection terms. However, the seasonal cycle of horizontal advection in selected regions shows that NEMO offline terms do not always compare well with PREFCLIM, sometimes less than online terms. Despite this difference, these results suggest the important role of small-scale variability in mixed-layer heat and salt budgets.

Global water level variability observed after the Hunga Tonga-Hunga Ha'apai volcanic tsunami of 2022 - Adam T. Devlin; David A. Jay; Stefan A. Talke; Jiayi Pan (Ocean Science, 2023)
The eruption of the Hunga Tonga-Hunga Ha'apai volcano on 15 January 2022 provided a rare opportunity to understand global tsunami impacts of explosive volcanism and to evaluate future hazards, including dangers from “volcanic meteotsunamis” (VMTs) induced by the atmospheric shock waves that followed the eruption. The propagation of the volcanic and marine tsunamis was analyzed using globally distributed 1 min measurements of air pressure and water level (WL) (from both tide gauges and deep-water buoys). The marine tsunami propagated primarily throughout the Pacific, reaching nearly 2 m at some locations, though most Pacific locations recorded maximums lower than 1 m. However, the VMT resulting from the atmospheric shock wave arrived before the marine tsunami and propagated globally, producing water level perturbations in the Indian Ocean, the Mediterranean, and the Caribbean. The resulting water level response of many Pacific Rim gauges was amplified, likely related to wave interaction with bathymetry. The meteotsunami repeatedly boosted tsunami wave energy as it circled the planet several times. In some locations, the VMT was amplified by as much as 35-fold relative to the inverse barometer due to near-Proudman resonance and topographic effects. Thus, a meteotsunami from a larger eruption (such as the Krakatoa eruption of 1883) could yield atmospheric pressure changes of 10 to 30 mb, yielding a 3–10 m near-field tsunami that would occur in advance of (usually) larger marine tsunami waves, posing additional hazards to local populations. Present tsunami warning systems do not consider this threat.

Improving statistical projections of ocean dynamic sea-level change using pattern recognition techniques - Víctor Malagón-Santos; Aimée B. A. Slangen; Tim H. J. Hermans; Sönke Dangendorf; Marta Marcos; Nicola Maher (Ocean Science, 2023)
Regional emulation tools based on statistical relationships, such as pattern scaling, provide a computationally inexpensive way of projecting ocean dynamic sea-level change for a broad range of climate change scenarios. Such approaches usually require a careful selection of one or more predictor variables of climate change so that the statistical model is properly optimized. Even when appropriate predictors have been selected, spatiotemporal oscillations driven by internal climate variability can be a large source of statistical model error. Using pattern recognition techniques that exploit spatial covariance information can effectively reduce internal variability in simulations of ocean dynamic sea level, significantly reducing random errors in regional emulation tools. Here, we test two pattern recognition methods based on empirical orthogonal functions (EOFs), namely signal-to-noise maximizing EOF pattern filtering and low-frequency component analysis, for their ability to reduce errors in pattern scaling of ocean dynamic sea-level change. We use the Max Planck Institute Grand Ensemble (MPI-GE) as a test bed for both methods, as it is a type of initial-condition large ensemble designed for an optimal characterization of the externally forced response. We show that the two methods tested here more efficiently reduce errors than conventional approaches such as a simple ensemble average. For instance, filtering only two realizations by characterizing their common response to external forcing reduces the random error by almost 60 %, a reduction that is only achieved by averaging at least 12 realizations. We further investigate the applicability of both methods to single-realization modeling experiments, including four CMIP5 simulations for comparison with previous regional emulation analyses. Pattern filtering leads to a varying degree of error reduction depending on the model and scenario, ranging from more than 20 % to about 70 % reduction in global-mean root mean squared error compared with unfiltered simulations. Our results highlight the relevance of pattern recognition methods as a tool to reduce errors in regional emulation tools of ocean dynamic sea-level change, especially when one or only a few realizations are available. Removing internal variability prior to tuning regional emulation tools can optimize the performance of the statistical model, leading to substantial differences in emulated dynamic sea level compared to unfiltered simulations.

An analogues-based forecasting system for Mediterranean marine-litter concentration - Gabriel Jordà; Javier Soto-Navarro (Ocean Science, 2023)
In this work, we explore the performance of a statistical forecasting system for marine-litter concentration in the Mediterranean Sea. In particular, we assess the potential skills of a system based on the analogues method. The system uses a historical database of marine-litter concentration simulated by a high-resolution realistic model and is trained to identify meteorological situations in the past that are similar to the forecasted ones. Then, the corresponding marine-litter concentrations of the past analogue days are used to construct the marine-litter concentration forecast. Due to the scarcity of observations, the forecasting system has been validated against a synthetic reality (i.e., the outputs from a marine-litter-modeling system). Different approaches have been tested to refine the system, and the results show that using integral definitions for the similarity function, based on the history of the meteorological situation, improves the system performance. We also find that the system accuracy depends on the domain of application being better for larger regions. Also, the method performs well in capturing the spatial patterns but performs worse in capturing the temporal variability, especially the extreme values. Despite the inherent limitations of using a synthetic reality to validate the system, the results are promising, and the approach has potential to become a suitable cost-effective forecasting method for marine-litter concentration.

Multiple mechanisms for chlorophyll a concentration variations in coastal upwelling regions: a case study east of Hainan Island in the South China Sea - Junyi Li; Min Li; Chao Wang; Quanan Zheng; Ying Xu; Tianyu Zhang; Lingling Xie (Ocean Science, 2023)
Using satellite observations from 2003 to 2020 and cruise observations from 2019 and 2021, this study reveals an unexpected minor role of upwelling in seasonal and interannual variations in chlorophyll a (Chl a) concentrations in the coastal upwelling region east of Hainan Island (UEH) in the northwestern South China Sea (NWSCS). The results show strong seasonal and interannual variability in the Chl a concentration in the core upwelling area of the UEH. Different from the strongest upwelling in summer, the Chl a concentration in the UEH area reaches a maximum of 1.18 mg m−3 in autumn and winter, with a minimum value of 0.74 mg m−3 in summer. The Chl a concentration in summer increases to as high as 1.0 mg m−3 with weak upwelling, whereas the maximum Chl a concentration in October increases to 2.5 mg m−3. The analysis of environmental factors shows that, compared to the limited effects of upwelling, the along-shelf coastal current from the northern shelf and the increased precipitation are crucially important to the Chl a concentration variation in the study area. These results provide new insights for predicting marine productivity in upwelling areas, i.e., multiple mechanisms, especially horizontal advection, should be considered in addition to the upwelling process.

Intraseasonal variability of the South Vietnam upwelling, South China Sea: influence of atmospheric forcing and ocean intrinsic variability - Marine Herrmann; Thai To Duy; Claude Estournel (Ocean Science, 2023)
The South Vietnam upwelling (SVU) develops off the Vietnamese coast (South China Sea, SCS) during the southwest summer monsoon over four main areas: the northern coastal upwelling (NCU), the southern coastal upwelling (SCU), the offshore upwelling (OFU) and the shelf off the Mekong River mouth (MKU). An ensemble of 10 simulations with perturbed initial conditions were run with the fine-resolution SYMPHONIE model (1 km inshore) to investigate the daily to intraseasonal variability of the SVU and the influence of the ocean intrinsic variability (OIV) during the strong SVU of summer 2018. The intraseasonal variability is similar for the SCU, MKU and OFU, driven to the first order by the wind variability. The MKU and SCU are induced by stable ocean dynamics (the northeastward then eastward boundary current) and have very little chaotic variability. The OIV has a stronger influence on OFU. In July, OFU mainly develops along the northern flank of the eastward jet. The influence of the OIV is strongest and related to the chaotic variability of the meridional position of the jet. In August, this position is stable and OFU develops mainly in the area of positive wind curl and cyclonic eddies north of the jet. The influence of the OIV, weaker than in July, is related to the organization of this mesoscale circulation. The NCU shows a behavior different from that observed in the other areas. In the heart of summer, a large-scale circulation formed by the eastward jet and eddy dipole is well established with an alongshore current preventing the NCU development. In early and late summer, this circulation is weaker, allowing a mesoscale circulation of strongly chaotic nature to develop in the NCU area. During those periods, the OIV influence on the NCU is very strong and related to the organization of this mesoscale circulation: the NCU is favored (annihilated) by offshore-oriented (alongshore) structures.

Current observed global mean sea level rise and acceleration estimated from satellite altimetry and the associated measurement uncertainty - Adrien Guérou; Benoit Meyssignac; Pierre Prandi; Michaël Ablain; Aurélien Ribes; François Bignalet-Cazalet (Ocean Science, 2023)
We present the latest release of the global mean sea level (GMSL) record produced by the French space agency Centre National d’Etudes Spatiales (CNES) and distributed on the AVISO+ website. This dataset is based on reprocessed along-track data, so-called L2P 21, of the reference missions TOPEX/Poseidon (TP) and Jason-1, Jason-2 and Jason-3. The L2P 21 CNES/AVISO+ GMSL record covers the period January 1993 to December 2021 and is now delivered with an estimate of its measurement uncertainties following the method presented in Ablain et al. (2019). Based on the latest calibration (Cal) and validation (Val) knowledge, we updated the uncertainty budget of the reference altimetry mission measurements and demonstrate that the CNES/AVISO+ GMSL record now achieves stability of performances of ± 0.3 mm yr−1 at the 90 % confidence level (C.L.) for its trend and ±0.05 mm yr−2 (90 % C.L.) for its acceleration over the 29 years of the altimetry record. Thanks to an analysis of the relative contribution of each measurement uncertainty budget contributor, i.e. the altimeter, the radiometer, the orbit determination and the geophysical corrections, we identified the current limiting factors to the GMSL monitoring stability and accuracy. We find that the radiometer wet troposphere correction (WTC) and the high-frequency errors with timescales shorter than 1 year are the major contributors to the GMSL measurement uncertainty over periods of 10 years (30 %–70 %), for both the trend and acceleration estimations. For longer periods of 20 years, the TP data quality is still a limitation, but more interestingly, the International Terrestrial Reference Frame (ITRF) realization uncertainties becomes dominant over all the other sources of uncertainty. Such a finding challenges the altimetry observing system as it is designed today and highlights clear topics of research to be explored in the future to help the altimetry community to improve the GMSL measurement accuracy and stability.

A simple diagnostic based on sea surface height with an application to central Pacific ENSO - Jufen Lai; Richard J. Greatbatch; Martin Claus (Ocean Science, 2023)
We use output from a freely running NEMO model simulation for the equatorial Pacific to investigate the utility of linearly removing the local influence of vertical displacements of the thermocline from variations in sea surface height. We show that the resulting time series of residual sea surface height, denoted ηnlti, measures variations in near-surface heat content that are independent of the local vertical displacement of the thermocline and can arise from horizontal advection, surface heat flux, and diapycnal mixing processes. We find that the variance of ηnlti and its correlation with sea surface temperature are focused on the Niño4 region. Furthermore, ηnlti averaged over the Niño4 region is highly correlated with indices of central Pacific El Niño–Southern Oscillation (CP ENSO), and its variance in 21-year running windows shows a strong upward trend over the past 50 years, corresponding to the emergence of CP ENSO following the 1976/77 climate shift. We show that ηnlti can be estimated from observations, using satellite altimeter data and a linear multi-mode model. The time series of ηnlti, especially when estimated using the linear model, show pronounced westward propagation in the western equatorial Pacific, arguing for an important role for zonal advective feedback in the dynamics of CP ENSO, in particular for cold events. We also present evidence that the role of the thermocline displacement in influencing sea surface height increased strongly after 2000 in the eastern part of the Niño4 region, at a time when CP ENSO was particularly active. Finally, the diagnostic is easy to compute and can be easily applied to mooring data or coupled climate models.

The sensitivity of primary productivity in Disko Bay, a coastal Arctic ecosystem, to changes in freshwater discharge and sea ice cover - Eva Friis Møller; Asbjørn Christensen; Janus Larsen; Kenneth D. Mankoff; Mads Hvid Ribergaard; Mikael Sejr; Philip Wallhead; Marie Maar (Ocean Science, 2023)
The Greenland ice sheet is melting, and the rate of ice loss has increased 6-fold since the 1980s. At the same time, the Arctic sea ice extent is decreasing. Meltwater runoff and sea ice reduction both influence light and nutrient availability in the coastal ocean, with implications for the timing, distribution, and magnitude of phytoplankton production. However, the integrated effect of both glacial and sea ice melt is highly variable in time and space, making it challenging to quantify. In this study, we evaluate the relative importance of these processes for the primary productivity of Disko Bay, west Greenland, one of the most important areas for biodiversity and fisheries around Greenland. We use a high-resolution 3D coupled hydrodynamic–biogeochemical model for 2004–2018 validated against in situ observations and remote sensing products. The model-estimated net primary production (NPP) varied between 90–147 gC m−2 yr−1 during 2004–2018, a period with variable freshwater discharges and sea ice cover. NPP correlated negatively with sea ice cover and positively with freshwater discharge. Freshwater discharge had a strong local effect within ∼ 25 km of the source-sustaining productive hot spots during summer. When considering the annual NPP at bay scale, sea ice cover was the most important controlling factor. In scenarios with no sea ice in spring, the model predicted a ∼ 30 % increase in annual production compared to a situation with high sea ice cover. Our study indicates that decreasing ice cover and more freshwater discharge can work synergistically and will likely increase primary productivity of the coastal ocean around Greenland.

Water mass transformation variability in the Weddell Sea in ocean reanalyses - Shanice T. Bailey; C. Spencer Jones; Ryan P. Abernathey; Arnold L. Gordon; Xiaojun Yuan (Ocean Science, 2023)
This study investigates the variability of water mass transformation (WMT) within the Weddell Gyre (WG). The WG serves as a pivotal site for the Meridional Overturning Circulation (MOC) and ocean ventilation because it is the primary origin of the largest volume of water mass in the global ocean: Antarctic Bottom Water (AABW). Recent mooring data suggest substantial seasonal and interannual variability of AABW properties exiting the WG, and studies have linked the variability to the large-scale climate forcings affecting wind stress in the WG region. However, the specific thermodynamic mechanisms that link variability in surface forcings to variability in water mass transformations and AABW export remain unclear. This study explores how current state-of-the-art data-assimilating ocean reanalyses can help fill the gaps in our understanding of the thermodynamic drivers of AABW variability in the WG via WMT volume budgets derived from Walin's classic WMT framework. The three ocean reanalyses used are the following: Estimating the Circulation and Climate of the Ocean state estimate (ECCOv4), Southern Ocean State Estimate (SOSE) and Simple Ocean Data Assimilation (SODA). From the model outputs, we diagnose a closed form of the water mass budget for AABW that explicitly accounts for transport across the WG boundary, surface forcing, interior mixing and numerical mixing. We examine the annual mean climatology of the WMT budget terms, the seasonal climatology and finally the interannual variability. Our finding suggests that the relatively coarse resolution of these models did not realistically capture AABW formation, export and variability. In ECCO and SOSE, we see strong interannual variability in AABW volume budget. In SOSE, we find an accelerating loss of AABW during 2005–2010, driven largely by interior mixing and changes in surface salt fluxes. ECCO shows a similar trend during a 4-year time period starting in late 2007 but also reveals such trends to be part of interannual variability over a much longer time period. Overall, ECCO provides the most useful time series for understanding the processes and mechanisms that drive WMT and export variability in the WG. SODA, in contrast, displays unphysically large variability in AABW volume, which we attribute to its data assimilation scheme. We also examine correlations between the WMT budgets and large-scale climate indices, including El Niño–Southern Oscillation (ENSO) and Southern Annular Mode (SAM), and find no strong relationships.

Global submesoscale diagnosis using along-track satellite altimetry - Oscar Vergara; Rosemary Morrow; Marie-Isabelle Pujol; Gérald Dibarboure; Clément Ubelmann (Ocean Science, 2023)
The ocean's sea surface height (SSH) field is a complex mix of motions in geostrophic balance and unbalanced motions including high-frequency tides, internal tides, and internal gravity waves. Barotropic tides are well estimated for altimetric SSH in the open ocean, but the SSH signals of internal tides remain. The transition scale, Lt, at which these unbalanced ageostrophic motions dominate balanced geostrophic motions is estimated for the first time using satellite altimetry. Lt is critical to define the spatial scales above which surface geostrophic currents can be inferred from SSH gradients. We use a statistical approach based on the analysis of 1 Hz altimetric SSH wavenumber spectra to obtain four geophysical parameters that vary regionally and seasonally: the background error, the spectral slope in the mesoscale range, a second spectral slope at smaller scales, and Lt. The mesoscale slope and error levels are similar to previous studies based on satellite altimetry. The break in the wavenumber spectra to a flatter spectral slope can only be estimated in midlatitude regions where the signal exceeds the altimetric noise level. Small values of Lt are observed in regions of energetic mesoscale activity, while larger values are observed towards low latitudes and regions of lower mesoscale activity. These results are consistent with recent analyses of in situ observations and high-resolution models. Limitations of our results and implications for reprocessed nadir and future swath altimetric missions are discussed.

The signature of NAO and EA climate patterns on the vertical structure of the Canary Current upwelling system - Tina Georg; Maria C. Neves; Paulo Relvas (Ocean Science, 2023)
The current study aims to analyse the vertical structure of the ocean during upwelling events using in situ and modelled data. Additionally, the influence of climate patterns, namely the North Atlantic Oscillation (NAO) and the East Atlantic (EA) pattern, on the vertical structure and their impact on the upwelling activity are assessed for a period of 25 years (1993–2017). The study focuses on the central part of the Canary Current (25–35∘ N) with persistent upwelling throughout the year, with an annual cycle and the strongest events from June to September. Upwelling is determined using two different approaches: one index is calculated based on temperature differences between the coastal and the offshore area, and the other is calculated based on wind data and the resulting Ekman transport. Different datasets were chosen according to the indices. Stable coastal upwelling can be observed in the study area for the analysed time span, with differences throughout the latitudes. A deepening of the isothermal layer depth and a cooling of temperatures are observed in the vertical structure of coastal waters, representing a deeper mixing of the ocean and the rise of cooler, denser water towards the surface. During years of a positive NAO, corresponding to a strengthening of the Azores High and the Icelandic Low, stronger winds lead to an intensification of the upwelling activity, an enhanced mixing of the upper ocean, and a deeper (shallower) isothermal layer along the coast (offshore). The opposite is observed in years of negative NAO. Both effects are enhanced in years with a coupled, opposite phase of the EA pattern and are mainly visible during winter months, where the effect of both indices is the greatest. The study therefore suggests that upwelling activities are stronger in winters of positive North Atlantic Oscillation coupled with a negative East Atlantic pattern and emphasizes the importance of interactions between the climate patterns and upwelling.

Surface circulation in the Gulf of Thailand from remotely sensed observations: seasonal and interannual timescales - Arachaporn Anutaliya (Ocean Science, 2023)
The Gulf of Thailand (GoT), a shallow semi-enclosed basin located in the western equatorial Pacific, undergoes much wind variabilities on both seasonal and interannual timescales that produce complex surface circulation. The local Ekman pumping modifies sea level in the northern GoT, while remote wind forcing influences sea level variability at the GoT western boundary, potentially through the coastal trapped Kelvin waves. The importance of the Ekman current on ageostrophic current is also important; the stronger influence of the Ekman current is found toward the southern part of the GoT. The GoT circulation reverses its direction seasonally following the monsoon wind reversal which is well-captured by the most dominant complex empirical orthogonal function explaining 28 % of the total circulation variance. During the monsoon transition, a strong meridional current along the western boundary that connects to the flow at the GoT southeastern entrance is observed. This implies high exchange between the GoT and the South China Sea and thus modification of the GoT water. On the interannual timescale, the GoT circulation is directly impacted by both the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD). Interestingly, the two climate modes have different spatial influences on the GoT circulation. The IOD dominates the interannual current along the GoT western boundary and the southern boundary of the observing domain (8∘ N), while the ENSO correlates with that in the interior. The results highlight the complex circulation pattern as being contributed by different dynamics over each region of the GoT.

Revisiting the global mean ocean mass budget over 2005–2020 - Anne Barnoud; Julia Pfeffer; Anny Cazenave; Robin Fraudeau; Victor Rousseau; Michaël Ablain (Ocean Science, 2023)
We investigate the performances of Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GRACE-FO) satellite gravimetry missions in assessing the ocean mass budget at the global scale over 2005–2020. For that purpose, we focus on the last years of the record (2015–2020) when GRACE and GRACE Follow-On faced instrumental problems. We compare the global mean ocean mass estimates from GRACE and GRACE Follow-On to the sum of its contributions from Greenland, Antarctica, land glaciers, terrestrial water storage and atmospheric water content estimated with independent observations. Significant residuals are observed in the global mean ocean mass budget at interannual timescales. Our analyses suggest that the terrestrial water storage variations based on global hydrological models likely contribute in large part to the misclosure of the global mean ocean mass budget at interannual timescales. We also compare the GRACE-based global mean ocean mass with the altimetry-based global mean sea level corrected for the Argo-based thermosteric contribution (an equivalent of global mean ocean mass). After correcting for the wet troposphere drift of the radiometer on board the Jason-3 altimeter satellite, we find that mass budget misclosure is reduced but still significant. However, replacing the Argo-based thermosteric component by the Ocean Reanalysis System 5 (ORAS5) or from the Clouds and the Earth's Radiant Energy System (CERES) top of the atmosphere observations significantly reduces the residuals of the mass budget over the 2015–2020 time span. We conclude that the two most likely sources of error in the global mean ocean mass budget are the thermosteric component based on Argo and the terrestrial water storage contribution based on global hydrological models. The GRACE and GRACE Follow-On data are unlikely to be responsible on their own for the non-closure of the global mean ocean mass budget.

Effects of including the adjoint sea ice rheology on estimating Arctic Ocean–sea ice state - Guokun Lyu; Armin Koehl; Xinrong Wu; Meng Zhou; Detlef Stammer (Ocean Science, 2023)
The adjoint assimilation method has been applied to coupled ocean and sea ice models for sensitivity studies and Arctic state estimations. However, the accuracy of the adjoint model is degraded by simplifications of the adjoint of the sea ice model, especially the adjoint sea ice rheologies. As part of ongoing developments in coupled ocean and sea ice estimation systems, we incorporate and approximate the adjoint of viscous-plastic sea ice dynamics (adjoint-VP) and compare it with the adjoint of free-drift sea ice dynamics (adjoint-FD) through assimilation experiments. Using the adjoint-VP results in a further cost reduction of 7.9 % in comparison to adjoint-FD, with noticeable improvements in the ocean temperature over the open water and the intermediate layers of the Arctic Ocean. Adjoint-VP adjusts the model input more efficiently than adjoint-FD does by involving different sea ice retreat processes. For instance, adjoint-FD melts the sea ice up to 1.0 m in the marginal seas from May to June by overadjusting air temperature (>8 ∘C); adjoint-VP reproduces the sea ice retreat with smaller adjustments to the atmospheric state within their prior uncertainty range. These developments of the adjoint model here lay the foundation for further improving Arctic Ocean and sea ice estimations by comprehensively adjusting the initial conditions, atmospheric forcings, and parameters of the model.

A numerical investigation on the energetics of a current along an ice-covered continental slope - Hengling Leng; Hailun He; Michael A. Spall (Ocean Science, 2023)
The Chukchi Slope Current is a westward-flowing current along the Chukchi slope, which carries Pacific-origin water from the Chukchi shelf into the Canada Basin and helps set the regional hydrographic structure and ecosystem. Using a set of experiments with an idealized primitive equation numerical model, we investigate the energetics of the slope current during the ice-covered period. Numerical calculations show that the growth of surface eddies is suppressed by the ice friction, while perturbations at mid-depths can grow into eddies, consistent with linear instability analysis. However, because the ice stress is spatially variable, it is able to drive Ekman pumping to decrease the available potential energy (APE) and kinetic energy of both the mean flow and mesoscale eddies over a vertical scale of 100 m, well outside the frictional Ekman layer. The rate at which the APE changes is determined by the vertical density flux, which is negative as the ice-induced Ekman pumping advects lighter (denser) water upward (downward). A scaling analysis shows that Ekman pumping will dominate the release of APE for large-scale flows, but the effect of baroclinic instability is also important when the horizontal scale of the mean flow is the baroclinic deformation radius and the eddy velocity is comparable to the mean flow velocity. Our numerical results highlight the importance of ice friction in the energetics of the slope current and eddies, and this may be relevant to other ice-covered regions.

Assimilation of sea surface salinities from SMOS in an Arctic coupled ocean and sea ice reanalysis - Jiping Xie; Roshin P. Raj; Laurent Bertino; Justino Martínez; Carolina Gabarró; Rafael Catany (Ocean Science, 2023)
In the Arctic, the sea surface salinity (SSS) plays a key role in processes related to water mixing and sea ice. However, the lack of salinity observations causes large uncertainties in Arctic Ocean forecasts and reanalysis. Recently the Soil Moisture and Ocean Salinity (SMOS) satellite mission was used by the Barcelona Expert Centre to develop an Arctic SSS product. In this study, we evaluate the impact of assimilating this data in a coupled ocean–ice data assimilation system. Using the deterministic ensemble Kalman filter from July to December 2016, two assimilation runs respectively assimilated two successive versions of the SMOS SSS product on top of a pre-existing reanalysis run. The runs were validated against independent in situ salinity profiles in the Arctic. The results show that the biases and the root-mean-squared differences (RMSD) of SSS are reduced by 10 % to 50 % depending on the area and highlight the importance of assimilating satellite salinity data. The time series of freshwater content (FWC) further shows that its seasonal cycle can be adjusted by assimilation of the SSS products, which is encouraging of the assimilation of SSS in a long-time reanalysis to better reproduce the Arctic water cycle.

Revisiting the tropical Atlantic western boundary circulation from a 25-year time series of satellite altimetry data - Djoirka Minto Dimoune; Florence Birol; Fabrice Hernandez; Fabien Léger; Moacyr Araujo (Ocean Science, 2023)
Geostrophic currents derived from altimetry are used to investigate the surface circulation in the western tropical Atlantic over the 1993–2017 period. Using six horizontal sections defined to capture the current branches of the study area, we investigate their respective variations at both seasonal and interannual timescales, as well as the spatial distribution of these variations, in order to highlight the characteristics of the currents on their route. Our results show that the central branch of the South Equatorial Current and its northern branch near the Brazilian coast, the North Brazil Current component located south of the Equator, and the Guyana Current have similar annual cycles, with maxima (minima) during late boreal winter (boreal fall) when the Intertropical Convergence Zone is at its southernmost (northernmost) location. In contrast, the seasonal cycles of the North Brazil Current branch located between the Equator and 7–8∘ N, its retroflected branch, the northern branch of the South Equatorial Current to the west of 35∘ W, and the North Equatorial Countercurrent show maxima (minima) during late boreal summer (boreal spring), following the remote wind stress curl strength variation. West of 32∘ W, an eastward current (the Equatorial Surface Current, ESC) is observed between 2–2∘ N, identified as the equatorial extension of the retroflected branch of the North Brazil Current. It is part of a large cyclonic circulation observed between 0–6∘ N and 35–45∘ W during boreal spring. We also observed a secondary North Brazil Current retroflection flow during the second half of the year, which leads to the two-core structure of the North Equatorial Countercurrent and might be related to the wind stress curl seasonal changes. To the east, the North Equatorial Countercurrent weakens and its two-core structure is underdeveloped due to the weakening of the wind stress. At interannual scales, depending on the side of the Equator examined, the North Brazil Current exhibits two opposite scenarios related to the phases of the tropical Atlantic Meridional Mode. At 32∘ W, the interannual variability of the North Equatorial Countercurrent and of the northern branch of the South Equatorial Current (in terms of both strength and/or latitudinal shift) are associated with the Atlantic Meridional Mode, whereas the variability of the Equatorial Surface Current intensity is associated with both the Atlantic Meridional Mode and Atlantic Zonal Mode phases.

How subsurface and double-core anticyclones intensify the winter mixed-layer deepening in the Mediterranean Sea - Alexandre Barboni; Solange Coadou-Chaventon; Alexandre Stegner; Briac Le Vu; Franck Dumas (Ocean Science, 2023)
The mixed layer is the uppermost layer of the ocean, connecting the atmosphere to the subsurface ocean through atmospheric fluxes. It is subject to pronounced seasonal variations: it deepens in winter due to buoyancy loss and shallows in spring while heat flux increases and restratifies the water column. A mixed-layer depth (MLD) modulation over this seasonal cycle has been observed within mesoscale eddies. Taking advantage of the numerous Argo floats deployed and trapped within large Mediterranean anticyclones over the last decades, we reveal for the first time this modulation at a 10 d temporal scale, free of the smoothing effect of composite approaches. The analysis of 16 continuous MLD time series inside 13 long-lived anticyclones at a fine temporal scale brings to light the importance of the eddy pre-existing vertical structure in setting the MLD modulation by mesoscale eddies. Extreme MLD anomalies of up to 330 m are observed when the winter mixed layer connects with a pre-existing subsurface anticyclonic core, greatly accelerating mixed-layer deepening. The winter MLD sometimes does not achieve such connection but homogenizes another subsurface layer, then forming a multi-core anticyclone with spring restratification. An MLD restratification delay is always observed, reaching more than 2 months in 3 out the 16 MLD time series. The water column starts to restratify outside anticyclones, while the mixed layer keeps deepening and cooling at the eddy core for a longer time. These new elements provide new keys for understanding anticyclone vertical-structure formation and evolution.

The effects of Hurricane Harvey on Texas coastal-zone chemistry - Piers Chapman; Steven F. DiMarco; Anthony H. Knap; Antonietta Quigg; Nan D. Walker (Ocean Science, 2023)
Hurricane Harvey deposited over 90×109 m3 of rainwater over central Texas, USA, during late August/early September 2017. During four cruises (June, August, September and November 2017) we observed changes in hydrography and nutrient and oxygen concentrations in Texas coastal waters. Despite intense terrestrial runoff, nutrient supply to the coastal ocean was transient, with little phytoplankton growth observed and no hypoxia. Observations suggest this was probably related to the retention of nutrients in the coastal bays and rapid uptake by phytoplankton of nutrients washed out of the bays, as well as dilution by the sheer volume of rainwater and the lack of significant carbon reserves in the sediments, despite the imposition of a strong pycnocline. By the November cruise conditions had apparently returned to normal, and no long-term effects were observed.

A turbulence data reduction scheme for autonomous and expendable profiling floats - Kenneth G. Hughes; James N. Moum; Daniel L. Rudnick (Ocean Science, 2023)
Autonomous and expendable profiling-float arrays such as those deployed in the Argo Program require the transmission of reliable data from remote sites. However, existing satellite data transfer rates preclude complete transmission of rapidly sampled turbulence measurements. It is therefore necessary to reduce turbulence data on board. Here we propose a scheme for onboard data reduction and test it with existing turbulence data obtained with a modified SOLO-II profiling float. First, voltage spectra are derived from shear probe and fast-thermistor signals. Then, we focus on a fixed-frequency band that we know to be unaffected by vibrations and that approximately corresponds to a wavenumber band of 5–25 cpm. Over the fixed-frequency band, we make simple power law fits that – after calibration and correction in post-processing – yield values for the turbulent kinetic energy dissipation rate ϵ and thermal-variance dissipation rate χ. With roughly 1 m vertical segments, this scheme reduces the necessary data transfer volume 300-fold to approximately 2.5 kB for every 100 m of a profile (when profiling at 0.2 m s−1). As a test, we apply our scheme to a dataset comprising 650 profiles and compare its output to that from our standard turbulence-processing algorithm. For ϵ, values from the two approaches agree within a factor of 2 87 % of the time; for χ, they agree 78 % of the time. These levels of agreement are greater than or comparable to that between the ϵ and χ values derived from two shear probes and two fast thermistors, respectively, on the same profiler.

Observation-based estimates of volume, heat, and freshwater exchanges between the subpolar North Atlantic interior, its boundary currents, and the atmosphere - Sam C. Jones; Neil J. Fraser; Stuart A. Cunningham; Alan D. Fox; Mark E. Inall (Ocean Science, 2023)
The Atlantic Meridional Overturning Circulation (AMOC) transports heat and salt between the tropical Atlantic and Arctic oceans. The interior of the North Atlantic subpolar gyre (SPG) is responsible for the much of the water mass transformation in the AMOC, and the export of this water to intensified boundary currents is crucial for projecting air–sea interaction onto the strength of the AMOC. However, the magnitude and location of exchange between the SPG and the boundary remains unclear. We present a novel climatology of the SPG boundary using quality-controlled CTD (conductivity–temperature–depth) and Argo hydrography, defining the SPG interior as the oceanic region bounded by 47∘ N and the 1000 m isobath. From this hydrography we find geostrophic flow out of the SPG around much of the boundary with minimal seasonality. The horizontal density gradient is reversed around western Greenland, where the geostrophic flow is into the SPG. Surface Ekman forcing drives net flow out of the SPG in all seasons with pronounced seasonality, varying between 2.45 ± 0.73 Sv in the summer and 7.70 ± 2.90 Sv in the winter. We estimate heat advected into the SPG to be between 0.14 ± 0.05 PW in the winter and 0.23 ± 0.05 PW in the spring, and freshwater advected out of the SPG to be between 0.07 ± 0.02 Sv in the summer and 0.15 ± 0.02 Sv in the autumn. These estimates approximately balance the surface heat and freshwater fluxes over the SPG domain. Overturning in the SPG varies seasonally, with a minimum of 6.20 ± 1.40 Sv in the autumn and a maximum of 10.17 ± 1.91 Sv in the spring, with surface Ekman the most likely mediator of this variability. The density of maximum overturning is at 27.30 kg m−3, with a second, smaller maximum at 27.54 kg m−3. Upper waters (σ0 27.54 kg m−3) exiting to the south. Our results support the present consensus that the formation and pre-conditioning of Subpolar Mode Water in the north-eastern Atlantic is a key determinant of AMOC strength.

On the ocean's response to enhanced Greenland runoff in model experiments: relevance of mesoscale dynamics and atmospheric coupling - Torge Martin; Arne Biastoch (Ocean Science, 2023)
Increasing Greenland Ice Sheet melting is anticipated to impact water mass transformation in the subpolar North Atlantic and ultimately the meridional overturning circulation. Complex ocean and climate models are widely applied to estimate magnitude and timing of related impacts under global warming. We discuss the role of the ocean mean state, subpolar water mass transformation, mesoscale eddies, and atmospheric coupling in shaping the response of the subpolar North Atlantic Ocean to enhanced Greenland runoff. In a suite of eight dedicated 60- to 100-year-long model experiments with and without atmospheric coupling, with eddy processes parameterized and explicitly simulated and with regular and significantly enlarged Greenland runoff, we find (1) a major impact by the interactive atmosphere in enabling a compensating temperature feedback, (2) a non-negligible influence by the ocean mean state biased towards greater stability in the coupled simulations, both of which make the Atlantic meridional overturning circulation less susceptible to the freshwater perturbation applied, and (3) a more even spreading and deeper mixing of the runoff tracer in the subpolar North Atlantic and enhanced inter-gyre exchange with the subtropics in the strongly eddying simulations. Overall, our experiments demonstrate the important role of mesoscale ocean dynamics and atmosphere feedback in projections of the climate system response to enhanced Greenland Ice Sheet melting and hence underline the necessity to advance scale-aware eddy parameterizations for next-generation climate models.

Seasonal cycle of sea surface temperature in the tropical Angolan Upwelling System - Mareike Körner; Peter Brandt; Marcus Dengler (Ocean Science, 2023)
The Angolan shelf system represents a highly productive ecosystem. Throughout the year the sea surface is cooler near the coast than further offshore. The lowest sea surface temperature (SST), strongest cross-shore temperature gradient, and maximum productivity occur in austral winter when seasonally prevailing upwelling-favourable winds are weakest. Here, we investigate the seasonal mixed layer heat budget to identify atmospheric and oceanic causes for heat content variability. By using different satellite and in situ data, we derive monthly estimates of surface heat fluxes, mean horizontal advection, and local heat content change. We calculate the heat budgets for the near-coastal and offshore regions separately to explore processes that lead to the observed SST differences. The results show that the net surface heat flux warms the coastal ocean stronger than further offshore, thus acting to damp spatial SST differences. Mean horizontal heat advection is dominated by meridional advection of warm water along the Angolan coast. However, its contribution to the heat budget is small. Ocean turbulence data suggest that the heat flux, due to turbulent mixing across the base of the mixed layer, is an important cooling term. This turbulent cooling, being strongest in shallow shelf regions, is capable of explaining the observed negative cross-shore temperature gradient. The residuum of the mixed layer heat budget and uncertainties of budget terms are discussed.

Unifying biological field observations to detect and compare ocean acidification impacts across marine species and ecosystems: what to monitor and why - Steve Widdicombe; Kirsten Isensee; Yuri Artioli; Juan Diego Gaitán-Espitia; Claudine Hauri; Janet A. Newton; Mark Wells; Sam Dupont (Ocean Science, 2023)
Approximately one-quarter of the CO2 emitted to the atmosphere annually from human activities is absorbed by the ocean, resulting in a reduction of seawater pH and shifts in seawater carbonate chemistry. This multi-decadal process, termed “anthropogenic ocean acidification” (OA), has been shown to have detrimental impacts on marine ecosystems. Recent years have seen a globally coordinated effort to measure the changes in seawater chemistry caused by OA, with best practices now available for these measurements. In contrast to these substantial advances in observing physicochemical changes due to OA, quantifying their biological consequences remains challenging, especially from in situ observations under real-world conditions. Results from 2 decades of controlled laboratory experiments on OA have given insight into the likely processes and mechanisms by which elevated CO2 levels affect biological process, but the manifestation of these process across a plethora of natural situations has yet to be fully explored. This challenge requires us to identify a set of fundamental biological and ecological indicators that are (i) relevant across all marine ecosystems, (ii) have a strongly demonstrated link to OA, and (iii) have implications for ocean health and the provision of ecosystem services with impacts on local marine management strategies and economies. This paper draws on the understanding of biological impacts provided by the wealth of previous experiments, as well as the findings of recent meta-analyses, to propose five broad classes of biological indicators that, when coupled with environmental observations including carbonate chemistry, would allow the rate and severity of biological change in response to OA to be observed and compared. These broad indicators are applicable to different ecological systems, and the methods for data analysis suggested here would allow researchers to combine biological response data across regional and global scales by correlating rates of biological change with the rate of change in carbonate chemistry parameters. Moreover, a method using laboratory observation to design an optimal observing strategy (frequency and duration) and observe meaningful biological rates of change highlights the factors that need to be considered when applying our proposed observation strategy. This innovative observing methodology allows inclusion of a wide diversity of marine ecosystems in regional and global assessments and has the potential to increase the contribution of OA observations from countries with developing OA science capacity.

Extension of Ekman (1905) wind-driven transport theory to the β plane - Nathan Paldor; Lazar Friedland (Ocean Science, 2023)
The seminal Ekman (1905) f-plane theory of wind-driven transport at the ocean surface is extended to the β plane by substituting the pseudo-angular momentum for the zonal velocity in the Lagrangian equation. When the β term is added, the equations become nonlinear, which greatly complicates the analysis. Though rotation relates the momentum equations in the zonal and the meridional directions, the transformation to pseudo-angular momentum greatly simplifies the longitudinal dynamics, which yields a clear description of the meridional dynamics in terms of a slow drift compounded by fast oscillations; this can then be applied to describe the motion in the zonal direction. Both analytical expressions and numerical calculations highlight the critical role of the Equator in determining the trajectories of water columns forced by eastward-directed (in the Northern Hemisphere) wind stress even when the water columns are initiated far from the Equator. Our results demonstrate that the averaged motion in the zonal direction depends on the amplitude of the meridional oscillations and is independent of the direction of the wind stress. The zonal drift is determined by a balance between the initial conditions and the magnitude of the wind stress, so it can be as large as the mean meridional motion; i.e., the averaged flow direction is not necessarily perpendicular to the wind direction.

Turbulent kinetic energy dissipation rate and associated fluxes in the western tropical Atlantic estimated from ocean glider observations - Peter M. F. Sheehan; Gillian M. Damerell; Philip J. Leadbitter; Karen J. Heywood; Rob A. Hall (Ocean Science, 2023)
Ocean gliders enable us to collect the high-resolution microstructure observations necessary to calculate the dissipation rate of turbulent kinetic energy, ε, on timescales of weeks to months: far longer than is normally possible using traditional ship-based platforms. Slocum gliders have previously been used to this end; here, we report the first detailed estimates of ε calculated using the Batchelor spectrum method on observations collected by a FP07 fast thermistor mounted on a Seaglider. We use these same fast thermistor observations to calculate ε following the Thorpe scale method and find very good agreement between the two methods. The Thorpe scale method yields larger values of ε, but the average difference, which is less than an order of magnitude, is smaller than reported elsewhere. The spatio-temporal distribution of ε is comparable for both methods. Maximum values of ε (10−7 W kg−1) are observed in the surface mixed layer; values of approximately 10−9 W kg−1 are observed between approximately 200 and 500 m depth. These two layers are separated by a 100 m thick layer of low ε (10−10 W kg−1), which is co-located with a high-salinity layer of Subtropical Underwater and a peak in the strength of stratification. We calculate the turbulent heat and salt fluxes associated with the observed turbulence. Between 200 and 500 m, ε induces downward fluxes of both properties that, if typical of the annual average, would have a very small influence on the heat and salt content of the overlying salinity-maximum layer. We compare these turbulent fluxes with two estimates of double-diffusive fluxes that occur in regions susceptible to salt fingers, such as the western tropical Atlantic. We find that the double-diffusive fluxes of both heat and salt are larger than the corresponding turbulent fluxes.

River effects on sea-level rise in the Río de la Plata estuary during the past century - Christopher G. Piecuch (Ocean Science, 2023)
Identifying the causes for historical sea-level changes in coastal tide-gauge records is important for constraining oceanographic, geologic, and climatic processes. The Río de la Plata estuary in South America features the longest tide-gauge records in the South Atlantic. Despite the relevance of these data for large-scale circulation and climate studies, the mechanisms underlying relative sea-level changes in this region during the past century have not been firmly established. I study annual data from tide gauges in the Río de la Plata and stream gauges along the Río Paraná and Río Uruguay to establish relationships between river streamflow and sea level over 1931–2014. Regression analysis suggests that streamflow explains 59 %±17 % of the total sea-level variance at Buenos Aires, Argentina, and 28 %±21 % at Montevideo, Uruguay (95 % confidence intervals). A long-term streamflow increase effected sea-level trends of 0.71±0.35 mm yr−1 at Buenos Aires and 0.48±0.38 mm yr−1 at Montevideo. More generally, sea level at Buenos Aires and Montevideo respectively rises by (7.3±1.8)×10-6 m and (4.7±2.6)×10-6 m per 1 m3 s−1 streamflow increase. These observational results are consistent with simple theories for the coastal sea-level response to streamflow forcing, suggesting a causal relationship between streamflow and sea level mediated by ocean dynamics. Findings advance understanding of local, regional, and global sea-level changes; clarify sea-level physics; inform future projections of coastal sea level and the interpretation of satellite data and proxy reconstructions; and highlight future research directions. Specifically, local and regional river effects should be accounted for in basin-scale and global mean sea-level budgets as well as reconstructions based on sparse tide-gauge records.

Fortnightly variability of Chl a in the Indonesian seas - Edward D. Zaron; Tonia A. Capuano; Ariane Koch-Larrouy (Ocean Science, 2023)
Twenty years of daily MODIS-Aqua ocean color observations (2002–2022) are used to identify periodic variability of near-surface chlorophyll (Chl a) in the Indonesian seas. The frequency spectrum of Chl a is dominated by the mean and low-frequency monsoonal variability; however, a prominent peak around the fortnightly tidal period, MSf, is present. Harmonic analysis is used to quantify and map the fortnightly Chl a signal, which is discovered to be significant along the continental shelves of NW Australia and at several sites associated with narrow passages between the Lesser Sunda Islands, within the Sulu Archipelago, and at a few other sites in the Philippines Archipelago. Fortnightly variability at the shallow coastal sites is attributed to the spring–neap cycle of barotropic ocean currents, while we hypothesize that the variability in deeper water near the island passages is due to the modulation of vertical nutrient fluxes by baroclinic tidal mixing. The results document the significance of tidal mixing and highlight the heterogeneous character of biophysical processes within the Indonesian seas.

Regionalizing the sea-level budget with machine learning techniques - Carolina M. L. Camargo; Riccardo E. M. Riva; Tim H. J. Hermans; Eike M. Schütt; Marta Marcos; Ismael Hernandez-Carrasco; Aimée B. A. Slangen (Ocean Science, 2023)
Attribution of sea-level change to its different drivers is typically done using a sea-level budget approach. While the global mean sea-level budget is considered closed, closing the budget on a finer spatial scale is more complicated due to, for instance, limitations in our observational system and the spatial processes contributing to regional sea-level change. Consequently, the regional budget has been mainly analysed on a basin-wide scale. Here we investigate the sea-level budget at sub-basin scales, using two machine learning techniques to extract domains of coherent sea-level variability: a neural network approach (self-organizing map, SOM) and a network detection approach (δ-MAPS). The extracted domains provide more spatial detail within the ocean basins and indicate how sea-level variability is connected among different regions. Using these domains we can close, within 1σ uncertainty, the sub-basin regional sea-level budget from 1993–2016 in 100 % and 76 % of the SOM and δ-MAPS regions, respectively. Steric variations dominate the temporal sea-level variability and determine a significant part of the total regional change. Sea-level change due to mass exchange between ocean and land has a relatively homogeneous contribution to all regions. In highly dynamic regions (e.g. the Gulf Stream region) the dynamic mass redistribution is significant. Regions where the budget cannot be closed highlight processes that are affecting sea level but are not well captured by the observations, such as the influence of western boundary currents. The use of the budget approach in combination with machine learning techniques leads to new insights into regional sea-level variability and its drivers.

Dimethyl sulfide cycling in the sea surface microlayer in the southwestern Pacific – Part 1: Enrichment potential determined using a novel sampler - Alexia D. Saint-Macary; Andrew Marriner; Theresa Barthelmeß; Stacy Deppeler; Karl Safi; Rafael Costa Santana; Mike Harvey; Cliff S. Law (Ocean Science, 2023)
Elevated dimethyl sulfide (DMS) concentrations in the sea surface microlayer (SML) have been previously related to DMS air–sea flux anomalies in the southwestern Pacific. To further address this, DMS, its precursor dimethylsulfoniopropionate (DMSP), and ancillary variables were sampled in the SML and also subsurface water at 0.5 m depth (SSW) in different water masses east of New Zealand. Despite high phytoplankton biomass at some stations, the SML chlorophyll a enrichment factor (EF) was low (

Imminent reversal of the residual flow through the Marsdiep tidal inlet into the Dutch Wadden Sea based on multiyear ferry-borne acoustic Doppler current profiler (ADCP) observations - Johan van der Molen; Sjoerd Groeskamp; Leo R. M. Maas (Ocean Science, 2022)
The Dutch Wadden Sea is a UN World Heritage Site connected to the North Sea by multiple tidal inlets. Although there are strong tidal currents flowing through these inlets, the magnitude and direction of the residual circulation in the western Dutch Wadden Sea is important for sediment, salinity and nutrient balances. We found that the direction of this residual flow is reversing. This residual circulation has been the subject of various studies since the 1970s, in which substantially different net volume fluxes were presented. Differences in tidal conditions in the main inlets, tidal rectification and meteorology were identified as driving mechanisms. Here we analysed almost 13 years of acoustic Doppler current profiler (ADCP) observations collected on the ferry crossing the Marsdiep tidal inlet in the Dutch Wadden Sea since 2009. The results are combined with earlier investigations covering the period 1998–2009. We find a significant trend in the magnitude of the residual volume flux, with decreasing export to the North Sea and with occasional imports observed in recent years. We hypothesise that this trend is related predominantly to changes in tides in the North Sea, which are caused by increased strength and duration of stratification in response to global warming. With warming projected to continue, we expect the residual flow in the Marsdiep to continue to reverse to full inflow within the current decade, with potential knock-on effects for the sediment budget and ecosystem of the western Wadden Sea.

Reconstruction of Mediterranean coastal sea level at different timescales based on tide gauge records - Jorge Ramos-Alcántara; Damià Gomis; Gabriel Jordà (Ocean Science, 2022)
A coastal sea level reconstruction based on tide gauge observations is developed and applied to the western basin of the Mediterranean sea. The reconstructions are carried out in four frequency bands and are based on an optimal interpolation method in which the correlation between tide gauge data and all coastal points has been determined from the outputs of a numerical model. The reconstructions for frequencies lower than 1 month use monthly observations from the Permanent Service for Mean Sea Level (PSMSL) database and cover the period from 1884 to 2019. For the reconstruction of higher frequencies, hourly observations from the Global Extreme Sea Level Analysis (GESLA-2) dataset are used and cover from 1980 to 2015. Total sea level is retrieved with high accuracy from the merging of the different frequency bands. Results of a cross-validation test show that independent tide gauge series are highly correlated with the reconstructions. Moreover, they correlate significantly better with the reconstructions than with altimetry data in all frequency bands, and therefore the reconstruction represents a valuable contribution to the attempts of recovering coastal sea level. The obtained reconstructions allow us to characterize the coastal sea level variability, estimate coastal sea level trends along the entire coastline, and examine the correlation between western Mediterranean coastal sea level and the main North Atlantic climate indices. The limitations and applicability of the method to other regions are also discussed.

Planktonic cnidarian responses to contrasting thermohaline and circulation seasonal scenarios in a tropical western boundary current system - Everton Giachini Tosetto; Arnaud Bertrand; Sigrid Neumann-Leitão; Alex Costa da Silva; Miodeli Nogueira Júnior (Ocean Science, 2022)
In western boundary current systems (WBCSs), strong currents flow coastward carrying oceanic water masses and their associated planktonic fauna. Variation in the intensity of these currents and in the continental runoff may affect the dynamic interplay between oceanic and coastal communities. In addition, changes in the continental runoff and the thermohaline structure modulate the primary production, adding complexity to the dynamics of these oligotrophic systems. These dynamics likely shape the planktonic cnidarian communities. To further understand such relationships, we used a comprehensive dataset encompassing samples collected above the shelf and slope and around oceanic seamounts and islands of the Fernando de Noronha Ridge in the western tropical South Atlantic, in two seasons characterised by distinct thermohaline structure and circulation patterns. Results show that in the tropical South Atlantic and, likely, other western boundary systems with narrow continental shelves, coastward currents spread oceanic waters and their associated cnidarian species over the continental shelf. However, while both coastal and oceanic communities co-occur when the continental runoff is notable, oceanic species dominate almost the entire shelf during the dry season characterised by a stronger boundary current intensity. We also conclude that when the mixed-layer depth and associated nutricline are shallower, the enhanced primary productivity supports larger populations of planktonic cnidarian species through a bottom–up control.

Interannual to decadal sea level variability in the subpolar North Atlantic: the role of propagating signals - Denis L. Volkov; Claudia Schmid; Leah Chomiak; Cyril Germineaud; Shenfu Dong; Marlos Goes (Ocean Science, 2022)
The gyre-scale, dynamic sea surface height (SSH) variability signifies the spatial redistribution of heat and freshwater in the ocean, influencing the ocean circulation, weather, climate, sea level, and ecosystems. It is known that the first empirical orthogonal function (EOF) mode of the interannual SSH variability in the North Atlantic exhibits a tripole gyre pattern, with the subtropical gyre varying out of phase with both the subpolar gyre and the tropics, influenced by the low-frequency North Atlantic Oscillation. Here, we show that the first EOF mode explains the majority (60 %–90 %) of the interannual SSH variance in the Labrador and Irminger Sea, whereas the second EOF mode is more influential in the northeastern part of the subpolar North Atlantic (SPNA), explaining up to 60 %–80 % of the regional interannual SSH variability. We find that the two leading modes do not represent physically independent phenomena. On the contrary, they evolve as a quadrature pair associated with a propagation of SSH anomalies from the eastern to the western SPNA. This is confirmed by the complex EOF analysis, which can detect propagating (as opposed to stationary) signals. The analysis shows that it takes about 2 years for sea level signals to propagate from the Iceland Basin to the Labrador Sea, and it takes 7–10 years for the entire cycle of the North Atlantic SSH tripole to complete. The observed westward propagation of SSH anomalies is linked to shifting wind forcing patterns and to the cyclonic pattern of the mean ocean circulation in the SPNA. The analysis of regional surface buoyancy fluxes in combination with the upper-ocean temperature and salinity changes suggests a time-dependent dominance of either air–sea heat fluxes or advection in driving the observed SSH tendencies, while the contribution of surface freshwater fluxes (precipitation and evaporation) is negligible. We demonstrate that the most recent cooling and freshening observed in the SPNA since about 2010 were mostly driven by advection associated with the North Atlantic Current. The results of this study indicate that signal propagation is an important component of the North Atlantic SSH tripole, as it applies to the SPNA.

Modelling the impact of anthropogenic measures on saltwater intrusion in the Weser estuary - Pia Kolb; Anna Zorndt; Hans Burchard; Ulf Gräwe; Frank Kösters (Ocean Science, 2022)
The Weser estuary has been subject to profound changes in topography in the past 100 years through natural variations and river engineering measures, leading to strong changes in hydrodynamics. These changes are also expected to have affected the dynamics of saltwater intrusion. Using numerical modelling, we examined saltwater intrusion in the Weser estuary in four different system states (1966, 1972, 1981, 2012). Models of each system state were set up with the respective topography and boundary values. We calibrated and validated each model individually to account for differences in sediments, bedforms, and the resolution of underlying bathymetric data between historical and recent system states. In simulations of 1 hydrological year, each with realistic forcing (hindcasting study), the influence of topography is overshadowed by the effects of other factors, particularly river discharge. At times of identical discharge, results indicate a landward shift of the salinity front between 1966 and 2012. Subsequent simulations with different topographies but identical boundary conditions (scenario study) confirm that topography changes in the Weser estuary affected saltwater intrusion. Solely through the topography changes, at a discharge of 300 m3 s−1, the position of the tidally averaged and depth-averaged salinity front shifted landwards by about 2.5 km between 1972 and 1981 and by another 1 km between 1981 and 2012. These changes are significant but comparatively small, since due to seasonal variations in run-off, the tidally averaged saltwater intrusion can vary by more than 20 km. An analysis of the salt flux through a characteristic cross section showed that saltwater intrusion in the Weser estuary is primarily driven by tidal pumping and only to a lesser degree due to estuarine circulation. However, results indicate that the contribution of individual processes has changed in response to anthropogenic measures.

Modelling floating riverine litter in the south-eastern Bay of Biscay: a regional distribution from a seasonal perspective - Irene Ruiz; Anna Rubio; Ana J. Abascal; Oihane C. Basurko (Ocean Science, 2022)
Although rivers contribute to the flux of litter to the marine environment, estimates of riverine litter amounts and detailed studies on floating riverine litter behaviour once it has reached the sea are still scarce. This paper provides an analysis of the seasonal behaviour of floating marine litter released by rivers within the south-eastern Bay of Biscay based on riverine litter characterizations, drifters, and high-frequency radar observations and Lagrangian simulations. Virtual particles were released in the coastal area as a proxy of the floating fraction of riverine litter entering from rivers and reaching the open waters. Particles were parameterized with a wind drag coefficient (Cd) to represent their trajectories and fate according to the buoyancy of the litter items. They were forced with numerical winds and measured currents provided by high-frequency radars covering selected seasonal week-long periods between 2009 and 2021. To gain a better insight into the type and buoyancy of the items, samples collected from a barrier placed at the Deba River (Spain) were characterized at the laboratory. Items were grouped into two categories: low-buoyancy items (objects not exposed to wind forcing, e.g. plastic bags) and highly buoyant items (objects highly exposed to wind forcing, e.g. bottles). Overall, low-buoyancy items encompassed almost 90 % by number and 68 % by weight. Weakly buoyant items were parameterized with Cd = 0 % and highly buoyant items with Cd = 4 %; this latter value is the result of the joint analysis of modelled and observed trajectories of four satellite drifting buoys released at the Adour (France), Deba (Spain), and Oria (Spain) river mouths. Particles parameterized with Cd = 4 % drifted faster towards the coast through the wind, notably during the first 24 h. In summer, over 97 % of particles beached after 1 week of simulation. In autumn this value fell to 54 %. In contrast, low-buoyancy items took longer to arrive at the shoreline, particularly during spring with fewer than 25 % of particles beached by the end of the simulations. The highest concentrations (>200 particles km−1) were recorded during summer for Cd = 4 % in the French region of Pyrénées-Atlantiques. Results showed that the regions in the study area were highly affected by rivers within or nearby the region itself. These results couple observations and a river-by-river modelling approach and can assist decision-makers on setting emergency responses to high fluxes of floating riverine litter and on defining future monitoring strategies for heavily polluted regions within the south-eastern Bay of Biscay.

Quantifying the impacts of the Three Gorges Dam on the spatial–temporal water level dynamics in the upper Yangtze River estuary - Huayang Cai; Hao Yang; Pascal Matte; Haidong Pan; Zhan Hu; Tongtiegang Zhao; Guangliang Liu (Ocean Science, 2022)
Understanding the alterations in spatial–temporal water level dynamics caused by natural and anthropogenic changes is essential for water resources management in estuaries, as this can directly impact the estuarine morphology, sediment transport, salinity intrusion, navigation conditions, and other factors. Here, we propose a simple triple linear regression model linking the water level variation on a daily timescale to the hydrodynamics at both ends of an estuary. The model was applied to the upper Yangtze River estuary (YRE) to examine the influence of the world's largest dam, the Three Gorges Dam (TGD), on the spatial–temporal water level dynamics within the estuary. It is shown that the regression model can accurately reproduce the water level dynamics in the upper YRE, with a root mean squared error (RMSE) of 0.061–0.150 m seen at five gauging stations for both the pre- and post-TGD periods. This confirms the hypothesis that the response of water level dynamics to hydrodynamics at both ends is mostly linear in the upper YRE. The regression model calibrated during the pre-TGD period was used to reconstruct the water level dynamics that would have occurred in the absence of the TGD's freshwater regulation. Results show that the spatial–temporal alterations in water levels during the post-TGD period are mainly driven by the variation in freshwater discharge due to the regulation of the TGD, which results in increased discharge during the dry season (from December to March) and a dramatic reduction in discharge during the wet-to-dry transitional period. The presented method to quantify the separate contributions made by changes in boundary conditions and geometry to spatial–temporal water level dynamics is particularly useful for determining scientific strategies for sustainable water resources management in dam-controlled or climate-driven estuaries worldwide.

Accuracy of numerical wave model results: application to the Atlantic coasts of Europe - Matias Alday; Fabrice Ardhuin; Guillaume Dodet; Mickael Accensi (Ocean Science, 2022)
Numerical wave models are generally less accurate in the coastal ocean than offshore. It is generally suspected that a number of factors specific to coastal environments can be blamed for these larger model errors: complex shoreline and topography, relatively short fetches, combination of remote swells and local wind seas, less accurate wind fields, presence of strong currents, bottom friction, etc. These factors generally have strong local variations, making it all the more difficult to adapt a particular model setup from one area to another. Here we investigate a wide range of modeling choices including forcing fields, spectral resolution, and parameterizations of physical processes in a regional model that covers most of the Atlantic and North Sea coasts. The effects of these choices on the model results are analyzed with buoy spectral data and wave parameter time series. Additionally, satellite altimeter data are employed to provide a more complete performance assessment of the modeled wave heights as a function of the distance to the coast and to identify areas where wave propagation is influenced by bottom friction. We show that the accurate propagation of waves from offshore is probably the most important factor on exposed shorelines, while other specific effects can be important locally, including winds, currents, and bottom friction.

Seasonal variation in eddy activity and associated heat/salt transport in the Bay of Bengal based on satellite, Argo, and 3D reprocessed data - Wei Cui; Jie Zhang; Jungang Yang (Ocean Science, 2022)
Based on satellite altimetry data spanning over 26 years in combination with Argo profile data or three-dimensional (3D) reprocessed thermohaline fields, the eddy synthesis method was used to construct vertical temperature and salinity structures of eddies in the Bay of Bengal, and the seasonal thermohaline properties of eddies and the heat and salt transport by eddies were analyzed. Analysis revealed that mesoscale eddy activities and the vertical thermohaline structures in the Bay of Bengal have evident seasonal variation. Temperature anomalies caused by eddies are usually between ±1 and ±3 ∘C (positive for anticyclonic eddies (AEs) and negative for cyclonic eddies (CEs)), and the magnitude varies seasonally. Salinity anomalies caused by eddies are small and disturbance signals in the southern bay due to the small vertical gradient of salinity there; salinity anomalies in the northern bay are generally between ±0.2 and ±0.3 psu, negative for AEs and positive for CEs. Owing to seasonal changes in both the eddy activity and the vertical thermohaline structure in the Bay of Bengal, the eddy-induced heat and salt transport in different seasons also changes substantially. Generally, high heat and salt transport is concentrated in eddy-rich regions, e.g., the western, northwestern, and eastern parts of the bay, the seas to the east of Sri Lanka, and the region to the southeast outside of the bay. The southern part of the bay shows weak salt transport owing to the inconsistent salinity signal within eddies. The result of the divergence of eddy heat transport illustrates that the 10–20 W m−2 value of the eddy-induced heat flux is comparable in magnitude with the annual mean air–sea net heat flux in the Bay of Bengal. Compared with the large-scale net heat flux and freshwater flux at the surface, the eddy-induced heat/freshwater transport can contribute substantially to regional and basin-scale heat/freshwater variability. This work provides data that could support further research on the heat and salt balance of the entire Bay of Bengal.

Ensemble quantification of short-term predictability of the ocean dynamics at a kilometric-scale resolution: a Western Mediterranean test case - Stephanie Leroux; Jean-Michel Brankart; Aurélie Albert; Laurent Brodeau; Jean-Marc Molines; Quentin Jamet; Julien Le Sommer; Thierry Penduff; Pierre Brasseur (Ocean Science, 2022)
We investigate the predictability properties of the ocean dynamics using an ensemble of short-term numerical regional ocean simulations forced by prescribed atmospheric conditions. In that purpose, we developed a kilometric-scale, regional model for the Western Mediterranean sea (MEDWEST60, at 1/60∘ horizontal resolution). A probabilistic approach is then followed, where a stochastic parameterization of model uncertainties is introduced in this setup to initialize ensemble predictability experiments. A set of three ensemble experiments (20 members and 2 months) are performed, one with the deterministic model initiated with perturbed initial conditions and two with the stochastic model, for two different amplitudes of stochastic model perturbations. In all three experiments, the spread of the ensemble is shown to emerge from the smallest scales (kilometric scale) and progressively upscales to the largest structures. After 2 months, the ensemble variance saturates over most of the spectrum, and the small scales (

Internal tides off the Amazon shelf during two contrasted seasons: interactions with background circulation and SSH imprints - Michel Tchilibou; Ariane Koch-Larrouy; Simon Barbot; Florent Lyard; Yves Morel; Julien Jouanno; Rosemary Morrow (Ocean Science, 2022)
The Amazon shelf break is a key region for internal tide (IT) generation. It also shows a large seasonal variation in circulation and associated stratification. This study, based on a high-resolution model (1/36∘) explicitly forced by tide, aims to better characterize how the ITs vary between two contrasted seasons. During the season from March to July (MAMJJ) the currents and mesoscale eddies are weak while the pycnocline is shallower and stronger. From August to December (ASOND) mean currents and mesoscale eddies are strong, and the pycnocline is deeper and weaker than in MAMJJ. For both seasons, semi-diurnal M2 ITs are generated on the shelf break mainly between the 100 and 1000 m isobath in the model. South of 2∘ N, the conversion from barotropic to baroclinic tide is more efficient in MAMJJ than in ASOND. Local dissipation of the coherent M2 at the generation sites is higher in MAMJJ (30 %) than in ASOND (22 %), because higher modes are favourably generated (modes 2 and 3), making the internal wave packet more dispersive. The remaining fraction (70 %–80 %) propagates away from the generation sites and mainly dissipates locally every ∼ 100 km, which corresponds to the mode 1 reflection beams. About 13 %, 30 %, and 40 % of the M2 coherent IT dissipates at the first, second, and third beams. M2 coherent baroclinic flux propagates more northward during MAMJJ while it seems to be blocked at 6∘ N during ASOND. There is no intensified dissipation of the coherent M2 that could explain the disappearance of the coherent flux. In fact, the flux at this location becomes more incoherent because of strong interaction with the currents. This has been shown in the paper using 25 h mean snapshots of the baroclinic flux that shows branching and stronger eastward deviation of the IT when interacting with mesoscale eddies and stratification during ASOND. Finally, we evaluated sea surface height (SSH) frequency and wavenumber spectra for subtidal (f 1/11 h−1) frequencies. Tidal frequencies explain most of the SSH variability for wavelengths between 250 and 70 km. Below 70 km, the SSH is mainly incoherent and supertidal. The length scale at which the SSH becomes dominated by unbalanced (non-geostrophic) IT was estimated to be around 250 km. Our results highlight the complexity of correctly predicting IT SSH in order to better observe mesoscale and sub-mesoscale eddies from existing and upcoming altimetric missions, notably the Surface Water Ocean Topography (SWOT) mission.

A numerical study of near-inertial motions in the Mid-Atlantic Bight area induced by Hurricane Irene (2011) - Peida Han; Xiping Yu (Ocean Science, 2022)
Hurricane Irene generated strong near-inertial currents (NICs) in ocean waters when passing over the Mid-Atlantic Bight (MAB) of the US East Coast in late August 2011. It is demonstrated that a combination of valuable field data and detailed model results can be taken advantage of to study the development and decay mechanism of this event. Numerical results obtained with the Regional Oceanic Modeling System (ROMS) are shown to agree well with the field data. Both computed and observed results show that the NICs were significant in most areas of the MAB region except in the nearshore area where the stratification was totally destroyed by the hurricane-induced strong mixing. Based on the energy budget, it is clarified that the near-inertial kinetic energy (NIKE) was mainly gained from the wind power during the hurricane event. In the deepwater region, NIKE was basically balanced by the vertical turbulence diffusion (40 %) and downward divergence (33 %), while in the continental shelf region, NIKE was mainly dissipated by the vertical turbulence diffusion (67 %) and partially by the bottom friction (24 %). Local dissipation of NIKE due to turbulence diffusion is much more closely related to the rate of the vertical shear rather than the intensity of turbulence. The strong vertical shear at the offshore side of the continental shelf led to a rapid dissipation of NIKE in this region.

Dimethyl sulfide cycling in the sea surface microlayer in the southwestern Pacific – Part 2: Processes and rates - Alexia D. Saint-Macary; Andrew Marriner; Stacy Deppeler; Karl A. Safi; Cliff S. Law (Ocean Science, 2022)
As the sea surface microlayer (SML) is the uppermost oceanic layer and differs in biogeochemical composition to the underlying subsurface water (SSW), it is important to determine whether processes in the SML modulate gas exchange, particularly for climate active gases. Enrichment of dimethyl sulfide (DMS) and its precursor dimethylsulfoniopropionate (DMSP) has been reported in the SML, but it remains unclear how this is maintained whilst DMS is lost to the atmosphere. To examine this, a comprehensive study of DMS source and sink processes, including production, consumption, and net response to irradiance, was carried out in deck-board incubations of SML water at five locations in different water masses in the southwestern Pacific east of New Zealand. Net consumption of DMSP and production of DMS in the light and dark occurred at all sites. The net response of DMS and DMSP to irradiance varied between stations but was always lower than conversion of DMSP to DMS in the dark. In addition, DMS photolytic turnover was slower than reported elsewhere, which was unexpected given the high light exposure in the SML incubations. Although no relationships were apparent between DMS process rates and biogeochemical variables, including chlorophyll a, bacteria, and phytoplankton groups, net bacterial DMSP consumption was correlated with DMSP and DMS concentrations and also dinoflagellate and Gymnodinium spp. biomass, supporting the findings of a companion study that dinoflagellates play an important role in DMS cycling in the SML. However, net DMS production rates and accumulation were low relative to regional air–sea DMS loss, indicating that DMS cycling within the SML is unlikely to influence regional DMS emissions.

The role of oceanic heat flux in reducing thermodynamic ice growth in Nares Strait and promoting earlier collapse of the ice bridge - Sergei Kirillov; Igor Dmitrenko; David G. Babb; Jens K. Ehn; Nikolay Koldunov; Søren Rysgaard; David Jensen; David G. Barber (Ocean Science, 2022)
The ice bridge in Nares Strait is a well-known phenomenon that affects the liquid and solid freshwater flux from the Arctic Ocean through the strait and controls the downstream North Open Water polynya in northern Baffin Bay. Recently, the ice bridge has been in a state of decline, either breaking up earlier in the year or not forming at all and thereby increasing the sea ice export out of the Arctic Ocean. The decline in the ice bridge has been ascribed to thinner and therefore weaker ice from the Arctic Ocean entering Nares Strait; however, local forcing also affects the state of the ice bridge and thereby influences when it breaks up. Using a variety of remotely sensed data we examine the spatial patterns of sea ice thickness within the ice bridge, highlighting the presence of negative ice thickness anomalies on both the eastern and western sides of the strait and identifying a recurrent sensible heat polynya that forms within the ice bridge near Cape Jackson in northwestern Greenland. Using the sea ice–ocean model FESOM2, we then attribute these ice thickness anomalies to the heat from warmer subsurface waters of Pacific and Atlantic origin that reduce thermodynamic ice growth throughout winter on the western and eastern sides, respectively. The consequently weaker and thinner areas within the ice bridge are then suggested to promote instability and earlier breakup. This work provides new insight into the structure of the Nares Strait ice bridge and highlights that warming of the modified Atlantic and/or Pacific waters that enter the strait may contribute to its further decline.

Exceptional freshening and cooling in the eastern subpolar North Atlantic caused by reduced Labrador Sea surface heat loss - Alan D. Fox; Patricia Handmann; Christina Schmidt; Neil Fraser; Siren Rühs; Alejandra Sanchez-Franks; Torge Martin; Marilena Oltmanns; Clare Johnson; Willi Rath; N. Penny Holliday; Arne Biastoch; et al. (Ocean Science, 2022)
Observations of the eastern subpolar North Atlantic in the 2010s show exceptional freshening and cooling of the upper ocean, peaking in 2016 with the lowest salinities recorded for 120 years. Published theories for the mechanisms driving the freshening include: reduced transport of saltier, warmer surface waters northwards from the subtropics associated with reduced meridional overturning; shifts in the pathways of fresher, cooler surface water from the Labrador Sea driven by changing patterns of wind stress; and the eastward expansion of the subpolar gyre. Using output from a high-resolution hindcast model simulation, we propose that the primary cause of the exceptional freshening and cooling is reduced surface heat loss in the Labrador Sea. Tracking virtual fluid particles in the model backwards from the eastern subpolar North Atlantic between 1990 and 2020 shows the major cause of the freshening and cooling to be an increased outflow of relatively fresh and cold surface waters from the Labrador Sea; with a minor contribution from reduced transport of warmer, saltier surface water northward from the subtropics. The cooling, but not the freshening, produced by these changing proportions of waters of subpolar and subtropical origin is mitigated by reduced along-track heat loss to the atmosphere in the North Atlantic Current. We analyse modelled boundary exchanges and water mass transformation in the Labrador Sea to show that since 2000, while inflows of lighter surface waters remain steady, the increasing output of these waters is due to reduced surface heat loss in the Labrador Sea beginning in the early 2000s. Tracking particles further upstream reveals that the primary source of the increased volume of lighter water transported out of the Labrador Sea is increased recirculation of water, and therefore longer residence times, in the upper 500–1000 m of the subpolar gyre.

Surface circulation properties in the eastern Mediterranean emphasized using machine learning methods - Georges Baaklini; Roy El Hourany; Milad Fakhri; Julien Brajard; Leila Issa; Gina Fifani; Laurent Mortier (Ocean Science, 2022)
The eastern Mediterranean surface circulation is highly energetic and composed of structures interacting stochastically. However, some main features are still debated, and the behavior of some fine-scale dynamics and their role in shaping the general circulation is yet unknown. In the following paper, we use an unsupervised neural network clustering method to analyze the long-term variability of the different mesoscale structures. We decompose 26 years of altimetric data into clusters reflecting different circulation patterns of weak and strong flows with either strain or vortex-dominated velocities. The vortex-dominated cluster is more persistent in the western part of the basin, which is more active than the eastern part due to the strong flow along the coast, interacting with the extended bathymetry and engendering continuous instabilities. The cluster that reflects a weak flow dominated the middle of the basin, including the Mid-Mediterranean Jet (MMJ) pathway. However, the temporal analysis shows a frequent and intermittent occurrence of a strong flow in the middle of the basin, which could explain the previous contradictory assessment of MMJ existence using in-situ observations. Moreover, we prove that the Levantine Sea is becoming more and more energetic as the activity of the main mesoscale features is showing a positive trend.

Detecting the most effective cleanup locations using network theory to reduce marine plastic debris: a case study in the Galapagos Marine Reserve - Stefanie L. Ypma; Quinten Bohte; Alexander Forryan; Alberto C. Naveira Garabato; Andy Donnelly; Erik van Sebille (Ocean Science, 2022)
The Galapagos Marine Reserve was established in 1986 to ensure protection of the islands' unique biodiversity. Unfortunately, the islands are polluted by marine plastic debris and the island authorities face the challenge to effectively remove plastic from its shorelines owing to limited resources. To optimize efforts, we have developed a methodology to identify the most effective cleanup locations on the Galapagos Islands using network theory. A network is constructed from a Lagrangian simulation describing the flow of macroplastic between the various islands within the Galapagos Marine Reserve, where the nodes represent locations along the coastline and the edges the likelihood of plastic leaving one location and beaching at another. We have found four network centralities that provide the best coastline ranking to optimize the cleanup effort based on various impact metrics. Locations with a high retention rate are particularly favorable for cleanup. The results indicate that using the most effective centrality for finding cleanup locations is a good strategy for heavily polluted regions if the distribution of marine plastic debris on the coastlines is unknown and limited cleanup resources are available.

A clustering approach to determine biophysical provinces and physical drivers of productivity dynamics in a complex coastal sea - Tereza Jarníková; Elise M. Olson; Susan E. Allen; Debby Ianson; Karyn D. Suchy (Ocean Science, 2022)
The balance between ocean mixing and stratification influences primary productivity through light limitation and nutrient supply in the euphotic ocean. Here, we apply a hierarchical clustering algorithm (Ward's method) to four factors relating to stratification (wind energy, freshwater index, water-column-averaged vertical eddy diffusivity, and halocline depth), as well as to depth-integrated phytoplankton biomass, extracted from a biophysical ocean model of the Salish Sea. Running the clustering algorithm on 4 years of model output, we identify distinct regions of the model domain that exhibit contrasting wind and freshwater input dynamics, as well as regions of varying water-column-averaged vertical eddy diffusivity and halocline depth regimes. The spatial regionalizations in physical variables are similar in all 4 analyzed years. We also find distinct interannually consistent biological zones. In the northern Strait of Georgia and Juan de Fuca Strait, a deeper winter halocline and episodic summer mixing coincide with higher summer diatom abundance, while in the Fraser River stratified central Strait of Georgia, shallower haloclines and stronger summer stratification coincide with summer flagellate abundance. Cluster-based model results and evaluation suggest that the Juan de Fuca Strait supports more biomass than previously thought. Our approach elucidates probable physical mechanisms controlling phytoplankton abundance and composition. It also demonstrates a simple, powerful technique for finding structure in large datasets and determining boundaries of biophysical provinces.

Major sources of North Atlantic Deep Water in the subpolar North Atlantic from Lagrangian analyses in an eddy-rich ocean model - Jörg Fröhle; Patricia V. K. Handmann; Arne Biastoch (Ocean Science, 2022)
The North Atlantic Deep Water (NADW) is a crucial component of the Atlantic meridional overturning circulation and is therefore an important factor of the climate system. In order to estimate the mean relative contributions, sources, and pathways of the NADW at the southern exit of the Labrador Sea, a Lagrangian particle experiment is performed. The particles were seeded according to the strength of the velocity field along the 53∘ N section and traced 40 years backward in time in the three-dimensional velocity and hydrography field. The resulting transport pathways, their sources and corresponding transit timescales were inferred. Our experiment shows that, of the 30.1 Sv of NADW passing 53∘ N on average, the majority of this water is associated with a diapycnal mass flux without contact to the atmosphere, accounting for 14.3 Sv (48 %), where 6.2 Sv originate from the Labrador Sea, compared to 4.7 Sv from the Irminger Sea. The second-largest contribution originates from the mixed layer with 7.2 Sv (24 %), where the Labrador Sea contribution (5.9 Sv) dominates over the Irminger Sea contribution (1.0 Sv). Another 5.7 Sv (19 %) of NADW crosses the Greenland–Scotland Ridge within the NADW density class, where about two-thirds pass the Denmark Strait, while one-third crosses the Iceland–Scotland Ridge. The NADW exported at 53∘ N is hence dominated by entrainment through the diapycnal mass flux and mixed-layer origin in the Labrador Sea.

Upper-ocean response to the passage of tropical cyclones in the Azores region - Miguel M. Lima; Célia M. Gouveia; Ricardo M. Trigo (Ocean Science, 2022)
Tropical cyclones (TCs) are extreme climate events that are known to strongly interact with the ocean through two mechanisms: dynamically through the associated intense wind stress and thermodynamically through moist enthalpy exchanges at the ocean surface. These interactions contribute to relevant oceanic responses during and after the passage of a TC, namely the induction of a cold wake and the production of chlorophyll (Chl a) blooms. This study aimed to understand these interactions in the Azores region, an area with relatively low cyclonic activity for the North Atlantic basin, since the area experiences much less intense events than the rest of the basin. Results for the 1998–2020 period showed that the averaged induced anomalies were on the order of +0.050 mg m−3 for Chl a and −1.615 ∘C for SST (sea surface temperature). Furthermore, looking at the role played by several TCs characteristics we found that the intensity of the TCs was the most important condition for the development of upper-ocean responses. Additionally, it was found that bigger TCs caused greater induced anomalies in both variables, while faster ones created greater Chl a responses, and TCs that occurred later in the season had greater TC-related anomalies. Two case studies (Ophelia in 2017 and Nadine in 2012) were conducted to better understand each upper-ocean response. Ophelia was shown to affect the SST at an earlier stage, while the biggest Chl a induced anomalies were registered at a later stage, allowing the conclusion that thermodynamic exchanges conditioned the SST more while dynamical mixing might have played a more important role in the later stage. Nadine showed the importance of the TC track geometry, revealing that the TC track observed in each event can impact a specific region for longer and therefore result in greater induced anomalies.

Import of Atlantic Water and sea ice controls the ocean environment in the northern Barents Sea - Øyvind Lundesgaard; Arild Sundfjord; Sigrid Lind; Frank Nilsen; Angelika H. H. Renner (Ocean Science, 2022)
The northern Barents Sea is a cold, seasonally ice-covered Arctic shelf sea region that has experienced major warming and sea ice loss in recent decades. Here, a 2-year observational record from two ocean moorings provides new knowledge about the seasonal hydrographic variability in the region and about the ocean exchange across its northern margin. The combined records of temperature, salinity, and currents show the advection of warmer and saltier waters of Atlantic origin into the Barents Sea from the north. The source of these warmer water masses is the Atlantic Water boundary current that flows along the continental slope north of Svalbard. Time-varying southward inflow through cross-shelf troughs was the main driver of the seasonal cycle in ocean temperature at the moorings. Inflows were intensified in autumn and early winter, in some cases occurring below the sea ice cover and halocline water. On shorter timescales, subtidal current variability was correlated with the large-scale meridional atmospheric pressure gradient, suggesting wind-driven modulation of the inflow. The mooring records also show that import of sea ice into the Barents Sea has a lasting impact on the upper ocean, where salinity and stratification are strongly affected by the amount of sea ice that has melted in the area. A fresh layer separated the ocean surface from the warm mid-depth waters following large sea ice imports in 2019, whereas diluted Atlantic Water was found close to the surface during episodes in autumn 2018 following a long ice-free period. Thus, the advective imports of ocean water and sea ice from surrounding areas are both key drivers of ocean variability in the region.

An update on dissolved methane distribution in the subtropical North Atlantic Ocean - Anna Kolomijeca; Lukas Marx; Sarah Reynolds; Thierry Cariou; Edward Mawji; Cedric Boulart (Ocean Science, 2022)
Methane (CH4) is a potent greenhouse gas and plays a significant role in recent increasing global temperatures. The oceans are a natural source of methane contributing to atmospheric methane concentrations, yet our understanding of the oceanic methane cycle is poorly constrained. Accumulating evidence indicates that a significant part of oceanic CH4 is produced in oxygenated surface waters as a by-product of phytoplanktonic activity. This study focused on the subtropical North Atlantic Ocean (26∘ N, 80′ W and 26∘ N, 18′ W) where the distribution of dissolved CH4 concentrations and associated air–sea fluxes during winter 2020 were investigated. Water samples from 64 stations were collected from the upper water column up to depths of 400 m. The upper oxic mixed layer was oversaturated in dissolved CH4 with concentrations ranging 3–7 nmol L−1, with the highest concentrations of 7–10 nmol L−1 found to the east of the transect, consistent with other subtropical regions of the world's oceans. The high anomalies of dissolved CH4 were found to be associated with phosphate-depleted waters and regions where the abundance of the ubiquitous picocyanobacteria Synechococcus and Prochlorococcus were elevated. Although other phytoplanktonic phyla cannot be excluded, this suggests that cyanobacteria contribute to the release of CH4 in this region. The calculation of air–sea fluxes further confirmed the subtropical North Atlantic Ocean as a source of CH4. This study provides evidence to corroborate the key role that picocyanobacteria play in helping to explain the oversaturation of CH4 found in surface mixed layer of the open ocean, otherwise known as the “ocean methane paradox”.

Global coarse-grained mesoscale eddy statistics based on integrated kinetic energy and enstrophy correlations - Imre M. Jánosi; Holger Kantz; Jason A. C. Gallas; Miklós Vincze (Ocean Science, 2022)
Recently, Jánosi et al. (2019) introduced the concept of a “vortex proxy” based on an observation of strong correlations between integrated kinetic energy and integrated enstrophy over a large enough surface area. When mesoscale vortices are assumed to exhibit a Gaussian shape, the two spatial integrals have particularly simple functional forms, and a ratio of them defines an effective radius of a “proxy vortex”. In the original work, the idea was tested over a restricted area in the Californian Current System. Here we extend the analysis to global scale by means of 25 years of AVISO altimetry data covering the (ice-free) global ocean. The results are compared with a global vortex database containing over 64 million mesoscale eddies. We demonstrate that the proxy vortex representation of surface flow fields also works globally and provides a quick and reliable way to obtain coarse-grained vortex statistics. Estimated mean eddy sizes (effective radii) are extracted in very good agreement with the data from the vortex census. Recorded eddy amplitudes are directly used to infer the kinetic energy transported by the mesoscale vortices. The ratio of total and eddy kinetic energies is somewhat higher than found in previous studies. The characteristic westward drift velocities are evaluated by a time-lagged cross-correlation analysis of the kinetic energy fields. While zonal mean drift speeds are in good agreement with vortex trajectory evaluation in the latitude bands 30–5∘ S and 5–30∘ N, discrepancies are exhibited mostly at higher latitudes on both hemispheres. A plausible reason for somewhat different drift velocities obtained by eddy tracking and cross-correlation analysis is the fact that the drift of mesoscale eddies is only one component of the surface flow fields. Rossby wave activities, coherent currents, and other propagating features on the ocean surface apparently contribute to the zonal transport of kinetic energy.

Hydrography, circulation, and response to atmospheric forcing in the vicinity of the central Getz Ice Shelf, Amundsen Sea, Antarctica - Vår Dundas; Elin Darelius; Kjersti Daae; Nadine Steiger; Yoshihiro Nakayama; Tae-Wan Kim (Ocean Science, 2022)
Ice shelves in the Amundsen Sea are thinning rapidly as ocean currents bring warm water into the cavities beneath the floating ice. Although the reported melt rates for the Getz Ice Shelf are comparatively low for the region, its size makes it one of the largest freshwater sources around Antarctica, with potential consequences for, bottom water formation downstream, for example. Here, we use a 2-year-long novel mooring record (2016–2018) and 16-year-long regional model simulations to describe, for the first time, the hydrography and circulation in the vicinity of the ice front between Siple and Carney Island. We find that, throughout the mooring record, temperatures in the trough remain below 0.15 ∘C, more than 1 ∘C lower than in the neighboring Siple and Dotson Trough, and we observe a mean current (0.03 m s−1) directed toward the ice shelf front. The variability in the heat transport toward the ice shelf appears to be governed by nonlocal ocean surface stress over the Amundsen Sea Polynya region, and northward to the continental shelf break, where strengthened westward ocean surface stress leads to increased southward flow at the mooring site. The model simulations suggest that the heat content in the trough during the observed period was lower than normal, possibly owing to anomalously low summertime sea ice concentration and weak winds.

Multi-platform study of the extreme bloom of the barrel jellyfish Rhizostoma pulmo (Cnidaria: Scyphozoa) in the northernmost gulf of the Mediterranean Sea (Gulf of Trieste) in April 2021 - Nydia Catalina Reyes Suárez; Valentina Tirelli; Laura Ursella; Matjaž Ličer; Massimo Celio; Vanessa Cardin (Ocean Science, 2022)
On 7 April 2021, an exceptional bloom of the scyphomedusa Rhizostoma pulmo was observed in the Gulf of Trieste (Italy). Blooms of this species in the northern Adriatic Sea have been reported since the late 1800s: the density of jellyfish observed in 2021 reached more than 10 specimens per square metre. We analyse the bloom from a multi-platform approach using observations and model data at different timescales. We attempt to explain the intensity of the bloom as a consequence of thermohaline and hydrodynamical conditions in the gulf. Meteo-oceanographic conditions that may have contributed to the exceptional aggregation of jellyfish observed along the northernmost coast of the Adriatic Sea are discussed in detail. Specifically, our results indicate that this bloom was enabled by (1) the presence of a high number of jellyfish in the gulf, probably linked to the anomalously warm sea conditions in spring 2020 and winter 2021, which may have favoured a longer reproductive period and enhanced survival of adult R. pulmo, respectively; and (2) strong wind events, such as the bora wind for the Gulf of Trieste, which enhanced upwelling and mixing processes in the gulf, bringing the jellyfish from the deeper waters to the surface and clustering them along the coast.

The influence of tides on the marine carbonate chemistry of a coastal polynya in the south-eastern Weddell Sea - Elise S. Droste; Mario Hoppema; Melchor González-Dávila; Juana Magdalena Santana-Casiano; Bastien Y. Queste; Giorgio Dall'Olmo; Hugh J. Venables; Gerd Rohardt; Sharyn Ossebaar; Daniel Schuller; Sunke Trace-Kleeberg; Dorothee C. E. Bakker (Ocean Science, 2022)
Tides significantly affect polar coastlines by modulating ice shelf melt and modifying shelf water properties through transport and mixing. However, the effect of tides on the marine carbonate chemistry in such regions, especially around Antarctica, remains largely unexplored. We address this topic with two case studies in a coastal polynya in the south-eastern Weddell Sea, neighbouring the Ekström Ice Shelf. The case studies were conducted in January 2015 (PS89) and January 2019 (PS117), capturing semi-diurnal oscillations in the water column. These are pronounced in both physical and biogeochemical variables for PS89. During rising tide, advection of sea ice meltwater from the north-east created a fresher, warmer, and more deeply mixed water column with lower dissolved inorganic carbon (DIC) and total alkalinity (TA) content. During ebbing tide, water from underneath the ice shelf decreased the polynya's temperature, increased the DIC and TA content, and created a more stratified water column. The variability during the PS117 case study was much smaller, as it had less sea ice meltwater input during rising tide and was better mixed with sub-ice shelf water. The contrasts in the variability between the two case studies could be wind and sea ice driven, and they underline the complexity and highly dynamic nature of the system. The variability in the polynya induced by the tides results in an air–sea CO2 flux that can range between a strong sink (−24 mmol m−2 d−1) and a small source (3 mmol m−2 d−1) on a semi-diurnal timescale. If the variability induced by tides is not taken into account, there is a potential risk of overestimating the polynya's CO2 uptake by 67 % or underestimating it by 73 %, compared to the average flux determined over several days. Depending on the timing of limited sampling, the polynya may appear to be a source or a sink of CO2. Given the disproportionate influence of polynyas on heat and carbon exchange in polar oceans, we recommend future studies around the Antarctic and Arctic coastlines to consider the timing of tidal currents in their sampling strategies and analyses. This will help constrain variability in oceanographic measurements and avoid potential biases in our understanding of these highly complex systems.

Variability of surface gravity wave field over a realistic cyclonic eddy - Gwendal Marechal; Charly de Marez (Ocean Science, 2022)
Recent remote sensing measurements and numerical studies have shown that surface gravity waves interact strongly with small-scale open ocean currents. Through these interactions, the significant wave height, the wave frequency, and the wave direction are modified. In the present paper, we investigate the interactions of surface gravity waves with a large and isolated realistic cyclonic eddy. This eddy is subject to instabilities, leading to the generation of specific features at both the mesoscale and submesoscale ranges. We use the WAVEWATCH III numerical framework to force surface gravity waves in the eddy before and after its destabilization. In the wave simulations the source terms are deactivated, and waves are initialized with different wave intrinsic frequencies. The study of these simulations illustrates how waves respond to the numerous kinds of instabilities in the large cyclonic eddy from a few hundred to a few tens of kilometres. Our findings show that the spatial variability of the wave direction, the mean period, and the significant wave height is very sensitive to the presence of submesoscale structures resulting from the eddy destabilization. The intrinsic frequency of the incident waves is key in the change of the wave direction resulting from the current-induced refraction and in the location, from the boundary where waves are generated, of the maximum values of significant wave height. However, for a given current forcing, the maximum values of the significant wave height are similar regardless of the frequency of the incident waves. In this idealized study it has been shown that the spatial gradients of wave parameters are sharper for simulations forced with the destabilized eddy. Because the signature of currents on waves encodes important information of currents, our findings suggest that the vertical vorticity of the current could be statistically estimated from the significant wave height gradients down to a very fine spatial scale. Furthermore, this paper shows the necessity to include currents in parametric models of sea-state bias; using a coarse-resolution eddy field may severely underestimate the sea-state-induced noise in radar altimeter measurements.

On the uncertainty associated with detecting global and local mean sea level drifts on Sentinel-3A and Sentinel-3B altimetry missions - Rémi Jugier; Michaël Ablain; Robin Fraudeau; Adrien Guerou; Pierre Féménias (Ocean Science, 2022)
An instrumental drift in the point target response (PTR) parameters has been detected on the Copernicus Sentinel-3A altimetry mission. It will affect the accuracy of sea level sensing, which could result in errors in sea level change estimates of a few tenths of a millimeter per year. In order to accurately evaluate this drift, a method for detecting global and regional mean sea level relative drifts between two altimetry missions is implemented. Associated uncertainties are also accurately calculated thanks to a detailed error budget analysis. A drift on both Sentinel-3A (S3A) and Sentinel-3B (S3B) global mean sea level (GMSL) is detected with values significantly higher than expected. For S3A, the relative GMSL drift detected is 1.0 mm yr−1 with Jason-3 and 1.3 mm yr−1 with SARAL/AltiKa. For S3B, the relative GMSL drift detected is −3.4 mm yr−1 with Jason-3 and −2.2 mm yr−1 with SARAL/AltiKa. The drift detected at global level does not show detectable regional variations above the uncertainty level of the proposed method. The investigations led by the altimeter experts can now explain the origin of this drift for S3A and S3B. The ability of the implemented method to detect a sea level drift with respect to the length of the common period is also analyzed. We find that the minimum detectable sea level drift over a 5-year period is 0.3 mm yr−1 at the global scale and 1.5 mm yr−1 at the 2400 km regional scale. However, these levels of uncertainty do not meet the sea level stability requirements for climate change studies.

Net community production in the northwestern Mediterranean Sea from glider and buoy measurements - Michael P. Hemming; Jan Kaiser; Jacqueline Boutin; Liliane Merlivat; Karen J. Heywood; Dorothee C. E. Bakker; Gareth A. Lee; Marcos Cobas García; David Antoine; Kiminori Shitashima (Ocean Science, 2022)
The Mediterranean Sea comprises just 0.8 % of the global oceanic surface, yet considering its size, it is regarded as a disproportionately large sink for anthropogenic carbon due to its physical and biogeochemical characteristics. An underwater glider mission was carried out in March–April 2016 close to the BOUSSOLE and DyFAMed time series moorings in the northwestern Mediterranean Sea. The glider deployment served as a test of a prototype ion-sensitive field-effect transistor pH sensor. Dissolved oxygen (O2) concentrations and optical backscatter were also observed by the glider and increased between 19 March and 1 April, along with pH. These changes indicated the start of a phytoplankton spring bloom, following a period of intense mixing. Concurrent measurements of CO2 fugacity and O2 concentrations at the BOUSSOLE mooring buoy showed fluctuations, in qualitative agreement with the pattern of glider measurements. Mean net community production rates (N) were estimated from glider and buoy measurements of dissolved O2 and inorganic carbon (DIC) concentrations, based on their mass budgets. Glider and buoy DIC concentrations were derived from a salinity-based total alkalinity parameterisation, glider pH and buoy CO2 fugacity. The spatial coverage of glider data allowed the calculation of advective O2 and DIC fluxes. Mean N estimates for the euphotic zone between 10 March and 3 April were (-17±36) for glider O2, (44±94) for glider DIC, (17±37) for buoy O2 and (49±86) mmolm-2d-1 for buoy DIC, all indicating net metabolic balance over these 25 d. However, these 25 d were actually split into a period of net DIC increase and O2 decrease between 10 and 19 March and a period of net DIC decrease and O2 increase between 19 March and 3 April. The latter period is interpreted as the onset of the spring bloom. The regression coefficients between O2 and DIC-based N estimates were 0.25 ± 0.08 for the glider data and 0.54 ± 0.06 for the buoy, significantly lower than the canonical metabolic quotient of 1.45±0.15. This study shows the added value of co-locating a profiling glider with moored time series buoys, but also demonstrates the difficulty in estimating N, and the limitations in achievable precision.

Four-dimensional temperature, salinity and mixed-layer depth in the Gulf Stream, reconstructed from remote-sensing and in situ observations with neural networks - Etienne Pauthenet; Loïc Bachelot; Kevin Balem; Guillaume Maze; Anne-Marie Tréguier; Fabien Roquet; Ronan Fablet; Pierre Tandeo (Ocean Science, 2022)
Despite the ever-growing number of ocean data, the interior of the ocean remains undersampled in regions of high variability such as the Gulf Stream. In this context, neural networks have been shown to be effective for interpolating properties and understanding ocean processes. We introduce OSnet (Ocean Stratification network), a new ocean reconstruction system aimed at providing a physically consistent analysis of the upper ocean stratification. The proposed scheme is a bootstrapped multilayer perceptron trained to predict simultaneously temperature and salinity (T−S) profiles down to 1000 m and the mixed-layer depth (MLD) from surface data covering 1993 to 2019. OSnet is trained to fit sea surface temperature and sea level anomalies onto all historical in situ profiles in the Gulf Stream region. To achieve vertical coherence of the profiles, the MLD prediction is used to adjust a posteriori the vertical gradients of predicted T−S profiles, thus increasing the accuracy of the solution and removing vertical density inversions. The prediction is generalized on a 1/4∘ daily grid, producing four-dimensional fields of temperature and salinity, with their associated confidence interval issued from the bootstrap. OSnet profiles have root mean square error comparable with the observation-based Armor3D weekly product and the physics-based ocean reanalysis Glorys12. The lowest confidence in the prediction is located north of the Gulf Stream, between the shelf and the current, where the thermohaline variability is large. The OSnet reconstructed field is coherent even in the pre-Argo years, demonstrating the good generalization properties of the network. It reproduces the warming trend of surface temperature, the seasonal cycle of surface salinity and mesoscale structures of temperature, salinity and MLD. While OSnet delivers an accurate interpolation of the ocean stratification, it is also a tool to study how the ocean stratification relates to surface data. We can compute the relative importance of each input for each T−S prediction and analyse how the network learns which surface feature influences most which property and at which depth. Our results demonstrate the potential of machine learning methods to improve predictions of ocean interior properties from observations of the ocean surface.

Compound flooding in convergent estuaries: insights from an analytical model - Ramin Familkhalili; Stefan A. Talke; David A. Jay (Ocean Science, 2022)
We investigate here the effects of geometric properties (channel depth and cross-sectional convergence length), storm surge characteristics, friction, and river flow on the spatial and temporal variability of compound flooding along an idealized, meso-tidal coastal-plain estuary. An analytical model is developed that includes exponentially convergent geometry, tidal forcing, constant river flow, and a representation of storm surge as a combination of two sinusoidal waves. Nonlinear bed friction is treated using Chebyshev polynomials and trigonometric functions, and a multi-segment approach is used to increase accuracy. Model results show that river discharge increases the damping of surge amplitudes in an estuary, while increasing channel depth has the opposite effect. Sensitivity studies indicate that the impact of river flow on peak water level decreases as channel depth increases, while the influence of tide and surge increases in the landward portion of an estuary. Moreover, model results show less surge damping in deeper configurations and even amplification in some cases, while increased convergence length scale increases damping of surge waves with periods of 12–72 h. For every modeled scenario, there is a point where river discharge effects on water level outweigh tide/surge effects. As a channel is deepened, this cross-over point moves progressively upstream. Thus, channel deepening may alter flood risk spatially along an estuary and reduce the length of a river estuary, within which fluvial flooding is dominant.

Kinematics of surface currents at the northern margin of the Gulf of Cádiz - Luciano de Oliveira Júnior; Paulo Relvas; Erwan Garel (Ocean Science, 2022)
The subtidal surface water circulation at the northern margin of the Gulf of Cádiz, at the southern extremity of the Iberian upwelling system, is described based on validated hourly high-frequency radar measurements from 2016 to 2020. Statistical analyses (mean, standard deviation, eccentricity and empirical orthogonal functions) are applied to the dataset, which is completed with ADCP time series from multiple moorings at five inner-shelf stations and ERA5 wind. Off the shelf, the main circulation pattern consists of a slope current, best developed in summer when north-westerlies dominate, in particular at the most exposed western region. Mechanisms other than upwelling must contribute to this flow in order to explain its seasonal persistence. The slope circulation reverses for regional wind events with an east component > 10 m s−1, approximately. On the shelf, currents are mainly alongshore and balanced. The circulation is generally continuous along the coast, except for weak (

Impact of a long-lived anticyclonic mesoscale eddy on seawater anomalies in the northeastern tropical Pacific Ocean: a composite analysis from hydrographic measurements, sea level altimetry data, and reanalysis model products - Kaveh Purkiani; Matthias Haeckel; Sabine Haalboom; Katja Schmidt; Peter Urban; Iason-Zois Gazis; Henko de Stigter; André Paul; Maren Walter; Annemiek Vink (Ocean Science, 2022)
Using observational data, satellite altimeters, and reanalysis model products, we have investigated eddy-induced seawater anomalies and heat and salt transport in the northeastern tropical Pacific Ocean. An eddy detection algorithm (EDA) was used to identify eddy formation at the Mexican Tehuantepec Gulf (TT) in July 2018 during an unusually strong summer wind event. The eddy separated from the coast with a mean translation velocity of 11 cm s−1 and a mean radius of 115 km and traveled 2050–2400 km westwards off the Central American coast, where it was followed at approx 114∘ W and 11∘ N for oceanographic observation between April and May 2019. The in situ observations show that the major eddy impacts are restricted to the upper 300 m of the water column and are traceable down to 1500 m water depth. In the eddy core at 92 m water depth an extreme positive temperature anomaly of 8.2 ∘C, a negative salinity anomaly of −0.78 psu, a positive fluorescence anomaly of +0.8 mg m−3, and a positive dissolved oxygen concentration anomaly of 137 µmol kg−1 are observed. Compared with annual climatological averages in 2018, the water trapped within the eddy is estimated to transport an average positive westward zonal heat anomaly of 85×1012 W and an average westward negative salt anomaly of -2.1×106 kg s−1. The heat transport is the equivalent of 1 % of the total annual zonal eddy-induced heat transport at this latitude in the Pacific Ocean. Understanding the dynamics of long-lived mesoscale eddies that may reach the seafloor in this region of the Pacific Ocean is especially important in light of potential deep-sea mining activities that are being targeted on this area.

The role of wind, mesoscale dynamics, and coastal circulation in the interannual variability of the South Vietnam Upwelling, South China Sea – answers from a high-resolution ocean model - Thai To Duy; Marine Herrmann; Claude Estournel; Patrick Marsaleix; Thomas Duhaut; Long Bui Hong; Ngoc Trinh Bich (Ocean Science, 2022)
The South Vietnam Upwelling (SVU) develops in the South China Sea (SCS) under the influence of southwest monsoon winds. To study the role of small spatiotemporal scales on the SVU functioning and variability, a simulation was performed over 2009–2018 with a high-resolution configuration (1 km at the coast) of the SYMPHONIE model implemented over the western region of the SCS. Its capability to represent ocean dynamics and water masses from daily to interannual scales and from coastal to regional areas is quantitatively demonstrated by comparison with available satellite data and four in situ datasets. The SVU interannual variability is examined for the three development areas already known: the southern (SCU) and northern (NCU) coastal upwelling areas and the offshore upwelling area (OFU). Our high-resolution model, together with in situ observations and high-resolution satellite data, moreover shows for the first time that upwelling develops over the Sunda Shelf off the Mekong Delta (MKU). Our results confirm for the SCU and OFU and show for the MKU the role of the mean summer intensity of wind and cyclonic circulation over the offshore area in driving the interannual variability of the upwelling intensity. They further reveal that other factors contribute to SCU and OFU variability. First, the intraseasonal wind chronology strengthens (in the case of regular wind peaks occurring throughout the summer for SCU or of stronger winds in July–August for OFU) or weakens (in the case of intermittent wind peaks for SCU) the summer average upwelling intensity. Second, the mesoscale circulation influences this intensity (multiple dipole eddies and associated eastward jets developing along the coast enhance the SCU intensity). The NCU interannual variability is less driven by the regional-scale wind (with weaker monsoon favoring stronger NCU) and more by the mesoscale circulation in the NCU area: the NCU is prevented (favored) when alongshore (offshore) currents prevail.

Evaluation of basal melting parameterisations using in situ ocean and melting observations from the Amery Ice Shelf, East Antarctica - Madelaine Rosevear; Benjamin Galton-Fenzi; Craig Stevens (Ocean Science, 2022)
Ocean-driven melting of Antarctic ice shelves is causing accelerating loss of grounded ice from the Antarctic continent. However, the ocean processes governing ice shelf melting are not well understood, contributing to uncertainty in projections of Antarctica's contribution to sea level. Here, we analyse oceanographic data and in situ measurements of ice shelf melt collected from an instrumented mooring beneath the centre of the Amery Ice Shelf, East Antarctica. This is the first direct measurement of basal melting from the Amery Ice Shelf and was made through the novel application of an upward-facing acoustic Doppler current profiler (ADCP). ADCP data were also used to map a region of the ice base, revealing a steep topographic feature or “scarp” in the ice with vertical and horizontal scales of ∼ 20 and ∼ 40 m, respectively. The annually averaged ADCP-derived melt rate of 0.51 ± 0.18 m yr−1 is consistent with previous modelling results and glaciological estimates. There is significant seasonal variation around the mean melt rate, with a 40 % increase in melting in May and a 60 % decrease in September. Melting is driven by temperatures ∼ 0.2 ∘C above the local freezing point and background and tidal currents, which have typical speeds of 3.0 and 10.0 cm s−1, respectively. We use the coincident measurements of ice shelf melt and oceanographic forcing to evaluate parameterisations of ice–ocean interactions and find that parameterisations in which there is an explicit dependence of the melt rate on current speed beneath the ice tend to overestimate the local melt rate at AM06 by between 200 % and 400 %, depending on the choice of drag coefficient. A convective parameterisation in which melting is a function of the slope of the ice base is also evaluated and is shown to underpredict melting by 20 % at this site. By combining these new estimates with available observations from other ice shelves, we show that the commonly used current speed-dependent parameterisation overestimates melting at all but the coldest and most energetic cavity conditions.

Attributing decadal climate variability in coastal sea-level trends - Sam Royston; Rory J. Bingham; Jonathan L. Bamber (Ocean Science, 2022)
Decadal sea-level variability masks longer-term changes due to natural and anthropogenic drivers in short-duration records and increases uncertainty in trend and acceleration estimates. When making regional coastal management and adaptation decisions, it is important to understand the drivers of these changes to account for periods of reduced or enhanced sea-level change. The variance in decadal sea-level trends about the global mean is quantified and mapped around the global coastlines of the Atlantic, Pacific, and Indian oceans from historical CMIP6 runs and a high-resolution ocean model forced by reanalysis data. We reconstruct coastal, sea-level trends via linear relationships with climate mode and oceanographic indices. Using this approach, more than one-third of the variability in decadal sea-level trends can be explained by climate indices at 24.6 % to 73.1 % of grid cells located within 25 km of a coast in the Atlantic, Pacific, and Indian oceans. At 10.9 % of the world's coastline, climate variability explains over two-thirds of the decadal sea-level trend. By investigating the steric, manometric, and gravitational components of sea-level trend independently, it is apparent that much of the coastal ocean variability is dominated by the manometric signal, the consequence of the open-ocean steric signal propagating onto the continental shelf. Additionally, decadal variability in the gravitational, rotational, and solid-Earth deformation (GRD) signal should not be ignored in the total. There are locations such as the Persian Gulf and African west coast where decadal sea-level variability is historically small that are susceptible to future changes in hydrology and/or ice mass changes that drive intensified regional GRD sea-level change above the global mean. The magnitude of variance explainable by climate modes quantified in this study indicates an enhanced uncertainty in projections of short- to mid-term regional sea-level trend.

Mechanism of generation and propagation characteristics of coastal trapped waves in the Black Sea - Müjdat Aydın; Şükrü Turan Beşiktepe (Ocean Science, 2022)
Coastal trapped waves (CTW) are a major mechanism to distribute energy from the atmosphere in the ocean and play a significant role in large-scale, low-frequency sea-level and current variability in continental shelf and slope areas. Despite their significance for coastal dynamics, observational evidence of the influence of CTWs on the large-scale circulation is rather limited. In this study, mode-1 coastal trapped waves that were captured at sea-level stations at five locations along the southern coast of the Black Sea are examined together with sea surface height reanalysis from the Copernicus Marine Service to reveal their generation mechanisms and role in the coastal dynamics. It is found that the CTWs were formed when water accumulated on the western shelf after gale-force alongshore winds in the western Black Sea. Excited waves propagate along the Black Sea coast from west to east with a speed of 2.3–2.6 m s−1 and transport the atmospherically induced energy all over the Black Sea. The coastal current generated by the CTWs (on the order of 1 m s−1 in magnitude) and the main Black Sea current merge and flow eastward as a single structure, resulting in intensification of the Black Sea circulation during winter. Hence, we present evidence of the influence of CTWs on the large-scale circulation.

Technical note: On seasonal variability of the M 2 tide - Richard D. Ray (Ocean Science, 2022)
Seasonal variability of the M2 ocean tide can be detected at many ports, perhaps most. Examination of the cluster of tidal constituents residing within the M2 tidal group can shed light on the physical mechanisms underlying seasonality. In the broadest terms these are astronomical, frictional–advective interactions, and climate processes; some induce annual modulations and some semiannual, in amplitude, phase, or both. This note reviews how this occurs and gives an example from each broad category. Phase conventions and their relationship with causal mechanisms, as well as nomenclature, are also addressed.

Deep through-flow in the Bight Fracture Zone - Tillys Petit; Virginie Thierry; Herlé Mercier (Ocean Science, 2022)
Iceland–Scotland Overflow Water (ISOW) is exported from the Nordic Seas into the Iceland Basin to feed the lower limb of the Meridional Overturning Circulation. The Bight Fracture Zone (BFZ) is known to be a major route for ISOW toward the Irminger Sea, but the role of this gateway in the evolution of ISOW properties over the subpolar gyre is unclear. A combination of ship-based and Deep-Argo data gathered between 2015 and 2018 allows us to investigate the pathways and hydrographic evolution of ISOW as it flows through the BFZ, as well as its influence on the North Atlantic Deep Water (NADW) properties in the Irminger Sea. The ISOW flow through the BFZ amounts to 0.8 ± 0.2 Sv and is mainly fed by the lighter part of the ISOW layer flowing west of 29–30∘ W as part of the East Reykjanes Ridge Current in the Iceland Basin. In the rift valley of the BFZ, between an eastern and a western sill, the bathymetry of the BFZ shapes a cyclonic circulation along which the ISOW layer is homogenized. The largest changes in ISOW properties are however observed downstream of the western sill, at the exit of the BFZ. There, ISOW is mixed isopycnally with comparatively fresher NADW circulating in the Irminger Sea. Hence, our analysis reveals the key role of the BFZ through-flow in the salinification of the NADW in the Irminger Current.

Coastal sea level monitoring in the Mediterranean and Black seas - Begoña Pérez Gómez; Ivica Vilibić; Jadranka Šepić; Iva Međugorac; Matjaž Ličer; Laurent Testut; Claire Fraboul; Marta Marcos; Hassen Abdellaoui; Enrique Álvarez Fanjul; Darko Barbalić; Benjamín Casas; et al. (Ocean Science, 2022)
Employed for over a century, the traditional way of monitoring sea level variability by tide gauges – in combination with modern observational techniques like satellite altimetry – is an inevitable ingredient in sea level studies over the climate scales and in coastal seas. The development of the instrumentation, remote data acquisition, processing, and archiving in the last decades has allowed the extension of the applications to a variety of users and coastal hazard managers. The Mediterranean and Black seas are examples of such a transition – while having a long tradition of sea level observations with several records spanning over a century, the number of modern tide gauge stations is growing rapidly, with data available both in real time and as a research product at different time resolutions. As no comprehensive survey of the tide gauge networks has been carried out recently in these basins, the aim of this paper is to map the existing coastal sea level monitoring infrastructures and the respective data availability. The survey encompasses a description of major monitoring networks in the Mediterranean and Black seas and their characteristics, including the type of sea level sensors, measuring resolutions, data availability, and existence of ancillary measurements, altogether collecting information about 240 presently operational tide gauge stations. The availability of the Mediterranean and Black seas sea level data in the global and European sea level repositories has been also screened and classified following their sampling interval and level of quality check, pointing to the necessity of harmonization of the data available with different metadata and series in different repositories. Finally, an assessment of the networks' capabilities for their use in different sea level applications has been done, with recommendations that might mitigate the bottlenecks and ensure further development of the networks in a coordinated way, a critical need in the era of human-induced climate changes and sea level rise.

Nonlocal and local wind forcing dependence of the Atlantic meridional overturning circulation and its depth scale - Tim Rohrschneider; Johanna Baehr; Veit Lüschow; Dian Putrasahan; Jochem Marotzke (Ocean Science, 2022)
We use wind sensitivity experiments to understand the wind forcing dependencies of the level of no motion as the depth of maximum overturning and the e-folding pycnocline scale, as well as their relationship to northward transport of the mid-depth Atlantic meridional overturning circulation (AMOC) which extends vertically to the depth of maximum overturning of the upper AMOC cell. In contrast to previous studies, we investigate the interplay of nonlocal and local wind effects on a decadal timescale. We use 30-year simulations with a high-resolution ocean general circulation model (OGCM) which is an eddy-resolving version of the Max Planck Institute for Meteorology ocean model (MPIOM). Our findings deviate from the common perspective that the AMOC is a nonlocal phenomenon only, because northward transport and its depth scales depend on both nonlocal Southern Ocean wind effects and local wind effects in the Northern Hemisphere downwelling region where Ekman pumping takes place. Southern Ocean wind forcing predominantly determines the magnitude of the pycnocline scale throughout the basin, whereas Northern Hemisphere winds additionally influence the level of no motion locally. In that respect, the level of no motion is a better proxy for northward transport and mid-depth velocity profiles than the pycnocline scale, since the wind forcing dependencies of the level of no motion and maximum overturning are equal. The changes in maximum overturning with wind forcing are explained by the changes in the level of no motion only. This is because wind-driven Ekman compensation is baroclinic and occurs above the level of no motion, and the internal vertical velocity shear that is not influenced by the external Ekman cells stays approximately constant. The analysis of the wind experiments suggests a hemisphere-dependent scaling of the strength of AMOC. We put forward the idea that the ability of numerical models to capture the spatial and temporal variations of the level of no motion is crucial to reproduce the mid-depth cell in an appropriate way both quantitatively and dynamically.

Causes of the 2015 North Atlantic cold anomaly in a global state estimate - Rachael N. C. Sanders; Daniel C. Jones; Simon A. Josey; Bablu Sinha; Gael Forget (Ocean Science, 2022)
The subpolar North Atlantic is an important part of the global ocean and climate system, with SST variability in the region influencing the climate of Europe and North America. While the majority of the global ocean exhibited higher than average surface temperatures in 2015, the subpolar North Atlantic experienced record low temperatures. This interannual cold anomaly is thought to have been driven by surface forcing, but detailed questions remain about how the anomaly was created and maintained. To better quantify and understand the processes responsible for the cold anomaly, we computed mixed-layer temperature budgets in the Estimating the Circulation and Climate of the Ocean (ECCO) Version 4 global ocean state estimate. State estimates have been brought into consistency with a large suite of observations without using artificial sources or sinks of heat, making them ideal for temperature budget studies. We found that strong surface forcing drove approximately 75 % of the initial anomalies in the cooling of the mixed layer in December 2013, while horizontal advection drove the remaining 25 %. The cold anomaly was then sequestered beneath the mixed layer. Re-emergence of the cold anomaly during the summer and autumn of 2014 was primarily the result of a strong temperature gradient across the base of the mixed layer, with vertical diffusion accounting for approximately 70 % of the re-emergence. Weaker surface warming of the mixed layer during the summer of 2015 enhanced the anomaly, causing a temperature minimum. Spatial patterns in the budgets also show large differences between the north and south of the anomaly region, with particularly strong initial surface cooling in the south related to the positive phase of the East Atlantic Pattern. It is important to note that this interannual cold anomaly, which is thought to be primarily driven by surface forcing, is distinct from the multi-decadal North Atlantic “warming hole”, which has been associated with changes in advection.

The Levantine Intermediate Water in the western Mediterranean and its interactions with the Algerian Gyres: insights from 60 years of observation - Katia Mallil; Pierre Testor; Anthony Bosse; Félix Margirier; Loic Houpert; Hervé Le Goff; Laurent Mortier; Ferial Louanchi (Ocean Science, 2022)
The presence of two large-scale cyclonic gyres in the Algerian Basin influences the general and eddy circulation, but their effects on water mass transfer remain poorly characterized. Our study has confirmed the presence of these gyres using the first direct current measurements of the whole water column collected during the SOMBA-GE2014 cruise, specifically designed to investigate these gyres. Using cruise sections and a climatology from 60 years of in situ measurements, we have also shown the effect of these gyres on the distribution at intermediate depth of Levantine Intermediate Water (LIW) with warmer (∼ 0.15 ∘C) and saltier (∼ 0.02) characteristics in the Algerian Basin than in the Provençal Basin. The Algerian Gyres, combined with the effect of anticyclonic Algerian Eddies, also impact horizontal density gradients with sinking of the isopycnals at the gyres' centers. Temporal cross-correlation of LIW potential temperature referenced to a signal observed southwest of Sardinia reveals a timescale of transit of 4 months to get to the center of the Algerian Basin. The LIW potential temperature and salinity trends, on average in the basin interior, are estimated to be +0.0022 ± 0.0002 ∘Cyr-1 and +0.0022 ± 0.0001 yr−1, respectively, over the 1968–2017 period and accelerating to +0.048 ± 0.003 ∘Cyr-1 and +0.0076 ± 0.0009 yr−1 over the 2013–2017 period.

Clustering analysis of the Sargassum transport process: application to beaching prediction in the Lesser Antilles - Didier Bernard; Emmanuel Biabiany; Raphaël Cécé; Romual Chery; Naoufal Sekkat (Ocean Science, 2022)
The massive Sargassum algae beachings observed over the past decade are a new natural hazard currently impacting the island states of the Caribbean region (human health, environmental damages, and economic losses). This study aims to improve the prediction of the surface current dynamic leading to beachings in the Lesser Antilles using clustering analysis methods. The input surface currents were derived from the Mercator model and the Hybrid Coordinate Ocean Model (HYCOM) outputs in which we integrated the windage effect. Past daily observations of Sargassum beaching on Guadeloupe coasts and satellite-based Sargassum offshore abundance were also integrated. Four representative current regimes were identified for both Mercator and HYCOM data. The analysis of the current sequences leading to beachings showed that the recurrence of two current regimes is related to the beaching peaks respectively observed in March and August. The performance score of the predictive model showed that the HYCOM data seem more suitable to assess coastal Sargassum hazard in the Lesser Antilles. For 1 year of tests (i.e., 2021), the decision tree accuracy respectively reached 70.1 % and 58.2 % for HYCOM and Mercator with a temporal uncertainty range ±3 d around the forecast date. The present clustering analysis predictive system, requiring lower computational resources compared to conventional forecast models, would help improve this risk management in the islands of the region.

Technical note: Tail behaviour of the statistical distribution of extreme storm surges - Tom Howard (Ocean Science, 2022)
The tail behaviour of the statistical distribution of extreme storm surges is conveniently described by a return level plot, consisting of water level (y axis) against average recurrence interval on a logarithmic scale (x axis). An average recurrence interval is often referred to as a “return period”. Hunter's allowance for sea-level rise gives a suggested amount by which to raise coastal defences in order to maintain the current level of flood risk, given an uncertain projection of future mean sea-level rise. The allowance is most readily evaluated by assuming that sea-level annual maxima follow a Gumbel distribution, and the evaluation is awkward if we use a generalized extreme value (GEV) fit. When we use a Gumbel fit, we are effectively assuming that the return level plot is a straight line. In other words, the shape parameter, which describes the curvature of the return level plot, is zero. On the other hand, coastal asset managers may need an estimate of the return period of unprecedented events even under current mean sea levels. For this purpose, curvature of the return level plot is usually accommodated by allowing a non-zero shape parameter whilst extrapolating the return level plot beyond the observations, using some kind of fit to observed extreme values (for example, a GEV fit to annual maxima). This might seem like a conflict: which approach is “correct”? Here I present evidence that the shape parameter varies around the coast of the UK and is consequently not zero. Despite this, I argue that there is no conflict: a suitably constrained non-zero-shape fit is appropriate for extrapolation and a Gumbel fit is appropriate for evaluation of Hunter's allowance.

Data-assimilation-based parameter estimation of bathymetry and bottom friction coefficient to improve coastal accuracy in a global tide model - Xiaohui Wang; Martin Verlaan; Jelmer Veenstra; Hai Xiang Lin (Ocean Science, 2022)
Global tide and surge models play a major role in forecasting coastal flooding due to extreme events or climate change. The model performance is strongly affected by parameters such as bathymetry and bottom friction. In this study, we propose a method that estimates bathymetry globally and the bottom friction coefficient in shallow waters for a global tide and surge model (GTSMv4.1). However, the estimation effect is limited by the scarcity of available tide gauges. We propose complementing sparse tide gauges with tide time series generated using FES2014. The FES2014 dataset outperforms the GTSM in most areas and is used as observations for the deep ocean and some coastal areas, such as Hudson Bay and Labrador, where tide gauges are scarce but energy dissipation is large. The experiment is performed with a computation- and memory-efficient iterative parameter estimation scheme (time–POD-based coarse incremental parameter estimation; POD: proper orthogonal decomposition) applied to the Global Tide and Surge Model (GTSMv4.1). Estimation results show that model performance is significantly improved for the deep ocean and shallow waters, especially in the European shelf, directly using the CMEMS tide gauge data in the estimation. The GTSM is also validated by comparing to tide gauges from UHSLC, CMEMS, and some Arctic stations in the year 2014.

Quasi-steady circulation regimes in the Baltic Sea - Taavi Liblik; Germo Väli; Kai Salm; Jaan Laanemets; Madis-Jaak Lilover; Urmas Lips (Ocean Science, 2022)
Circulation plays an essential role in the creation of physical and biogeochemical fluxes in the Baltic Sea. The main aim of the work was to study the quasi-steady circulation patterns under prevailing forcing conditions. A total of 6 months of continuous vertical profiling and fixed-point measurements of currents, two month-long underwater glider surveys, and numerical modeling were applied in the central Baltic Sea. The vertical structure of currents was strongly linked to the location of the two pycnoclines: the seasonal thermocline and the halocline. The vertical movements of pycnoclines and velocity shear maxima were synchronous. The quasi-steady circulation patterns were in geostrophic balance and highly persistent. The persistent patterns included circulation features such as upwelling, downwelling, and boundary currents, as well as a sub-halocline gravity current. The patterns had a prevailing zonal scale of 5–60 km as well as considerably higher magnitude and different direction than the long-term mean circulation pattern. A northward (southward) geostrophic boundary current in the upper layer was observed along the eastern coast of the central Baltic in the case of southwesterly (northerly) wind. The geostrophic current at the boundary was often a consequence of wind-driven, across-shore advection. The sub-halocline quasi-permanent gravity current with a width of 10–30 km from the Gotland Deep to the north over the narrow sill separating the Fårö Deep and Nothern Deep was detected in the simulation, and it was confirmed by an Argo float trajectory. According to the simulation, a strong flow, mostly to the north, with a zonal scale of 5 km occurred at the sill. This current is an important deeper limb of the overturning circulation of the Baltic Sea. The current was stronger with northerly winds and restricted by the southwesterly winds. The circulation regime had an annual cycle due to seasonality in the forcing. The boundary current was stronger and more frequent northward during the winter period. The sub-halocline current towards the north was strongest in March–May and weakest in November–December.

Seasonal extrema of sea surface temperature in CMIP6 models - Yanxin Wang; Karen J. Heywood; David P. Stevens; Gillian M. Damerell (Ocean Science, 2022)
CMIP6 model sea surface temperature (SST) seasonal extrema averaged over 1981–2010 are assessed against the World Ocean Atlas (WOA18) observational climatology. We propose a mask to identify and exclude regions of large differences between three commonly used climatologies (WOA18, WOCE-Argo Global Hydrographic climatology (WAGHC) and the Hadley Centre Sea Ice and Sea Surface Temperature data set (HadISST)). The biases in SST seasonal extrema are largely consistent with the annual mean SST biases. However, the amplitude and spatial pattern of SST bias vary seasonally in the 20 CMIP6 models assessed. Large seasonal variations in the SST bias occur in eastern boundary upwelling regions, polar regions, the North Pacific and the eastern equatorial Atlantic. These results demonstrate the importance of evaluating model performance not simply against annual mean properties. Models with greater vertical resolution in their ocean component typically demonstrate better representation of SST extrema, particularly seasonal maximum SST. No significant relationship of SST seasonal extrema with horizontal ocean model resolution is found.

Coastal high-frequency radars in the Mediterranean – Part 2: Applications in support of science priorities and societal needs - Emma Reyes; Eva Aguiar; Michele Bendoni; Maristella Berta; Carlo Brandini; Alejandro Cáceres-Euse; Fulvio Capodici; Vanessa Cardin; Daniela Cianelli; Giuseppe Ciraolo; Lorenzo Corgnati; Vlado Dadić; et al. (Ocean Science, 2022)
The Mediterranean Sea is a prominent climate-change hot spot, with many socioeconomically vital coastal areas being the most vulnerable targets for maritime safety, diverse met-ocean hazards and marine pollution. Providing an unprecedented spatial and temporal resolution at wide coastal areas, high-frequency radars (HFRs) have been steadily gaining recognition as an effective land-based remote sensing technology for continuous monitoring of the surface circulation, increasingly waves and occasionally winds. HFR measurements have boosted the thorough scientific knowledge of coastal processes, also fostering a broad range of applications, which has promoted their integration in coastal ocean observing systems worldwide, with more than half of the European sites located in the Mediterranean coastal areas. In this work, we present a review of existing HFR data multidisciplinary science-based applications in the Mediterranean Sea, primarily focused on meeting end-user and science-driven requirements, addressing regional challenges in three main topics: (i) maritime safety, (ii) extreme hazards and (iii) environmental transport process. Additionally, the HFR observing and monitoring regional capabilities in the Mediterranean coastal areas required to underpin the underlying science and the further development of applications are also analyzed. The outcome of this assessment has allowed us to provide a set of recommendations for future improvement prospects to maximize the contribution to extending science-based HFR products into societally relevant downstream services to support blue growth in the Mediterranean coastal areas, helping to meet the UN's Decade of Ocean Science for Sustainable Development and the EU's Green Deal goals.

Coastal high-frequency radars in the Mediterranean – Part 1: Status of operations and a framework for future development - Pablo Lorente; Eva Aguiar; Michele Bendoni; Maristella Berta; Carlo Brandini; Alejandro Cáceres-Euse; Fulvio Capodici; Daniela Cianelli; Giuseppe Ciraolo; Lorenzo Corgnati; Vlado Dadić; Bartolomeo Doronzo; et al. (Ocean Science, 2022)
Due to the semi-enclosed nature of the Mediterranean Sea, natural disasters and anthropogenic activities impose stronger pressures on its coastal ecosystems than in any other sea of the world. With the aim of responding adequately to science priorities and societal challenges, littoral waters must be effectively monitored with high-frequency radar (HFR) systems. This land-based remote sensing technology can provide, in near-real time, fine-resolution maps of the surface circulation over broad coastal areas, along with reliable directional wave and wind information. The main goal of this work is to showcase the current status of the Mediterranean HFR network and the future roadmap for orchestrated actions. Ongoing collaborative efforts and recent progress of this regional alliance are not only described but also connected with other European initiatives and global frameworks, highlighting the advantages of this cost-effective instrument for the multi-parameter monitoring of the sea state. Coordinated endeavors between HFR operators from different multi-disciplinary institutions are mandatory to reach a mature stage at both national and regional levels, striving to do the following: (i) harmonize deployment and maintenance practices; (ii) standardize data, metadata, and quality control procedures; (iii) centralize data management, visualization, and access platforms; and (iv) develop practical applications of societal benefit that can be used for strategic planning and informed decision-making in the Mediterranean marine environment. Such fit-for-purpose applications can serve for search and rescue operations, safe vessel navigation, tracking of marine pollutants, the monitoring of extreme events, the investigation of transport processes, and the connectivity between offshore waters and coastal ecosystems. Finally, future prospects within the Mediterranean framework are discussed along with a wealth of socioeconomic, technical, and scientific challenges to be faced during the implementation of this integrated HFR regional network.

Tracer and observationally derived constraints on diapycnal diffusivities in an ocean state estimate - David S. Trossman; Caitlin B. Whalen; Thomas W. N. Haine; Amy F. Waterhouse; An T. Nguyen; Arash Bigdeli; Matthew Mazloff; Patrick Heimbach (Ocean Science, 2022)
Use of an ocean parameter and state estimation framework – such as the Estimating the Circulation and Climate of the Ocean (ECCO) framework – could provide an opportunity to learn about the spatial distribution of the diapycnal diffusivity parameter (κρ) that observations alone cannot due to gaps in coverage. However, we show that the inclusion of misfits to observed physical variables – such as in situ temperature, salinity, and pressure – currently accounted for in ECCO is not sufficient, as κρ from ECCO does not agree closely with any observationally derived product. These observationally derived κρ products were inferred from microstructure measurements, derived from Argo and conductivity–temperature–depth (CTD) data using a strain-based parameterization of fine-scale hydrographic structure, or calculated from climatological and seafloor data using a parameterization of tidal mixing. The κρ products are in close agreement with one another but have both measurement and structural uncertainties, whereas tracers can have relatively small measurement uncertainties. With the ultimate goal being to jointly improve the ECCO state estimate and representation of κρ in ECCO, we investigate whether adjustments in κρ due to inclusion of misfits to a tracer – dissolved oxygen concentrations from an annual climatology – would be similar to those due to inclusion of misfits to observationally derived κρ products. We do this by performing sensitivity analyses with ECCO. We compare multiple adjoint sensitivity calculations: one configuration uses misfits to observationally derived κρ, and the other uses misfits to observed dissolved oxygen concentrations. We show that adjoint sensitivities of dissolved oxygen concentration misfits to the state estimate's control space typically direct κρ to improve relative to the observationally derived values. These results suggest that the inclusion of oxygen in ECCO's misfits will improve κρ in ECCO, particularly in (sub)tropical regions.

Counter-rotating eddy pair in the Luzon Strait - Ruili Sun; Peiliang Li; Yanzhen Gu; Fangguo Zhai; Yunwei Yan; Bo Li; Yang Zhang (Ocean Science, 2022)
Based on satellite remote-sensing observation data and Hybrid Coordinate Ocean Model (HYCOM) re-analysis data, we studied the counter-rotating eddy pair in the Luzon Strait (LS). Statistical analysis reveals that when an anti-cyclonic mesoscale eddy (AE) (cyclonic mesoscale eddy [CE]) in the Northwest Pacific (NWP) gradually approaches the east side of the LS, a CE (an AE) gradually forms on the west side of the LS, and it is defined as the AE (CE) mode of the counter-rotating eddy pair in the LS. The counter-rotating eddy pair exhibits obvious seasonal variation: the AE mode mainly occurs in the summer half of the year, while the CE mode mainly occurs in the winter half of the year. The mean durations of the AE and CE modes are both approximately 70 d. Based on the vorticity budget equation and energy analysis, the dynamic mechanism of counter-rotating eddy-pair occurrence is determined to be as follows: the AE (CE) on the east side of the LS causes a positive (negative) vorticity anomaly through horizontal velocity shear on the west side of the LS, and the positive (negative) vorticity anomaly is transported westward by the zonal advection of the vorticity, finally leading to the formation of the CE (AE) on the west side of the LS.

Influence of cyclonic and anticyclonic eddies on plankton in the southeastern Mediterranean Sea during late summertime - Natalia Belkin; Tamar Guy-Haim; Maxim Rubin-Blum; Ayah Lazar; Guy Sisma-Ventura; Rainer Kiko; Arseniy R. Morov; Tal Ozer; Isaac Gertman; Barak Herut; Eyal Rahav (Ocean Science, 2022)
Planktonic food webs were studied contemporaneously in a mesoscale cyclonic (upwelling, ∼ 13 months old) and an anticyclonic (downwelling, ∼ 2 months old) eddy as well as in an uninfluenced background situation in the oligotrophic southeastern Mediterranean Sea (SEMS) during late summer 2018. We show that integrated nutrient concentrations were higher in the cyclone compared to the anticyclone or the background stations by 2–13-fold. Concurrently, Synechococcus and Prochlorococcus were the dominant autotrophs abundance-wise in the oligotrophic anticyclone (∼ 300 × 1010 cells m−2). In the cyclone, functional groups such as dinoflagellates, Prymnesiophyceae and Ochrophyta contributed substantially to the total phytoplankton abundance (∼ 14 × 1010 cells m−2), which was ∼ 65 % lower at the anticyclone and background stations (∼ 5 × 1010 cells m−2). Primary production was highest in the cyclonic eddy (191 mgCm-2d-1) and 2–5-fold lower outside the eddy area. Heterotrophic prokaryotic cell-specific activity was highest in the cyclone (∼ 10 fgCcell-1d-1), while the least productive cells were found in the anticyclone (4 fgCcell-1d-1). Total zooplankton biomass in the upper 300 m was 10-fold higher in the cyclone compared with the anticyclone or background stations (1337 vs. 112–133 mg C m−2, respectively). Copepod diversity was much higher in the cyclone (44 species), compared to the anticyclone (6 small-size species). Our results highlight that cyclonic and anticyclonic eddies show significantly different community structure and food-web dynamics in oligotrophic environments, with cyclones representing productive oases in the marine desert of the SEMS.

Currents generated by the sea breeze in the southern Caspian Sea - Mina Masoud; Rich Pawlowicz (Ocean Science, 2022)
The sea breeze system is the dominant atmospheric forcing at high frequency in the southern Caspian Sea. Here, we describe and interpret current meter observations on the continental margins of the southern Caspian from 2012 to 2013 to identify and characterize the water column's response to the sea breeze system. Time series analysis provides evidence for diurnal baroclinic current signals of O(0.02 m s−1) and surface height changes of O(0.03 m). A two-layer model, including interfacial and bottom friction, is developed to further investigate the sea breeze response. This model is able to reproduce the structure, amplitudes, and phases of observed diurnal current fluctuations, explaining half of the variance in observational current response at frequencies at 1 cpd and higher. The sea breeze response thus results in a “tide-like” daily cycle, which is actually linked to the local forcing all along the southern Caspian coast.

Interannual variation of a bottom cold water mass in the Seto Inland Sea, Japan - Junying Zhu; Jie Shi; Xinyu Guo (Ocean Science, 2022)
A bottom cold water mass (BCWM) is a widespread physical oceanographic phenomenon in coastal seas, and its temperature variability has an important effect on the marine ecological environment. In this study, the interannual variation of the BCWM in Iyo-Nada (INCWM), a semi-enclosed bay in the Seto Inland Sea, Japan, from 1994 to 2015 and its response to air–sea heat flux change were investigated using monthly observational data and a hydrodynamic model. Surrounded by the isotherm of 18 ∘C, the observed multi-year average water temperature inside the INCWM was 17.58 ∘C with a standard deviation of 0.27 ∘C, while the mean area of INCWM was 5.73 × 105 m2 with a standard deviation of 4.35 × 105 m2. The interannual variation of average water temperature of INCWM showed a negative correlation with its area that indicates a low temperature corresponds to a big area. In addition, the interannual variation of the average temperature inside INCWM showed positive correlations with the local water temperature from April to July and with remote water temperature below 10 m in an adjacent strait in July. Differing from previously studied BCWMs, which had interannual variations depending closely on the water temperature before the warming season, the interannual variation of INCWM is more sensitive to the air–sea heat flux during the warming season than that in the previous winter. Further, by comparing several BCWMs, we found that the BCWM size is a key factor in understanding the heat transfer process responsible for the interannual variation of BCWMs in coastal seas. These findings will help us to understand the response of bottom cold water mass in coastal seas to sea surface forcing change.

Marine heatwaves in the Arabian Sea - Abhisek Chatterjee; Gouri Anil; Lakshmi R. Shenoy (Ocean Science, 2022)
Marine heatwaves (MHWs) are prolonged warm sea condition events that can have a destructive impact on marine ecosystems. The documentation of MHWs and assessment of their impacts is largely confined to a few regional seas or to global mean studies. The north Indian Ocean received almost no attention in this regard despite the fact that this ocean basin, particularly the Arabian Sea, has been warming at the most rapid pace among the other tropical basins in recent decades. This study shows the characteristics of MHW events for the Arabian Sea during 1982–2019. Our analysis shows that the duration (frequency) of MHWs exhibits a rapidly increasing trend of ∼20 d per decade (1.5–2 events per decade) in the northern Arabian Sea and the southeastern Arabian Sea close to the west coast of India, which is a multifold increase in MHW days (frequency) from the 80s. Notably, since the beginning of the satellite record, the years 2010 and 2016 have exhibited the maximum number of heatwave days when more than 75 % of days of the pre-monsoon and summer monsoon season experience heatwaves. The accelerated trend of the heatwave days is found to be driven by the rapid rise in the mean sea surface temperature (SST) of the Arabian Sea in the recent decade. Moreover, longer heatwave days are also associated with dominant climate modes. Among them, the Indian Ocean Basin Mode via the decaying phase of El Niño is the most influential mode contributing to more than 70 %–80 % of observed heatwave days in this basin. Further analysis of the most prolonged observed heatwave during April–June 2010 indicates that surface heat flux associated with the weaker latent heat loss and the shallow mixed layer was the primary cause of this event. Further, we note that the pre-monsoon cyclonic storms in the Arabian Sea often contribute to the waning of such heatwaves in the basin.

Technical note: TEOS-10 Excel – implementation of the Thermodynamic Equation Of Seawater – 2010 in Excel - Carlos Gil Martins; Jaimie Cross (Ocean Science, 2022)
This paper and associated software implement the Thermodynamic Equation Of Seawater – 2010 (TEOS-10) in Excel for an efficient estimation of Absolute Salinity (SA), Conservative Temperature (Θ), and derived thermodynamic properties of seawater – potential density (σΘ), in situ density (ρSA,Θ,p), and sound speed (c). Vertical profile template plots for these parameters are included as is an SA–Θ diagram template, which includes plotting of the density field (computation of user-selected σΘ lines is included). Absolute Salinity can be directly measured with the aid of a densimeter (IOC, SCOR and IAPSO, 2010: p. 82), but in TEOS-10 its estimation relies on the interpolation of data from casts of seawater from the world ocean (IOC, SCOR and IAPSO, 2010), and the Excel workbook introduced here (TEOS-10 Excel, available at https://doi.org/10.5281/zenodo.4748829) includes a subset of the TEOS-10 look-up tables necessary for this estimation, namely the Absolute Salinity Anomaly [deltaSA_ref] and the Absolute Salinity Anomaly Ratio [SAAR_ref] look-up tables. As the user simply needs to paste new data into the spreadsheet to automatically compute the oceanographic parameters referred above, this tool may prove to be useful for all who are not comfortable using the full-featured TEOS-10 programming language environments (e.g. MATLAB, FORTRAN, C) but rather need a simpler way of computing fundamental properties of seawater (e.g. density, sound speed) while adhering to current standards. Returned values are the same (up to 15 decimal places; i.e. difference = 0.000000000000000) as the ones obtained with the MATLAB version of the GSW (Gibbs Sea Water) toolbox (McDougall and Barker, 2011) available at the TEOS-10 website (https://www.teos-10.org, last access: 13 April 2022). This paper describes the Excel workbook, its use, and the included VBA (Visual Basic for Applications) functions. Quality control against the GSW toolbox is also addressed, namely issues detected with the interpolated values returned by the toolbox when there are missing values in the reference look-up table. In these situations, the GSW toolbox replaces missing values with a level pressure horizontal interpolation of neighbour points, while it is clear from the testing results that vertical interpolation, which was then implemented in TEOS-10 Excel, returns a more robust solution.

Contribution of a constellation of two wide-swath altimetry missions to global ocean analysis and forecasting - Mounir Benkiran; Pierre-Yves Le Traon; Gérald Dibarboure (Ocean Science, 2022)
Swath altimetry is likely to revolutionize our ability to monitor and forecast ocean dynamics. To meet the requirements of the EU Copernicus Marine Service, a constellation of two wide-swath altimeters is envisioned for the long-term (post-2030) evolution of the Copernicus Sentinel 3 topography mission. A series of observing system simulation experiments (OSSEs) is carried out to quantify the expected performances. The OSSEs use a state-of-the-art high-resolution (1/12∘) global ocean data assimilation system similar to the one used operationally by the Copernicus Marine Service. Flying a constellation of two wide-swath altimeters will provide a major improvement of our capabilities to monitor and forecast the oceans. Compared to the present situation with three nadir altimeters flying simultaneously, the sea surface height (SSH) analysis and 7 d forecast error are globally reduced by about 50 % in the OSSEs. With two wide-swath altimeters, the quality of SSH 7 d forecasts is equivalent to the quality of SSH analysis errors from three nadir altimeters. Our understanding of ocean currents is also greatly improved (30 % improvements at the surface and 50 % at 300 m depth). The resolution capabilities will be drastically improved and will be closer to 100 km wavelength compared to about 250 km today. Flying a constellation of two wide-swath altimeters thus looks to be a very promising solution for the long-term evolution of the Sentinel 3 constellation and the Copernicus Marine Service.

Ocean bubbles under high wind conditions – Part 2: Bubble size distributions and implications for models of bubble dynamics - Helen Czerski; Ian M. Brooks; Steve Gunn; Robin Pascal; Adrian Matei; Byron Blomquist (Ocean Science, 2022)
Bubbles formed by breaking waves in the open ocean influence many surface processes but are poorly understood. We report here on detailed bubble size distributions measured during the High Wind Speed Gas Exchange Study (HiWinGS) in the North Atlantic, during four separate storms with hourly averaged wind speeds from 10–27 m s−1. The measurements focus on the deeper plumes formed by advection downwards (at 2 m depth and below), rather than the initial surface distributions. Our results suggest that bubbles reaching a depth of 2 m have already evolved to form a heterogeneous but statistically stable population in the top 1–2 m of the ocean. These shallow bubble populations are carried downwards by coherent near-surface circulations; bubble evolution at greater depths is consistent with control by local gas saturation, surfactant coatings and pressure. We find that at 2 m the maximum bubble radius observed has a very weak wind speed dependence and is too small to be explained by simple buoyancy arguments. For void fractions greater than 10−6, bubble size distributions at 2 m can be fitted by a two-slope power law (with slopes of −0.3 for bubbles of radius

Ocean bubbles under high wind conditions – Part 1: Bubble distribution and development - Helen Czerski; Ian M. Brooks; Steve Gunn; Robin Pascal; Adrian Matei; Byron Blomquist (Ocean Science, 2022)
The bubbles generated by breaking waves are of considerable scientific interest due to their influence on air–sea gas transfer, aerosol production, and upper ocean optics and acoustics. However, a detailed understanding of the processes creating deeper bubble plumes (extending 2–10 m below the ocean surface) and their significance for air–sea gas exchange is still lacking. Here, we present bubble measurements from the HiWinGS expedition in the North Atlantic in 2013, collected during several storms with wind speeds of 10–27 m s−1. A suite of instruments was used to measure bubbles from a self-orienting free-floating spar buoy: a specialised bubble camera, acoustical resonators, and an upward-pointing sonar. The focus in this paper is on bubble void fractions and plume structure. The results are consistent with the presence of a heterogeneous shallow bubble layer occupying the top 1–2 m of the ocean, which is regularly replenished by breaking waves, and deeper plumes which are only formed from the shallow layer at the convergence zones of Langmuir circulation. These advection events are not directly connected to surface breaking. The void fraction distributions at 2 m depth show a sharp cut-off at a void fraction of 10−4.5 even in the highest winds, implying the existence of mechanisms limiting the void fractions close to the surface. Below wind speeds of 16 m s−1 or a wind-wave Reynolds number of RHw=2×106, the probability distribution of void fraction at 2 m depth is very similar in all conditions but increases significantly above either threshold. Void fractions are significantly different during periods of rising and falling winds, but there is no distinction with wave age. There is a complex near-surface flow structure due to Langmuir circulation, Stokes drift, and wind-induced current shear which influences the spatial distribution of bubbles within the top few metres. We do not see evidence for slow bubble dissolution as bubbles are carried downwards, implying that collapse is the more likely termination process. We conclude that the shallow and deeper bubble layers need to be studied simultaneously to link them to the 3D flow patterns in the top few metres of the ocean. Many open questions remain about the extent to which deep bubble plumes contribute to air–sea gas transfer. A companion paper (Czerski et al., 2022) addresses the observed bubble size distributions and the processes responsible for them.

Weakening and warming of the European Slope Current since the late 1990s attributed to basin-scale density changes - Matthew Clark; Robert Marsh; James Harle (Ocean Science, 2022)
Oceanic influences on shelf seas are mediated by flow along and across continental slopes, with consequences for regional hydrography and ecosystems. Here we present evidence for the variable North Atlantic influence on European shelf seas over the last 4 decades using ocean analysis and reanalysis products, as well as an eddy-resolving ocean model hindcast. To first order, flows oriented along isobaths at the continental slope are related to the poleward increase in density in the adjacent deep ocean that supports a geostrophic inflow towards the slope. In the North Atlantic, this density gradient and associated inflow have undergone substantial, sometimes abrupt, changes in recent decades. Inflow in the range 10–15 Sv is identified with eastward transport in temperature classes at 30∘ W in the latitude range 45–60∘ N. Associated with major subpolar warming around 1997, a cool and fresh branch of the Atlantic inflow was substantially reduced, while a warm and more saline inflow branch strengthened, with respective changes of the order of 5 Sv. Total inflow fell from ∼ 15 Sv pre-1997 to ∼ 10 Sv post-1997. In the model hindcast, particle tracking is used to trace the origins of poleward flows along the continental slope to the west of Ireland and Scotland before and after 1997. Backtracking particles up to 4 years, a range of subtropical and subpolar pathways is identified from a statistical perspective. In broad terms, cold, fresh waters of subpolar provenance were replaced by warm, saline waters of subtropical provenance. These changes have major implications for the downstream shelf regions that are strongly influenced by Atlantic inflow, in particular the northern North Sea, where “subtropicalization” of ecosystems has already been observed since the late 1990s.

Decomposing oceanic temperature and salinity change using ocean carbon change - Charles E. Turner; Peter J. Brown; Kevin I. C. Oliver; Elaine L. McDonagh (Ocean Science, 2022)
As the planet warms due to the accumulation of anthropogenic CO2 in the atmosphere, the interaction of surface ocean carbonate chemistry and the radiative forcing of atmospheric CO2 leads to the global ocean sequestering heat and carbon in a ratio that is nearly constant in time. This ratio has been approximated as globally uniform, enabling the intimately linked patterns of ocean heat and carbon uptake to be derived. Patterns of ocean salinity also change as the Earth system warms due to hydrological cycle intensification and perturbations to air–sea freshwater fluxes. Local temperature and salinity change in the ocean may result from perturbed air–sea fluxes of heat and fresh water (excess temperature, salinity) or from reorganisation of the preindustrial temperature and salinity fields (redistributed temperature, salinity), which are largely due to circulation changes. Here, we present a novel method in which the redistribution of preindustrial carbon is diagnosed and the redistribution of temperature and salinity is estimated using only local spatial information. We demonstrate this technique in the NEMO ocean general circulation model (OGCM) coupled to the MEDUSA-2 biogeochemistry model under an RCP8.5 scenario over 1860–2099. The excess changes (difference between total and redistributed property changes) are thus calculated. We demonstrate that a global ratio between excess heat and temperature is largely appropriate regionally with key regional differences consistent with reduced efficiency in the transport of carbon through the mixed layer base at high latitudes. On centennial timescales, excess heat increases everywhere, with the North Atlantic being a key site of excess heat uptake over the 21st century, accounting for 25 % of the total. Excess salinity meanwhile increases in the Atlantic but is generally negative in other basins, consistent with increasing atmospheric transport of fresh water out of the Atlantic. In the North Atlantic, changes in the inventory of excess salinity are detectable in the late 19th century, whereas increases in the inventory of excess heat do not become significant until the early 21st century. This is consistent with previous studies which find salinification of the subtropical North Atlantic to be an early fingerprint of anthropogenic climate change. Over the full simulation, we also find the imprint of Atlantic meridional overturning circulation (AMOC) slowdown through significant redistribution of heat away from the North Atlantic and of salinity to the South Atlantic. Globally, temperature change at 2000 m is accounted for by both redistributed and excess heat, but for salinity the excess component accounts for the majority of changes at the surface and at depth. This indicates that the circulation variability contributes significantly less to changes in ocean salinity than to heat content. By the end of the simulation excess heat is the largest contribution to density change and steric sea level rise, while excess salinity greatly reduces spatial variability in steric sea level rise through density compensation of excess temperature patterns, particularly in the Atlantic. In the Atlantic, redistribution of the preindustrial heat and salinity fields also produces generally compensating changes in sea level, though this compensation is less clear elsewhere. The regional strength of excess heat and salinity signals grows through the model run in response to the evolving forcing. In addition, the regional strength of the redistributed temperature and salinity signals also grows, indicating increasing circulation variability or systematic circulation change on timescales of at least the model run.

A newly reconciled dataset for identifying sea level rise and variability in Dublin Bay - Amin Shoari Nejad; Andrew C. Parnell; Alice Greene; Peter Thorne; Brian P. Kelleher; Robert J. N. Devoy; Gerard McCarthy (Ocean Science, 2022)
We provide an updated sea level dataset for Dublin for the period 1938–2016 at yearly resolution. Using a newly collated sea level record for Dublin Port, as well as two nearby tide gauges at Arklow and Howth Harbour, we perform data quality checks and calibration of the Dublin Port record by adjusting the biased high water level measurements that affect the overall calculation of mean sea level (MSL). To correct these MSL values, we use a novel Bayesian linear regression that includes the mean low water values as a predictor in the model. We validate the re-created MSL dataset and show its consistency with other nearby tide gauge datasets. Using our new corrected dataset, we estimate a rate of sea level rise of 1.1 mm yr−1 during 1953–2016 (95 % credible interval from 0.6 to 1.6 mm yr−1), and a rate of 7 mm yr−1 during 1997–2016 (95 % credible interval from 5 to 8.8 mm yr−1). The overall sea level rise is in line with expected trends, but large multidecadal variability has led to higher rates of rise in recent years.

Untangling the mistral and seasonal atmospheric forcing driving deep convection in the Gulf of Lion: 2012–2013 - Douglas Keller Jr.; Yonatan Givon; Romain Pennel; Shira Raveh-Rubin; Philippe Drobinski (Ocean Science, 2022)
Deep convection in the Gulf of Lion is believed to be primarily driven by the mistral winds. However, our findings show that the seasonal atmospheric change provides roughly two-thirds of the buoyancy loss required for deep convection to occur for the year 2012 to 2013, with the mistral supplying the final third. Two NEMOMED12 ocean simulations of the Mediterranean Sea were run from 1 August 2012 to 31 July 2013, forced with two sets of atmospheric-forcing data from a RegIPSL coupled run within the Med-CORDEX framework. One set of atmospheric-forcing data was left unmodified, while the other was filtered to remove the signal of the mistral. The control simulation featured deep convection, while the seasonal simulation did not. A simple model was derived by relating the anomaly scale forcing (the difference between the control and seasonal runs) and the seasonal scale forcing to the ocean response through the stratification index. This simple model revealed that the mistral's effect on buoyancy loss depends more on its strength rather than its frequency or duration. The simple model also revealed that the seasonal cycle of the stratification index is equal to the net surface heat flux over the course of the year, with the stratification maximum and minimum occurring roughly at the fall and spring equinoxes.

Simultaneous estimation of ocean mesoscale and coherent internal tide sea surface height signatures from the global altimetry record - Clément Ubelmann; Loren Carrere; Chloé Durand; Gérald Dibarboure; Yannice Faugère; Maxime Ballarotta; Frédéric Briol; Florent Lyard (Ocean Science, 2022)
This study proposes an approach to estimate the ocean sea surface height signature of coherent internal tides from a 25-year along-track altimetry record, with a single inversion over time, resolving both internal tide contributions and mesoscale eddy variability. The inversion is performed on a reduced-order basis of topography and practically achieved with a conjugate gradient. The particularity of this approach is to mitigate the potential aliasing effects between mesoscales and internal tide estimation from the uneven altimetry sampling (observing the sum of these components) by accounting for their statistics simultaneously, while other methods generally use a prior mesoscale. The four major tidal components are considered (M2, K1, S2, O1) over the period 1992–2017 on a global configuration. From the solution, we use altimetry data after 2017 for independent validation in order to evaluate the performance of the simultaneous inversion and compare it with an existing model.

Inherent optical properties of dissolved and particulate matter in an Arctic fjord (Storfjorden, Svalbard) in early summer - Tristan Petit; Børge Hamre; Håkon Sandven; Rüdiger Röttgers; Piotr Kowalczuk; Monika Zablocka; Mats A. Granskog (Ocean Science, 2022)
There have been considerable efforts to understand the hydrography of Storfjorden (Svalbard). A recurring winter polynya with large sea ice production makes it an important region of dense water formation at the scale of the Arctic Ocean. In addition, this fjord is seasonally influenced by freshwater inputs from sea ice melt and the surrounding islands of the Svalbard archipelago, which impacts the hydrography. However, the understanding of factors controlling the optical properties of the waters in Storfjorden are lacking and are crucial for the development of more accurate regional bio-optical models. Here, we present results from the first detailed optical field survey of Storfjorden conducted in early summer of 2020. Our observations are based on spectrometric analysis of water samples and in situ vertical profiles with an absorption and attenuation meter, a fluorometer, and a conductivity, temperature, and depth (CTD) sensor. In addition to the expected seasonal contribution from phytoplankton, we find that in early summer waters in Storfjorden are optically complex with a significant contribution from coloured dissolved organic matter (CDOM, 33 %–64 % of the non-water absorption at 443 nm) despite relatively low CDOM concentrations and in the nearshore or near the seabed from non-algal particles (up to 61 % of the non-water absorption at 550 nm). In surface waters, the spatial variability of light attenuation was mainly controlled by inorganic suspended matter originating from river runoff. A distinct subsurface maximum of light attenuation was largely driven by a subsurface phytoplankton bloom, controlled by stratification resulting from sea ice melt. Lastly, the cold dense bottom waters of Storfjorden from winter sea ice production periodically overflows the sill at the mouth of the fjord and can thus reach the Fram Strait. It contained elevated levels of both non-algal particles and dissolved organic matter, which are likely caused by the dense flows of the nepheloid layer interacting with the sea bed.

There and back again, a journey of many pathways: conceptualising the marine organic carbon cycle - Maike Iris Esther Scheffold; Inga Hense (Ocean Science, 2022)
Understanding and determining the pathways that organic carbon (OC) takes in the ocean is one of the pressing tasks of our time, as the fate of OC in the ocean is linked to the climate system and the functionality of marine ecosystems. The multitude and complexity of these pathways are typically investigated with sophisticated, mainly quantitative, methods focused on individual pathways in order to resolve their interactions and processes as realistically as possible. In addition to these approaches towards understanding and recreating complexity, there is a need to identify commonalities and differences between individual OC pathways and define their overarching structures. Such structures can provide a framework for the growing number of partly overlapping concepts, which conceptualise selected OC pathways, and promote more systematic comparisons and consistent communication, especially between different disciplines. In response, we propose a (visual) concept in which we define such higher-level “structures” by comparing and condensing marine OC pathways based on their sequences of processes and the layers of the marine system in which they operate. The resulting structures comprise “closed loops”, three remineralisation and two recalcitrant dissolved organic carbon (rDOC) loops that close within marine systems, and “open loops”, condensing pathways leaving the marine system for the atmosphere or deeper sediment layers. In addition, we provide a synthesis of embedded processes, OC pools, and process-performing organisms (agents) embedded in these loops. By translating a definition of the biological carbon pump (BCP) into our concept, we show how the application and discussion of our defined structures facilitate a consistent visualisation, a systematic comparison of differently resolved concepts and studies, and integration of these into the larger picture of the marine OC cycle. As a complement to quantitative studies and descriptions of individual pathways, our concept decomposes the complexity of OC pathways by defining new universal structures. These structures provide a skeleton that can be adapted to different systems and filled with life by the users.

Forecasting hurricane-forced significant wave heights using a long short-term memory network in the Caribbean Sea - Brandon J. Bethel; Wenjin Sun; Changming Dong; Dongxia Wang (Ocean Science, 2022)
A long short-term memory (LSTM) neural network is proposed to predict hurricane-forced significant wave heights (SWHs) in the Caribbean Sea (CS) based on a dataset of 20 CS, Gulf of Mexico, and western Atlantic hurricane events collected from 10 buoys from 2010–2020. SWH nowcasting and forecasting are initiated using LSTM on 0, 3, 6, 9, and 12 h horizons. Through examining study cases Hurricanes Dorian (2019), Sandy (2012), and Igor (2010), results illustrate that the model is well suited to forecast hurricane-forced wave heights much more rapidly at a significantly cheaper computational cost compared to numerical wave models, with much less required expertise. Forecasts are highly accurate with regards to observations. For example, Hurricane Dorian nowcasts had correlation (R), root mean square error (RMSE), and mean absolute percentage error (MAPE) values of 0.99, 0.16 m, and 2.6 %, respectively. Similarly, on the 3, 6, 9, and 12 h forecasts, results produced R (RMSE; MAPE) values of 0.95 (0.51 m; 7.99 %), 0.92 (0.74 m; 10.83 %), 0.85 (1 m; 13.13 %), and 0.84 (1.24 m; 14.82 %), respectively. In general, the model can provide accurate predictions within 12 h (R≥0.8) and errors can be maintained at under 1 m within 6 h of forecast lead time. However, the model also consistently overpredicted the maximum observed SWHs. From a comparison of LSTM with a third-generation wave model, Simulating Waves Nearshore (SWAN), it was determined that when using Hurricane Dorian as a case example, nowcasts were far more accurate with regards to the observations. This demonstrates that LSTM can be used to supplement, but perhaps not replace, computationally expensive numerical wave models for forecasting extreme wave heights. As such, addressing the fundamental problem of phase shifting and other errors in LSTM or other data-driven forecasting should receive greater scrutiny from Small Island Developing States. To improve models results, additional research should be geared towards improving single-point LSTM neural network training datasets by considering hurricane track and identifying the hurricane quadrant in which buoy observations are made.

Salinity as a key control on the diazotrophic community composition in the southern Baltic Sea - Christian Furbo Reeder; Ina Stoltenberg; Jamileh Javidpour; Carolin Regina Löscher (Ocean Science, 2022)
Over the next decade, the Baltic Sea is predicted to undergo severe changes including decreased salinity due to altered precipitation related to climate changes. This will likely impact the distribution and community composition of Baltic Sea dinitrogen-fixing (N2-fixing) microbes, among which heterocystous cyanobacteria are especially adapted to low salinities and may expand to waters with currently higher salinity, including the Danish Strait and Kattegat, while other high-salinity-adapted N2 fixers might decrease in abundance. In order to explore the impact of salinity on the distribution and activity of different diazotrophic clades, we followed the natural salinity gradient from the eastern Gotland and Bornholm basins through the Arkona Basin to the Kiel Bight and combined N2 fixation rate measurements with a molecular analysis of the diazotrophic community using the key functional marker gene for N2 fixation nifH, as well as the key functional marker genes anfD and vnfD, encoding for the two alternative nitrogenases. We detected N2 fixation rates between 0.7 and 6 nmol N L−1 d−1, and the diazotrophic community was dominated by the cyanobacterium related to Nodularia spumigena and the small unicellular, cosmopolitan cyanobacterium UCYN-A. Nodularia was present in gene abundances between 8.07 × 105 and 1.6 × 107 copies L−1 in waters with salinities of 10 and below, while UCYN-A reached gene abundances of up to 4.5 × 107 copies L−1 in waters with salinity above 10. Besides those two cyanobacterial diazotrophs, we found several clades of proteobacterial N2 fixers and alternative nitrogenase genes associated with Rhodopseudomonas palustris, a purple non-sulfur bacterium. Based on principal component analysis (PCA), salinity was identified as the primary parameter describing the diazotrophic distribution, while pH and temperature did not have a significant influence on the diazotrophic distribution. While this statistical analysis will need to be explored in direct experiments, it gives an indication for the future development of diazotrophy in a freshening Baltic Sea with UCYN-A retracting to more saline North Sea waters and heterocystous cyanobacteria expanding as salinity decreases.

The capabilities of Sentinel-MSI (2A/2B) and Landsat-OLI (8/9) in seagrass and algae species differentiation using spectral reflectance - Abderrazak Bannari; Thamer Salim Ali; Asma Abahussain (Ocean Science, 2022)
This paper assesses the reflectance difference values between the respective spectral bands in the visible and near-infrared (VNIR) of Sentinel 2A/2B Multi-Spectral Instrument (MSI) and Landsat 8/9 Operational Land Imager (OLI) sensors for seagrass, algae, and mixed species discrimination and monitoring in a shallow marine environment southeast of Bahrain Island in the Arabian Gulf. To achieve these, a field survey was conducted to collect samples of seawater, underwater sediments, seagrass (Halodule uninervis and Halophila stipulacea), and algae (green and brown). In addition, an experimental mode was established in a goniometric laboratory to simulate the marine environment, and spectral measurements were performed using an Analytical Spectral Devices (ASD) spectroradiometer. Measured spectra and their transformation using the continuum-removed reflectance spectral (CRRS) approach were analyzed to assess spectral separability among separate or mixed species at varying coverage rates. Afterward, the spectra were resampled and convolved in the solar-reflective spectral bands of MSI and OLI sensors and converted into water vegetation indices (WVIs) to investigate the potential of red, green, and blue bands for seagrass and algae species discrimination. The results of spectral and CRRS analyses highlighted the importance of the blue, green, and near-infrared (NIR) wavelengths for seagrass and algae detection and likely discrimination based on hyperspectral measurements. However, when resampled and convolved in MSI and OLI bands, spectral information loses the specific and unique absorption features and becomes more generalized and less precise. Therefore, relying on the multispectral bandwidth of MSI and OLI sensors, it is difficult or even impossible to differentiate or to map seagrass and algae individually at the species level. Instead of the red band, the integration of the blue or the green band in WVI increases their power to discriminate submerged aquatic vegetation (SAV), particularly the water adjusted vegetation index (WAVI), water enhanced vegetation index (WEVI), and water transformed difference vegetation index (WTDVI). These results corroborate the spectral and the CRRS analyses. However, despite the power of blue wavelength to penetrate deeper into the water, it also leads to a relative overestimation of dense SAV coverage due to more scattering in this part of the spectrum. Furthermore, statistical fits (p

Technical note: Turbulence measurements from a light autonomous underwater vehicle - Eivind H. Kolås; Tore Mo-Bjørkelund; Ilker Fer (Ocean Science, 2022)
A self-contained turbulence instrument from Rockland Scientific was installed on a light autonomous underwater vehicle (AUV) from OceanScan Marine Systems and Technology Lda. We report on the data quality and discuss limitations of dissipation estimated from two shear probes during a deployment in the Barents Sea in February 2021. The AUV mission lasted for 5 h, operating at a typical horizontal speed of 1.1 m s−1. The AUV was programmed to find and cross the maximum along-path thermal gradient at 10, 20 and 30 m depths along 4 km transects. Although the AUV vibrations contaminate the shear probe records, the noise is mitigated by removing vibration-induced components from shear spectra using the accelerometer signal measured in multiple directions. Dissipation rate estimates in the observed transects varied in the range 1×10-8 and 6×10-6 W kg−1, with the values from the two orthogonal probes typically in agreement to within a factor of 2. Dissipation estimates from the AUV show good agreement with nearby vertical microstructure profiles obtained from the ship during the transects, indicating that the turbulence measurements from the AUV are reliable for this relatively turbulent environment. However, the lowest reliable dissipation rates are limited to 5×10-8 W kg−1, making this setup unfit for use in quiescent environments.

Sea-level variability and change along the Norwegian coast between 2003 and 2018 from satellite altimetry, tide gauges, and hydrography - Fabio Mangini; Léon Chafik; Antonio Bonaduce; Laurent Bertino; Jan Even Ø. Nilsen (Ocean Science, 2022)
Sea-level variations in coastal areas can differ significantly from those in the nearby open ocean. Monitoring coastal sea-level variations is therefore crucial to understand how climate variability can affect the densely populated coastal regions of the globe. In this paper, we study the sea-level variability along the coast of Norway by means of in situ records, satellite altimetry data, and a network of eight hydrographic stations over a period spanning 16 years (from 2003 to 2018). At first, we evaluate the performance of the ALES-reprocessed coastal altimetry dataset (1 Hz posting rate) by comparing it with the sea-level anomaly from tide gauges over a range of timescales, which include the long-term trend, the annual cycle, and the detrended and deseasoned sea-level anomaly. We find that coastal altimetry and conventional altimetry products perform similarly along the Norwegian coast. However, the agreement with tide gauges in terms of trends is on average 6 % better when we use the ALES coastal altimetry data. We later assess the steric contribution to the sea level along the Norwegian coast. While longer time series are necessary to evaluate the steric contribution to the sea-level trends, we find that the sea-level annual cycle is more affected by variations in temperature than in salinity and that both temperature and salinity give a comparable contribution to the detrended and deseasoned sea-level variability along the entire Norwegian coast. A conclusion from our study is that coastal regions poorly covered by tide gauges can benefit from our satellite-based approach to study and monitor sea-level change and variability.

Occurrence of structural aluminium (Al) in marine diatom biological silica: visible evidence from microscopic analysis - Qian Tian; Dong Liu; Peng Yuan; Mengyuan Li; Weifeng Yang; Jieyu Zhou; Huihuang Wei; Junming Zhou; Haozhe Guo (Ocean Science, 2022)
The global marine biogeochemical cycle of aluminium (Al) is believed to be driven by marine diatoms, due to the uptake of dissolved Al (DAl) by living diatoms from surface seawater. The occurrence of Al in diatom biogenic silica (BSi) can inhibit the dissolution of BSi, thus benefiting the effects of the ballast role of diatoms in the biological pump and forming a coupled Si–Al biogeochemical cycle. However, the occurrence characteristic of Al in marine diatoms is still unclear. In particular, whether or not Al is incorporated into the structure of BSi of living diatoms is unrevealed, resulting in difficulties in understanding the biogeochemical behaviours of Al. In this study, Thalassiosira weissflogii, a widely distributed marine diatom in marginal seas, was selected as the model to evaluate the occurrence of structural Al in BSi based on culturing experiments with the addition of DAl. The structural Al in BSi was detected by combining focused ion beam (FIB) scanning electron microscopy and energy-dispersive X-ray spectroscopy (EDS) mapping analysis. Visible, direct evidence of structural Al in living BSi was obtained, and the distribution and content of this Al were revealed by the EDS-mapping analysis. The effects of structural Al on BSi dissolution–inhibition are discussed based on the content of this Al. The fundamental results indicate the significant contribution of marine diatoms to the biogeochemical migration of marine Al.

Passive tracer advection in the equatorial Pacific region: statistics, correlations and a model of fractional Brownian motion - Imre M. Jánosi; Amin Padash; Jason A. C. Gallas; Holger Kantz (Ocean Science, 2022)
Evaluating passive tracer advection is a common tool to study flow structures, particularly Lagrangian trajectories ranging from molecular scales up to the atmosphere and oceans. Here we report on numerical experiments in the region of the tropical Pacific (20∘ S–20∘ N), where 6600 tracer parcels are advected from a regular initial configuration (along a meridional line at 110∘ W between 15∘ S and 15∘ N) during periods of 1 year for 25 years altogether. We exploit AVISO surface flow fields and solve the kinematic equation for passive tracer movement in the 2D advection tests. We demonstrate that the strength of the advection defined by mean monthly westward displacements of the tracer clouds exhibit surprisingly large inter- and intra-annular variabilities. Furthermore, an analysis of cross-correlations between advection strength and the El-Niño and Southern Oscillation (SOI) indices reveal a significant anticorrelation between advection intensity and ONI (the Oceanic Niño Index) and a weaker positive correlation with SOI, both with a time lag of about 3 months (the two indices are strongly anticorrelated near real time). The statistical properties of advection (time-dependent mean squared displacement and first passage time distribution) suggest that the westward-moving tracers can be mapped into a simple 1D stochastic process, namely fractional Brownian motion. We fit the model parameters and show by numerical simulations of the fractional Brownian motion model that it is able to reproduce the observed statistical properties of the tracers' trajectories well. We argue that a traditional explanation based on the superposition of ballistic drift and a diffusion term yields different statistics and is incompatible with our observations.

Carbon and nitrogen dynamics in the coastal Sea of Japan inferred from 15 years of measurements of stable isotope ratios of Calanus sinicus - Ken-ichi Nakamura; Atsushi Nishimoto; Saori Yasui-Tamura; Yoichi Kogure; Misato Nakae; Naoki Iguchi; Haruyuki Morimoto; Taketoshi Kodama (Ocean Science, 2022)
Both nitrogen and carbon dynamics have changed in the Sea of Japan. We hypothesized that the carbon and nitrogen stable isotope ratios (δ13C and δ15N) of the copepod Calanus sinicus could record changes in the coastal environment of the Sea of Japan. Consequently, these isotope ratios were monitored during the spring at four stations from 2006 to 2020 to identify the changes in carbon and nitrogen dynamics. The δ13C values ranged from −24.7 ‰ to −15.0 ‰ and decreased from the spring bloom (February–March) to the post-bloom (June–July) seasons. These variations were attributed to changes in the physiology of both C. sinicus and phytoplankton δ13C contents. The δ15N values range from 2.8 ‰ to 8.8 ‰, indicating that C. sinicus is a secondary producer; the tendency of the δ15N values to increase from the bloom to the post-bloom seasons was attributable to an increase in the δ15N of phytoplankton. A generalized linear model (GLM) approach indicated that >70 % of the variations in δ13C can be explained by sea surface temperature (SST), sea surface chlorophyll a concentration (SSC), carbon:nitrogen ratio of C. sinicus (C/N ratio), and geographic differences. The residuals of δ13C in the GLM decreased yearly (−0.035 ‰ yr−1). The GLM for δ15N of C. sinicus indicated that δ15N varies with the stage or sex in addition to SST, SSC, C/N ratio, and geographic differences. The δ15N values of female C. sinicus and stage V copepodites were the lowest and highest, respectively. The residuals of δ15N in the GLM did not exhibit a significant interannual trend. These results suggest that the carbon isotope ratio in the secondary producer has linearly changed in the coastal Sea of Japan over the past 15 years. Moreover, the changes in carbon dynamics of this area have been recorded and observed to impact the marine ecosystem, while the nitrogen dynamics have not been recorded despite the increasing nitrogenous nutrient inputs in this sea.

Using machine learning and beach cleanup data to explain litter quantities along the Dutch North Sea coast - Mikael L. A. Kaandorp; Stefanie L. Ypma; Marijke Boonstra; Henk A. Dijkstra; Erik van Sebille (Ocean Science, 2022)
Coastlines potentially harbor a large part of litter entering the oceans, such as plastic waste. The relative importance of the physical processes that influence the beaching of litter is still relatively unknown. Here, we investigate the beaching of litter by analyzing a data set of litter gathered along the Dutch North Sea coast during extensive beach cleanup efforts between the years 2014 and 2019. This data set is unique in the sense that data are gathered consistently over various years by many volunteers (a total of 14 000) on beaches that are quite similar in substrate (sandy). This makes the data set valuable to identify which environmental variables play an important role in the beaching process and to explore the variability of beach litter concentrations. We investigate this by fitting a random forest machine learning regression model to the observed litter concentrations. We find that tides play an especially important role, where an increasing tidal variability and tidal height leads to less litter found on beaches. Relatively straight and exposed coastlines appear to accumulate more litter. The regression model indicates that transport of litter through the marine environment is also important in explaining beach litter variability. By understanding which processes cause the accumulation of litter on the coast, recommendations can be given for more effective removal of litter from the marine environment, such as organizing beach cleanups during low tides at exposed coastlines. We estimate that 16 500–31 200 kg (95 % confidence interval) of litter is located along the 365 km of Dutch North Sea coastline.

Global distribution and variability of subsurface chlorophyll a concentrations - Sayaka Yasunaka; Tsuneo Ono; Kosei Sasaoka; Kanako Sato (Ocean Science, 2022)
Chlorophyll a (Chl a) often exhibits a maximum concentration in the subsurface layer rather that at the surface. The depth of the Chl a maximum primarily depends on the balance between light penetration from the surface and the nutrient supply from the deep ocean. However, a global map of subsurface Chl a concentrations based on observations has not been presented yet. In this study, we integrate Chl a concentration data from recent biogeochemical floats and historical ship-based (and other) observations and present global maps of subsurface Chl a concentrations with related variables. The subsurface Chl a maximum was observed globally throughout the oceans: at depths greater than 80 m in the subtropics and tropics (30∘ S to 30∘ N); in the 40–80 m depth range in the tropics, in the Southern Ocean (south of 40∘ S), and at the midlatitudes (30–40∘ N/S) in the North Pacific; and at depths of less than 40 m in the northern subarctic (north of 40∘ N). The observed maxima all lie below the mixed-layer depth for the entire year in the subtropics and tropics and during summer in the midlatitudes and the northern subarctic. The depths of the subsurface Chl a maxima are greater than those of the photosynthetically active layer in the subtropics but shallower in the tropics and midlatitudes. In the subtropics, a seasonal increase in oxygen below the mixed layer implies substantial new biological production, which corresponds to 10 % of the net primary production in that region. During El Niño, subsurface Chl a concentrations are higher in the middle and eastern equatorial Pacific but lower to the west in comparison with La Niña, a pattern which is opposite to that on the surface. The spatiotemporal variability of the Chl a concentrations described here has implications to not only for the biogeochemical cycling in the ocean but also for understanding the thermal structure and dynamics of the ocean via absorption of shortwave radiation.

Autonomous methane seep site monitoring offshore western Svalbard: hourly to seasonal variability and associated oceanographic parameters - Knut Ola Dølven; Bénédicte Ferré; Anna Silyakova; Pär Jansson; Peter Linke; Manuel Moser (Ocean Science, 2022)
Improved quantification techniques of natural sources are needed to explain variations in atmospheric methane. In polar regions, high uncertainties in current estimates of methane release from the seabed remain. We present unique 10- and 3-month time series of bottom water measurements of physical and chemical parameters from two autonomous ocean observatories deployed at separate intense seabed methane seep sites (91 and 246 m depth) offshore western Svalbard from 2015 to 2016. Results show high short-term (100–1000 nmol L−1 within hours) and seasonal variation, as well as higher (2–7 times) methane concentrations compared to previous measurements. Rapid variability is explained by uneven distribution of seepage and changing ocean current directions. No overt influence of tidal hydrostatic pressure or water temperature variations on methane concentration was observed, but an observed negative correlation with temperature at the 246 m site fits with hypothesized seasonal blocking of lateral methane pathways in the sediments. Negative correlation between bottom water methane concentration (and variability) and wind forcing, concomitant with signs of weaker water column stratification, indicates increased potential for methane release to the atmosphere in fall and winter. We present new information about short- and long-term methane variability and provide a preliminary constraint on the uncertainties that arise in methane inventory estimates from this variability.

Responses of estuarine circulation to the morphological evolution in a convergent, microtidal estuary - Rui Zhang; Bo Hong; Lei Zhu; Wenping Gong; Heng Zhang (Ocean Science, 2022)
The Huangmaohai estuary (HE) is a funnel-shaped microtidal estuary in the west of the Pearl River Delta (PRD) in southern China. Since China reformed and opened up in 1978, extensive human activities have occurred and greatly changed the estuary's topography and modified its hydrodynamics. In this study, we examined the morphological evolution by analyzing remote sensing data with ArcGIS tools and studied the responses of hydrodynamics to the changes in topography from 1977 to 2010 by using the Delft3D model. We took the changes in estuarine circulation during neap tides in dry seasons as an example. The results show that human reclamation caused a narrowing of the estuary, and channel dredging deepened the estuary. These human activities changed both the longitudinal and lateral estuarine circulations. The longitudinal circulation was observed to increase with the deepening and narrowing of the estuary. The lateral circulation experienced changes in both the magnitude and pattern. The momentum balance analysis shows that when the depth and width changed simultaneously, the longitudinal estuarine circulation was modulated by both the channel deepening and width reduction, in which the friction, pressure gradient force, and advection terms were altered. The analysis of the longitudinal vortex dynamics indicates that the changes in the vertical shear of the longitudinal flow, lateral salinity gradient, and vertical mixing were responsible for the change in the lateral circulation. The changes in water depth are the dominant factors affecting lateral circulation intensity. This study has implications for sediment transport and morphological evolution in estuaries heavily impacted by human interventions.

Interannual variability of sea level in the southern Indian Ocean: local vs. remote forcing mechanisms - Marion Kersalé; Denis L. Volkov; Kandaga Pujiana; Hong Zhang (Ocean Science, 2022)
The subtropical southern Indian Ocean (SIO) has been described as one of the world's largest heat accumulators due to its remarkable warming during the past 2 decades. However, the relative contributions of remote (of Pacific origin) forcing and local wind forcing to the variability of heat content and sea level in the SIO have not been fully attributed. Here, we combine a general circulation model, an analytic linear reduced-gravity model, and observations to disentangle the spatial and temporal inputs of each forcing component on interannual to decadal timescales. A sensitivity experiment is conducted with artificially closed Indonesian straits to physically isolate the Indian Ocean and Pacific Ocean, intentionally removing the Indonesian Throughflow (ITF) influence on the Indian Ocean heat content and sea level variability. We show that the relative contribution of the signals originating in the equatorial Pacific vs. signals caused by local wind forcing to the interannual variability of sea level and heat content in the SIO is dependent on location within the basin (low latitude vs. midlatitude and western side vs. eastern side of the basin). The closure of the ITF in the numerical experiment reduces the amplitude of interannual-to-decadal sea level changes compared to the simulation with a realistic ITF. However, the spatial and temporal evolution of sea level patterns in the two simulations remain similar and correlated with El Niño–Southern Oscillation (ENSO). This suggests that these patterns are mostly determined by local wind forcing and oceanic processes, linked to ENSO via the “atmospheric bridge” effect. We conclude that local wind forcing is an important driver for the interannual changes of sea level, heat content, and meridional transports in the SIO subtropical gyre, while oceanic signals originating in the Pacific amplify locally forced signals.

Impact of acoustic Doppler current profiler (ADCP) motion on structure function estimates of turbulent kinetic energy dissipation rate - Brian D. Scannell; Yueng-Djern Lenn; Tom P. Rippeth (Ocean Science, 2022)
Turbulent mixing is a key process in the transport of heat, salt, and nutrients in the marine environment, with fluxes commonly derived directly from estimates of the turbulent kinetic energy dissipation rate, ε. Time series of ε estimates are therefore useful in helping to identify and quantify key biogeochemical processes. The velocity structure function method can be used to determine time series of ε estimates using along-beam velocity measurements from suitably configured acoustic Doppler current profilers (ADCPs). Shear in the background current can bias such estimates; therefore, standard practice is to deduct the mean or linear trend from the along-beam velocity over the period of an observation burst. This procedure is effective if the orientation of the ADCP to the current remains constant over the burst period. However, if the orientation of the ADCP varies, a proportion of the velocity difference between bins is retained in the structure function and the resulting ε estimates will be biased. Long-term observations from a mooring with three inline ADCPs show the heading oscillating with an angular range that depends on the flow speed: from large, slow oscillations at low flow speeds to smaller, higher-frequency oscillations at higher flow speeds. The mean tilt was also determined by the flow speed, whilst the tilt oscillation range was primarily determined by surface wave height. Synthesised along-beam velocity data for an ADCP subject to sinusoidal oscillation in a sheared flow indicate that the retained proportion of the potential bias is primarily determined by the angular range of the oscillation, with the impact varying between beams depending on the mean heading relative to the flow. Since the heading is typically unconstrained in a tethered mooring, heading oscillation is likely to be the most significant influence on the retained bias for a given level of shear. Use of an instrument housing designed to reduce oscillation would mitigate the impact, whilst if the shear is linear over the observation depth range, the bias can be corrected using a modified structure function method designed to correct for bias due to surface waves.

Wind-driven upwelling and surface nutrient delivery in a semi-enclosed coastal sea - Ben Moore-Maley; Susan E. Allen (Ocean Science, 2022)
Wind-driven upwelling is an important control on surface nutrients and water properties in stratified lakes and seas. In this study, a high-resolution biophysical coupled model is used to investigate upwelling in the Strait of Georgia on the Canadian Pacific coast. The model is forced with surface winds from a high-resolution atmospheric forecast and has been tuned in previous studies to reproduce extensive observations of water level, temperature, salinity, nutrients and chlorophyll with competitive skill relative to similar models of the study region. A total of 5 years of hourly surface nitrate and temperature fields are analyzed in order to characterize the dominant upwelling patterns of the basin. A prevailing along-axis wind pattern steered by mountainous topography produces episodic upwelling along the western shore during the spring and fall southeasterlies and along the eastern shore during the summer northwesterlies, as indicated by positive nitrate anomalies. Principal component analysis reveals that these cross-axis upwelling patterns account for nearly one-third of the surface nitrate variance during the summer productive season. By contrast, nearly half of the surface temperature variance over the same period is dominated by a single, combined mixing and diurnal heating–cooling pattern. The principal components associated with these patterns correlate with along-axis wind stress in a manner consistent with these physical interpretations. The cross-axis upwelling response to wind is similar to other dynamically wide basins where the baroclinic Rossby deformation radius is smaller than the basin width. However, the nitrate anomaly during upwelling along the eastern shore is stronger in the northern basin, which may be indicative of an along-axis pycnocline tilt or an effect of the background along-axis stratification gradient due to the Fraser River. Our findings highlight an important spatiotemporal consideration for future ecosystem monitoring.

Characterization of the Atlantic Water and Levantine Intermediate Water in the Mediterranean Sea using 20 years of Argo data - Giusy Fedele; Elena Mauri; Giulio Notarstefano; Pierre Marie Poulain (Ocean Science, 2022)
Atlantic Water (AW) and Levantine Intermediate Water (LIW) are important water masses that play a crucial role in the internal variability of the Mediterranean thermohaline circulation. To be more specific, their variability and interaction, along with other water masses that characterize the Mediterranean basin, such as the Western Mediterranean Deep Water (WMDW), contribute to modify the Mediterranean Outflow through the Strait of Gibraltar, and hence they may influence the stability of the global thermohaline circulation. This work aims to characterize AW and LIW in the Mediterranean Sea, taking advantage of the large observational dataset (freely available on https://argo.ucsd.edu, https://www.ocean-ops.org, last access: 17 January 2022; Wong et al., 2020) provided by Argo floats from 2001 to 2019. AW and LIW were identified using different diagnostic methods, highlighting the inter-basin variability and the strong zonal gradient that both denote the two water masses in this marginal sea. Their temporal variability was also investigated over the last 2 decades, providing a more robust view of AW and LIW characteristics, which have only been investigated using very short periods in previous studies due to a lack of data. A clear salinification and warming trend characterize AW and LIW over the last 2 decades (∼ 0.007 ± 0.140 and 0.006 ± 0.038 yr−1; 0.026 ± 0.715 and 0.022 ± 0.232 ∘C yr−1, respectively). The salinity and temperature trends found at sub-basin scale are in good agreement with previous results. The strongest trends are found in the Adriatic basin in the properties of both AW and LIW.

Components of 21 years (1995–2015) of absolute sea level trends in the Arctic - Carsten Bjerre Ludwigsen; Ole Baltazar Andersen; Stine Kildegaard Rose (Ocean Science, 2022)
The Arctic Ocean is at the frontier of the fast-changing climate in the northern latitudes, and sea level trends are a bulk measure of ongoing processes related to climate change. Observations of sea level in the Arctic Ocean are nonetheless difficult to validate with independent measurements, and this is globally the region where the sea level trend (SLT) is most uncertain. The aim of this study is to create a satellite-independent reconstruction of Arctic SLT, as it is observed by altimetry and tide gauges (TGs). Previous studies use Gravity Recovery and Climate Experiment (GRACE) observations to estimate the manometric (mass component of) SLT. GRACE estimates, however, are challenged by large mass changes on land, which are difficult to separate from much smaller ocean mass changes. Furthermore, GRACE is not available before 2003, which significantly limits the period and makes the trend more vulnerable to short-term changes. As an alternative approach, this study estimates the climate-change-driven Arctic manometric SLT from the Arctic sea level fingerprints of glaciers, Greenland, Antarctica and glacial isostatic adjustment (GIA) with the addition of the long-term inverse barometer (IB) effect. The halosteric and thermosteric components complete the reconstructed Arctic SLT and are estimated by interpolating 300 000 temperature (T) and salinity (S) in situ observations. The SLT from 1995–2015 is compared to the observed SLT from altimetry and 12 selected tide gauges (TGs) corrected for vertical land movement (VLM). The reconstructed estimate manifests the salinity-driven halosteric component as dominating the spatial SLT pattern with variations between −7 and 10 mm yr−1. The manometric SLT in comparison is estimated to be 1–2 mm yr−1 for most of the Arctic Ocean. The reconstructed SLT shows a larger sea level rise in the Beaufort Sea compared to altimetry, an issue that is also identified by previous studies. There is a TG-observed sea level rise in the Siberian Arctic in contrast to the sea level fall from the reconstructed and altimetric estimate. From 1995–2015 the reconstructed SLT agrees within the 68 % confidence interval with the SLT from observed altimetry in 87 % of the Arctic between 65∘ N and 82∘ N (R=0.50) and with 5 of 12 TG-derived (VLM-corrected) SLT estimates. The residuals are seemingly smaller than results from previous studies using GRACE estimates and modeled T–S data. The spatial correlation of the reconstructed SLT to altimetric SLT during the GRACE period (2003–2015) is R=0.38 and R=0.34/R=0.37 if GRACE estimates are used instead of the constructed manometric component. Thus, the reconstructed manometric component is suggested as a legitimate alternative to GRACE that can be projected into the past and future.

Refined estimates of water transport through the Åland Sea in the Baltic Sea - Antti Westerlund; Elina Miettunen; Laura Tuomi; Pekka Alenius (Ocean Science, 2022)
Water exchange through the Åland Sea (in the Baltic Sea) greatly affects the environmental conditions in the neighbouring Gulf of Bothnia. Recently observed changes in the eutrophication status of the Gulf of Bothnia may be connected to changing nutrient fluxes through the Åland Sea. Pathways and variability of sub-halocline northward-bound flows towards the Bothnian Sea are important for these studies. While the general nature of the water exchange is known, that knowledge is based on only a few studies that are somewhat limited in detail. Notably, no high-resolution modelling studies of water exchange in the Åland Sea area have been published. In this study, we present a configuration of the NEMO 3D hydrodynamic model for the Åland Sea–Archipelago Sea area at around 500 m horizontal resolution. We then use it to study the water exchange in the Åland Sea and volume transports through the area. We first ran the model for the years 2013–2017 and validated the results, with a focus on the simulated current fields. We found that the model reproduced current direction distributions and layered structure of currents in the water column with reasonably good accuracy. Next, we used the model to calculate volume transports across several transects in the Åland Sea. These calculations provided new details about water transport in the area. Time series of monthly mean volume transports showed consistent northward transport in the deep layer. In the surface layer there was more variability: while net transport was towards the south, in several years some months in late summer or early autumn showed net transport to the north. Furthermore, based on our model calculations, it seems that dynamics in the Lågskär Deep are more complex than has been previously understood. While Lågskär Deep is the primary route of deep-water exchange, a significant volume of deep water still enters the Åland Sea through the depression west of the Lågskär Deep. Better spatial and temporal coverage of current measurements is needed to further refine the understanding of water exchange in the area. Future studies of transport and nutrient dynamics will eventually enable a deeper understanding of eutrophication changes in the Gulf of Bothnia.

A framework to evaluate and elucidate the driving mechanisms of coastal sea surface p CO 2 seasonality using an ocean general circulation model (MOM6-COBALT) - Alizée Roobaert; Laure Resplandy; Goulven G. Laruelle; Enhui Liao; Pierre Regnier (Ocean Science, 2022)
The temporal variability of the sea surface partial pressure of CO2 (pCO2) and the underlying processes driving this variability are poorly understood in the coastal ocean. In this study, we tailor an existing method that quantifies the effects of thermal changes, biological activity, ocean circulation and freshwater fluxes to examine seasonal pCO2 changes in highly variable coastal environments. We first use the Modular Ocean Model version 6 (MOM6) and biogeochemical module Carbon Ocean Biogeochemistry And Lower Trophics version 2 (COBALTv2) at a half-degree resolution to simulate coastal CO2 dynamics and evaluate them against pCO2 from the Surface Ocean CO2 Atlas database (SOCAT) and from the continuous coastal pCO2 product generated from SOCAT by a two-step neuronal network interpolation method (coastal Self-Organizing Map Feed-Forward neural Network SOM-FFN, Laruelle et al., 2017). The MOM6-COBALT model reproduces the observed spatiotemporal variability not only in pCO2 but also in sea surface temperature, salinity and nutrients in most coastal environments, except in a few specific regions such as marginal seas. Based on this evaluation, we identify coastal regions of “high” and “medium” agreement between model and coastal SOM-FFN where the drivers of coastal pCO2 seasonal changes can be examined with reasonable confidence. Second, we apply our decomposition method in three contrasted coastal regions: an eastern (US East Coast) and a western (the Californian Current) boundary current and a polar coastal region (the Norwegian Basin). Results show that differences in pCO2 seasonality in the three regions are controlled by the balance between ocean circulation and biological and thermal changes. Circulation controls the pCO2 seasonality in the Californian Current; biological activity controls pCO2 in the Norwegian Basin; and the interplay between biological processes and thermal and circulation changes is key on the US East Coast. The refined approach presented here allows the attribution of pCO2 changes with small residual biases in the coastal ocean, allowing for future work on the mechanisms controlling coastal air–sea CO2 exchanges and how they are likely to be affected by future changes in sea surface temperature, hydrodynamics and biological dynamics.

Arctic sea level variability from high-resolution model simulations and implications for the Arctic observing system - Guokun Lyu; Nuno Serra; Meng Zhou; Detlef Stammer (Ocean Science, 2022)
Two high-resolution model simulations are used to investigate the spatiotemporal variability of the Arctic Ocean sea level. The model simulations reveal barotropic sea level variability at periods of

On the circulation, water mass distribution, and nutrient concentrations of the western Chukchi Sea - Jaclyn Clement Kinney; Karen M. Assmann; Wieslaw Maslowski; Göran Björk; Martin Jakobsson; Sara Jutterström; Younjoo J. Lee; Robert Osinski; Igor Semiletov; Adam Ulfsbo; Irene Wåhlström; Leif G. Anderson (Ocean Science, 2022)
Substantial amounts of nutrients and carbon enter the Arctic Ocean from the Pacific Ocean through the Bering Strait, distributed over three main pathways. Water with low salinities and nutrient concentrations takes an eastern route along the Alaskan coast, as Alaskan Coastal Water. A central pathway exhibits intermediate salinity and nutrient concentrations, while the most nutrient-rich water enters the Bering Strait on its western side. Towards the Arctic Ocean, the flow of these water masses is subject to strong topographic steering within the Chukchi Sea with volume transport modulated by the wind field. In this contribution, we use data from several sections crossing Herald Canyon collected in 2008 and 2014 together with numerical modelling to investigate the circulation and transport in the western part of the Chukchi Sea. We find that a substantial fraction of water from the Chukchi Sea enters the East Siberian Sea south of Wrangel Island and circulates in an anticyclonic direction around the island. This water then contributes to the high-nutrient waters of Herald Canyon. The bottom of the canyon has the highest nutrient concentrations, likely as a result of addition from the degradation of organic matter at the sediment surface in the East Siberian Sea. The flux of nutrients (nitrate, phosphate, and silicate) and dissolved inorganic carbon in Bering Summer Water and Winter Water is computed by combining hydrographic and nutrient observations with geostrophic transport referenced to lowered acoustic Doppler current profiler (LADCP) and surface drift data. Even if there are some general similarities between the years, there are differences in both the temperature–salinity and nutrient characteristics. To assess these differences, and also to get a wider temporal and spatial view, numerical modelling results are applied. According to model results, high-frequency variability dominates the flow in Herald Canyon. This leads us to conclude that this region needs to be monitored over a longer time frame to deduce the temporal variability and potential trends.

Plastics in the Indian Ocean – sources, transport, distribution, and impacts - Charitha Pattiaratchi; Mirjam van der Mheen; Cathleen Schlundt; Bhavani E. Narayanaswamy; Appalanaidu Sura; Sara Hajbane; Rachel White; Nimit Kumar; Michelle Fernandes; Sarath Wijeratne (Ocean Science, 2022)
Plastic debris is the most common and exponentially increasing human pollutant in the world's ocean. The distribution and impact of plastic in the Pacific and Atlantic oceans have been the subject of many publications but not so the Indian Ocean (IO). Some of the IO rim countries have the highest population densities globally and mismanagement of plastic waste is of concern in many of these rim states. Some of the most plastic-polluted rivers empty into the IO, with all this suggesting that the IO receives a tremendous amount of plastic debris each year. However, the concentration, distribution, and impacts of plastics in the IO are poorly understood as the region is under-sampled compared to other oceans. In this review, we discuss sources and sinks, which are specific to the IO. We also discuss unique atmospheric, oceanographic, and topographic features of the IO that control plastic distribution, such as reversing wind directions due to the monsoon, fronts, and upwelling regions. We identify hotspots of possible plastic accumulation in the IO, which differ between the two hemispheres. In the southern IO, plastics accumulate in a garbage patch in the subtropical gyre. However, this garbage patch is not well defined, and plastics may leak into the southern Atlantic or the Pacific Ocean. There is no subtropical gyre and associated garbage in the northern IO due to the presence of landmasses. Instead, the majority of buoyant plastics most likely end up on coastlines. Finally, we identify the vast knowledge gaps concerning plastics in the IO and point to the most pressing topics for future investigation.

Swell hindcast statistics for the Baltic Sea - Jan-Victor Björkqvist; Siim Pärt; Victor Alari; Sander Rikka; Elisa Lindgren; Laura Tuomi (Ocean Science, 2021)
The classic characterisation of swell as regular, almost monochromatic, wave trains does not necessarily accurately describe swell in water bodies shielded from the oceanic wave climate. In such enclosed areas the locally generated swell waves still contribute to processes at the air and seabed interfaces, and their presence can be quantified by partitioning wave components based on their speed relative to the wind. We present swell statistics for the semi-enclosed Baltic Sea using 20 years of swell-partitioned model data. The swell significant wave height was mostly under 2 m, and in the winter (DJF) the mean significant swell height was typically less than 0.4 m; higher swell was found in limited nearshore areas. Swell waves were typically short (under 5 s), with mean periods over 8 s being rare. In open-sea areas the average ratio of swell energy (to total energy) was mostly below 0.4 – significantly less than in the World Ocean. Certain coastal areas were swell dominated over half the time, mostly because of weak winds (U

Assimilating realistically simulated wide-swath altimeter observations in a high-resolution shelf-seas forecasting system - Robert R. King; Matthew J. Martin (Ocean Science, 2021)
The impact of assimilating simulated wide-swath altimetry observations from the upcoming Surface Water and Ocean Topography (SWOT) mission is assessed using observing system simulation experiments (OSSEs). These experiments use the Met Office 1.5 km resolution North West European Shelf analysis and forecasting system. In an effort to understand the importance of future work to account for correlated errors in the data assimilation scheme, we simulate SWOT observations with and without realistic correlated errors. These are assimilated in OSSEs along with simulated observations of the standard observing network, also with realistic errors added. It was found that while the assimilation of SWOT observations without correlated errors reduced the RMSE (root mean squared error) in sea surface height (SSH) and surface current speeds by up to 20 %, the inclusion of correlated errors in the observations degraded both the SSH and surface currents, introduced an erroneous increase in the mean surface currents and degraded the subsurface temperature and salinity. While restricting the SWOT data to the inner half of the swath and applying observation averaging with a 5 km radius negated most of the negative impacts, it also severely limited the positive impacts. To realise the full benefits in the prediction of the ocean mesoscale offered by wide-swath altimetry missions, it is crucial that methods to ameliorate the effects of correlated errors in the processing of the SWOT observations and account for the correlated errors in the assimilation are implemented.

Nutrient ratios driven by vertical stratification regulate phytoplankton community structure in the oligotrophic western Pacific Ocean - Zhuo Chen; Jun Sun; Ting Gu; Guicheng Zhang; Yuqiu Wei (Ocean Science, 2021)
The stratification of the upper oligotrophic ocean has a direct impact on biogeochemistry by regulating the components of the upper-ocean environment that are critical to biological productivity, such as light availability for photosynthesis and nutrient supply from the deep ocean. We investigated the spatial distribution pattern and diversity of phytoplankton communities in the western Pacific Ocean (WPO) in the autumn of 2016, 2017, and 2018. Our results showed the phytoplankton community structure mainly consisted of cyanobacteria, diatoms, and dinoflagellates, while the abundance of Chrysophyceae was negligible. Phytoplankton abundance was high from the equatorial region to 10∘ N and decreased with increasing latitude in spatial distribution. Phytoplankton also showed a strong variation in the vertical distribution. The potential influences of physicochemical parameters on phytoplankton abundance were analyzed by a structural equation model (SEM) to determine nutrient ratios driven by vertical stratification to regulate phytoplankton community structure in the typical oligotrophic ocean. Regions with strong vertical stratification were more favorable for cyanobacteria, whereas weak vertical stratification was more conducive to diatoms and dinoflagellates. Our study shows that stratification is a major determinant of phytoplankton community structure and highlights that physical processes in the ocean control phytoplankton community structure by driving the balance of chemical elements, providing a database to better predict models of changes in phytoplankton community structure under future ocean scenarios.

Rectified tidal transport in Lofoten–Vesterålen, northern Norway - Eli Børve; Pål Erik Isachsen; Ole Anders Nøst (Ocean Science, 2021)
Vestfjorden in northern Norway, a major spawning ground for the northeast Arctic cod, is sheltered from the continental shelf and open ocean by the Lofoten–Vesterålen archipelago. The archipelago, however, is well known for hosting strong and vigorous tidal currents in its many straits, currents that can produce significant time-mean tracer transport from Vestfjorden to the shelf outside. We use a purely tidally driven unstructured-grid ocean model to look into non-linear tidal dynamics and the associated tracer transport through the archipelago. Of particular interest are two processes: tidal pumping through the straits and tidal rectification around islands. The most prominent tracer transport is caused by tidal pumping through the short and strongly non-linear straits Nordlandsflaget and Moskstraumen near the southern tip of the archipelago. Here, tracers from Vestfjorden are transported tens of kilometers westward out on the outer shelf. Further north, weaker yet notable tidal pumping also takes place through the longer straits Nappstraumen and Gimsøystraumen. The other main transport route out of Vestfjorden is south of the island of Røst. Here, the transport is primarily due to tracer advection by rectified anticyclonic currents around the island. There is also an anticyclonic circulation cell around the island group Mosken–Værøy, and both cells have flow speeds up to 0.2 m s−1, magnitudes similar to the observed background currents in the region. These high-resolution simulations thus emphasize the importance of non-linear tidal dynamics for transport of floating particles, like cod eggs and larvae, in the region.

Progress in understanding of Indian Ocean circulation, variability, air–sea exchange, and impacts on biogeochemistry - Helen E. Phillips; Amit Tandon; Ryo Furue; Raleigh Hood; Caroline C. Ummenhofer; Jessica A. Benthuysen; Viviane Menezes; Shijian Hu; Ben Webber; Alejandra Sanchez-Franks; Deepak Cherian; Emily Shroyer; et al. (Ocean Science, 2021)
Over the past decade, our understanding of the Indian Ocean has advanced through concerted efforts toward measuring the ocean circulation and air–sea exchanges, detecting changes in water masses, and linking physical processes to ecologically important variables. New circulation pathways and mechanisms have been discovered that control atmospheric and oceanic mean state and variability. This review brings together new understanding of the ocean–atmosphere system in the Indian Ocean since the last comprehensive review, describing the Indian Ocean circulation patterns, air–sea interactions, and climate variability. Coordinated international focus on the Indian Ocean has motivated the application of new technologies to deliver higher-resolution observations and models of Indian Ocean processes. As a result we are discovering the importance of small-scale processes in setting the large-scale gradients and circulation, interactions between physical and biogeochemical processes, interactions between boundary currents and the interior, and interactions between the surface and the deep ocean. A newly discovered regional climate mode in the southeast Indian Ocean, the Ningaloo Niño, has instigated more regional air–sea coupling and marine heatwave research in the global oceans. In the last decade, we have seen rapid warming of the Indian Ocean overlaid with extremes in the form of marine heatwaves. These events have motivated studies that have delivered new insight into the variability in ocean heat content and exchanges in the Indian Ocean and have highlighted the critical role of the Indian Ocean as a clearing house for anthropogenic heat. This synthesis paper reviews the advances in these areas in the last decade.

The diurnal cycle of p CO 2 in the coastal region of the Baltic Sea - Martti Honkanen; Jens Daniel Müller; Jukka Seppälä; Gregor Rehder; Sami Kielosto; Pasi Ylöstalo; Timo Mäkelä; Juha Hatakka; Lauri Laakso (Ocean Science, 2021)
The direction and magnitude of carbon dioxide fluxes between the atmosphere and the sea are regulated by the gradient in the partial pressure of carbon dioxide (pCO2) across the air–sea interface. Typically, observations of pCO2 at the sea surface are carried out by using research vessels and ships of opportunity, which usually do not resolve the diurnal cycle of pCO2 at a given location. This study evaluates the magnitude and driving processes of the diurnal cycle of pCO2 in a coastal region of the Baltic Sea. We present pCO2 data from July 2018 to June 2019 measured in the vicinity of the island of Utö at the outer edge of the Archipelago Sea, and quantify the relevant physical, biological, and chemical processes controlling pCO2. The highest monthly median of diurnal pCO2 variability (31 µatm) was observed in August and predominantly driven by biological processes. Biological fixation and mineralization of carbon led to sinusoidal diurnal pCO2 variations, with a maximum in the morning and a minimum in the afternoon. Compared with the biological carbon transformations, the impacts of air–sea fluxes and temperature changes on pCO2 were small, with their contributions to the monthly medians of diurnal pCO2 variability being up to 12 and 5 µatm, respectively. During upwelling events, short-term pCO2 variability (up to 500 µatm within a day) largely exceeded the usual diurnal cycle. If the net annual air–sea flux of carbon dioxide at our study site and for the sampled period is calculated based on a data subset that consists of only one regular measurement per day, the bias in the net exchange depends on the sampling time and can amount up to ±12 %. This finding highlights the importance of continuous surface pCO2 measurements at fixed locations for the assessment of the short-term variability of the carbonate system and the correct determination of air–sea CO2 fluxes.

A tidally driven fjord-like strait close to an amphidromic region - Sissal Vágsheyg Erenbjerg; Jon Albretsen; Knud Simonsen; Erna Lava Olsen; Eigil Kaas; Bogi Hansen (Ocean Science, 2021)
The strait studied in this paper, “Sundalagið Norður”, is the northern part of a narrow body of seawater separating the two largest islands in the Faroe Islands (Faroes). It has shallow sills in both ends and considerably deeper waters in between. South of the southern end of the strait there is an amphidromic region for the semidiurnal tides so that the tidal range is much lower south of the strait than north of it. The resulting tidal forcing generates periodically varying inflow of seawater across the northern sill, but only a part of that manages to cross the narrow and shallow southern sill. Combined with a large input of freshwater, this gives the strait a fjord-like character. To investigate how this fjord-like character affects the circulation within the strait and its exchanges with outside waters, a pilot project was initiated to simulate the dynamics of the strait with a high-resolution ocean model for a month. The model simulations show clearly the dominance of tidal forcing over freshwater (estuarine) and wind on timescales up to a day. On longer timescales, the simulations indicate systematic variations in the net flows (averaged over a diurnal tidal period) through both the upper and deeper layers. These long-period variations of net flow in the model simulations are forced by sea level differences between both ends of the strait generated by the dominant fortnightly and monthly tidal constituents (Mf, MSf, Mm, MSm). Harmonic analysis of sea level records from two tide gauges located off each end of the strait demonstrates that this behaviour is not a model artefact and it has pronounced effects on the strait. Not only does it induce long-period (mainly fortnightly) variations in the net flow through the strait, but it also generates variations in the estuarine characteristics. According to the model simulations, periods with net southward flow, typically lasting a week, have a strait-like character with net southward flow almost everywhere. Periods with net northward flow, in contrast, have a more fjord-like character with stronger salinity stratification and a southward counter-flow in the deep layer. This also induces a large difference in renewal rate of the deep water between the two periods, which is important to consider for human utilization of the strait, especially the local aquaculture plant. The combination of topographic, freshwater, and tidal characteristics creating these long-period variations is rather unusual, and it is not known whether similar systems exist elsewhere, but the long-period variations tend to be masked by the stronger semidiurnal and diurnal variations and may easily be overlooked.

Mean sea level and tidal change in Ireland since 1842: a case study of Cork - David T. Pugh; Edmund Bridge; Robin Edwards; Peter Hogarth; Guy Westbrook; Philip L. Woodworth; Gerard D. McCarthy (Ocean Science, 2021)
Knowledge of regional changes in mean sea level and local changes in tides are crucial to inform effective climate adaptation. An essential element is the availability of accurate observations of sea level. Sea level data in the Republic of Ireland, prior to the establishment of the National Tide Gauge Network in the mid-2000s, are limited but belie a wealth of historical data available in archival form. In this study, we digitize records located in Cork Harbour, Ireland, from 1842 and show how short-duration (6–8 weeks), high-quality data with a large interval to the present can accurately inform tidal and mean sea level changes. We consider error sources in detail. We estimate for the main M2 tidal constituent that the accuracy of these historical measurements is 1 % and 2 min for amplitude and phase, respectively, once adjustments for seasonal and nodal effects are made. Our mean sea level estimates are accurate to the 2 cm level, once adjustments for atmospheric and seasonal effects are made. Our results show tidal stability with a 2 % change in the amplitude of the M2 component, 4 min change in the phase over a period of 177 years, and mean sea level rise of 40 cm in the Cork Harbour area from 1842 to 2019, approximately in line with global mean sea level trends plus local glacial isostatic adjustment. More broadly, we show that with careful seasonal, nodal, and atmospheric corrections, together with knowledge of benchmark provenance, these historic, survey-oriented data can accurately inform of sea level changes.

Dynamics of fortnightly water level variations along a tide-dominated estuary with negligible river discharge - Erwan Garel; Ping Zhang; Huayang Cai (Ocean Science, 2021)
Observations indicate that the fortnightly fluctuations in the mean amplitude of water level increase in the upstream direction along the lower half of a tide-dominated estuary (the Guadiana Estuary), with negligible river discharge, but remain constant upstream. Analytical solutions reproducing the semi-diurnal wave propagation shows that this pattern results from reflection effects at the estuary head. The phase difference between velocity and elevation increases from the mouth to the head (where the wave has a standing nature) as the timing of high and low water levels come progressively closer to slack water. Thus, the tidal (flood–ebb) asymmetry in discharge is reduced in the upstream direction. It becomes negligible along the upper estuary half as the mean sea level remains constant despite increased friction due to wave shoaling. Observations of a flat mean water level along a significant portion of an upper estuary suggest a standing wave character and, thus, indicate significant reflection of the propagating semi-diurnal wave at the head. Details of the analytical model show that changes in the mean depth or length of semi-arid estuaries, in particular for macrotidal locations, affect the fortnightly tide amplitude and, thus, the upstream mass transport and inundation regime. This has significant potential impacts on the estuarine environment in terms of ecosystem management.

On the low western Pacific sea levels observed prior to strong East Pacific El Niños - David J. Webb (Ocean Science, 2021)
Research, based on both observations and ocean models, has indicated that strong East Pacific El Niños are preceded by the development of unusually low sea levels at the western end of the North Equatorial Trough early in the year. This results in an increase in the strength of the North Equatorial Counter Current which, aided by low sea levels due to the annual Rossby wave, then transports West Pacific Warm Pool water to the central and eastern Pacific. Here, an ocean model is used to investigate the factors affecting sea level prior to the 1982–1983 East Pacific El Niño. The results indicate that during this period the reduction in sea level was caused by Ekman pumping, due to the local winds, acting over a period of many months. The north–south distances involved mean that such upwelling can be connected to the westerly wind phase of Madden–Julian Oscillations. Since the 1980s, the quality and quantity of global wind measurement has steadily improved. So, if the Ekman pumping hypothesis is correct, the better-quality wind data available prior to the 1997–1998 and 2015–2016 East Pacific El Niños should also show unusually large Ekman pumping in the same region, over a similar long period. This is tested and found to be correct. However, a study of the integrated Ekman pumping for the period 1976 to 2015 indicates that in some years similar events did not develop despite a comparable amount of pumping during the first part of the year. The results imply that significant Ekman pumping early in the year is a necessary but not sufficient requirement for the development of a strong East Pacific El Niño.

Background stratification impacts on internal tide generation and abyssal propagation in the western equatorial Atlantic and the Bay of Biscay - Simon Barbot; Florent Lyard; Michel Tchilibou; Loren Carrere (Ocean Science, 2021)
The forthcoming SWOT altimetric missions aim to resolve the mesoscale with an unprecedented spatial resolution and swath. However, high-frequency processes, such as tides, are undersampled in time and aliased onto lower frequencies, so they need to be corrected properly. Unlike barotropic tides, internal tides (ITs) are not completely stationary and have significant temporal variability due to their interactions with the ocean circulation and the stratification variability. Stratification changes impact both the generation and the propagation of ITs. The present study proposes a methodology to quantify the impacts of background stratification using a clustering method for the classification of a broad range of stratification and idealized modeling of ITs in the frequency domain. The methodology is successfully tested in the western equatorial Atlantic and in the Bay of Biscay. For the western equatorial Atlantic, a single pycnocline is observed and only the two first vertical modes of ITs have significant amplitudes. With no variation in the stratification intensity, the variation in the depth of this single pycnocline linearly impacts the elevation amplitude, energy fluxes and surface wavelength of the two modes. In the Bay of Biscay, there is a permanent deep pycnocline and secondary seasonal pycnoclines near the surface. No proxy have been found to describe the changes in ITs, so a seasonal climatology is explored. The seasonality of the stratification strongly affects the elevation amplitudes as well as the energy fluxes of modes 1, 2 and 3. The distribution of the modes vary with the background stratification, changing the horizontal scales of the ITs.

Defining Southern Ocean fronts using unsupervised classification - Simon D. A. Thomas; Daniel C. Jones; Anita Faul; Erik Mackie; Etienne Pauthenet (Ocean Science, 2021)
Oceanographic fronts are transitions between thermohaline structures with different characteristics. Such transitions are ubiquitous, and their locations and properties affect how the ocean operates as part of the global climate system. In the Southern Ocean, fronts have classically been defined using a small number of continuous, circumpolar features in sea surface height or dynamic height. Modern observational and theoretical developments are challenging and expanding this traditional framework to accommodate a more complex view of fronts. Here, we present a complementary new approach for calculating fronts using an unsupervised classification method called Gaussian mixture modelling (GMM) and a novel inter-class parameter called the I-metric. The I-metric approach produces a probabilistic view of front location, emphasising the fact that the boundaries between water masses are not uniformly sharp across the entire Southern Ocean. The I-metric approach uses thermohaline information from a range of depth levels, making it more general than approaches that only use near-surface properties. We train the GMM using an observationally constrained state estimate in order to have more uniform spatial and temporal data coverage. The probabilistic boundaries defined by the I-metric roughly coincide with several classically defined fronts, offering a novel view of this structure. The I-metric fronts appear to be relatively sharp in the open ocean and somewhat diffuse near large topographic features, possibly highlighting the importance of topographically induced mixing. For comparison with a more localised method, we also use an edge detection approach for identifying fronts. We find a strong correlation between the edge field of the leading principal component and the zonal velocity; the edge detection method highlights the presence of jets, which are supported by thermal wind balance. This more localised method highlights the complex, multiscale structure of Southern Ocean fronts, complementing and contrasting with the more domain-wide view offered by the I-metric. The Sobel edge detection method may be useful for defining and tracking smaller-scale fronts and jets in model or reanalysis data. The I-metric approach may prove to be a useful method for inter-model comparison, as it uses the thermohaline structure of those models instead of tracking somewhat ad hoc values of sea surface height and/or dynamic height, which can vary considerably between models. In addition, the general I-metric approach allows front definitions to shift with changing temperature and salinity structures, which may be useful for characterising fronts in a changing climate.

Using feature-based verification methods to explore the spatial and temporal characteristics of the 2019 chlorophyll- a bloom season in a model of the European Northwest Shelf - Marion Mittermaier; Rachel North; Jan Maksymczuk; Christine Pequignet; David Ford (Ocean Science, 2021)
Two feature-based verification methods, thus far only used for the diagnostic evaluation of atmospheric models, have been applied to compare ∼7 km resolution pre-operational analyses of chlorophyll-a (Chl-a) concentrations to a 1 km gridded satellite-derived Chl-a concentration product. The aim of this study was to assess the value of applying such methods to ocean models. Chl-a bloom objects were identified in both data sets for the 2019 bloom season (1 March to 31 July). These bloom objects were analysed as discrete (2-D) spatial features, but also as space–time (3-D) features, providing the means of defining the onset, duration and demise of distinct bloom episodes and the season as a whole. The new feature-based verification methods help reveal that the model analyses are not able to represent small coastal bloom objects, given the coarser definition of the coastline, also wrongly producing more bloom objects in deeper Atlantic waters. Model analyses' concentrations are somewhat higher overall. The bias manifests itself in the size of the model analysis bloom objects, which tend to be larger than the satellite-derived bloom objects. The onset of the bloom season is delayed by 26 d in the model analyses, but the season also persists for another month beyond the diagnosed end. The season was diagnosed to be 119 d long in the model analyses, compared to 117 d from the satellite product. Geographically, the model analyses and satellite-derived bloom objects do not necessarily exist in a specific location at the same time and only overlap occasionally.

Winter observations alter the seasonal perspectives of the nutrient transport pathways into the lower St. Lawrence Estuary - Cynthia Evelyn Bluteau; Peter S. Galbraith; Daniel Bourgault; Vincent Villeneuve; Jean-Éric Tremblay (Ocean Science, 2021)
The St. Lawrence Estuary connects the Great Lakes with the Atlantic Ocean. The accepted view, based on summer conditions, is that the estuary's surface layer receives its nutrient supply from vertical mixing processes. This mixing is caused by the estuarine circulation and tides interacting with the topography at the head of the Laurentian Channel. During winter when ice forms, historical process-based studies have been limited in scope. Winter monitoring has been typically confined to vertical profiles of salinity and temperature as well as near-surface water samples collected from a helicopter for nutrient analysis. In 2018, however, the Canadian Coast Guard approved a science team to sample in tandem with its ice-breaking and ship escorting operations. This opportunistic sampling provided the first winter turbulence observations, which covered the largest spatial extent ever measured during any season within the St. Lawrence Estuary and the Gulf of St. Lawrence. The nitrate enrichment from tidal mixing resulted in an upward nitrate flux of about 30 nmol m−2 s−1, comparable to summer values obtained at the same tidal phase. Further downstream, deep nutrient-rich water from the gulf was mixed into the subsurface nutrient-poor layer at a rate more than an order of magnitude smaller than at the head. These fluxes were compared to the nutrient load of the upstream St. Lawrence River. Contrary to previous assumptions, fluvial nitrate inputs are the most significant source of nitrate in the estuary. Nitrate loads from vertical mixing processes would only exceed those from fluvial sources at the end of summer when fluvial inputs reach their annual minimum.

The riddle of eastern tropical Pacific Ocean oxygen levels: the role of the supply by intermediate-depth waters - Olaf Duteil; Ivy Frenger; Julia Getzlaff (Ocean Science, 2021)
Observed oxygen minimum zones (OMZs) in the tropical Pacific Ocean are located above intermediate-depth waters (IDWs), defined here as the 500–1500 m water layer. Typical climate models do not represent IDW properties well and are characterized by OMZs that are too deep-reaching. We analyze the role of the IDW in the misrepresentation of oxygen levels in a heterogeneous subset of ocean models characterized by a horizontal resolution ranging from 0.1 to 2.8∘. First, we show that forcing the extratropical boundaries (30∘ S and N) to observed oxygen values results in a significant increase in oxygen levels in the intermediate eastern tropical region. Second, we highlight the fact that the Equatorial Intermediate Current System (EICS) is a key feature connecting the western and eastern part of the basin. Typical climate models lack in representing crucial aspects of this supply at intermediate depth, as the EICS is basically absent in models characterized by a resolution lower than 0.25∘. These two aspects add up to a “cascade of biases” that hampers the correct representation of oxygen levels at intermediate depth in the eastern tropical Pacific Ocean and potentially future OMZ projections.

Decadal sea-level variability in the Australasian Mediterranean Sea - Patrick Wagner; Claus W. Böning (Ocean Science, 2021)
Strong regional sea-level trends, mainly related to basin-wide wind stress anomalies, have been observed in the western tropical Pacific over the last 3 decades. Analyses of regional sea level in the densely populated regions of the neighbouring Australasian Mediterranean Sea (AMS; also called tropical Asian seas) are hindered by its complex topography and respective studies are sparse. We used a series of global eddy-permitting ocean models, including a high-resolution configuration that resolves the AMS with 120∘ horizontal resolution, forced by a comprehensive atmospheric forcing product over 1958–2016 to characterize the patterns and magnitude of decadal sea-level variability in the AMS. The nature of this variability is elucidated further by sensitivity experiments with interannual variability restricted to either the momentum or buoyancy fluxes, building on an experiment employing a repeated-year forcing without interannual variability in all forcing components. Our results suggest that decadal fluctuations of the El Niño–Southern Oscillation (ENSO) account for over 80 % of the variability in all deep basins of the region, except for the central South China Sea (SCS). Changes related to the Pacific Decadal Oscillation (PDO) are most pronounced in the shallow Arafura and Timor seas and in the central SCS. On average, buoyancy fluxes account for less than 10 % of decadal SSH variability, but this ratio is highly variable over time and can reach values of up to 50 %. In particular, our results suggest that buoyancy flux forcing amplifies the dominant wind-stress-driven anomalies related to ENSO cycles. Intrinsic variability is mostly negligible except in the SCS, where it accounts for 25 % of the total decadal SSH variability.

Western boundary circulation and coastal sea-level variability in Northern Hemisphere oceans - Samuel Tiéfolo Diabaté; Didier Swingedouw; Joël Jean-Marie Hirschi; Aurélie Duchez; Philip J. Leadbitter; Ivan D. Haigh; Gerard D. McCarthy (Ocean Science, 2021)
The northwest basins of the Atlantic and Pacific oceans are regions of intense western boundary currents (WBCs): the Gulf Stream and the Kuroshio. The variability of these poleward currents and their extensions in the open ocean is of major importance to the climate system. It is largely dominated by in-phase meridional shifts downstream of the points at which they separate from the coast. Tide gauges on the adjacent coastlines have measured the inshore sea level for many decades and provide a unique window on the past of the oceanic circulation. The relationship between coastal sea level and the variability of the western boundary currents has been previously studied in each basin separately, but comparison between the two basins is missing. Here we show for each basin that the inshore sea level upstream of the separation points is in sustained agreement with the meridional shifts of the western boundary current extension over the period studied, i.e. the past 7 (5) decades in the Atlantic (Pacific). Decomposition of the coastal sea level into principal components allows us to discriminate this variability in the upstream sea level from other sources of variability such as the influence of large meanders in the Pacific. Our result extends previous findings limited to the altimetry era and suggests that prediction of inshore sea-level changes could be improved by the inclusion of meridional shifts of the western boundary current extensions as predictors. Long-duration tide gauges, such as Key West, Fernandina Beach or Hosojima, could be used as proxies for the past meridional shifts of the western boundary current extensions.

Sea surface salinity short-term variability in the tropics - Frederick M. Bingham; Susannah Brodnitz (Ocean Science, 2021)
Using data from the Global Tropical Moored Buoy Array, we study the validation process for satellite measurement of sea surface salinity (SSS). We compute short-term variability (STV) of SSS, variability on timescales of 2–17 d. It is a proxy for subfootprint variability over a 100 km footprint as seen by a satellite measuring SSS. We also compute representation error, which is meant to mimic the SSS satellite validation process where footprint averages are compared to pointwise in situ values. We present maps of these quantities over the tropical array. We also look at seasonality in the variability of SSS and find which months have maximum and minimum amounts. STV is driven at least partly by rainfall. Moorings exhibit larger STV during rainy periods than during non-rainy ones. The same computations are also done using output from a high-resolution global ocean model to see how it might be used to study the validation process. The model gives good estimates of STV, in line with the moorings, although tending to have smaller values.

Distribution of suspended particulate matter at the equatorial transect in the Atlantic Ocean - Vadim Sivkov; Ekaterina Bubnova (Ocean Science, 2021)
A suspended particulate matter distribution against a hydrographical background was studied at the oceanographic transect across the equatorial Atlantic in the year 2000. An area of abnormally high suspended matter volume concentrations was found above the Sierra Leone Rise in the entire water column (eastern part of the transect). The suggested explanation for the anomaly is based on the ballast hypothesis whereby solid particles are incorporated as ballast into suspended biogenic aggregates, leading to increased velocities of sinking. This occurs within the Northwest African upwelling area, where the plankton exposed to the Saharan dust abundance form a significant number of aggregates, which are later transported equatorward via the Canary Current. An intermediate nepheloid layer associated with the Deep Western Boundary Current was recorded from the South American Slope at depths of 3200–3700 to 4300 m above the Para Abyssal Plain. Antarctic Bottom Water enriched in suspended matter was found mostly in the troughs at 40–41∘ W. It was detached from the bottom, coinciding with the core of the flow due to the bottom rise “dam” located up-stream. The grain size of particles along the entire transect has a polymodal distribution with 2–4 and 8–13 µm modes. The registered rise in percentage in some parts of the transect of the 7–21 µm sized particles suggests the presence of the well-known coarse mode (20–60 µm) formed by aggregation of transparent exopolymer particles (mucus).

Flow separation, dipole formation, and water exchange through tidal straits - Ole Anders Nøst; Eli Børve (Ocean Science, 2021)
We investigate the formation and evolution of dipole vortices and their contribution to water exchange through idealized tidal straits. Self-propagating dipoles are important for transporting and exchanging water properties through straits and inlets in coastal regions. In order to obtain a robust dataset to evaluate flow separation, dipole formation and evolution, and the effect on water exchange, we conduct 164 numerical simulations, varying the width and length of the straits as well as the tidal forcing. We show that dipoles form and start propagating at the time of flow separation, and their vorticity originates in the velocity front formed by the separation. We find that the dipole propagation velocity is proportional to the tidal velocity amplitude and twice as large as the dipole velocity derived for a dipole consisting of two point vortices. We analyze the processes creating a net water exchange through the straits and derive a kinematic model dependent on dimensionless parameters representing strait length, dipole travel distance, and dipole size. The net tracer transport resulting from the kinematic model agrees closely with the numerical simulations and provides an understanding of the processes controlling net water exchange.

Interannual variability in contributions of the Equatorial Undercurrent (EUC) to Peruvian upwelling source water - Gandy Maria Rosales Quintana; Robert Marsh; Luis Alfredo Icochea Salas (Ocean Science, 2021)
Time-varying sources of upwelling waters off the coast of northern Peru are analyzed in a Lagrangian framework, tracking virtual particles backwards in time for 12 months. Particle trajectories are calculated with temperature, salinity and velocity fields from a hindcast spanning 1988–2007, obtained with an eddy-resolving (1/12∘) global configuration of the Nucleus for European Modelling of the Ocean (NEMO) ocean model. At 30 and 100 m, where coastal upwelling rates exceed 50 m month−1, particles are seeded at monthly intervals in proportion to the upwelling rate. Ensemble maps of particle concentration, age, depth, temperature, salinity and density reveal that a substantial but variable fraction of the particles upwelling off Peru arrives via the Equatorial Undercurrent (EUC). Particles follow the EUC core within the depth range 125–175 m, characterized by temperatures

Atmospherically forced sea-level variability in western Hudson Bay, Canada - Igor A. Dmitrenko; Denis L. Volkov; Tricia A. Stadnyk; Andrew Tefs; David G. Babb; Sergey A. Kirillov; Alex Crawford; Kevin Sydor; David G. Barber (Ocean Science, 2021)
In recent years, significant trends toward earlier breakup and later freeze-up of sea ice in Hudson Bay have led to a considerable increase in shipping activity through the Port of Churchill, which is located in western Hudson Bay and is the only deep-water ocean port in the province of Manitoba. Therefore, understanding sea-level variability at the port is an urgent issue crucial for safe navigation and coastal infrastructure. Using tidal gauge data from the port along with an atmospheric reanalysis and Churchill River discharge, we assess environmental factors impacting synoptic to seasonal variability of sea level at Churchill. An atmospheric vorticity index used to describe the wind forcing was found to correlate with sea level at Churchill. Statistical analyses show that, in contrast to earlier studies, local discharge from the Churchill River can only explain up to 5 % of the sea-level variability. The cyclonic wind forcing contributes from 22 % during the ice-covered winter–spring season to 30 % during the ice-free summer–fall season due to cyclone-induced storm surges generated along the coast. Multiple regression analysis revealed that wind forcing and local river discharge combined can explain up to 32 % of the sea-level variability at Churchill. Our analysis further revealed that the seasonal cycle of sea level at Churchill appears to be impacted by the seasonal cycle in atmospheric circulation rather than by the seasonal cycle in local discharge from the Churchill River, particularly post-construction of the Churchill River diversion in 1977. Sea level at Churchill shows positive anomalies for September–November compared to June–August. This seasonal difference was also revealed for the entire Hudson Bay coast using satellite-derived sea-level altimetry. This anomaly was associated with enhanced cyclonic atmospheric circulation during fall, reaching a maximum in November, which forced storm surges along the coast. Complete sea-ice cover during winter impedes momentum transfer from wind stress to the water column, reducing the impact of wind forcing on sea-level variability. Expanding our observations to the bay-wide scale, we confirmed the process of wind-driven sea-level variability with (i) tidal-gauge data from eastern Hudson Bay and (ii) satellite altimetry measurements. Ultimately, we find that cyclonic winds generate sea-level rise along the western and eastern coasts of Hudson Bay at the synoptic and seasonal timescales, suggesting an amplification of the bay-wide cyclonic geostrophic circulation in fall (October–November), when cyclonic vorticity is enhanced, and Hudson Bay is ice-free.

Role of air–sea fluxes and ocean surface density in the production of deep waters in the eastern subpolar gyre of the North Atlantic - Tillys Petit; M. Susan Lozier; Simon A. Josey; Stuart A. Cunningham (Ocean Science, 2021)
Wintertime convection in the North Atlantic Ocean is a key component of the global climate as it produces dense waters at high latitudes that flow equatorward as part of the Atlantic Meridional Overturning Circulation (AMOC). Recent work has highlighted the dominant role of the Irminger and Iceland basins in the production of North Atlantic Deep Water. Dense water formation in these basins is mainly explained by buoyancy forcing that transforms surface waters to the deep waters of the AMOC lower limb. Air–sea fluxes and the ocean surface density field are both key determinants of the buoyancy-driven transformation. We analyze these contributions to the transformation in order to better understand the connection between atmospheric forcing and the densification of surface water. More precisely, we study the impact of air–sea fluxes and the ocean surface density field on the transformation of subpolar mode water (SPMW) in the Iceland Basin, a water mass that “pre-conditions” dense water formation downstream. Analyses using 40 years of observations (1980–2019) reveal that the variance in SPMW transformation is mainly influenced by the variance in density at the ocean surface. This surface density is set by a combination of advection, wind-driven upwelling and surface fluxes. Our study shows that the latter explains ∼ 30 % of the variance in outcrop area as expressed by the surface area between the outcropped SPMW isopycnals. The key role of the surface density in SPMW transformation partly explains the unusually large SPMW transformation in winter 2014–2015 over the Iceland Basin.

Drifting dynamics of the bluebottle ( Physalia physalis ) - Daniel Lee; Amandine Schaeffer; Sjoerd Groeskamp (Ocean Science, 2021)
Physalia physalis, also called the bluebottle in Australia, is a colonial animal resembling a jellyfish that is well known to beachgoers for the painful stings delivered by its tentacles. Despite being a common occurrence, the origin of the bluebottle before reaching the coastline is not well understood, and neither is the way it drifts at the surface of the ocean. Previous studies used numerical models in combination with simple assumptions to calculate the drift of this species, excluding complex drifting dynamics. In this study, we provide a new parameterization for Lagrangian modelling of the bluebottle by considering the similarities between the bluebottle and a sailboat. This allows us to compute the hydrodynamic and aerodynamic forces acting on the bluebottle and use an equilibrium condition to create a generalized model for calculating the drifting speed and course of the bluebottle under any wind and ocean current conditions. The generalized model shows that the velocity of the bluebottle is a linear combination of the ocean current velocity and the wind velocity scaled by a coefficient (“shape parameter”) and multiplied by a rotation matrix. Adding assumptions to this generalized model allows us to retrieve models used in previous literature. We discuss the sensitivity of the model to different parameters (shape, angle of attack and sail camber) and explore different cases of wind and current conditions to provide new insights into the drifting dynamics of the bluebottle.

A dynamically based method for estimating the Atlantic meridional overturning circulation at 26° N from satellite altimetry - Alejandra Sanchez-Franks; Eleanor Frajka-Williams; Ben I. Moat; David A. Smeed (Ocean Science, 2021)
The large-scale system of ocean currents that transport warm waters in the upper 1000 m northward and return deeper cooler waters southward is known as the Atlantic meridional overturning circulation (AMOC). Variations in the AMOC have significant repercussions for the climate system; hence, there is a need for long-term monitoring of AMOC fluctuations. Currently the longest record of continuous directly measured AMOC changes is from the RAPID-MOCHA-WBTS programme, initiated in 2004. The RAPID programme and other mooring programmes have revolutionised our understanding of large-scale circulation; however, by design they are constrained to measurements at a single latitude and cannot tell us anything pre-2004. Nearly global coverage of surface ocean data from satellite altimetry has been available since the launch of the TOPEX/Poseidon satellite in 1992 and has been shown to provide reliable estimates of surface ocean transports on interannual timescales including previous studies that have investigated empirical correlations between sea surface height variability and the overturning circulation. Here we show a direct calculation of ocean circulation from satellite altimetry of the upper mid-ocean transport (UMO), the Gulf Stream transport through the Florida Straits (GS), and the AMOC using a dynamically based method that combines geostrophy with a time mean of the vertical structure of the flow from the 26∘ N RAPID moorings. The satellite-based transport captures 56 %, 49 %, and 69 % of the UMO, GS, and AMOC transport variability, respectively, from the 26∘ N RAPID array on interannual (18-month) timescales. Further investigation into the vertical structure of the horizontal transport shows that the first baroclinic mode accounts for 83 % of the interior geostrophic variability, and the combined barotropic and first baroclinic mode representation of dynamic height accounts for 98 % of the variability. Finally, the methods developed here are used to reconstruct the UMO and the AMOC for the time period pre-dating RAPID, 1993 to 2003. The effective implementation of satellite-based method for monitoring the AMOC at 26∘ N lays down the starting point for monitoring large-scale circulation at all latitudes.

Causes of uncertainties in the representation of the Arabian Sea oxygen minimum zone in CMIP5 models - Henrike Schmidt; Julia Getzlaff; Ulrike Löptien; Andreas Oschlies (Ocean Science, 2021)
Open-ocean oxygen minimum zones (OMZs) occur in regions with high biological productivity and weak ventilation. They restrict marine habitats and alter biogeochemical cycles. Global models generally show a large model–data misfit with regard to oxygen. Reliable statements about the future development of OMZs and the quantification of their interaction with climate change are currently not possible. One of the most intense OMZs worldwide is located in the Arabian Sea (AS). We give an overview of the main model deficiencies with a detailed comparison of the historical state of 10 climate models from the 5th Coupled Model Intercomparison Project (CMIP5) that present our present-day understanding of physical and biogeochemical processes. Most of the models show a general underestimation of the OMZ volume in the AS compared to observations that is caused by an overly shallow layer of oxygen-poor water in the models. The deviation of oxygen values in the deep AS is the result of oxygen levels that are too high simulated in the Southern Ocean formation regions of Indian Ocean Deep Water in the models compared to observations and uncertainties in the deepwater mass transport from the Southern Ocean northward into the AS. Differences in simulated water mass properties and ventilation rates of Red Sea Water and Persian Gulf Water cause different mixing in the AS and thus influence the intensity of the OMZ. These differences in ventilation rates also point towards variations in the parameterizations of the overflow from the marginal seas among the models. The results of this study are intended to foster future model improvements regarding the OMZ in the AS.

In situ observations of turbulent ship wakes and their spatiotemporal extent - Amanda T. Nylund; Lars Arneborg; Anders Tengberg; Ulf Mallast; Ida-Maja Hassellöv (Ocean Science, 2021)
In areas of intensive ship traffic, ships pass every 10 min. Considering the amount of ship traffic and the predicted increase in global maritime trade, there is a need to consider all types of impacts shipping has on the marine environment. While the awareness about, and efforts to reduce, chemical pollution from ships is increasing, less is known about physical disturbances, and ship-induced turbulence has so far been completely neglected. To address the potential importance of ship-induced turbulence on, e.g., gas exchange, dispersion of pollutants, and biogeochemical processes, a characterisation of the temporal and spatial scales of the turbulent wake is needed. Currently, field measurements of turbulent wakes of real-size ships are lacking. This study addresses that gap by using two different methodological approaches: in situ and ex situ observations. For the in situ observations, a bottom-mounted acoustic Doppler current profiler (ADCP) was placed at 32 m depth below the shipping lane outside Gothenburg harbour. Both the acoustic backscatter from the air bubbles in the wake and the dissipation rate of turbulent kinetic energy were used to quantify the turbulent wake depth, intensity, and temporal longevity for 38 ship passages of differently sized ships. The results from the ADCP measurements show median wake depths of 13 m and several occasions of wakes reaching depths > 18 m, which is in the same depth range as the seasonal thermocline in the Baltic Sea. The temporal longevity of the observable part of the wakes had a median of around 10 min and several passages of > 20 min. In the ex situ approach, sea surface temperature was used as a proxy for the water mass affected by the turbulent wake (thermal wake), as lowered temperature in the ship wake indicates vertical mixing in a thermally stratified water column. Satellite images of the thermal infrared sensor (TIRS) onboard Landsat-8 were used to measure thermal wake width and length, in the highly frequented and thus major shipping lane north of Bornholm, Baltic Sea. Automatic information system (AIS) records from both the investigated areas were used to identify the ships inducing the wakes. The satellite analysis showed a median thermal wake length of 13.7 km (n=144), and the longest wake extended over 60 km, which would correspond to a temporal longevity of 1 h 42 min (for a ship speed of 20 kn). The median thermal wake width was 157.5 m. The measurements of the spatial and temporal scales are in line with previous studies, but the maximum turbulent wake depth (30.5 m) is deeper than previously reported. The results from this study, combined with the knowledge of regional high traffic densities, show that ship-induced turbulence occurs at temporal and spatial scales large enough to imply that this process should be considered when estimating environmental impacts from shipping in areas with intense ship traffic.

Filtering method based on cluster analysis to avoid salinity drifts and recover Argo data in less time - Emmanuel Romero; Leonardo Tenorio-Fernandez; Iliana Castro; Marco Castro (Ocean Science, 2021)
Currently there is a huge amount of freely available hydrographic data, and it is increasingly important to have easy access to it and to be provided with as much information as possible. Argo is a global collection of around 4000 active autonomous hydrographic profilers. Argo data go through two quality processes, real time and delayed mode. This work shows a methodology to filter profiles within a given polygon using the odd–even algorithm; this allows analysis of a study area, regardless of size, shape or location. The aim is to offer two filtering methods and to discard only the real-time quality control data that present salinity drifts. This takes advantage of the largest possible amount of valid data within a given polygon. In the study area selected as an example, it was possible to recover around 80 % in the case of the first filter that uses cluster analysis and 30 % in the case of the second, which discards profilers with salinity drifts, of the total real-time quality control data that are usually discarded by the users due to problems such as salinity drifts. This allows users to use any of the filters or a combination of both to have a greater amount of data within the study area of their interest in a matter of minutes, rather than waiting for the delayed-mode quality control that takes up to 12 months to be completed. This methodology has been tested for its replicability in five selected areas around the world and has obtained good results.

Effects of strongly eddying oceans on multidecadal climate variability in the Community Earth System Model - André Jüling; Anna von der Heydt; Henk A. Dijkstra (Ocean Science, 2021)
Climate variability on multidecadal timescales appears to be organized in pronounced patterns with clear expressions in sea surface temperature, such as the Atlantic Multidecadal Variability and the Pacific Decadal Oscillation. These patterns are now well studied both in observations and global climate models and are important in the attribution of climate change. Results from CMIP5 models have indicated large biases in these patterns with consequences for ocean heat storage variability and the global mean surface temperature. In this paper, we use two multi-century Community Earth System Model simulations at coarse (1∘) and fine (0.1∘) ocean model horizontal grid spacing to study the effects of the representation of mesoscale ocean flows on major patterns of multidecadal variability. We find that resolving mesoscale ocean flows both improves the characteristics of the modes of variability with respect to observations and increases the amplitude of the heat content variability in the individual ocean basins. In the strongly eddying model, multidecadal variability increases compared to sub-decadal variability. This shift of spectral power is seen in sea surface temperature indices, basin-scale surface heat fluxes, and the global mean surface temperature. This implies that the current CMIP6 model generation, which predominantly does not resolve the ocean mesoscale, may systematically underestimate multidecadal variability.

Lagrangian eddy tracking reveals the Eratosthenes anticyclonic attractor in the eastern Levantine Basin - Alexandre Barboni; Ayah Lazar; Alexandre Stegner; Evangelos Moschos (Ocean Science, 2021)
Statistics of anticyclonic eddy activity and eddy trajectories in the Levantine Basin over the 2000–2018 period are analyzed using the DYNED-Atlas database, which links automated mesoscale eddy detection by the Angular Momentum Eddy Detection and Tracking Algorithm (AMEDA) algorithm to in situ oceanographic observations. This easternmost region of the Mediterranean Sea, delimited by the Levantine coast and Cyprus, has a complex eddying activity, which has not yet been fully characterized. In this paper, we use Lagrangian tracking to investigate the eddy fluxes and interactions between different subregions in this area. The anticyclonic structure above the Eratosthenes Seamount is identified as hosting an anticyclone attractor, constituted by a succession of long-lived anticyclones. It has a larger radius and is more persistent (staying in the same position for up to 4 years with successive merging events) than other eddies in this region. Quantification of anticyclone flux shows that anticyclones that drift towards the Eratosthenes Seamount are mainly formed along the Israeli coast or in a neighboring area west of the seamount. The southeastern Levantine area is isolated, with no anticyclone transfers to or from the western part of the basin, defining the effective attraction basin for the Eratosthenes anticyclone attractor. Co-localized in situ profiles inside eddies provide quantitative information on their subsurface physical anomaly signature, whose intensity can vary greatly with respect to the dynamical surface signature intensity. Despite interannual variability, the so-called Eratosthenes anticyclone attractor stores a larger amount of heat and salt than neighboring anticyclones, in a deeper subsurface anomaly that usually extends down to 500 m. This suggests that this attractor could concentrate heat and salt from this subbasin, which will impact the properties of intermediate water masses created there.

Control of oceanic circulation on sediment distribution in the southwestern Atlantic margin (23 to 55° S) - Michel Michaelovitch de Mahiques; Roberto Violante; Paula Franco-Fraguas; Leticia Burone; Cesar Barbedo Rocha; Leonardo Ortega; Rosangela Felicio dos Santos; Bianca Sung Mi Kim; Rubens Cesar Lopes Figueira; Marcia Caruso Bícego (Ocean Science, 2021)
In this study, we interpret the role played by ocean circulation in sediment distribution on the southwestern Atlantic margin using radiogenic Nd and Pb isotopes. The latitudinal trends for Pb and Nd isotopes reflect the different current systems acting on the margin. The utilization of the sediment fingerprinting method allowed us to associate the isotopic signatures with the main oceanographic features in the area. We recognized differences between Nd and Pb sources to the Argentinean shelf (carried by the flow of Subantarctic Shelf Water) and slopes (transported by deeper flows). Sediments from Antarctica extend up to the Uruguayan margin, carried by the Upper and Lower Circumpolar Deep Water. Our data confirm that, for shelf and intermediate areas (the upper 1200 m), the transfer of sediments from the Argentinean margin to the north of 35∘ S is limited by the Subtropical Shelf Front and the basin-wide recirculated Antarctic Intermediate Water. On the southern Brazilian inner and middle shelf, it is possible to recognize the northward influence of the Río de la Plata sediments carried by the Plata Plume Water. Another flow responsible for sediment transport and deposition on the outer shelf and slope is the southward flow of the Brazil Current. Finally, we propose that the Brazil–Malvinas Confluence and the Santos Bifurcation act as boundaries of geochemical provinces in the area. A conceptual model of sediment sources and transport is provided for the southwestern Atlantic margin.

Regional imprints of changes in the Atlantic Meridional Overturning Circulation in the eddy-rich ocean model VIKING20X - Arne Biastoch; Franziska U. Schwarzkopf; Klaus Getzlaff; Siren Rühs; Torge Martin; Markus Scheinert; Tobias Schulzki; Patricia Handmann; Rebecca Hummels; Claus W. Böning (Ocean Science, 2021)
A hierarchy of global 1/4∘ (ORCA025) and Atlantic Ocean 1/20∘ nested (VIKING20X) ocean–sea-ice models is described. It is shown that the eddy-rich configurations performed in hindcasts of the past 50–60 years under CORE and JRA55-do atmospheric forcings realistically simulate the large-scale horizontal circulation, the distribution of the mesoscale, overflow and convective processes, and the representation of regional current systems in the North and South Atlantic. The representation of the Atlantic Meridional Overturning Circulation (AMOC), and in particular the long-term temporal evolution, strongly depends on numerical choices for the application of freshwater fluxes. The interannual variability of the AMOC instead is highly correlated among the model experiments and also with observations, including the 2010 minimum observed by RAPID at 26.5∘ N. This points to a dominant role of the wind forcing. The ability of the model to represent regional observations in western boundary current (WBC) systems at 53∘ N, 26.5∘ N and 11∘ S is explored. The question is investigated of whether WBC systems are able to represent the AMOC, and in particular whether these WBC systems exhibit similar temporal evolution to that of the zonally integrated AMOC. Apart from the basin-scale measurements at 26.5∘ N, it is shown that in particular the outflow of North Atlantic Deepwater at 53∘ N is a good indicator of the subpolar AMOC trend during the recent decades, once provided in density coordinates. The good reproduction of observed AMOC and WBC trends in the most reasonable simulations indicate that the eddy-rich VIKING20X is capable of representing realistic forcing-related and ocean-intrinsic trends.

Evaluating high-frequency radar data assimilation impact in coastal ocean operational modelling - Jaime Hernandez-Lasheras; Baptiste Mourre; Alejandro Orfila; Alex Santana; Emma Reyes; Joaquín Tintoré (Ocean Science, 2021)
The impact of the assimilation of HFR (high-frequency radar) observations in a high-resolution regional model is evaluated, focusing on the improvement of the mesoscale dynamics. The study area is the Ibiza Channel, located in the western Mediterranean Sea. The resulting fields are tested against trajectories from 13 drifters. Six different assimilation experiments are compared to a control run (no assimilation). The experiments consist of assimilating (i) sea surface temperature, sea level anomaly, and Argo profiles (generic observation dataset); the generic observation dataset plus (ii) HFR total velocities and (iii) HFR radial velocities. Moreover, for each dataset, two different initialization methods are assessed: (a) restarting directly from the analysis after the assimilation or (b) using an intermediate initialization step applying a strong nudging towards the analysis fields. The experiments assimilating generic observations plus HFR total velocities with the direct restart provide the best results, reducing by 53 % the average separation distance between drifters and virtual particles after the first 48 h of simulation in comparison to the control run. When using the nudging initialization step, the best results are found when assimilating HFR radial velocities with a reduction of the mean separation distance by around 48 %. Results show that the integration of HFR observations in the data assimilation system enhances the prediction of surface currents inside the area covered by both antennas, while not degrading the correction achieved thanks to the assimilation of generic data sources beyond it. The assimilation of radial observations benefits from the smoothing effect associated with the application of the intermediate nudging step.

Can assimilation of satellite observations improve subsurface biological properties in a numerical model? A case study for the Gulf of Mexico - Bin Wang; Katja Fennel; Liuqian Yu (Ocean Science, 2021)
Given current threats to ocean ecosystem health, there is a growing demand for accurate biogeochemical hindcasts, nowcasts, and predictions. Provision of such products requires data assimilation, i.e., a comprehensive strategy for incorporating observations into biogeochemical models, but current data streams of biogeochemical observations are generally considered insufficient for the operational provision of such products. This study investigates to what degree the assimilation of satellite observations in combination with a priori model calibration by sparse BGC-Argo profiles can improve subsurface biogeochemical properties. The multivariate deterministic ensemble Kalman filter (DEnKF) has been implemented to assimilate physical and biological observations into a three-dimensional coupled physical–biogeochemical model, the biogeochemical component of which has been calibrated by BGC-Argo float data for the Gulf of Mexico. Specifically, observations of sea surface height, sea surface temperature, and surface chlorophyll were assimilated, and profiles of both physical and biological variables were updated based on the surface information. We assessed whether this leads to improved subsurface distributions, especially of biological properties, using observations from five BGC-Argo floats that were not assimilated. An alternative light parameterization that was tuned a priori using BGC-Argo observations was also applied to test the sensitivity of data assimilation impact on subsurface biological properties. Results show that assimilation of the satellite data improves model representation of major circulation features, which translate into improved three-dimensional distributions of temperature and salinity. The multivariate assimilation also improves the agreement of subsurface nitrate through its tight correlation with temperature, but the improvements in subsurface chlorophyll were modest initially due to suboptimal choices of the model's optical module. Repeating the assimilation run by using the alternative light parameterization greatly improved the subsurface distribution of chlorophyll. Therefore, even sparse BGC-Argo observations can provide substantial benefits for biogeochemical prediction by enabling a priori model tuning. Given that, so far, the abundance of BGC-Argo profiles in the Gulf of Mexico and elsewhere has been insufficient for sequential assimilation, updating 3D biological properties in a model that has been well calibrated is an intermediate step toward full assimilation of the new data types.

Simulated zonal current characteristics in the southeastern tropical Indian Ocean (SETIO) - Nining Sari Ningsih; Sholihati Lathifa Sakina; Raden Dwi Susanto; Farrah Hanifah (Ocean Science, 2021)
Detailed ocean currents in the southeastern tropical Indian Ocean adjacent to southern Sumatran and Javan coasts have not been fully explained because of limited observations. In this study, zonal current characteristics in the region have been studied using simulation results of a 1/8∘ global hybrid coordinate ocean model from 1950 to 2013. The simulated zonal currents across three meridional sections were then investigated using an empirical orthogonal function (EOF), where the first three modes account for 75 %–98 % of the total variance. The first temporal mode of EOF is then investigated using ensemble empirical mode decomposition (EEMD) to distinguish the signals. This study has revealed distinctive features of currents in the South Java Current (SJC) region, the Indonesian Throughflow (ITF)–South Equatorial Current (SEC) region, and the transition zone between these regions. The vertical structures of zonal currents in southern Java and offshore Sumatra are characterized by a one-layer flow. Conversely, a two-layer flow is observed in the nearshore and transition regions of Sumatra. Current variation in the SJC region has peak energies that are sequentially dominated by semiannual, intraseasonal, and annual timescales. Meanwhile, the transition zone is characterized by semiannual and intraseasonal periods with pronounced interannual variations. In contrast, interannual variability associated with El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) modulates the prominent intraseasonal variability of current in the ITF–SEC region. ENSO has the strongest influence at the outflow ITF, while the IOD's strongest influence is in southwestern Sumatra, with the ENSO (IOD) leading the current by 4 months (1 month). Moreover, the contributions (largest to smallest) of each EEMD mode at the nearshore of Java and offshore Sumatra are intraseasonal, semiannual, annual, interannual, and long-term fluctuations. The contribution of long-term variation (19.2 %) in the far offshore eastern Indian Ocean is larger than the interannual (16.3 %) and annual (14.7 %) variations. Future studies should be conducted to investigate this long-term variation.