Search

Accurate prediction of ocean surface currents is important for marine safety, ship routing, tracking of pollutants and in coupled forecasting. Presently, velocity observations are not routinely assimilated in global ocean forecasting systems, largely due to the sparsity of the observation network. Several satellite missions are now being proposed with the capability to measure Total Surface Current Velocities (TSCV). If successful, these would substantially increase the coverage of ocean current observations and could improve accuracy of ocean current forecasts through data assimilation. In this paper, Observing System Simulation Experiments (OSSEs) are used to assess the impact of assimilating TSCV in the Met Office's global ocean forecasting system. Synthetic observations are generated from a high-resolution model run for all standard observation types (sea surface temperature, profiles of temperature and salinity, sea level anomaly and sea ice concentration) as well as TSCV observations from a Sea surface KInematics Multiscale monitoring (SKIM) like satellite. The assimilation of SKIM like TSCV observations is tested over an 11 month period. Preliminary experiments assimilating idealised single TSCV observations demonstrate that ageostrophic velocity corrections are not well retained in the model. We propose a method for improving ageostrophic currents through TSCV assimilation by initialising Near Inertial Oscillations with a rotated incremental analysis update (IAU) scheme. The OSSEs show that TSCV assimilation has the potential to significantly improve the prediction of velocities, particularly in the Western Boundary Currents, Antarctic Circumpolar Current and in the near surface equatorial currents. For global surface velocity the analysis root-mean-squareerrors (RMSEs) are reduced by 23% and there is a 4-day gain in forecast RMSE. There are some degradations to the subsurface in the tropics, generally in regions with complex vertical salinity structures. However, outside of the tropics, improvements are seen to velocities throughout the water column. Globally there are also improvements to temperature and sea surface height when TSCV are assimilated. The TSCV assimilation largely corrects the geostrophic ocean currents, but results using the rotated IAU method show that the energy at inertial frequencies can be improved with this method. Overall, the experiments demonstrate significant potential benefit of assimilating TSCV observations in a global ocean forecasting system. [ABSTRACT FROM AUTHOR]

This study presents a generalised characterisation of the surface circulation over the northern shelf of the Gulf of Cadiz, based on 4 years of high-resolution satellite altimetry data from Sentinel-3A and wind model data. The altimetry-based surface zonal currents, adjusted for bottom-drag and wind effects, are compared with a generic CMEMS product and validated against in-situ ADCP measurements. The proposed altimetry product demonstrates superior performance than the CMEMS product, accurately reflecting surface circulation direction compared to in-situ measurements (r = 0.77, RMSE = 0.10 m/s, bias = 0.01 m/s). The use of the bottom-drag and wind-corrected/uncorrected altimetry product for spatiotemporal analysis of the shelf circulation revealed the distinct contributions of wind-driven and geostrophic components in different basin sectors. The results show that over the western basin, positive (eastward) surface currents were predominantly driven by westerly winds, while only occasionally, westward flows coincided with easterly winds, suggesting a higher control of the geostrophic component over the westward flows. In contrast, over the eastern basin, both eastward and westward flows were found to be primarily driven by favourable winds. Additionally, the analysis of Absolute Dynamic Topography (ADT) values along the whole basin showed the presence of ADT gradients both along-shore and cross-shore over the shelf, contributing to geostrophic flows. Finally, the seasonal analysis showed that eastward circulation tends to dominate during the spring and summer months, related to the upwelling season in the Gulf of Cadiz and associated westerly winds. Westward flows prevail during the winter months, related to easterly winds and the rebalancing of the along-shore sea level gradient during relaxed upwelling conditions. The findings demonstrate a significant improvement in the use of satellite altimetry data to study complex oceanographic dynamics in coastal areas, where both spatial and temporal variability are high. Moreover, the similarity of our results to those obtained from in-situ systems supports the use of altimetry data and publicly available wind models to support oceanographic studies in remote or resource-limited areas. [ABSTRACT FROM AUTHOR]

The ocean regulates the Earth's climate by transporting heat from the equator to the poles. Here, we use satellite‐based sea surface observations of air‐sea heat fluxes and eddy detection to investigate the mesoscale heat transport. "Mesoscale" refers to both the Eulerian perspective as the spatio‐temporal scales of ∼100 km and ∼1 month, as well as the Lagrangian aspect as isolated vortices identified from the dynamic topography. Paradoxically, there are a considerable number of mesoscale eddies inconsistent between their surface thermal and dynamic signals, that is, cold‐core anticyclones and warm‐core cyclones are globally prevalent. On account of such inconsistency, we show that the mesoscale meridional heat transport carried by geostrophic components is 10 times larger than (and opposite in direction to) that of the wind‐driven Ekman components. An offset between SSH‐SST coherent and incoherent eddies in the Ekman heat transport is apparent, whereas the geostrophic heat transport is contained within coherent eddies. Plain Language Summary: Rotating water masses in the ocean, termed as mesoscale eddies, are dynamically important features because they can transport nutrients and displace heat in the Earth's climate system. In the same manner as the atmospheric weather systems that are associated with high or low pressures, mesoscale eddies are accompanied by sea level anomalies and therefore can be observed by satellite images. Conventionally, anticyclonic eddies detected as clockwise vortices in the Northern Hemisphere display warm temperature anomalies at the center. However, recent studies have found that the sea surface temperature signal of an eddy might be uncertain and even contradictory to their dynamical circulation patterns. In this study, we combine sea surface observations of temperature and velocity fields to determine the patterns of global ocean heat transport associated with mesoscale eddies. Key Points: The geostrophic component of meridional heat transport (MHT) at the mesoscale is 10 times larger than the Ekman componentSSH‐SST coherent eddies dominate the spatial patterns of MHT at the mesoscaleThough incomparable in magnitude with the large scale MHT, mesoscale eddies can still transport 30 TW of heat near sea surface [ABSTRACT FROM AUTHOR]

Realistic representation of monthly sea level anomalies in coastal regions has been a challenge for global ocean reanalyses. This is especially the case in coastal regions where sea levels are influenced by western boundary currents such as near the U.S. Atlantic Coast and the Gulf of Mexico. For these regions, most ocean reanalyses compare poorly to observations. Problems in reanalyses include errors in data assimilation and horizontal resolutions that are too coarse to simulate energetic currents like the Gulf Stream and Loop Current System. However, model capabilities are advancing with improved data assimilation and higher resolution. Here, we show that some current-generation ocean reanalyses produce monthly sea level anomalies with improved skill when compared to satellite altimetry observations of sea surface heights. Using tide gauge observations for coastal verification, we find the highest skill associated with the GLORYS12 and HYCOM ocean reanalyses. Both systems assimilate altimetry observations and have eddy-resolving horizontal resolutions (1/12°). We found less skill in three other ocean reanalyses (ACCESS-S2, ORAS5, and ORAP6) with coarser, though still eddy-permitting, resolutions (1/4°). The operational reanalysis from ECMWF (ORAS5) and their pilot reanalysis (ORAP6) provide an interesting comparison because the latter assimilates altimetry globally and with more weight, as well as assimilating ocean observations over continental shelves. We find these attributes associated with improved skill near many tide gauges. We also assessed an older reanalysis (CFSR), which has the lowest skill likely due to its lower resolution (1/2°) and lack of altimetry assimilation. ACCESS-S2 likewise does not assimilate altimetry, although its skill is much better than CFSR and only somewhat lower than ORAS5. Since coastal flooding is influenced by sea level anomalies, the recent development of skilful ocean reanalyses on monthly timescales may be useful for better understanding the physical processes associated with flood risks. [ABSTRACT FROM AUTHOR]

Global ocean surface currents estimated from satellite derived data based on a regular global grid are affected by the grid’s shape and placement. Due to different neighbourhood relationships, the rectangular lat/lon grids lose accuracy when interpolating and fitting elevation data. Hexagonal grids have shown to be advantageous due to their isotropic, uniform neighbourhood. Considering these merits, this paper aims to estimate global ocean surface current using a global isotropic hexagonal grid from satellite remote sensing data. First, gridded satellite altimeter data and wind data with different resolutions are interpolated into the centre of the global isotropic hexagonal grid. Then, geostrophic and Ekman currents components are estimated according to the Lagerlof Ocean currents theory. Finally, the inversion results are verified. By analyzing the results, we conclude that the ocean surface currents estimated based on the global isotropic hexagonal grid have considerable accuracy, with improvement over rectangular lat/lon grids.

We present a study on the ocean surface currents reconstruction by merging Level-4 (L4, gap-free) altimeter-derived geostrophic currents and satellite sea surface temperature. Building upon past studies on the multi-variate reconstruction of geostrophic currents from satellite observations, we regionalized and optimized an algorithm to improve the altimeter-derived surface circulation estimates in the North Atlantic Ocean. A ten-year-long time series (2010–2019) is presented and validated by means of in situ observations. The newly optimized algorithm allowed us to improve the currents estimate along the main axis of the Gulf Stream and in correspondence of well-known upwelling areas in the North Eastern Atlantic, with percentage improvements of around 15% compared to standard operational altimetry products. [ABSTRACT FROM AUTHOR]

Based on the relation between volume transport and sea level difference across the Korea/Tsushima Strait (KTS), in this study, we estimated the volume transport of the Tsushima Warm Current over the past four decades and examined its long-term trend using an ensemble empirical model decomposition (EEMD) method. The estimated geostrophic transports showed very good agreement with those measured by acoustic Doppler current profiler and those estimated from satellite altimetry-based sea level data. This corroborates the reliability of the long-term estimation. EEMD analysis revealed an acceleration in volume transport in the KTS, with a persistent increase in the trend rate occurring after the mid-1980s. Further analyses indicate that the trend shift in the KTS closely coincided with an upper-ocean warming trend in the East/Japan Sea (EJS) after the mid-1980s. Analyzing large-scale wind patterns suggests that long-term trends of wind stress curls over the North Pacific are likely a key driver of the trend shifts in the inflow through the KTS and the upper-ocean warming in the EJS over the past four decades. [ABSTRACT FROM AUTHOR]

In this study, we investigate the record-breaking intensification and abrupt weakening of the East Korea Warm Current (EKWC) in the summer of 2021. We analyzed the ocean data assimilation products resolving this event to examine the association between the abrupt changes in the EKWC and various oceanic/atmospheric factors. The results indicate that during the summer of 2021, the EKWC extended northward beyond its climatology, reaching up to 40°N with the maximum speed of 1.16 m s-1 on August 1. In mid-August, the EKWC underwent a rapid weakening, returning to its climatological level. We could attribute the temporal variability in the anomalous EKWC in 2021 to the distinct temporal variability in the dynamic height anomalies between coastal and offshore regions. The offshore variability in the dynamic height anomaly, which is related to warm eddy variability, led to an anomalously increased EKWC velocity (up to 0.59 m s-1) during the EKWC peak velocity period in 2021. However, anomalous coastal downwelling induced by a weak northerly wind anomaly decelerated the EKWC by -0.06 m s-1 in the same period. In mid-August, a typhoon-related northerly wind induced a sudden rise in the coastal dynamic height anomaly, resulting in a rapid weakening of the EKWC. Our findings suggest that changes in geostrophic current related to warm eddies and typhoons have substantially contributed to the temporal variability in the EKWC, improving our understanding of the temporal variability in the western boundary currents. [ABSTRACT FROM AUTHOR]

The offshore transport of Greenland coastal waters influenced by freshwater input from ice sheet melting during summer plays an important role in ocean circulation and biological processes in the Labrador Sea. Many previous studies over the last decade have investigated shelfbreak transport processes in the region, primarily using ocean model simulations. Here, we use 27 years of surface geostrophic velocity observations from satellite altimetry, modified to include Ekman dynamics based on atmospheric reanalysis, and virtual particle releases to investigate seasonal and interannual variability in transport of coastal water in the Labrador Sea. Two sets of tracking experiments were pursued, one using geostrophic velocities only, and another using total velocities including the wind effect. Our analysis revealed substantial seasonal variability, even when only geostrophic velocities were considered. Water from coastal southwest Greenland is generally transported northward into Baffin Bay, although westward transport off the west Greenland shelf increases in fall and winter due to winds. Westward offshore transport is increased for water from southeast Greenland so that, in some years, water originating near the east Greenland coast during summer can be transported into the central Labrador Sea and the convection region. When wind forcing is considered, long-term trends suggest decreasing transport of Greenland coastal water during the melting season toward Baffin Bay, and increasing transport into the interior of the Labrador Sea for water originating from southeast Greenland during summer, where it could potentially influence water column stability. Future studies using higher-resolution velocity observations are needed to capture the role of submesoscale variability in transport pathways in the Labrador Sea. [ABSTRACT FROM AUTHOR]

The East Kamchatka and East Sakhalin Currents (EKC and ESC) are the western boundary currents of the subarctic North Pacific and Okhotsk Sea. Variability in the EKC and ESC velocities could exert a substantial effect on ecosystems and fish stocks in the southwestern Bering Sea and Okhotsk Sea. Using satellite-derived data (sea surface heights, geostrophic current velocities, and sea surface temperatures, 2002–2020), we demonstrate that changes in zonal wind generate sea level variations on the shelf in the southwestern Bering Sea over a period of 18–29 days and with an amplitude of 5–20 cm. The ebb/flood events on the shelf lead to changes in the velocity, direction, and position of the EKC. The sea level anomalies propagate along the western Kamchatka, northern Kuril Islands and the northern and western Okhotsk Sea and result in the variability of geostrophic current velocities in the ESC zone. The strengthening (weakening) of ESC leads to an increase (a decrease) in SST in the southern part of the Okhotsk Sea by 1–3 °C. In the northwestern Okhotsk Sea, in addition to wind-induced variability, there are temporary changes in the geostrophic currents with a period of 14 days caused by fortnightly tides. [ABSTRACT FROM AUTHOR]

Mesoscale eddies play a crucial role in the transport of mass, heat, salt and nutrients, exerting significant influence on ocean circulation patterns, biogeochemical processes and the global climate system. Based on Lagrangian-Averaged Vorticity Deviation (LAVD) method, this study applies 27 years (1993–2019) of geostrophic current velocity data to detect Rotationally Coherent Lagrangian Vortices (RCLVs) in the Northwest Pacific (NWP; 10°N–30°N, 115°E–155°E), with the spatiotemporal characteristics of Eulerian Sea Surface Height Eddies (SSH eddies) and RCLVs being compared. A higher number of SSH eddies and RCLVs can be observed in spring and winter, and their inter-annual variations are similar. SSH eddies show higher generation number and larger radius in the Subtropical Countercurrent region, while RCLVs occur more favorably in the ocean basin. The propagation speed distributions of both eddy types are nearly identical and decrease with increasing latitude. Due to the material coherent transport maintained by RCLVs within a finite time interval, the coherent cores of RCLVs are considerably smaller in scale as compared to those of SSH eddies. The average zonal transports induced by SSH eddies and RCLVs are estimated to be −0.82 Sv and −0.51 Sv (1 Sv = 106 m3/s), respectively. For non-overlapping SSH eddies with RCLVs, approximately 80% of the water within the eddy leaks out during the eddy's lifespan. In the case of overlapping SSH eddies, the ratio of coherent water inside the eddy decreases with increasing radius, and the leakage rate is around 58%. Finally, an examination of 36 shedding RCLVs events from the Kuroshio near the Luzon Strait, which induce an average zonal transport of −0.14 Sv, reveals that 54% of the water within the shedding RCLVs originates from the Kuroshio. [ABSTRACT FROM AUTHOR]

In this study, we present an extension to existing numerical retrackers of synthetic-aperture radar (SAR) altimetry signals. To our knowledge at the time of writing this manuscript, it offers the most consistent retrieval of geophysical parameters compared to low-resolution mode (LRM) retracking results. We achieve this by additionally estimating the standard deviation of vertical wave-particle velocities σ v and a new parameter u x , linked to a residual Doppler in the returned radar echoes, which can be related to wind speed and direction. Including this new parameter into the SAR stack retracker mitigates sea surface height estimation errors by up to two centimeters for Sentinel-6MF SAR mode results. Additionally, we found a closed-form equation to describe u x as a function of eastward and northward wind variables, which allows mitigating the effects of this parameter on a SAR stack within level 1B processing and generating a lookup table to correct sea surface height estimates in SAR mode. This additionally opens up the door to estimating the wind speed and direction from SAR altimetry stacks. Additionally, we discuss how this new retracker performs with respect to different planned future baseline processor changes of Sentinel-6MF, namely F09 and F10, by attempting to imitate their level 2 processing. This is achieved by processing cycles 017 to 051 (nearly a full year) of Sentinel-6MF level 1A data on a global scale. We observe that the new retracking method is, on average, more accurate with respect to LRM. However, there is a slight increase in measurement noise due to the introduction of an additional parameter. To ensure that the results of the new retracker are not biased, we retrack using both the new method and the SINCS-OV ZSK retracker on Sentinel-6MF stack data produced in a Monte Carlo simulation. We analyze the simulation results with respect to accuracy, precision, and correlations between estimated parameters. We show that the accuracy of the new retracker is better than SINCS-OV ZSK but less precise, which could be related to higher correlation coefficients—especially with respect to the new parameter u x —between estimated parameters. [ABSTRACT FROM AUTHOR]

For over 25 years, satellite altimetry has provided invaluable information about the ocean dynamics at many scales. In particular, gridded sea surface height (SSH) maps allow us to estimate the mesoscale geostrophic circulation in the ocean. However, conventional interpolation techniques rely on static optimal interpolation schemes, hence limiting the estimation of non-linear dynamics at scales not well sampled by altimetry (i.e., below 150–200 km at mid-latitudes). To overcome this limitation in the resolution of small-scale SSH structures (and thus small-scale geostrophic currents), a back-and-forth nudging algorithm combined with a quasi-geostrophic model, a technique called BFN-QG, has been successfully applied on simulated SSH data in observing system simulation experiments (OSSEs). The result is a significant reduction in interpolation error and an improvement in the space–time resolutions of the experimental gridded product compared to those of operational products. In this study, we propose that the BFN-QG be applied to real altimetric SSH data in a highly turbulent region spanning a part of the Agulhas Current. The performances are evaluated within observing system experiments (OSEs) that use independent data (such as independent SSH, sea surface temperature and drifter data) as ground truth. By comparing the mapping performances to the ones obtained with operational products, we show that the BFN-QG improves the mapping of short, energetic mesoscale structures and associated geostrophic currents both in space and time. In particular, the BFN-QG improves (i) the spatial effective resolution of the SSH maps by a factor of 20 %, (ii) the zonal and (especially) the meridional geostrophic currents, and (iii) the prediction of Lagrangian transport for lead times up to 10 d. Unlike the results obtained in the OSSEs, the OSEs reveal more contrasting performances in low-variability regions, which are discussed in the paper. [ABSTRACT FROM AUTHOR]

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 (<150 km) have a strong signature on the horizontal geostrophic strain rate and for all eddy diagnostics in the Western Boundary Current and along the meandering Agulhas Extension. We investigate the horizontal cascade of energy using a coarse-graining technique, and we observe that the wavelength range where the inverse cascade occurs is biased towards larger mesoscale wavelengths with today's altimetric sampling. We also calculate the projected sampling of the eddy diagnostics under the SWOT swaths built with the NASA–CNES simulator to include the satellite position and realistic noise. For the swaths, a neural network noise mitigation method is implemented to reduce the residual SWOT random error before calculating eddy diagnostics. In terms of SSH, observable wavelengths of 15 to 20 km are retrieved after neural network noise mitigation, as opposed to wavelengths larger than 40 km before the noise reduction. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Although El Niño‐Southern Oscillation (ENSO) and its global impacts through teleconnection have been known for decades, if and how the mean currents and mesoscale eddies in the Caribbean Sea are linked to ENSO remains an open question. Here, by analyzing satellite observations and an ocean reanalysis product, we found a close connection between mean currents, eddies in the Caribbean Sea and ENSO on interannual timescales. Strong El Niño events result in enhanced north‐south sea surface height differences and consequently stronger mean currents in the Caribbean Sea, and the opposite happens during La Niña events. The eddy kinetic energy responds to ENSO via eddy‐mean flow interaction, primarily through baroclinic instability, which releases the available potential energy stored in the mean currents to mesoscale eddies. Our results suggest some predictability of the mean currents and eddies in the Caribbean Sea, particularly during strong El Niño and La Niña events. Plain Language Summary: We explored the potential impacts of El Niño‐Southern Oscillation (ENSO) on the circulation and mesoscale eddies in the Caribbean Sea. We found ENSO‐related synchronized changes in mean currents and eddies across the entire Caribbean Sea. The connection between mean currents and ENSO is established through ENSO's impact on the north‐south sea surface height (SSH) difference in the Caribbean Sea, which determines the strength of the geostrophic jet. During strong El Niño events, the easterly wind anomalies will increase the north‐south SSH difference through Ekman transport, and consequently generate stronger mean currents. During strong La Niña events, the opposite happens. Through baroclinic instability, available potential energy stored in the mean currents will be transferred to eddies and results in ENSO‐modified interannual variations of eddy kinetic energy. Our results suggest that interannual variations of mean currents and eddies in the Caribbean Sea might be predictable, particularly during strong El Niño and La Niña events. Key Points: Interannual variations of mean currents and eddies in the Caribbean Sea are linked to El Niño‐Southern Oscillation (ENSO)ENSO‐induced wind pattern changes modulate the north‐south sea surface height differences and hence the mean currents in the Caribbean SeaInterannual variation of eddy kinetic energy in the Caribbean Sea is controlled by baroclinic instability [ABSTRACT FROM AUTHOR]

The influences of sea surface wind on the oceanic mesoscale eddy are complex. By integrating our self-developed surface drifters with satellite observations, we examined the influence of sea surface wind on the distribution of water masses and biomass within the interior of an anticyclonic eddy. Ten drifters were deployed in the northern South China Sea in the spring of 2021. Eventually, six were trapped in an anticyclonic mesoscale eddy for an extended period. Interestingly, the drifters' trajectories were not symmetric around the eddy center, displaying a significant offset of the distance from the wind turns to the southerly wind. Particle tracking experiments demonstrated that this departure could mainly be attributed to wind-driven ageostrophic currents. This is due to the strength of wind-driven ageostrophic currents being more comparable to geostrophic currents when accompanied by a deflection between the directions of the wind-driven current and the eddy's translation. The drifters' derived data indicated that sub-mesoscale ageostrophic currents within the eddy contributed to this asymmetric trajectory, with Ekman and non-Ekman components playing a role. Furthermore, the evolution of ocean color data provided corroborating evidence of these dynamic processes, highlighting the importance of ageostrophic processes within mesoscale eddies. [ABSTRACT FROM AUTHOR]

Ocean bottom pressure-gauge (OBP) records play an essential role in seafloor geodesy. Oceanographic fluctuations in OBP data, however, pose as a significant noise source in seafloor transient crustal deformation observations, including slow slip events (SSEs), making it crucial to evaluate them quantitatively. To extract the significant fluctuation phenomena common to multiple observation networks, including oceanographic fluctuations and tectonic signals, we applied principal component analysis (PCA) to the 3-year Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) OBP time series for 40 stations during 2016–2019. PCA could separate several oceanographic signals based on the characteristics of their spatial distributions, although evident transient tectonic signals could not be confirmed from the observed pressure records during this observed period. The spatial distribution of the first four principal components (PCs) reflected the common component, inclined component along sea depth, longitude component, and parabola-like pattern, respectively. By subtracting each PC (in particular, PC-2 and PC-4) from the time series, we could significantly reduce the sea depth dependence of OBP records, which has been highlighted in several previous studies and is also evident in this region. We interpreted PCs 2–4 as the reflection of the strength and meandering of ocean geostrophic currents based on a comparison with the PC spatial distribution of the numerical oceanographic models. In addition, to evaluate the ability of PCA to separate transient tectonic signal from OBP time series, including oceanographic fluctuations, we conducted a synthetic ramp assuming an SSE by rectangular fault and then applied PCA. The assumed synthetic tectonic signal could be separated from the oceanographic signals and included in the principal component independently depending on its amplitude, suggesting that the spatial distribution of each PC would change if the amplitude of the synthetic signal were sufficiently large. We propose a transient event-detection method based on the spatial distribution difference of a specific PC with or without a tectonic signal. We used the normalized inner product (NIP) between these PCs as the indicator of their similarities. This method can detect transient tectonic signals more significantly than the moment-magnitude scale of 5.9 from OBP records. [ABSTRACT FROM AUTHOR]

Globally available and highly vertical resolved wind fields are crucial for the analysis of atmospheric dynamics for the benefit of climate studies. Most observation techniques have problems to fulfill both requirements. Especially in the tropics and in the southern hemisphere more wind data availability is required. In this study we investigate the potential of radio occultation (RO) data for climate-oriented wind field monitoring in the tropics, with a specific focus on the equatorial area between ±5° latitude. In this region, the geostrophic balance breaks down, due to the Coriolis force term approaching zero. One further aim is to understand how the individual wind components of the geostrophic balance and equatorial balance approximations bridge across the equator and where each component breaks down. We analyze the equatorial balance equation within this latitude band. In a wider range over the tropics, we derive the RO wind fields also using the geostrophic approximation and we compared the RO winds with ERA5 data. From analyzing first the zonal and meridional wind component, we find that the meridional wind component is more volatile in its derivation, however the total wind speed benefits from a computation of both wind components. Investigating next the bias between the RO and ERA5 computed winds, we find that the systematic data bias is smaller than the bias resulting from the approximation itself. As a final aspect we inspected the monthly-mean RO wind data over the full example year 2009. The bias in the upper troposphere and lower stratosphere is mainly smaller than ±2ms-1, which is in line with the wind field requirements of the World Meteorological Organization. This is encouraging for the use of RO wind fields in climate monitoring over the entire globe including the equatorial region. [ABSTRACT FROM AUTHOR]

By decomposing observed El Niño and La Niña events into a strong group and a weak group, respectively, we discovered that the strong La Niña group has its peak center more towards the west compared to the weak La Niña group, whereas the strong El Niño group has its peak center more towards the east compared to the weak El Niño group. The cause of this structure asymmetry is investigated through an ocean mixed-layer heat budget analysis. It was found that the asymmetry is closely linked to the longitudinal distribution of SST anomaly (SSTA) skewness along the equator, and is fundamentally caused by nonlinear dynamic heating, especially nonlinear horizontal temperature advection. It was demonstrated that near the equatorial central Pacific, the anomalous zonal and meridional currents generate negative nonlinear zonal and meridional temperature advection anomalies for both the El Niño and La Niña events, thus favoring a stronger La Niña and a weaker El Niño. Over the eastern Pacific, due to the dominant geostrophic zonal current anomalies and the southward shift of SSTA centers, nonlinear horizontal temperature advection anomalies tend to be positive for both the El Niño and La Niña, thus favoring a stronger growth of El Niño than La Niña. Nonlinear vertical temperature advection anomalies play minor roles in the central Pacific and tend to partially offset the nonlinear horizontal advection effect in the equatorial eastern Pacific. [ABSTRACT FROM AUTHOR]

Global long-term stable wind fields are valuable information for climate analyses of atmospheric dynamics. Their monitoring remains a challenging task, given shortcomings of available observations. One promising option for progress is the use of radio occultation (RO) satellite data, based on which the winds are estimated using the geostrophic approximation. Hence, in this study we focus on two goals, explored through European Re-Analysis ERA5 and RO datasets, using monthlymean January and July data over 2007-2020 with a 2.5° × 2.5° resolution. First, we compare actual and geostrophic ERA5 wind speeds to evaluate the validity of the geostrophic approximation. Second, we test how well ERA5 and RO geostrophic winds agree. We find the geostrophic approximation to work well within 2 m s-1 accuracy almost globally (5°-85° latitude), especially over the summer hemisphere; larger differences (more than 5 m s-1) may occur in the winter stratosphere. We noticed the effect of large mountain ranges on the wind flow as a wave-like pattern, also in the difference between RO and ERA5 geostrophic winds, pointing to effects of different geopotential height estimations. Generally, RO and ERA5 geostrophic winds showed very good agreement. In the long-term, systematic differences in decadal trends of higher than 0.5 m s-1 per decade were found at subtropical latitudes, mainly related to observing system changes in the year 2016 that influenced ERA5. Together with the validity of the geostrophic approximation, this indicates that the long-term stability of RO-derived wind field monitoring can provide added value to reanalysis winds, for the benefit of climate monitoring and analyses. [ABSTRACT FROM AUTHOR]

It is not fully understood how near-inertial kinetic energy (NIKE) is spatially distributed near Tsushima oceanic front (TOF) as a typhoon travels across the region. Underneath TOF, a year-round mooring covering a major part of water column was implemented in 2019. During summer, three massive typhoons (Krosa, Tapah, and Mitag) consecutively traversed the frontal area and delivered a substantial amount of NIKE into surface mixed layer. According to a mixed-layer slab model, NIKE was widely distributed near the cyclone's track. The mooring observation exhibited the vertical distribution and pathways of surface-generated NIKE in response to the successive typhoon events. According to the modal decomposition, first three modes mostly explain the NIKE's elevations following the typhoon events. According to ray-tracing experiments based on the internal-wave theory, large-scale near-inertial waves (NIWs) rapidly descend to a depth greater than 1000 m, while mesoscale NIWs slowly descend and rarely reached beyond the main pycnocline. Following the passage of Tapah, a profound energy mass was found nearly stationary at shallow depths coincident with vertical shear of geostrophic current. We infer that the descending rate of NIWs fell and then they were amplified through the energy conservation when the waves came from the north side of TOF. [ABSTRACT FROM AUTHOR]

We expand on a recent determination of the first global energy spectrum of the ocean's surface geostrophic circulation (Storer et al., 2022, https://doi.org/10.1038/s41467-022-33031-3) using a coarse‐graining (CG) method. We compare spectra from CG to those from spherical harmonics by treating land in a manner consistent with the boundary conditions. While the two methods yield qualitatively consistent domain‐averaged results, spherical harmonics spectra are too noisy at gyre‐scales (>1,000 km). More importantly, spherical harmonics are inherently global and cannot provide local information connecting scales with currents geographically. CG shows that the extra‐tropics mesoscales (100–500 km) have a root‐mean‐square (rms) velocity of ∼15 cm/s, which increases to ∼30–40 cm/s locally in the Gulf Stream and Kuroshio and to ∼16–28 cm/s in the ACC. There is notable hemispheric asymmetry in mesoscale energy‐per‐area, which is higher in the north due to continental boundaries. We estimate that ≈25%–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)‐resolved by pre‐SWOT satellite products. Spectra of the time‐mean circulation show that most of its energy (up to 70%) resides in stationary eddies with characteristic scales smaller than (<500 km). This highlights the preponderance of "standing" small‐scale structures in the global ocean due to the temporally coherent forcing by boundaries. By coarse‐graining in space and time, we compute the first spatio‐temporal global spectrum of geostrophic circulation from AVISO and NEMO. These spectra show that every length‐scale evolves over a wide range of time‐scales with a consistent peak at ≈200 km and ≈2–3 weeks. Plain Language Summary: Traditionally, "eddies" are identified as time‐varying features relative to a background time‐mean flow. As such, "mean" does not imply large length‐scale. Standing eddies or meanders due to topography have little time‐variation, but can have significant energy at small length‐scales that are unresolved and need to be parameterized in coarse climate simulations. Similarly, "eddy" or "time‐varying" do not imply small length‐scale, such as large‐scale motions from Rossby waves or fluctuations of the Kuroshio. Another common method is Fourier analysis in "representative" ocean boxes that cannot capture the circulation's planetary scales. We overcome these limitations thanks to recent advances: (a) a method for calculating spectra by coarse‐graining, (b) properly defining convolutions on the sphere, which "blur" oceanic flow in a way that preserves its underlying symmetries, opening the door for global "wavelet" analysis and, more generally, spatial coarse‐graining, and (c) FlowSieve: an efficient parallel code. We employ coarse‐graining in space‐time to gain new insights into the global oceanic circulation, including how much energy resides in its different spatial structures and how they vary in time. Key Points: Coarse‐graining, which disentangles flow concurrently in scale and space, reveals hemispheric asymmetry in mesoscale energy‐per‐area due to boundariesCoarse‐graining spectra of the time‐mean velocity show that most (up to 70%) of its energy resides in "standing" small‐scale eddies <500 kmWe estimate that ≈25%–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)‐resolved by pre‐SWOT satellite products [ABSTRACT FROM AUTHOR]

We present a study on the ocean surface currents reconstruction by merging Level-4 (L4, gap-free) altimeter-derived geostrophic currents and satellite sea surface temperature. Building upon past studies on the multi-variate reconstruction of geostrophic currents from satellite observations, we regionalized and optimized an algorithm to improve the altimeter-derived surface circulation estimates in the North Atlantic Ocean. A ten-year-long time series (2010–2019) is presented and validated by means of in situ observations. The newly optimized algorithm allowed us to improve the currents estimate along the main axis of the Gulf Stream and in correspondence of well-known upwelling areas in the North Eastern Atlantic, with percentage improvements of around 15% compared to standard operational altimetry products.

Gridded sea surface height (SSH) maps estimated from satellite altimetry are widely used for estimating surface ocean geostrophic currents. Satellite altimeters observe SSH along one‐dimensional tracks widely spaced in space and time, making accurately reconstructing the two‐dimensional (2D) SSH field challenging. Traditionally, SSH is mapped using optimal interpolation (OI). However, OI artificially smooths the SSH field leading to high mapping errors in regions with rapidly‐evolving mesoscale features such as western boundary currents. Motivated by the dynamical relation between SSH and sea surface temperature (SST) and the notion that even the chaotic evolution of mesoscale ocean turbulence may contain repeating patterns, we outline a deep learning (DL) approach where a neural network is trained to reconstruct 2D SSH by synthesizing altimetry and SST observations. In the Gulf Stream Extension region, dominated by mesoscale variability, our DL method substantially improves the SSH reconstruction compared to existing methods. Our SSH map has 17% lower root‐mean‐square error and resolves spatial scales 30% smaller than OI compared against independent altimeter observations. Surface geostrophic currents calculated from our map are closer to surface drifter observations and appear qualitatively more realistic, with stronger currents, a clearer separation between the Gulf Stream and neighboring eddies, and the appearance of smaller coherent eddies missed by other methods. Our map yields significant re‐estimations of important dynamical quantities such as eddy kinetic energy, vorticity, and strain rate. Applying our DL method to produce a global SSH product may provide a more accurate and higher resolution product for studying mesoscale ocean turbulence. Plain Language Summary: Satellites observe small variations in the height of the sea surface but with large gaps in the observations. Having an estimate of the two‐dimensional sea surface height field allows one to estimate surface ocean currents, so filling in the gaps between the observations is an important problem. The traditionally‐used method for filling in the gaps between sea surface height observations struggles when there are lots of small‐scale, rapidly‐interacting ocean currents. We developed a deep learning (DL) model to estimate the sea surface height field more accurately. We achieved this by combining the sea surface height observations with satellite observations of sea surface temperature. The relationship between sea surface temperature and height is non‐trivial, but our DL model learned to use information from the sea surface temperature observations in the places where sea surface height wasn't observed to improve the accuracy of the sea surface height estimate. We applied and tested our method in the Gulf Stream and demonstrated that our sea surface height map is more accurate than that from the traditional method. Our method provides a more accurate sea surface height map which could allow us in future to learn new lessons about small‐scale surface currents in the ocean. Key Points: We developed a deep learning method that significantly improves the accuracy and resolution of gridded sea surface height anomaliesThis data‐driven method takes advantage of combining sea surface temperature and altimetry observationsInferred surface geostrophic currents are quantitatively and qualitatively more realistic than those from existing sea surface height maps [ABSTRACT FROM AUTHOR]

The ocean circulation dynamics can be represented as a high‐dimensional multi‐scale nonlinear system with inhomogenous meanfields and topography. Numerical simulations of the ocean dynamics resolve the large scales of motion on a computational grid, with the unresolved subgrid interactions parameterized. These simulations are highly dependent upon the grid resolution, unless the subgrid terms are appropriately parameterized. There are five fundamental classes of subgrid interactions: eddy‐eddy; eddy‐topographic; eddy‐meanfield; meanfield‐meanfield; and meanfield‐topographic. Scale dependent parameterizations representing each of these interaction classes are presented here in oceanic flows for the first time. Subgrid parameterizations are calculated and validated in baroclinic quasi‐geostrophic simulations of idealized Antarctic Circumpolar Current flows with representative mean currents and ocean floor topography. The parameterization coefficients are derived from the coarse grained statistics of high resolution reference simulations in spectral space. Stochastic and deterministic parameterizations are developed for the eddy‐eddy interactions, and deterministic forms for the remaining classes. The kinetic energy spectra and meanfield resulting from large eddy simulations (LES) adopting these coefficients accurately replicate those of the reference simulation. The eddy‐eddy interactions are dominant, but all classes need to be parameterized for best results. For LES where baroclinic instability is explicitly resolved the stochastic variants out‐perform the deterministic ones across all scales. When baroclinic instability is not explicitly resolved, the stochastic variants out‐perform the deterministic cases at the large scales, but introduce some distortions at the smallest resolved scales. This study provides an assessment of the relevant strengths of the subgrid interaction classes in an idealized yet representative ocean. Plain Language Summary: Numerical simulation of oceanic circulations is incredibly challenging due in part to the large range of scales. There are oscillations on the size of entire ocean basins, through to fine scale centimeter sized turbulence, and everything in between. Global simulations of the ocean typically capture the large scale features on a discrete grid. To ensure these simulations are physically representative, the influence of the missing smaller scales needs to be captured accurately. Here we present a first calculation of parameterizations governing the interactions between the resolved dynamics and the missing small scale eddies, meanfield and topographic features. Reduced resolution simulations incorporating these parameterizations without tuning have statistics in very good agreement with the high resolution simulations from which they were derived. Key Points: Subgrid models of all interaction classes between eddies, meanfield and topography are calculated from reference global ocean simulationsInter‐eddy interactions are modeled both stochastically and deterministically, with deterministic models used for the other interactionsCoarse simulations using these subgrid coefficients without tuning are in very good agreement with those of the high‐resolution simulation [ABSTRACT FROM AUTHOR]

The Kuroshio Current in the North Pacific displays path changes on an interannual‐to‐decadal time scale. In an idealized barotropic quasi‐geostrophic model of the double‐gyre wind‐driven circulation under stochastic wind‐stress forcing, such variability can occur due to transitions between different equilibrium states. The high‐dimensionality of the problem makes it challenging to determine the probability of these transitions under the influence of stochastic noise. Here we present a new method to estimate these transition probabilities, using a Dynamical Orthogonal (DO) field approach. In the DO approach, the solution of the stochastic partial differential equations system is decomposed using a Karhunen–Loève expansion and separate problems arise for the ensemble mean state and the so‐called time‐dependent DO modes. The original method is first reformulated in a matrix approach which has much broader application potential to various (geophysical) problems. Using this matrix‐DO approach, we are able to determine transition probabilities in the double‐gyre problem and to identify transition paths between the different states. This analysis also leads to the understanding which conditions are most favorable for transition. Plain Language Summary: In an idealized wind‐driven double‐gyre circulation model of the midlatitude ocean circulation, a northward displaced jet (jet‐up state) and a southward displaced jet (jet‐down state) exist under the same forcing conditions. The existence of these states can help explain behavior of, for example, the Kuroshio Current in the North Pacific, which undergoes path variations on an interannual‐to‐decadal time scale. Under the influence of wind‐stress noise, the model circulation may undergo a transition between the jet‐up state to the jet‐down state and vice versa. Determining transition probabilities between the two states is challenging as the model has many degrees of freedom and transition probabilities are very small when the noise is low. In this paper, we use a new formulation of a model order reduction technique, the so‐called Dynamical Orthogonal (DO) field approach, to determine both transition probabilities and transition paths. Key Points: A matrix version of the Dynamical Orthogonal (DO) field method is presentedDO is used to determine transition probabilities in multi‐stable stochastic dynamical systemsThe method is applied to study transition behavior of ocean western boundary currents [ABSTRACT FROM AUTHOR]

Residual current analysis of multiple stations' periodic observational data for sea currents, and multiple voyages, multiple seasons from 2018 to 2022, revealed the existence of a strong southwest current zone, 20–30 m underwater and within the coastal current area of Yangjiang, western Guangdong, in the northwest region of the South China Sea. The velocity of the residual currents in the strong–current zone was 38.2–100.3 cm/s. Observational data for wind, sea currents, salinity and tide from multiple coastal stations in the spring (from March 2019 to May 2019) and summer (from June 2019 to August 2019) of 2019 demonstrated the existence of abrupt strong currents in the coastal current sea area of Yangjiang, western Guangdong. Analyses of continuous sequence hydro meteorological data for the Yangjiang area indicated that wind stress was the main factor determining the direction of sea currents; increases in the near-shore water level produce a westward geostrophic current and a wind current in the shallow sea area; their joint effect was the key factor determining the velocity of the coastal current in this sea area. In spring and summer, when Pearl River runoff into the sea reaches a peak, and under the action of the northeast wind, the water level on the west coast of Guangdong rose. This created a barotropic pressure–gradient force from the shore to the outer sea, generating a southwestward geostrophic current in the same direction as the wind-driven current. The joint action of the wind-driven and geostrophic currents then generated a sudden southwestward coastal current. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

A high precision and fine resolution reconstruction of sea surface current is beneficial to the exploration of complicated ocean dynamic processes. Existing studies mainly use satellite sea level and wind stress fields to derive sea surface geostrophic and Ekman currents, and build physical inversion models for global or regional oceans. Despite the obvious success, there are a variety of typical dynamic processes in the ocean such as mesoscale eddies and small-scale waves, and any product of surface current that neglects the contribution of wave motion would be, at best, incomplete. In this context, we present a precise sea surface current product at 15 m depth named GEST (Geostrophic-Ekman-Stokes-Tide) by analyzing the coupling relationship between ocean surface components that correspond to different physical processes and the actual currents as observed by drifting buoys. The GEST was generated based on Ekman, geostrophic currents, and waved-induced Stokes drift and TPXO9 tidal currents. Specifically, in the calculation of Ekman currents, local applicability is taken into account to ensure that the friction layer of wind-driven current can reach where drifters could operate as normal. A comparison proves that combining multi-scale theory and the condition of local applicability improves estimation results by up to 3.6 cm/s compared with OSCAR product, and 0.3 cm/s with GlobCurrent. Furthermore, by comparing the reconstructed products with 1° and 0.25° resolution, we find that the higher resolution not only reveals more details of ocean currents especially mesoscale eddy energy associated with geostrophic currents, but also improves the accuracy by up to 5.62 cm/s. [ABSTRACT FROM AUTHOR]

This study assesses the capability of Surface Quasi-Geostrophy (SQG) to reconstruct the three-dimensional (3D) dynamics in four critical areas of the Arctic Ocean: the Nordic, Barents, East Siberian, and Beaufort Seas. We first reconstruct the upper ocean dynamics from TOPAZ4 reanalysis of sea surface height (SSH), surface buoyancy (SSB), and surface velocities (SSV) and validate the results with the geostrophic and total TOPAZ4 velocities. The reconstruction of upper ocean dynamics using SSH fields is in high agreement with the geostrophic velocities, with correlation coefficients greater than 0.8 for the upper 400 m. SSH reconstructions outperform surface buoyancy reconstructions, even in places near freshwater inputs from river discharges, melting sea ice, and glaciers. Surface buoyancy fails due to the uncorrelation of SSB and subsurface potential vorticity (PV). Reconstruction from surface currents correlates to the total TOPAZ4 velocities with correlation coefficients greater than 0.6 up to 200 m. In the second part, we apply the SQG approach validated with the reanalysis outputs to satellite-derived sea level anomalies and validate the results against in-situ measurements. Due to lower water column stratification, the SQG approach's performance is better in fall and winter than in spring and summer. Our results demonstrate that using surface information from SSH or surface velocities, combined with information on the stratification of the water column, it is possible to effectively reconstruct the upper ocean dynamics in the Arctic and Subarctic Seas up to 400 m. Future remote sensing missions in the Arctic Ocean, such as SWOT, Seastar, WaCM, CIMR, and CRISTAL, will produce enhanced SSH and surface velocity observations, allowing SQG schemes to characterize upper ocean 3D mesoscale dynamics up to 400 m with higher resolutions and lower uncertainties. [ABSTRACT FROM AUTHOR]

The atmosphere of Venus exhibits equatorial planetary-scale waves that are suspected to play an important role in its complex atmospheric circulation. Due to its particularly long sidereal day (243 terrestrial days against 24 h for the Earth), the Venusian waves must be described with the momentum equations for a cyclostrophic regime, but efforts to derive analytical wave solutions have been scarce. Following a classic approach for the terrestrial quasi-geostrophic regime, I present analytical solutions for equatorial waves in the atmosphere of Venus, assuming a single layer of a homogeneous incompressible fluid with a free surface and focusing on two asymptotic cases described by the ratio of their non-dimensional frequency and zonal wavenumber. One of the dispersion relations that has been obtained describes waves on a small spatial scale propagating upstream relative to the zonal flow, which is associated with a Rossby-type wave called "centrifugal". The solutions for the other asymptotic case were interpreted as inertio-surface waves, which describe planetary-scale waves that can propagate "upstream" and "downstream" relative to the zonal winds and have null group velocity. These new wave solutions stress relevant differences between waves in geostrophic and cyclostrophic regimes and may be applicable to Saturn's moon, Titan, and Venus-like exoplanets. [ABSTRACT FROM AUTHOR]

Coastal polynyas in the Ross Sea are important source regions of high-salinity shelf water (HSSW) – the precursor of Antarctic Bottom Water that supplies the lower limb of the thermohaline circulation. Here, the response of sea ice production and HSSW formation to synoptic-scale and mesoscale cyclones was investigated for the Ross Ice Shelf Polynya (RISP) using a coupled ocean–sea ice–ice shelf model targeted on the Ross Sea. When synoptic-scale cyclones prevailed over RISP, sea ice production (SIP) increased rapidly by 20 %–30 % over the entire RISP. During the passage of mesoscale cyclones, SIP increased by about 2 times over the western RISP but decreased over the eastern RISP, resulting respectively from enhancement in the offshore and onshore winds. HSSW formation mainly occurred in the western RISP and was enhanced responding to the SIP increase under both types of cyclones. Promoted HSSW formation could persist for 12–60 h after the decay of the cyclones. The HSSW exports across the Drygalski Trough and the Glomar Challenger Trough were positively correlated with the meridional wind. Such correlations are mainly controlled by variations in geostrophic ocean currents that result from sea surface elevation change and density differences. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We assess the contribution of flow over the Northeast Greenland Shelf (NEGS) to the total freshwater transport (FWT) through the Fram Strait. We analyze available observations since 2013 consisting of hydrographic sections, new mooring data, and surface‐geostrophic velocity and ERA5 winds. We provide first seasonal estimates of the FWT over the entire shelf and find a net southward FWT between 47 and 56 mSv (reference salinity 34.9), or ∼40%–45% of the total FWT through the Fram Strait. In summer, the long‐term mean FWT east of 8°W cumulates to 27 mSv southward, there is no trend toward increasing FWT, and new mooring data show a large seasonal cycle. We show that the NEGS contributes significantly to the Arctic freshwater outflow, however, uncertainty is high. To close the Arctic freshwater budget and provide better FWT estimates to the North Atlantic, sustained year‐round observations on the NEGS are essential. Plain Language Summary: We use new observations since 2013 to calculate transport over the Northeaster Greenland Shelf for the first time, and investigate how much it adds to the transport of the East Greenland Current over the shelf break which is monitored by moorings. Over most of the shelf, direct observations of ocean properties are limited to summer. Thus, we use remote‐sensing and wind reanalysis data, with shelf‐wide and year‐round coverage, but higher uncertainty, to estimate the year‐round transport of low‐salinity water (freshwater, less than 34.9) over the entire shelf. We find that on the shelf, freshwater transport is at least 40% of the total in Fram Strait. Still, uncertainty is high and year‐round, direct observations in the western part of the shelf are required to provide accurate estimates for the Arctic freshwater outflow. Key Points: Summer freshwater transport (FWT) over the shelf cumulates to −27 mSv and shows no trend in timeNew mooring data on the shelf at 10.1°W show a very large seasonal cycleYear‐round FWT over the shelf is estimated to ∼−52 mSv or 43% of the total freshwater outflow from the Fram Strait [ABSTRACT FROM AUTHOR]

Barotropic and baroclinic tides were simulated for the Coral and Tasman Seas off eastern Australia and compared to a simulation without tides. Both simulations included geostrophic currents and replicated a current analogous to the East Australian Current (EAC). Tides and tidal currents in most of this region are weak, generally 1–2 cm s−1, with the exception of the far northwest portion of the Coral Sea. Even these weak tides were found to impact the mean EAC-like current transport, enhancing it by 1–4 Sv in some areas. Southward flow increased over the continental shelf. Tides did not appear to impact eddy formation, size, or rotational speed; however, they affected eddy propagation. Cyclonic eddies propagated northward faster with tides than without tides. Tides impacted cross-shelf transport of colder water, with significantly more on-shore transport occurring with tides, particularly equatorward of the diurnal critical latitudes. Cross-shelf transport of nutrient rich water onto the shelf is important in this oligotrophic region. Although the prime source of vertical shear and mixing were mean currents and eddies, tides played a secondary role. Tides influenced mixing by increasing vertical temperature diffusivities to 10−4 to 10−3 m2 s−1 over portions of the continental slope and over rough topography, particularly in regions near the diurnal critical latitudes (27°–30°). In conclusion, even small tides can significantly impact the circulation through their effects on the mean currents, eddy rotation velocities, eddy propagation, and mixing. [ABSTRACT FROM AUTHOR]

Global ocean surface currents estimated from satellite derived data based on a regular global grid are affected by the grid's shape and placement. Due to different neighbourhood relationships, the rectangular lat/lon grids lose accuracy when interpolating and fitting elevation data. Hexagonal grids have shown to be advantageous due to their isotropic, uniform neighbourhood. Considering these merits, this paper aims to estimate global ocean surface current using a global isotropic hexagonal grid from satellite remote sensing data. First, gridded satellite altimeter data and wind data with different resolutions are interpolated into the centre of the global isotropic hexagonal grid. Then, geostrophic and Ekman currents components are estimated according to the Lagerlof Ocean currents theory. Finally, the inversion results are verified. By analyzing the results, we conclude that the ocean surface currents estimated based on the global isotropic hexagonal grid have considerable accuracy, with improvement over rectangular lat/lon grids. [ABSTRACT FROM AUTHOR]

Satellite altimetry missions flying over the ice-covered Arctic Ocean have opened the possibility of further understanding changes in the ocean beneath the sea ice. This requires complex processing of satellite signals emerging from the sea surface in leads within the sea ice, with efforts to generate consistent Arctic-wide datasets of sea surface height ongoing. The aim of this paper is to provide and assess a novel gridded dataset of sea surface height anomaly and geostrophic velocity, which incorporates both the ice-covered and open ocean areas of the Arctic. Data from the CryoSat-2 mission in the period 2011–2020 were gridded at monthly intervals, up to 88 ∘ N, using the Data-Interpolating Variational Analysis (DIVA) method. To examine the robustness of our results, we compare our dataset to independent satellite data, mooring time series and Arctic-wide hydrographic observations. We find that our dataset is well correlated with independent satellite data at monthly timescales. Comparisons to in situ ocean observations show that our dataset provides reliable information on the variability of sea surface height and surface geostrophic currents over geographically diverse regions of the Arctic Ocean and different dynamical regimes and sea ice states. At all comparison sites we find agreement with in situ observed variability at seasonal to interannual timescales. Furthermore, we find that our geostrophic velocity fields can resolve the variability of boundary currents wider than about 50 km, a result relevant for studies of Arctic Ocean circulation. Additionally, large-scale seasonal features emerge. Sea surface height exhibits a wintertime Arctic-wide maximum, with the highest amplitude over the shelves. Also, we find a basin-wide seasonal acceleration of Arctic slope currents in winter. We suggest that this dataset can be used to study not only the large-scale sea surface height and circulation, but also the regionally confined boundary currents. The dataset is available in netCDF format from PANGAEA at 10.1594/PANGAEA.931869. [ABSTRACT FROM AUTHOR]

We present a new gridded sea surface height and current dataset produced by combining observations from nadir altimeters and drifting buoys. This product is based on a multiscale and multivariate mapping approach that offers the possibility to improve the physical content of gridded products by combining the data from various platforms and resolving a broader spectrum of ocean surface dynamic than in the current operational mapping system. The dataset covers the entire global ocean and spans from 1 July 2016 to 30 June 2020. The multiscale approach decomposes the observed signal into different physical contributions. In the present study, we simultaneously estimate the mesoscale ocean circulations as well as part of the equatorial wave dynamics (e.g. tropical instability and Poincaré waves). The multivariate approach is able to exploit the geostrophic signature resulting from the synergy of altimetry and drifter observations. Sea-level observations in Arctic leads are also used in the merging to improve the surface circulation in this poorly mapped region. A quality assessment of this new product is proposed with regard to an operational product distributed in the Copernicus Marine Service. We show that the multiscale and multivariate mapping approach offers promising perspectives for reconstructing the ocean surface circulation: observations of leads contribute to improvement of the coverage in delivering gap-free maps in the Arctic and observations of drifters help to refine the mapping in regions of intense dynamics where the temporal sampling must be accurate enough to properly map the rapid mesoscale dynamics. Overall, the geostrophic circulation is better mapped in the new product, with mapping errors significantly reduced in regions of high variability and in the equatorial band. The resolved scales of this new product are therefore between 5 % and 10 % finer than the Copernicus product (https://doi.org/10.48670/moi-00148 , Pujol et al., 2022b). [ABSTRACT FROM AUTHOR]

This study focuses on the determination and validation of a new global mean sea surface (MSS) model, named the Shandong University of Science and Technology 2020 (SDUST2020), with a grid size of 1′×1′. This new model was established with a 19-year moving average method and fused multi-satellite altimetry data over a 27-year period (from January 1993 to December 2019). The data of HaiYang-2A, Jason-3, and Sentinel-3A were first ingested in the SDUST2020 MSS but not in any other global MSS model, such as the CLS15 and DTU18 MSS models. Validations, including comparisons with the CLS15 and DTU18 MSS models, GPS-leveled tide gauges, and altimeter data, were performed to evaluate the quality of the SDUST2020 MSS model, all of which showed that the SDUST2020 MSS model is accurate and reliable. The SDUST2020 MSS dataset is freely available at the site (data DOI: 10.5281/zenodo.6555990, Yuan et al., 2022). [ABSTRACT FROM AUTHOR]

• Sea surface height and SSH derived geostrophic current from S3A are evaluated. • Measurable accuracy improvements are observed in DD/SAR compared to LRM. • DD/SAR altimetry can resolve small scale current features in the NW Atlantic shelf. In this study, we evaluate Sentinel-3A satellite synthetic aperture radar (SAR) altimeter observations along the Northwest Atlantic coast, spanning the Nova Scotian Shelf, Gulf of Maine, and Mid-Atlantic Bight. Comparisons are made of altimeter sea surface height (SSH) measurements from three different altimeter data processing approaches: fully-focused synthetic aperture radar (FFSAR), un-focused SAR (UFSAR), and conventional low-resolution mode (LRM). Results show that fully-focused SAR data always outperform LRM data and are comparable or slightly better than the nominal un-focused SAR product. SSH measurement noise in both SAR-mode datasets is significantly reduced compared to LRM. FFSAR SSH 20-Hz noise levels, derived from 80-Hz FFSAR data, are lower than 20-Hz UFSAR SSH with 25% noise reduction offshore of 5 km, and 55–70% within 5 km of the coast. The offshore noise improvement is most likely due to the higher native along-track data posting rate (80 Hz for FFSAR, and 20 Hz for UFSAR), while the large coastal improvement indicates an apparent FFSAR data processing advantage approaching the coastlines. FFSAR-derived geostrophic ocean current estimates exhibit the lowest bias and noise when compared to in situ buoy-measured currents. Assessment at short spatial scales of 5–20 km reveals that Sentinel-3A SAR data provide sharper and more realistic measurement of small-scale sea surface slopes associated with expected nearshore coastal currents and small-scale gyre features that are much less well resolved in conventional altimetric LRM data. [ABSTRACT FROM AUTHOR]

Canada ratified the United Nations Convention on the Law of the Sea (UNCLOS) in 2003. With that ratification is an obligation to submit data and information to the U.N. pertaining to the limits of the country's extended continental shelf (ECS); the portion of the juridical continental shelf that extends beyond 200 nautical miles. A team of Canadian scientists, managers, and legal experts that included representation from three Federal Departments (Natural Resources Canada, Fisheries and Oceans Canada, and Global Affairs Canada) with additional support from other departments, spent 13 years compiling and acquiring data to provide the scientific evidence to support delineation of Canada's seaward most maritime limit. The submission has the potential to provide Canada with 2.4 million km2 of additional submarine landmass in the Atlantic and the Arctic oceans over which Canada exercises sovereign rights for the purpose of exploring and exploiting its natural resources. Specific information such as the tectonic framework of the continental margin, the geomorphology of the margin and in particular the continental slope, the geologic nature of adjoined ridges, rises, and plateaux, and sediment thickness within adjacent basins are examples of fundamental pieces of geoscientific information needed to substantiate Canada's outermost maritime limits. This paper highlights a number of segments of Canada's continental margins to showcase this scientific evidence and how it is applied in the UNCLOS context. In doing so, the paper demonstrates the geologic complexity of Canada's margins as illustrated in scientific publications that have resulted from these new data collections, while at the same time presenting new scientific evidence and interpretations. This collection of data and information provides a wealth of new knowledge in Canada's offshore regions. The massive data compilation in the Atlantic led to conception of continental margins, in a source-to-sink scenario, as having an equilibrium base level or graded form, comparable to river systems. Departures from this shape relate to the interplay of sedimentary processes and in particular to those processes that do not fit the source-to-sink paradigm. For example, a significant part of the Atlantic margin is shown to be heavily influenced by along-slope geostrophic currents that generated massive contourite drift deposits. These deposits reflect lateral transport of sediment that had a significant impact on the morphology of the margin. The role of mass transport processes in shaping continental margins is also highlighted, and in particular the collapses of entire segments of the margin were observed. The prominent role mass failure processes play in delivering sediment to the adjacent abyssal plain is also critical in the ECS context. These observations challenge the entrenched notion of a continental margin comprising a shelf, slope, and rise and in particular the concept of the "continental rise". Prior to 2006, regions of the Arctic Ocean seaward of the Canadian landmass had fewer than 5000 km of seismic reflection data. The massive efforts of Arctic coastal States to map their margins for ECS purposes have led to a leap in technological advances to acquire data in ice-covered seas and have led to a wealth of new geoscientific knowledge. Perhaps foremost amongst this knowledge is demonstration that Canada Basin is indeed a fully developed ocean basin, albeit significantly infilled with sediment. Based on this knowledge and identification of related structures, new realistic tectonic scenarios for opening of the Amerasia Basin are proposed that include a significant component of transform or strike-slip motions. With seismic velocity and rock sample information, the continental nature of Alpha and Mendeleev ridges has been substantiated. Even bathymetric data were lacking in the Arctic and new editions of seafloor maps now support grids of 500 m spacing; although some regions remain sparse. Once thought to be relatively stagnant, sedimentary processes such as found in many ocean basins were discovered in the Arctic Ocean. Evidence of geostrophic currents, sediment mass failures, and deep-sea turbidity current channels were found to be ubiquitous, even in the deepest parts of the Arctic's basins. [ABSTRACT FROM AUTHOR]

We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable. The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography. The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing. We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence. [ABSTRACT FROM AUTHOR]

The variation of boundary layer circulation caused by the influence of complex underlying surface is one of the reasons why it is difficult to forecast hourly heavy rainfall (HHR) in the middle Yangtze River Valley (YRV). Based on the statistics of high-resolution observation data, it is found that the low resolution data underestimate the frequency of HHR in the mountain that are between the twain-lake basins in the middle YRV (TLB-YRV). The HHR frequency of mountainous area in the TLB-YRV is much higher than that of Dongting Lake on its left and is equivalent to the HHR frequency of Poyang Lake on its right. The hourly reanalysis data of ERA5 were used to study the variation of boundary layer circulation when HHR occurred. It can be found that the boundary layer circulation corresponding to different underlying surfaces changed under the influence of the weather system. Firstly, the strengthening of the weather system in the early morning resulted in the strengthening of the southwest low-level air flow, which intensified the uplift of the windward slope air flow on the west and south slopes of the mountainous areas in the TLB-YRV. As a result, the sunrise HHR gradually increases from the foot of the mountain. The high-frequency HHR period of sunrise occurs when the supergeostrophic effect is weakened, the low-level vorticity and frontal forcing are strengthened, and the water vapor flux convergence begins to weaken. Secondly, the high-frequency HHR period of the sunset is caused by stronger local uplift and more unstable atmospheric stratification, but the enhanced local uplift is caused by the coupling of the terrain forcing of the underlying surface and the enhanced northern subgeostrophic flow, which causes the HHR to start closer to the mountain top at sunset than at sunrise. [ABSTRACT FROM AUTHOR]

The newly developed fast Earth system model CLIMBER-X is presented. The climate component of CLIMBER-X consists of a 2.5-D semi-empirical statistical–dynamical atmosphere model, a 3-D frictional–geostrophic ocean model, a dynamic–thermodynamic sea ice model and a land surface model. All the model components are discretized on a regular lat–long grid with a horizontal resolution of 5∘×5∘. The model has a throughput of ∼ 10 000 simulation years per day on a single node with 16 CPUs on a high-performance computer and is designed to simulate the evolution of the Earth system on temporal scales ranging from decades to >100000 years. A comprehensive evaluation of the model performance for the present day and the historical period shows that CLIMBER-X is capable of realistically reproducing many observed climate characteristics, with results that generally lie within the range of state-of-the-art general circulation models. The analysis of model performance is complemented by a thorough assessment of climate feedbacks and model sensitivities to changes in external forcings and boundary conditions. Limitations and applicability of the model are critically discussed. CLIMBER-X also includes a detailed representation of the global carbon cycle and is coupled to an ice sheet model, which will be described in separate papers. CLIMBER-X is available as open-source code and is expected to be a useful tool for studying past climate changes and for the investigation of the long-term future evolution of the climate. [ABSTRACT FROM AUTHOR]

Purpose. The aim of the work is to clarify the spatial structure of the climatic dynamic fronts (geostrophic current jets) and to estimate the relationship between their position and that of the large-scale temperature fronts on the surface of the Scotia Sea and the adjacent water areas in the southwestern part of the Atlantic sector of Antarctica. Methods and Results. The daily averaged data arrays of the CMEMS (1993–2017) and NOAA OI SST (1982–2017) reanalysis at the regular 0.25° grid were used. The CMEMS reanalysis contains the sea surface geostrophic velocity values, the NOAA OI SST reanalysis – the sea surface temperature ones which were reduced to the climatic form through their averaging for each month of the corresponding periods. Position of the current jets and the temperature fronts was determined using the maximums of the geostrophic velocity components and the extremes of the temperature horizontal gradients. The updated scheme of the average long-term position of dynamic fronts was constructed. It shows that in the areas of the most pronounced bottom topography inhomogeneities (the northern boundary of the Falkland Plateau and the Tierra del Fuego shelf, the boundaries of the Falkland Islands shelf and the Birdwood Bank, the Shackleton Ridge and the South Shetland Islands shelf), the fronts do not change their latitudinal position during a year. It is revealed that in most of the water area, the temperature horizontal gradient extremes (temperature fronts) correspond to the geostrophic velocity maximums (dynamic fronts). The Northern and Central Branches of the Antarctic Circumpolar Current are most clearly manifested in the temperature field. In general, in the water area under study, the average annual latitudinal position of the Subantarctic and Antarctic Polar Fronts is displaced to the south relative to the position of the Northern and Central branches jets of the Antarctic Circumpolar Current by 0.25–0.5° and 0.25–1°, respectively. Conclusions. It is shown that, being influenced by the bottom topography, the large-scale jets of geostrophic currents form intense topographic meanders and recirculation branches which are stably manifested on the climatic scale. The Antarctic Circumpolar Current branches being affected by the bottom topography, can merge forming the joint flows, and then diverge forming a system of separate jets again. It is found that the main spatial features of frontal structure in the geostrophic velocities field persist throughout the whole year and are conditioned mainly by the bottom topography. The most of the dynamic fronts are shown to be clearly pronounced in the temperature field on the surface during a year. A high level of linear correlation between the positions of temperature fronts and current jets was revealed; the correlation coefficient values are 0.6–0.97.

Filtering methods are usually used to combine the mean sea surface (MSS) and geoid (computable by global geopotential model (GGM)) into a common subspace, to model mean dynamic topography (MDT), which may lead to signal leakage and distortion problems. The use of the rigorous least squares (LS) method and multivariate objective analysis (MOA) alleviates these problems, and the derived MDTs from these two methods show better performance than MDTs derived from filtering methods. However, the advantages and disadvantages of these two methods have not been evaluated, and no direct comparison has yet been conducted between these two approaches regarding the performances in MDT recovery. In this study, we compare the performances of the MOA method with the LS method, providing information with respect to the usability of different methods in MDT modeling over regions with heterogeneous ocean states and hydrological conditions. We combined a recently published mean sea surface called DTU21MSS, and a satellite-only GGM named GO_CONS_GCF_2_DIR_R6, for MDT computation over four typical study areas. The results showed that the MDTs derived from the LS method outperformed the MOA method, especially over coastal regions and ocean current areas. The root mean square (RMS) of the discrepancies between the LS-derived MDT and the ocean reanalysis data was lower than the RMS of the discrepancies computed from the MOA method, by a magnitude of 1–2 cm. The formal error of the MDT estimated by the LS method was more reasonable than that derived from the MOA method. Moreover, the geostrophic velocities calculated by the LS-derived MDT were more consistent with buoy data than those calculated by the MOA-derived solution, by a magnitude of approximately 1 cm/s. The reason can be attributed to the fact that the LS method forms the design matrix segmentally, based on the error characteristics of the GGM, and suppresses high-frequency noise by applying constraints in different frequency bands, which improves the quality of the computed MDT. Our studies highlight the superiority of the LS-derived method versus the MOA method in MDT modeling. [ABSTRACT FROM AUTHOR]

The article presents the first evidence of a dual kinetic energy (KE) cascade in the ocean by analyzing in situ velocity measurements from surface drifters. Topics discussed include seasonal modulation of surface KE, seasonality of interscale KE transfers, and statistical metrics and error estimates.

In the Arabian Sea, southeast of the Arabian peninsula, an oceanic dipole, named the Ras Al Hadd (RAH) dipole, is formed each year, lying near the Ras Al Hadd cape. The RAH dipole is the association of a cyclonic eddy (CE) to the northeast, with an anticyclonic eddy (AE) to the southwest. This dipole intensifies in the summer monsoon and disappears during the winter monsoon. This dipole has been described previously, but mostly for its surface expression, and for short time intervals. Here, we describe the 3D structure of this dipole over the 2000–2015 period, by combining colocalized ARGO float profiler data (a total of 7552 profiles inside and outside the RAH dipole) with angular momentum eddy detection and tracking algorithm (AMEDA) surface data. We show first the different water masses in and near the RAH dipole. The presence of the Persian Gulf water (PGW) below 200 m depth is confirmed in both eddies. Arabian Sea high salinity water (ASHSW) is found exclusively in the AE; a layer of fresh and cold water is observed above 100 m depth in both eddies. By analyzing the potential density structures, we show that the CE has a surface-intensified structure while the AE is subsurface-intensified. The sea level anomaly shows a 0.04 m elevation above the AE and a 0.2 m depression over the CE. The CE has a faster geostrophic velocity, (vertical velocity, respectively) 0.6 m s−1 than the AE, 0.15 m s−1 (respectively, 3 m day−1 for the CE and 0.6 m day−1 for the AE). After presenting the vertical structure of the dipole, we show the dominance of the nonlinear Ekman pumping in the CE over the linear pumping affecting the dipole. As a consequence, we explain the CE's longer lifetime by its intensity and shallowness, and by its sensitivity to the interaction with the atmosphere (in particular the wind stress) and with neighboring eddies. We examined the possible (co)existence of symmetric, barotropic, and baroclinic instabilities in both eddies. These instabilities coexist near the surface in both eddies. They are intensified for the CE, which suggests that the CE is unstable and the AE is rather stable or may need a long time to be unstable. [ABSTRACT FROM AUTHOR]