Publications

In Preparation
Li, X. , A. Solodoch, E. Rahav, and Y. Lehahn. Wind-relative sampling orientation drives systematic bias in sea-surface microlayer enrichment estimates (In Preparation).
Wei, H. , A. Stewart, A. Medvedev, K. Srinivasan, A. Solodoch, G. Manucharyan, and McC. A. Hogg. Machine-learning-enabled satellite monitoring of ocean overturning circulation (In Preparation).
Solodoch, A. , H. Gildor, Y. Toledo, R. Barkan, V. Verma, Y. Fadida, and Y. Lehahn. Observations of a Deep-Reaching Submesoscale Eddy (In Preparation).
Solodoch, A. , A.L. Stewart, A. K Morrison, McC. A. Hogg, and A. N Garabato. Ocean Ventilation through the Antarctic Margins (In Preparation).
Submitted
Medvedev, A. , H. Wei, K. Armour, K. Srinivasan, A. Solodoch, A. Stewart, and . Basin-wide sea-surface observations reveal post-2000 emergence of AMOC weakening (Submitted). . Publisher's Version
2026
Obeid, J.P. , Y. Lehahn, A. Solodoch, E. Rahav, T. Treibitz, and A Shalabney. Optical Measurement of Oil Slick Thickness Using Active Interferometry. Frontiers in Marine Science 13 (2026). . Publisher's VersionAbstract

 

Spills of liquid petroleum hydrocarbons are a growing concern worldwide, posing great risks to marine life and community services. Identifying and treating oil spills is operationally and scientifically challenging and compounded by the difficulty in accurately obtaining real-time measurements of the oil thickness slicks. Here, we present a method that allows precise real-time measurement of oil slick thickness, based on active optical interferometry. A series of laboratory experiments with common hydrocarbon pollutant types, namely crude oil and gas condensate, showed that our method yields precise thickness measurements for slicks in the thickness range 0.382 - 23.3 (μm), with an accuracy of 95%. The proposed spectral-domain active interferometric system enables direct and physically grounded retrieval of oil film thickness without mechanical scanning and without reliance on ambient illumination. In principle, the system can be adapted for deployment at sea, opening the way for real-time, in situ thickness measurements that will improve oil-spill mitigation efforts and contribute to a deeper understanding of processes at the ocean-atmosphere interface.

 

obeid2026_oilthicknessinterferometry.pdf
Fadida, Y. , V. Verma, R. Barkan, Biton. E., A. Solodoch, and Y. Lehahn. Coastal-to-offshore submesoscale horizontal stirring enhances wintertime phytoplankton blooms in the ultra-oligotrophic Eastern Mediterranean Sea. Ocean Science 22, no. 1 (2026): 329-343. . Publisher's VersionAbstract

The large seasonal increases in marine photosynthetic organisms – i.e., phytoplankton blooms – are a ubiquitous oceanic phenomenon that contributes to the removal of carbon dioxide from the atmosphere and supports the growth of larger marine organisms. The underlying mechanisms controlling the intensity and timing of these blooms have been proposed to be dominated by vertical transport and mixing processes that are enhanced at fine-scale frontal and filamental circulations, commonly known as submesoscale currents. Here we show that the winter blooms characteristic of the ultra-oligotrophic waters of the Eastern Mediterranean Sea, which manifest as a seasonal increase in satellite-derived levels of surface chlorophyll, are intensified by enhanced horizontal stirring induced by the submesoscale currents. Using ocean color remote sensing data and high-resolution numerical simulations, we demonstrate that the intensification of submesoscale currents in winter efficiently connect the coastal waters and the ultra-oligotrophic open-sea waters, thereby enriching the latter with chlorophyll-rich waters. A climatological chlorophyll time series comparison between two different regions equidistant to the Nile River Delta indicates that this submesoscale horizontal stirring mechanism accounts for the ∼ 24.8 % larger wintertime increase in surface chlorophyll observed downstream of the Nile Delta. These results shed new light on the processes governing phytoplankton bloom intensity and emphasize the important role of submesoscale horizontal stirring in modulating the marine ecosystem.

2025
Ashkenazy, Y. , H. Gildor, and A. Solodoch. Statistical analysis of ocean currents in the Eastern Mediterranean. EGUsphere (2025). . Publisher's Version
H., Wei. , K. Srinivasan, A.L. Stewart, A. Solodoch, G.E. Manucharyan, and McC. A Hogg. Full-depth Reconstruction of Long-Term Meridional Overturning Circulation Variability from Satellite-Measurable Quantities via Machine Learning . Journal of Advances in Modeling Earth Systems 17, no. 7 (2025). . Publisher's VersionAbstract
The meridional overturning circulation (MOC) plays a crucial role in the global distribution of heat, carbon, and other climate‐relevant tracers. Monitoring the evolution of MOC is essential for understanding climate variability, yet direct MOC observations are sparse and geographically limited. Although satellite measurements have shown potential for short‐term monitoring of the MOC, it remains unclear whether MOC variability on decadal and longer timescales can be detected remotely. In this study, we leverage machine learning to reconstruct long‐term MOC variability from satellite‐measurable quantities, using climate simulations under pre‐industrial conditions. We demonstrate that our proposed non‐local dual‐branch neural network (DBNN) effectively reconstructs both the strength and vertical structure of the Atlantic MOC (AMOC) and the Southern Ocean MOCs across sub‐annual to multi‐decadal timescales. Using a neural network interpretation technique, we identify ocean bottom pressure near the western boundary and along dense‐water export pathways as the dominant input features for MOC reconstruction. This indicatesthat DBNN's predictions can be interpreted as an approximation of geostrophic balance. The DBNN also effectively reconstructs the AMOC in the equatorial region, where geostrophy breaks down. This success is attributed to the capability of DBNN in utilizing latitudinally non‐local ocean bottom pressure information and the meridional coherence of AMOC variability. Additionally, the DBNN accurately reconstructs Southern Ocean MOCs using only sea surface height and zonal wind stress as inputs, thereby avoiding reliance on ocean bottom pressure, which is subject to considerable measurement uncertainty in practice. This work demonstrates the possibility of continuous, long‐term MOC monitoring using satellite measurements
wei-2025-fulldepth-reconstruction-of-longter.pdf
Solodoch, A. , R. Barkan, and H. Gildor. The East Mediterranean boundary current and eddy spectrum. In Book chapter - in "Hod Hayam 2". (in Hebrew), 2025.
Solodoch, A. , Y. Toledo, V. Grigorieva, and Y. Lehahn. Retrieval of Surface Waves Spectrum from UAV Nadir Video. IEEE Transactions on Geoscience and Remote Sensing (2025). . Publisher's VersionAbstract
Sea surface wave spectrum measurements are necessary for a host of basic research questions as well as for engineering and societal needs. However, most measurement techniques require great investment in infrastructure and time-intensive deployment techniques. We propose a new approach of wave measurement from standard video footage recorded by low-cost Unmanned Aerial Vehicles (UAV). We address UAV nadir imagery, which are particularly simple to obtain operationally. The method relies on the fact that optical contrast of surface gravity waves is proportional to their steepness. We present a robust methodology of regularized inversion of the optical imagery spectra, resulting in retrieval of the three-dimensional wavenumber-frequency sea surface height spectrum. The system was tested in several sea trials and in different bathymetric depths and sea state conditions. The resulting wave bulk parameters and spectral characteristics are in good agreement with collocated measurements from wave buoys and bottom-mounted acoustic sensors. Simple deployment, mobility, and flexibility in spatial coverage show a great potential of UAVs to significantly enhance the availability of wave measurements.
2024
Jeffree, J. , A.M.C. Hogg, A.K Morrison, A. Solodoch, A. Stewart, and R. McGirr. GRACE Satellite Observations of Antarctic Bottom Water Transport Variability. Journal of Geophysical Research - Oceans 129 (2024). . Publisher's VersionAbstract

Antarctic Bottom Water (AABW) formation and transport constitute a key component of the global ocean circulation. Direct observations suggest that AABW volumes and transport rates may be decreasing, but these observations are too temporally or spatially sparse to determine the cause. To address this problem, we develop a new method to reconstruct AABW transport variability using data from the GRACE (Gravity Recovery and Climate Experiment) satellite mission. We use an ocean general circulation model to investigate the relationship between ocean bottom pressure and AABW: we calculate both of these quantities in the model, and link them using a regularized linear regression. Our reconstruction from modeled ocean bottom pressure can capture 65%–90% of modeled AABW transport variability, depending on the ocean basin. When realistic observational uncertainty values are added to the modeled ocean bottom pressure, the reconstruction can still capture 30%–80% of AABW transport variability. Using the same regression values, the reconstruction skill is within the same range in a second, independent, general circulation model. We conclude that our reconstruction method is not unique to the model in which it was developed and can be applied to GRACE satellite observations of ocean bottom pressure. These advances allow us to create the first global reconstruction of AABW transport variability over the satellite era. Our reconstruction provides information on the interannual variability of AABW transport, but more accurate observations are needed to discern AABW transport trends.

grace_aabw_paper.pdf
Malul, D. , H. Berman, A. Solodoch, Tal. O., N. Barak, G. Mizrahi, I. Berenshtein, et al. Directional swimming patterns in jellyfish aggregations. In revision.. Current Biology 34, no. 17 (2024): 4033-4038. . Publisher's VersionAbstract
Having a profound influence on marine and coastal environments worldwide, jellyfish hold significant scientific, economic, and public interest. The predictability of outbreaks and dispersion of jellyfish is limited by a fundamental gap in our understanding of their movement. Although there is evidence that jellyfish may actively affect their position, the role of active swimming in controlling jellyfish movement, and the characteristics of jellyfish swimming behavior, are not well understood. Consequently, jellyfish are often regarded as passively drifting or randomly moving organisms, both conceptually and in process studies. Here we show that the movement of jellyfish is controlled by distinctly directional swimming patterns, which are oriented against the direction of surface gravity waves. Taking a Lagrangian viewpoint from drone videos that allows the tracking of multiple adjacent jellyfish, and focusing the scyphozoan jellyfish Rhopilema nomadica as a model organism, we show that the behavior of individual jellyfish translates into a synchronized directional swimming of the aggregation as a whole. Numerical simulations show that this counter-wave swimming behavior results in biased correlated random-walk movement patterns that reduce the risk of stranding, thus providing jellyfish with an adaptive advantage critical to their survival. Our results emphasize the importance of active swimming in regulating jellyfish movement, and open the way for a more accurate representation in model studies, thus improving the predictability of jellyfish outbreaks and their dispersion, and contributing to our ability to mitigate their possible impact on coastal infrastructure and populations.
directionalswimminginjellyfish.pdf
Verma, V. , R. Barkan, A. Solodoch, H. Gildor, and Y. Toledo. The Eastern Mediterranean boundary current: seasonality, stability, and spiral formation. Journal of Physical Oceanography 54, no. 9 (2024): 1971-1989. 10.1175/JPO-D-23-0207.1. Publisher's VersionAbstract
Seasonal variability and the effect of bottom interaction on the dynamics of the along-slope boundary current flowing around the Levantine basin are investigated using nested high-resolution simulations of the Eastern Mediterranean Sea. The numerical solutions show a persistent boundary current year-round that is ≈ 60 km wide and ≈ 200 m deep. An enstrophy balance diagnostic reveals significant bottom-drag influence on the boundary current, leading to anticyclonic vorticity generation in thin regions along the coast, which in turn become unstable and roll into surface intensified anticyclonic spirals characterized by (1) Rossby numbers. An eddy kinetic energy generation analysis suggests that a mix of baroclinic and barotropic instabilities are likely responsible for the spiral formation. The boundary current and spirals play a crucial role in the cross-shore transport of materials. In winter, the anticyclonic spirals frequently interact and exchange material with the energetic offshore submesoscale flow field. In summer, when the offshore flow structures are relatively less energetic, the spirals remain confined to the boundary current region as they are advected by the boundary current and undergo an upscale kinetic energy (KE) cascade that is manifested in spiral merging, and growth up to 100 km in diameter. In both seasons, a coarse-graining analysis demonstrates that the cross-scale KE fluxes are spatially localized in coherent structures. The upscale KE fluxes typically occur within the spirals, while the downscale KE fluxes are confined to fronts and filaments at spiral peripheries.
Stewart, A.L. , Y. Wang, A. Solodoch, R. Chen, and J.C. McWilliams. Formation of eastern boundary undercurrents via mesoscale eddy rectification. Journal of Physical Oceanography (2024).Abstract

Eastern Boundary Upwelling Systems (EBUSs) host equatorward wind-driven near-surface currents overlying poleward subsurface undercurrents. Various previous theories for these undercurrents have emphasized the role of poleward alongshore pressure gradient forces (APF). Energetic mesoscale variability may also serve to accelerate undercurrents via mesoscale stirring of the potential vorticity gradient imposed by the continental slope. However, it remains unclear whether this eddy rectification mechanism contributes substantially to driving poleward undercurrents in EBUS. This study isolates the influence of eddy rectification on undercurrents via a suite of idealized simulations forced either by alongshore winds, with or without an APF, or by randomly-generated mesoscale eddies. It is found that the simulations develop undercurrents with strengths comparable to those found in nature in both wind-forced and randomly forced experiments. Analysis of the momentum budget reveals that the along-isobath undercurrent flow is accelerated by isopycnal advective eddy momentum fluxes and the APF, and retarded by frictional drag. The undercurrent acceleration may manifest as eddy momentum fluxes or as topographic form stress depending on the coordinate system used to compute the momentum budget, which reconciles these findings with previous work that linked eddy acceleration of the undercurrent to topographic form stress. The leading-order momentum balance motivates a scaling for the strength of the undercurrent that explains most of the variance across the simulations. These findings indicate that eddy rectification is of comparable importance to the APF in driving poleward undercurrents in EBUSs, and motivate further work to diagnose this effect in high-resolution models and observations, and to parameterize it in coarse-resolution ocean/climate models.

stewart-2024-formation-of-eastern-boundary-under.pdf
2023
Solodoch, Aviv , Roy Barkan, Vicky Verma, Hezi Gildor, Yaron Toledo, Pavel Khain, and Yoav Levi. Basin scale to Submesoscale variability of the East-Mediterranean Sea upper circulation. Journal of Physical Oceanography (2023).Abstract

 

The East Mediterranean Sea (EMS) circulation has previously been characterized as dominated by gyres, mesoscale eddies, and disjoint boundary currents. We develop nested high resolution EMS numerical simulations to examine circulation variability spectrum with emphasis on the yet unexplored EMS submesoscale. It is identified in model and altimetry data that there is a continuous cyclonic boundary current (BC) encircling the Levantine basin, rather than several disjoint currents. This EMS BC advects eddy chains downstream and is identified as a principle source of regional mesoscale and submesoscale variability. During the seasonal fall to winter mixed layer deepening, energetic submesoscale (O(10 km)) eddies, fronts, and filaments emerge throughout the basin, characterized by O(1) Rossby numbers. 

We identify the EMS submesoscale temporal scales as ~1-5 days through spatio-temoral spectral analysis. Mooring data confirms the EMS winter energization of submesoscales with the model-identified times scales. The submesoscale is associated with a ~ k-2 kinetic energy (KE) wavenumber (k) spectral slope, shallower than the quasigeostrophic-like  ~ k-3 slope diagnosed at summer. The shallowness of the winter spectrum is shown to be due to divergent subinertial motions, consistent with the Boyd 1992 theoretical model rather than surface quasigeostrophy. Using coarse graining, we diagnose a seasonal inverse (forward) KE cascade above (below) 30 km scales due to rotational (divergent) motions, and show that these commence after completion of the fall submesosacle energization. We also show via coarse graining that at scales larger than several 100 kms, the spectrum becomes near constant and a weak forward cascade occurs (from gyre scales) to mesoscales.

 

Stewart, Andrew L, Nicole K Neumann, and Aviv Solodoch. ``Eddy'' Saturation of the Antarctic Circumpolar Current by Standing Waves. Journal of Physical Oceanography 53 (2023): 1161–1181.Abstract
It is now well established that changes in the zonal wind stress over the ACC do not lead to changes in its baroclinicity nor baroclinic transport, a phenomenon referred to as “eddy saturation”. Previous studies provide contrasting dynamical mechanisms for this phenomenon: on one extreme, changes in the winds lead to changes the efficiency with which transient eddies transfer momentum to the sea floor; on the other, structural adjustments of the ACC’s standing meanders increase the efficiency of momentum transfer. In this study the authors investigate the relative importance of these mechanisms using an idealized, isopycnal channel model of the ACC. Via separate diagnoses of the model’s time-mean flow and eddy diffusivity, the authors decompose the model’s response to changes in wind stress into contributions from transient eddies and the mean flow. A key result is that holding the transient eddy diffusivity constant while varying the mean flow very closely compensates changes in the wind stress, whereas holding the mean flow constant and varying the eddy diffusivity does not. This implies that “eddy saturation” primarily occurs due to adjustments in the ACC’s standing waves/meanders, rather than due to adjustments of transient eddy behavior. The authors derive a quasi-geostrophic theory for ACC transport saturation by standing waves, in which the transient eddy diffusivity is held fixed, and thus provides dynamical insights into standing wave adjustment to wind changes. These findings imply that representing eddy saturation in global models requires adequate resolution of the ACC’s standing meanders, with wind-responsive parameterizations of the transient eddies being of secondary importance.
stewart2023_topoeddysatacc.pdf
Solodoch, Aviv , Andrew L Stewart, Andrew McC Hogg, and Georgy E Manucharyan. Machine Learning-derived inference of the Meridional Overturning Circulation drom satellite-observable variables in an ocean state estimate. Journal of Advances in Modeling Earth Systems 15 (2023): e2022MS003370.Abstract

The oceanic Meridional Overturning Circulation (MOC) plays a key role in the climate system, and monitoring its evolution is a scientific priority. Monitoring arrays have been established at several latitudes in the Atlantic Ocean, but other latitudes and oceans remain unmonitored for logistical reasons. This study explores the possibility of inferring the MOC from globally-available satellite measurements via machine learning (ML) tecchniques, using the ECCOV4 state estimate as a test bed. The methodological advantages of the present approach include the use purely of available satellite data, its applicability to multiple basins within a single ML framework, and the ML model simplicity (a feed-forward fully connected neural network with small number of neurons). The ML model exhibits high skill in MOC reconstruction in the Atlantic, Indo-Pacific, and Southern Oceans. The approach achieves a higher skill in predicting the model Southern Ocean abyssal MOC than has previously been achieved via a dynamically-based approach. The skill of the model is quantified as a function of latitude in each ocean basin, as well as a function of the time scale of MOC variability. We find that ocean bottom pressure generally has the highest reconstruction skill potential, followed by zonal wind stress. We addition ally test which combinations of variables are optimal. Furthermore, ML interpretabil ity techniques are used to show that high reconstruction skill in the Southern Ocean is mainly due to (NN processing of) bottom pressure variability at a few prominent bathy metric ridges. Finally, the potential for reconstructing MOC strength estimates from real satellite measurements is discussed.

moc_ml.pdf
2022
Solodoch, A , AL Stewart, McC A Hogg, AK Morrison, AE Kiss, AF Thompson, SG Purkey, and L Cimoli. How does Antarctic Bottom Water cross the Southern Ocean?. Geophysical Research Letters 49 (2022): e2021GL097211.Abstract
Antarctic Bottom Water (AABW), which fills the global ocean abyss, is derived from dense water that forms in several distinct Antarctic shelf regions. Previous modeling stud ies have reached conflicting conclusions regarding export pathways of AABW across the Southern Ocean and the degree to which AABW originating from distinct source regions are blended during their export. This study addresses these questions using passive tracer deployments in a 61-year global high-resolution (0.1 degree) ocean/sea-ice simulation. Two distinct export “conduits” are identified: Weddell Sea- and Prydz Bay-sourced AABW are blended together and exported mainly to the Atlantic and Indian Oceans, while Ross Sea and Adelie Land-sourced AABW are exported mainly to the Pacific Ocean. Northward transport of each tracer occurs almost exclusively (>90%) within a single conduit. These findings imply that regional changes in AABW production may impact the three-dimensional structure of the global overturning circulation.
aabw_tracerpathways.pdf
2021
Solodoch, Aviv , Andrew L Stewart, and James C McWilliams. Formation of anticyclones above topographic depressions. Journal of Physical Oceanography 51 (2021): 207–228.Abstract
Long-lived anticyclonic eddies (ACs) have been repeatedly observed over several North Atlantic basins characterized by bowl-like topographic depressions. Motivated by these previous findings, the authors conduct numerical simulations of the spin-down of eddies initialized in idealized topographic bowls. In experiments with 1 or 2 isopycnal layers, it is found that a bowl-trapped AC is an emergent circulation pattern under a wide range of parameters. The trapped AC, often formed by repeated mergers of ACs over the bowl interior, is characterized by anomalously low potential vorticity (PV). Several PV segregation mechanisms that can contribute to the AC formation are examined. In one-layer experiments, the dynamics of the AC are largely determined by a nonlinearity parameter () that quantifies the vorticity of the AC relative to the bowl’s topographic PV gradient. The AC is trapped in the bowl for low . 1, but for moderate values (0.5 . . 1) partial PV segregation allows the AC to reside at finite distances from the center of the bowl. For higher & 1, eddies freely cross the topography and the AC is not confined to the bowl. These regimes are characterized across a suite of model experiments using and a PV homogenization parameter. Two-layer experiments show that the trapped AC can be top- or bottom-intensified, as determined by the domain-mean initial vertical energy distribution. These findings contrast with previous theories of mesoscale turbulence over topography that predict the formation of a prograde slope current, but do not predict a trapped AC.
anticyclonesabovetopo.pdf