Search

Ocean eddies affect large‐scale circulation and induce a kinetic energy cascade through their non‐linear interactions. However, since global observations of eddy dynamics come from satellite altimetry maps that smooth eddies and distort their geometry, the strength of this cascade is underestimated. Here, we use deep learning to improve observational estimates of global surface geostrophic currents and explore the implications for the cascade. By synthesizing multi‐modal satellite observations of sea surface height (SSH) and temperature, we achieve up to a 30% improvement in spatial resolution over the community‐standard SSH product. This reveals numerous strongly interacting eddies that were previously obscured by smoothing. In many regions, these newly resolved eddies lead to nearly an order‐of‐magnitude increase in the upscale kinetic energy cascade that peaks in spring and is strong enough to drive the seasonality of large mesoscale eddies. Our study suggests that deep learning can be a powerful paradigm for satellite oceanography. Plain Language Summary: We developed a deep learning method to estimate global maps of surface ocean currents from satellite observations with significantly improved resolution and accuracy compared to existing methods. These maps dramatically improve our ability to observe eddy dynamics and the impact of eddies on the transfer of energy between scales in the ocean. Our study suggests that deep learning can be a powerful paradigm for satellite oceanography. Key Points: We develop the first deep learning global estimates of surface ocean currents from multi‐modal satellite observationsOur deep learning method is able to map surface currents with state‐of‐the‐art resolution and accuracyThe diagnosed kinetic energy cascade is an order of magnitude higher compared to conventional altimetry products [ABSTRACT FROM AUTHOR]

Harnessing ocean kinetic energy has emerged as a promising renewable energy solution in recent years. However, identifying optimal locations for extracting this energy remains a significant challenge. This study presents a novel scheme to estimate the total surface current (TSC) as permanent surface current by integrating geodetic data and in-situ measurements. The TSC is typically a combination of the geostrophic current, derived from dynamic topography, and the Ekman current. We utilize NOAA's Ekman current data to complement the geostrophic current and obtain the TSC. To further enhance the accuracy of the TSC estimates, we employ a 3DVAR data assimilation method, incorporating local current meter observations. The results are verified against two control current meter stations. The data-assimilation process resulted in an improvement of 4 to 15 cm/s in the precision of calculated TSC. Using the assimilated TSC data, we then assess the kinetic energy potential and identify six regions with the most significant promise for marine kinetic energy extraction. This innovative approach can assist researchers and policymakers in targeting the most suitable locations for harnessing renewable ocean energy. [ABSTRACT FROM AUTHOR]

The recurring occurrence of green tides as an ecological disaster has been reported annually in the Yellow Sea. While remote sensing technology effectively tracks the scale, extent, and duration of green tide outbreaks, there is limited research on the underlying driving mechanisms of green tide drift transport and the determination of the leeway coefficient. This study investigates the green tide transport mechanism and evaluates the feasibility of estimating the leeway coefficient by analyzing green tide drift velocities obtained from Geostationary Ocean Color Imager-II (GOCI-II) images using the maximum cross-correlation (MCC) technique and leeway method across various time intervals alongside ocean current and wind speed data. The results reveal the following: (1) Significant spatial variations in green tide movement, with a distinct boundary at 34°40′N. (2) Short-term green tide transport is primarily influenced by tidal forces, while wind and ocean currents, especially the combined Ekman and geostrophic current component, predominantly govern net transport. (3) Compared to 1, 3, and 7 h intervals, estimating the leeway coefficient with a 25 h interval is feasible for moderate-resolution geostationary images, yielding values consistent with previous studies. This study offers new insights into exploring the transport mechanisms of green tides through remote sensing-driven velocity. [ABSTRACT FROM AUTHOR]

We analyze the feasibility of so-called radial surface currents as derived from Sentinel-1 Wave mode Doppler shift observations to support the geodetic estimation of the mean dynamic topography (MDT). Broadly defined as subtracting the geoid from the mean sea surface, this separation problem suffers from data inhomogeneities and usually requires strong prior assumptions about the MDT's smoothness. Recent advancements in calibration and current retrieval methodology for Sentinel-1 data give reason to assess potential gains from this observation type, which can be linked to the gradient of the MDT under geostrophic approximation. Due to current lack of long time-series, we synthesize 10 years of observations from daily surface currents grids and include these in a parametric least-squares adjustment of the MDT. Our results utilizing the synthetic data in place of smoothness constraints are promising, showing 2 cm RMS agreement with a state-of-the-art MDT model. [ABSTRACT FROM AUTHOR]

By decomposing the total sea surface current into its geostrophic and ageostrophic components, we examined the contribution of each to the long-term variability of the total sea surface current. Our findings demonstrate the importance of geostrophic currents in Baltic Sea gyre formations. Additionally, ageostrophic currents contribute significantly to the flow across the region. Quantifying the difference between total sea surface current fields has revealed two dominant general sea surface circulation patterns in the Baltic Sea, whose characteristics correspond to the monthly mean climatology of sea surface current fields in May and December. Subsequently, a machine learning technique was employed to effectively detect the types of sea surface circulation patterns using wind vectors and sea level anomaly fields. This underscored the combined influence of sea level anomaly-driven and wind-driven components in shaping surface current vectors in the Baltic Sea, consistent with geostrophic and ageostrophic decompositions. [ABSTRACT FROM AUTHOR]

In order to address the issue of limited common coverage and high cost in mapping ocean surface vector current by two (or more) high-frequency surface wave radars, this paper proposes a single-station surface wave radar vector current inversion algorithm. The feasibility of this algorithm has been validated in the Taiwan Strait. Based on the ocean dynamic characteristics of the Taiwan Strait, the algorithm utilizes the radial current obtained from a high-frequency surface wave radar (HFSWR) in Fujian Province to invert the ocean surface vector current. The surface vector current can be decomposed into three primary dynamic components: tidal currents, wind-driven currents, and geostrophic currents. Firstly, tidal current forecasting models and Ekman and Stokes theories are used to calculate the tidal and wind-driven currents in the Taiwan Strait, respectively. Subsequently, the directions of geostrophic currents in the Taiwan Strait are determined with sea surface height data, and the magnitudes of the geostrophic currents are constrained using the radial current from the single HFSWR. Finally, the three components are added together to obtain the vector current. Comparative results demonstrate that the efficacy of the algorithm has been validated through field experiments (with two HFSWRs and two drifting buoys) conducted in the southwestern of the Taiwan Strait. Further research is needed on the applicability of this algorithm to other sea areas and monitoring systems. [ABSTRACT FROM AUTHOR]

Exceptional precision and excellent resolution reconstruction of sea surface currents are beneficial for exploring complex oceanic dynamic processes. Normally, this required physical inversion models for global or regional oceans are constructed to reconstruct oceanic currents. These models are based on the analysis of sea surface geostrophic and Ekman currents derived from satellite observations of sea level and wind stress fields. Nevertheless, the presence of various typical dynamic processes in marine environments, such as mesoscale eddies and small-scale waves, continues to pose challenges in accurately reconstructing oceanic currents. Meanwhile, any product of surface current that neglects the contribution of wave motion would, at best, be incomplete. Therefore, in this paper, we introduce an accurate sea surface current product at a depth of 15 m, named GEST (Geostrophic-Ekman-Stokes-Tide). This product is produced by a multi-scale dynamics-informed neural network that learns the intricate representation of concealed characteristics in Ekman, geostrophic currents, wave-induced Stokes drift, and TPXO9 tidal currents. Its structure design is predicated upon the intricate coupling relationships between various ocean surface components and the veritable currents discerned by the deployment of drift buoys, with each ocean surface component correlating to discrete physical processes. Compared with the prevailing product, the GEST confers an elevation in precision by approximately 9.2 cm/s over the traditional multinomial fitting method, 10.4 cm/s beyond the OSCAR, and 8.81 cm/s surpassing GlobCurrent. [ABSTRACT FROM AUTHOR]

Satellite altimetry offers a unique approach for direct sea surface current observation, but it is limited to measuring the surface‐constrained geostrophic component. Ageostrophic dynamics, prevalent at horizontal scales below 100 km and time scales below 10 days, are often underestimated by ocean reanalyzes employing data assimilation schemes. To address this limitation, we introduce a novel deep learning scheme, rooted in a variational data assimilation formulation with trainable observations and a priori terms, that harnesses the synergies between satellite‐derived sea surface observations, namely sea surface height (SSH) and sea surface temperature (SST), to enhance sea surface current reconstruction. Numerical experiments, conducted using realistic simulations, in a case study area of the Gulf Stream, demonstrate the potential of the proposed scheme to capture ageostrophic dynamics at time scales of 2.5–3.0 days and horizontal scales of 0.5°–0.7°. The analysis of diverse observation configurations, encompassing nadir along‐track altimetry, wide‐swath SWOT (Surface Water and Ocean Topography) altimetry, and SST data, highlights the pivotal role of SST features in retrieving a significant portion of the ageostrophic dynamics (approximately 47%). These findings underscore the potential of deep learning and 4DVarNet schemes in improving ocean reanalyzes and enhancing our understanding of ocean dynamics. Plain Language Summary: Satellite altimetry provides a unique means for direct observation of sea surface currents, but it is confined to the geostrophic component, limiting the recovery of a substantial portion of mesoscale sea surface currents in operational products. To address this limitation, we present a novel deep learning framework, rooted in a variational data assimilation paradigm, that unlocks new avenues for leveraging the synergistic relationships between satellite‐derived sea surface observations, namely sea surface height and sea surface temperature. This innovative scheme demonstrates its remarkable potential to enhance sea surface current reconstruction and recover a substantial portion of the elusive ageostrophic dynamics. Numerical experiments, employing realistic simulations, in a case study area along the Gulf Stream, underscore the efficacy of our proposed approach. These findings support the pivotal role of physics‐informed deep learning in maximizing the utilization of available multimodal observation data sets and numerical simulations to elucidate partially observed sea surface dynamics. Key Points: We present end‐to‐end deep learning schemes to improve the reconstruction of total sea surface currents from satellite‐derived observationsExperiments in a region of the Gulf Stream support the synergistic analysis of sea surface temperature and sea surface height dataThe strain of sea surface dynamics is a proxy of the uncertainty of the retrieved estimation [ABSTRACT FROM AUTHOR]

For the past 20 years, Global Navigation Satellite System reflectometry (GNSS-R) technology has successfully shown its potential for remote sensing of the Earth's surface, including ocean and land surfaces. It is a multistatic radar that uses the GNSS signals reflected from the Earth's surface to extract land and ocean characteristics. Because of its numerous advantages such as low cost, multiple signal sources, and all-day/weather and high-spatiotemporal-resolution observations, this new technology has attracted the attention of many researchers. One of its most promising applications is GNSS-R ocean altimetry, which can complement existing techniques such as tide gauging and radar satellite altimetry. Since this technology for ocean altimetry was first proposed in 1993, increasing progress has been made including diverse methods for processing reflected signals (such as GNSS interferometric reflectometry, conventional GNSS-R, and interferometric GNSS-R), different instruments (such as an RHCP antenna with one geodetic receiver, a linearly polarized antenna, and a system of simultaneously used RHCP and LHCP antennas with a dedicated receiver), and different platform applications (such as ground-based, air-borne, or space-borne). The development of multi-mode and multi-frequency GNSS, especially for constructing the Chinese BeiDou Global Navigation Satellite System (BDS-3), has enabled more free signals to be used to further promote GNSS-R applications. The GNSS has evolved from its initial use of GPS L1 and L2 signals to include other GNSS bands and multi-GNSS signals. Using more advanced, multi-frequency, and multi-mode signals will bring new opportunities to develop GNSS-R technology. In this paper, studies of GNSS-R altimetry are reviewed from four perspectives: (1) classifications according to different data processing methods, (2) different platforms, (3) development of different receivers, and (4) our work. We overview the current status of GNSS-R altimetry and describe its fundamental principles, experiments, recent applications to ocean altimetry, and future directions. [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 coastal ocean off western Patagonia is one of the main coastal regions with high freshwater inputs from rivers, rain, and glaciers in the Southern Hemisphere. This study conducts an analysis of the seasonal and interannual variations in sea surface salinity and meridional geostrophic transports, specifically focusing on the Cape Horn Current, using improved satellite-derived data of sea surface salinity (SSS) and geostrophic velocities spanning an ∼11-year period (September 2011–August 2022). Our results reveal a clear salinity minimum in a coastal band between 42–54°S associated with the highest freshwater content. The average geostrophic currents are stronger south of 49°S, in line with the location of the Cape Horn Current. The average salinity minimum tends to disappear south of 54°S, with salinity values increasing slightly southward. The seasonal cycle of salinity shows the most pronounced minimum in summer (∼33.2–33.4). The greatest variability in salinity (standard deviation of salinity fields) occurs in the southern region of the Cape Horn Current. Hovmöller plots reveal two cores of minimum salinity observed in spring and summer (∼33.3–33.4). The freshwater off the Gulf of Penas contributes to the northern core. The meridional geostrophic transport differs between the northern and southern sections, with transports predominantly towards the Equator (Pole) north (south) of about 47–48°S during spring–summer. There is a marked seasonal variability in the magnitude and northern limit of the southward-flowing Cape Horn Current, being extended further north during winter and with a maximum average magnitude during summer–fall (about − 2 × 10 4 m2 s−1). On the interannual scale, a major drop in surface salinity occurred off northern and central Patagonia during 2018–2019. Finally, a potential long-term freshening trend is observed in the coastal area off southern Patagonia (south of 52°S), although prolonged data records are essential to confirm this pattern. [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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

• 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]

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]

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 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]

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]

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 & multivariate mapping approach that offers the possibility to improve the physical content of gridded products by combining the data from various platforms and in 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 2016-07-01 to 2020-06-30. 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 against the DUACS operational product distributed in the Copernicus Marine Service. We show that the multiscale & multivariate mapping approach offers promising perspectives for reconstructing the ocean surface circulation: leads observations contribute to improve the coverage in delivering gap free maps in the Arctic; drifters observations 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 effective resolution of this new product is hence between 5% and 10% finer than the Copernicus product. [ABSTRACT FROM AUTHOR]

Six oceanographic cruises in a NE-SW transect were made nearshore of southern Sinaloa and Nayarit from March 2006 through May 2008, where no in situ hydrographic data are available. Applying the Thermodynamic Equation of Seawater 2010 (TEOS-10) to the observations, the hydrography and geostrophic currents of the region were characterized. Results indicate that surface variability (0-50 m) emerged mainly from seasonal atmospheric forcing. A relative salinity maximum was present during all cruises below this surface layer, which is attributed to a water mass intrusion of Subtropical Subsurface Water that could be associated with the Mexican Coastal Current. Another water mass intrusion is from the California Current. Samples from the 2007-2008 La Niña produced an uncommon circulation, where water flowing from the Gulf of California along the coast of Sinaloa was observed, opposite to what is commonly known as a mean circulation. This uncommon circulation matches the generation of anticyclonic eddies around the Islas Marias archipelago. [ABSTRACT FROM AUTHOR]

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

Ocean spiciness is referred to as density-compensated temperature and salinity variations with warm (cool) and salty (fresh) waters having high (low) spiciness, respectively. The structure and evolution of density-compensated (warm/salty or cool/fresh) spiciness anomalies are investigated in the North Pacific thermocline using Argo observations during the period 2004–20. Two well-organized decadal spiciness events are identified through isopycnal surface analyses. One warm/salty spiciness anomaly of about 0.15°C and 0.05 g kg−1 temperature and salinity perturbations on the 25 kg m−3 isopycnal surface appeared in the eastern subtropical North Pacific at (18°–30°N, 120°–150°W) in 2007, which then migrated southwestward along the mean circulation and arrived in the western tropics at (∼15°N, 145°E–175°W) in 2012–13, with the reduced salinity perturbation of about 0.043 g kg−1. Another cool/fresh spiciness anomaly of about −0.2°C and −0.07 g kg−1 temperature and salinity perturbations originated from the eastern subtropics at (120°–150°W, ∼30°N) in 2014 and followed a similar advective pathway during the period from 2014–15 to 2019–20. The subduction pathway can be adequately determined by the mean Montgomery stream function on the isopycnal surface; the propagation direction and speed are in good agreement with the expectation for the role played by advection due to the mean geostrophic current. Moreover, the subducted decadal spiciness anomalies can have negative feedback on sea surface temperature (SST) variability in the western tropical Pacific through the diapycnal processes. The identifications of these density-compensated spiciness anomalies and their propagation pathways provide a clear illustration of the oceanic extratropics-tropics interactions in the North Pacific Ocean. [ABSTRACT FROM AUTHOR]

This data set was established from two open sources: (a) geoid undulation (N) from EIGEN‐6C4, which is a static global combined gravity field model up to degree and order 2,190, and (b) 'World ocean geostrophic velocity inverted from World Ocean Atlas 2013 with the P‐vector method' (NCEI accession 0,121,576). With the given non‐positive values of N in the oceans, absolute geostrophic currents (u, v) are easily obtained on N with 1°×1° resolution from the second data set except the equatorial zone (5oS – 5oN) due to the non‐existence of the geostrophic balance. Altogether, the data set contains 15,868 (u, v) data pairs. The data set shows that the hypothetical situation that there would be no flow on the geoid does not exist. [ABSTRACT FROM AUTHOR]