It is generally agreed that the resolution of a regular quadrilateral mesh is the side length of quadrilateral cells. There is less agreement on what is the resolution of triangular meshes, exacerbated by the fact that the numbers of edges or cells on triangular meshes are approximately three or two times larger than that of vertices. However, the geometrical resolution of triangular meshes, i.e. maximum wavenumbers that can be represented on such meshes, is a well defined quantity, known from solid state physics. These wavenumbers are related to a smallest common mesh cell (primitive unit cell), and the set of mesh translations that map it into itself. The wavenumbers do not depend on whether discrete degrees of freedom are placed on vertices, cells or edges. The resolution is defined by the height of triangles.

Preservation of organic carbon (OC) in marine and terrestrial deposits is enhanced by bonding with reactive iron (FeR) phases. The association of OC with FeR (OC-FeR) provides physical protection and hinders microbiological degradation. Roughly 20% of all OC stored in unconsolidated marine sediments and 40% of all OC present in Quaternary terrestrial deposits is preserved as OC-FeR, but this value varies from 10 to 80% across depositional environments. In this work, we provide a new assessment of global OC-FeR burial rates in both marine and terrestrial environments, using published estimates of the fraction of OC associated with FeR, carbon burial, and probabilistic modelling. We estimate the marine OC-FeR sink at between 31 – 70 Mt C yr-1 (mean 52 Mt C yr-1), and the terrestrial OC-FeR sink at between 171 - 946 Mt C yr-1 (mean 472 Mt C yr-1). In marine environments, continental shelves (mean 17 Mt C yr-1) and deltaic/estuarine environments (mean 11 Mg C yr-1) are the primary locations of OC-FeR burial. On land, croplands (279 Mt C yr-1) and grasslands (121 Mt C yr-1) dominate the OC-FeR burial budget. Changes in the Earth system through geological time likely alter the OC-FeR pools, particularly in marine locations. For example, periods of intense explosive volcanism may lead to increased net OC-FeR burial in marine sediments. Our work highlights the importance of OC-FeR in marine carbon burial and demonstrates how OC-FeR burial rates may be an order of magnitude greater in terrestrial environments, those potentially most sensitive to anthropogenic impacts.

Ocean circulation is dominated by turbulent geostrophic eddy fields with typical scales ranging from 10 km to 300 km. At mesoscales (> 50 km), the size of eddy structures varies regionally following the Rossby radius of deformation. The variability of the scale of smaller eddies is not well known due to the limitations in existing numerical simulations and satellite capability. But it is well established that oceanic flows (< 50km) generally exhibit strong seasonality. In this study, we present a basin-scale analysis of coherent structures down to 10\,km in the North Atlantic Ocean using two submesoscale-permitting ocean models, a NEMO-based North Atlantic simulation with a horizontal resolution of 1/60 (NATL60) and an HYCOM-based Atlantic simulation with a horizontal resolution of 1/50 (HYCOM50). We investigate the spatial and temporal variability of the scale of eddy structures with a particular focus on eddies with scales of 10 to 100\,km, and examine the impact of the seasonality of submesoscale energy on the seasonality and distribution of coherent structures in the North Atlantic. Our results show an overall good agreement between the two models in terms of surface wavenumber spectra and seasonal variability. The key findings of the paper are that (i) the mean size of ocean eddies show strong seasonality; (ii) this seasonality is associated with an increased population of submesoscale eddies (10\,–\,50\,km) in winter; and (iii) the net release of available potential energy associated with mixed layer instability is responsible for the emergence of the increased population of submesoscale eddies in wintertime.

The western tropical Pacific (WTP) exhibits large interannual sea level anomalies (SLAs), and the sea level falling in El Niño is evidently stronger than the rising in La Niña. The asymmetry is most prominent near 160°E with the response to El Niño larger by three times and becomes less obvious near the western boundary. Sensitivity experiments of a simplified ocean model suggest that the asymmetry in surface wind forcing structure between El Niño and La Niña is critical. The El Niño’s westerly wind anomaly patch locates more east than the La Niña’s easterly wind patch during the mature stage, and its upwelling effects are accumulated over a wider longitude range and cause stronger negative SLAs in the WTP. Near the western boundary, however, upwelling effects are attenuated by easterly wind anomalies during El Niño conditions. The asymmetric ocean responses to ENSO winds may participate in the asymmetry of ENSO cycle.

The traditional ocean color remote sensing usually focuses on using optical inversion models to estimate the properties of in-water components from the above-surface spectra, so we call it the spectrum-concentration (SC) scheme. Unlike the SC scheme, this study proposed a new research scheme, distribution-distribution (DD) scheme, which uses statistical inference models to estimate the possibility distribution of these in-water components, based on the possibility distribution of the observed spectra. The DD scheme has the advantages that (1) it can rapidly give the key and overview information of the interest water, instead of using the SC scheme to compute each image pixel, (2) it can assist the SC scheme to improve their models and parameters, and (3) it can provide more valuable information for better understanding and indicating the features and dynamics of aquatic environment. In this study, based on Landsat-8 images, we analyzed the spectral possibility distributions (SPD) of 688 global water and found many of them were normal, lognormal, and exponential distributions, but with diverse patterns in distribution parameters such as the mean, standard deviation, skewness and kurtosis. Furthermore, we used Monte-Carlo and Hydrolight simulations to study the theoretical and statistical connections between the possibility distributions of in-water components and SPDs. The simulation results were basically consistent with the observations on the real water. Then by using the simulation and field measured data, we proposed a bootstrap-based DD scheme and developed some simple statistical inference models to estimate the distribution parameters of yellow substance in lakes. Since DD scheme is still on its early stage, we also suggested some potential and useful topics for the future work.

One third of all coastlines worldwide consist of permafrost. Many of these permafrost coasts are presently exposed to greater environmental forcing as a consequence of climate change, such as a lengthening of the open water season, intensified storms, and higher water and air temperatures. As a result, increasing erosion rates are currently reported from various sites across the Arctic. It is crucial to synthetize these data on Arctic shoreline dynamics in order to improve our understanding on present coastal dynamics on the pan-Arctic scale. A first synthesis product was released in form of the Arctic Coastal Dynamics databse in 2012, which included data published until 2009 (Lantuit et al., 2012). Since then, numerous publications and data products were published on short and long term changes of Arctic coasts across a wide range of study sites. We made an extensive literature review of publications released within the last 10 years and updated the shoreline change data section in the Arctic Coastal Dynamics database. While in 2009 for one percent of the Arctic shoreline data on coastal dynamics was available, the addition of new data leads to a broader data coverage, which is mainly the effect of the greater availability of remotely sensed products for analyses conducted in these remote regions. Further, the additional data allow us to update the current mean rate of Arctic shoreline change.

Some of the Earth system data products such as those from NASA airborne and field investigations (a.k.a. campaigns), are highly heterogeneous and cross-disciplinary, making the data extremely challenging to manage. For example, airborne and field campaign measurements tend to be sporadic over a period of time, with large gaps. Data products generated are of various processing levels and utilized for a wide range of inter- and cross-disciplinary research and applications. Data and derived products have been historically stored in a variety of domain-specific standard (and some non-standard) formats and in various locations such as NASA Distributed Active Archive Centers (DAACs), NASA airborne science facilities, field archives, or even individual scientists’ computer hard drives. As a result, airborne and field campaign data products have often been managed and represented differently, making it onerous for data users to find, access, and utilize campaign data. Some difficulties in discovering and accessing the campaign data originate from the incomplete data product and contextual metadata that may contain details relevant to the campaign (e.g. campaign acronym and instrument deployment locations), but tend to lack other significant information needed to understand conditions surrounding the data. Such details can be burdensome to locate after the conclusion of a campaign. Utilizing consistent terminology, essential for improved discovery and reuse, is also challenging due to the variety of involved disciplines. To help address the aforementioned challenges faced by many repositories and data managers handling airborne and field data, this presentation will describe stewardship practices developed by the Airborne Data Management Group (ADMG) within the Interagency Implementation and Advanced Concepts Team (IMPACT) under the NASA’s Earth Science Data systems (ESDS) Program.

Mangrove forests with complex root systems contribute to increased coastal protection through drag effects. Previous flume studies proposed a predictive model of drag in Rhizophora mangrove forests based on quadratic drag law. However, its general applicability on mangrove forests in the field has not been tested. To fill this knowledge gap, this study quantified drag in a 17-year-old planted Rhizophora mangrove forest using a comprehensive measurement of hydrodynamics and vegetation morphology. The vegetation projected area density, a, showed an approximate exponential increase towards the bed, mainly due to root branching. This vertical variation led to enhanced vegetation drag per unit water volume relative to velocity with decreasing water depth. Alternatively, the drag per vegetation projected area solely depended on the square of velocity, indicating association with the quadratic drag law. The derived drag coefficient (CD) was 1.0 ± 0.2 for tide-driven currents, consistent with previous flume studies. By using the mean value of derived CD (1.0), it was confirmed that the quadratic drag model expresses well the field-measured drag. We also presented a method for predicting a value for a, another unknown parameter in the drag model, using an empirical Rhizophora root model, and confirmed a successful prediction of a and drag. Therefore, the drag in a Rhizophora mangrove forest can be accurately predicted only using the input parameters of the Rhizophora root model – stem diameter and tree density. This provides insights into effectively implementing the drag model in hydrodynamic models for better representation of mangroves’ coastal protection function.

Time series of shipboard observations in the southern Arafura Sea near the Tiwi Islands indicated that the water column dynamics differed between the east and west sides of the islands. On the west side, the water column, characterized by temperature, salinity, and velocity, was barotropic and tidal advection dominated. On the east side, the water column was baroclinic and internal tides were present along with tidal advection. These conditions affected the distribution of the turbidity and fluorescence in the water column. Likewise, the influence of the daily solar radiation cycle reached the bottom on the western side, but was limited to the upper layer above the thermocline on the eastern side. The fluorescence peaks also differed between the east and west sides, with the eastern side dominated by the semidiurnal tides and the western side by the daily solar cycle. Fluorescence integrated over the water column was much higher on the eastern side than the western side. Also on the eastern side, fluorescence was limited to the lower layer, while on the western side, it encompassed the entire water column at times and peaked below the warmer, higher oxygenated water generated by solar radiation and surface mixing. These dynamics have distinct implications for biological productivity and also may affect a proposed tidal power system in the region.

This study investigates trends in global tropical cyclone (TC) activity from 1990–2020, a period where observational platforms are mostly consistent. Several global TC metrics have decreased during this period, with significant decreases in hurricanes and Accumulated Cyclone Energy (ACE). Most of this decrease has been driven by significant downward trends in the western North Pacific. Globally, short-lived named storms, 24-hr intensification periods of >=50 kt day-1 and TC-related damage have increased significantly. The increase in short-lived named storms is likely due to technological improvements, while rapidly intensifying TC increases may be fueled by higher potential intensity. Damage increases are largely due to increased coastal assets. The decreasing trends in hurricane numbers and global ACE are likely due to the trend towards a more La Niña-like base state from 1990–2020, favoring TC activity in the North Atlantic and suppressing TC activity in the eastern and western North Pacific.

Eastern Boundary Currents Systems are typically studied as a whole due to their dynamical similarities, mainly because Ekman pumping is predominant at these currents, and they typically have low kinetic energy. In this study, we used the output of a high-resolution global simulation to make a dynamical comparison among the California, Canary, Peru, and Benguela currents during the winter and summer months, focusing on submesoscale motions (Ro ~ 1) in both the frequency-wavenumber and space-time domains. After we confirmed the presence of submesoscale activity and isolated it from mesoscale motions, we found that their divergence and vorticity fields follow similar seasonal patterns in the near-diurnal frequency range, despite regional differences. The results showed that heat fluxes at the ocean surface, along with weak to moderate wind stresses, significantly impact the modulation of submesoscale vorticity and divergence fields at diurnal frequencies.

Understanding and predicting sea ice dynamics and ice-ocean feedback processes requires accurate descriptions of momentum fluxes across the ice-ocean interface. In this study, we present observations from an array of moorings in the Beaufort Sea. Using a force-balance approach, we determine ice-ocean drag coefficient values over an annual cycle and a range of ice conditions. Statistics from high resolution ice draft measurements are used to calculate expected drag coefficient values from morphology-based parameterization schemes. With both approaches, drag coefficient values ranged from approximately 1-10×10^-3, with a minimum in fall and a maximum at the end of spring, consistent with previous observations. The parameterizations do a reasonable job of predicting the observed drag values if the under ice geometry is known, and reveal that keel drag is the primary contributor to the total ice-ocean drag coefficient. When translations of bulk model outputs to ice geometry are included in the parameterizations, they overpredict drag on floe edges, leading to the inverted seasonal cycle seen in prior models. Using these results to investigate the efficiency of total momentum flux across the atmosphere-ice-ocean interface suggests an inter-annual trend of increasing coupling between the atmosphere and the ocean.

The Reykjanes Ridge is a major topographic feature that lies south of Iceland in the North- Atlantic Ocean and strongly influences the Subpolar Gyre (SPG) circulation. Based on velocity and hydrographic measurements carried out along the crest of the Ridge from the Icelandic continental shelf to 50°N during the RREX cruise in June-July 2015, we derived the first direct estimates of volume and water masses transports over the Ridge. The circulation was mainly westward north of 53.35°N and eastward south of it. The westward transport was estimated at 21.9 ± 2.5 Sv (Sv = 10 6 m 3 s -1 ) and represents the SPG intensity. The westward flows followed two main pathways at 57°N near Bight Fracture Zone and at 59 – 62°N. We argue that those pathways were respectively connected to the northern branch of the North Atlantic Current and to the Sub-Arctic Front that were intersected by the southern part of the section. In addition to this horizontal circulation, mixing and bathymetry shaped the water mass distribution. Water mass transformations in the Iceland Basin lead to the formation of weakly stratified SubPolar Mode Water (SPMW). We explain why SPMW, which was the main contributor in terms of water mass to the westward flow, was denser at 57°N than at 59–62°N along the Ridge. At higher densities, both Intermediate Water, defined by a dissolved oxygen minimum, and Icelandic Slope Water contributed as much to the westward transport across the Ridge as the sum of Labrador Sea Water and Iceland-Scotland Overflow Water.

High performance liquid chromatography (HPLC) remains one of the most widely applied methods currently available for estimation of phytoplankton taxonomy from ocean samples. This method measures the concentrations of phytoplankton pigments, some of which are useful chemotaxonomic markers that can be used to diagnose the relative abundance of phytoplankton groups. Here, we use HPLC phytoplankton pigment concentrations measured on surface water samples from 38 field surveys for a total of over 3,000 distinct samples that cover every major ocean basin and represent a wide range of ecological regimes. The data compilation has been quality controlled to remove measurements below pigment detection limits and outliers from the linear regression of total chlorophyll-a concentration with total accessory pigment concentrations and only samples from labs that have participated in round-robin quality assurance experiments (e.g. NASA SeaHARRE) have been included. We assess the environmental and spatial drivers controlling the global distribution and co-variability of individual phytoplankton pigments. Preliminary results of hierarchical clustering show strong differentiation in phytoplankton pigments following known relationships between phytoplankton size class and relative pigment concentration, partitioning their contributions by micro-, nano-, and pico-phytoplankton size classes. However, the exact clusters relationships change when the data are divided by ocean basin or latitude. We also use statistical techniques, including EOFs and network-based exploration, to examine the associations between groups of pigments over a range of environmental conditions on local to global scales and diagnose the main controls on these associations.

Natural methane gas release from the seafloor is a widespread phenomenon that occurs at cold seeps along most continental margins. Since their discovery in the early 1980s, seeps have been the focus of intensive research, partly aimed to refine the global carbon budget. However, deep-sea research is challenging and expensive and, to date, few works have successfully monitored the variability of methane gas release over long time periods (> 1 yr). Long-term monitoring is necessary to study the mechanisms that control seabed gas release. The M³ project, funded by the German Ministry of Education and Research, aims to study the temporal and spatial variability of gas emissions at the Southern Hydrate Ridge (SHR) by acoustically monitoring and quantifying gas effluxes over several years. Located 850 m deep on the Cascadia accretionary prism offshore Oregon, the SHR is one of the most studied seep sites and persistent but variable gas release has been observed for more than 20 years. Since 2015, the Ocean Observatories Initiative’s (OOI) Cabled Array observatory, provides power supply and two-way communication to the SHR, making it an ideal site for continuous long-term monitoring work. In this work, we present how we will take advantage of the OOI infrastructure and deploy several instruments on the seabed for at least 1.5 year. A multi-beam “overview” sonar mounted on a rotor will identify every gas bubble stream located within 200 m from the sonar location. A scanning “quantification” sonar will be used to estimate the amount of gas that is released from discrete gas streams. A camera system and a CTD probe will help process and analyze the hydro-acoustic data. All instruments will be powered and controlled from land through the OOI infrastructure. We present the instrument design, the operation protocol, as well as the data processing steps and expected results.

The Amazon River has the largest volume on earth, making up 15–20% of the annual fluvial discharge into oceans. The neighboring Pará River mixes with the Amazon River waters in the Amazon Estuary before forming a plume that extends into the Atlantic. Despite the global importance of these rivers, dissolved trace metal fluxes from this estuary remain unknown. Here we present data for dissolved (<0.2 µm) trace metals (Al, Mn, Fe, Co, Ni, Cu, Zn, Cd, Pb and U) in the Amazon Estuary during the high discharge season (April–May 2018). We observed distinct trace metal signatures for the Amazon and Pará Rivers, reflecting different catchment areas. Concentrations of the particle-reactive elements (Mn, Fe and Pb) decreased rapidly at low-salinity (S≤2), resulting in the highest estuarine removal (86–94% in the Amazon; 61–70% for the Pará). Co, Ni and Cu removal was comparatively low in both river transects (6–39%), while Cd was the only element with a consistent net input. Chemical fluxes were estimated using (a) endmember concentrations and estuarine removal and (b) combining trace element concentrations with 228Ra fluxes. Relative to global total river fluxes, the Amazon and Pará Rivers combined contribute 21% of dissolved Cu and 18% of dissolved Ni during the high discharge season, but account for comparatively low fractions of Mn, Fe, Co and Zn. These data quantify, for the first time, the trace metal output from the world’s largest and 5th largest river into the Atlantic Ocean, filling a critical gap in knowledge of this globally-important region.

A supermodel connects different models interactively so that their systematic errors compensate and achieve a model with superior performance. It differs from the standard non-interactive multi-model ensembles (NI), which combines model outputs a-posteriori. We formulate the first supermodel framework for Earth System Models (ESMs) and use data assimilation to synchronise models. The ocean of three ESMs is synchronised every month by assimilating pseudo sea surface temperature (SST) observations generated from them. Discrepancies in grid and resolution are handled by constructing the synthetic pseudo-observations on a common grid. We compare the performance of two supermodel approaches to that of the NI for 1980—2006. In the first (EW), the models are connected to the equal-weight multi-model mean, while in the second (SINGLE), they are connected to a single model. Both versions achieve synchronisation in locations where the ocean drives the climate variability. The time variability of the supermodel multi-model mean SST is reduced compared to the observed variability; most where synchronisation is not achieved and is bounded by NI. The damping is larger in EW than in SINGLE because EW yields additional damping of the variability in the individual models. Hence, under partial synchronisation, the part of variability that is not synchronised gets damped in the multi-model average pseudo-observations, causing a deflation during the assimilation. The SST bias in individual models of EW is reduced compared to that of NI, and so is its multi-model mean in the synchronised regions. The performance of a trained supermodel remains to be tested.

Interaction between ocean waves and sea ice may play an important role in sea ice retreat in the Arctic. However, it is difficult to quantify the change of ocean waves propagating in ice as nearly no available measurements. Although SAR has shown the capability of imaging ocean waves in ice-covered areas, there are few attempts to retrieve two-dimensional ocean wave spectra (OWS) by SAR. In this study, we applied the previously developed nonlinear inversion scheme, i.e., the MPI scheme, to retrieve OWS by the Sentinel-1 SAR data acquired in the Barents Sea, where waves penetrate deeply in ice. We compared the retrieved spectra by different combinations of modulation transfer functions (MTFs) used in the MPI scheme, i.e., the same MTFs as those used in retrievals in open water, neglecting both the hydrodynamic and tilt modulations in the MTFs, and neglecting the hydrodynamic modulation but remaining the tilt modulation (a new one fitted in this study for HH-polarized SAR data over ice) in the MTFs. As no in situ measurements (e.g., by directional buoys) available, we compared the simulated SAR image spectra based on the retrievals with the observational SAR image spectra to quantify their respective performances. The comparisons suggest that neglecting hydrodynamic modulation can significantly improve the retrievals. Remaining tilt modulation can further improve the retrievals, particularly for range-travelling waves. The study enhances the understanding of principles of SAR imaging waves in ice and provides basics for retrievals of ocean wave spectra by SAR data in ice-covered areas.

We argue that ocean general circulation models and observations based on Ekman or geostrophic balance provide estimates of the Lagrangian-mean ocean velocity field averaged over surface waves — the total time-averaged velocity that advects oceanic tracers, particles, and water parcels. This interpretation contradicts an assumption often made in ocean transport studies that numerical models and observations based on dynamical balances estimate the Eulerian-mean velocity — the velocity time-averaged at a fixed position and only _part_ of the total ocean velocity. Our argument uses the similarity between the wave-averaged Lagrangian-mean momentum equations appropriate at large oceanic scales, and the momentum equations solved by “wave-agnostic” general circulation models that neglect surface wave effects. We further our case by comparing a realistic, global, “wave-agnostic” general circulation ocean model to a wave-averaged Lagrangian-mean general circulation ocean model at eddy-permitting 1/4-degree resolution, and find that the wave-agnostic velocity field is almost identical to the wave-averaged Lagrangian-mean velocity.

The water mass assembly of Nares Strait is variable, owing to fluctuating wind forcings over the Arctic Basins, and irregular northward flows from the West Greenland Current (WGC) in Baffin Bay. Here we characterize the physico-chemical properties of the water masses entering Nares Strait in August 2014, and we employ an extended optimum multi-parameter (OMP) water mass analysis to estimate the mixing fractions of predefined source water masses, and to distinguish the role of physical and biological processes in governing the distribution of dissolved inorganic carbon (DIC) in Nares Strait. We show the first documented evidence of Siberian shelf waters in Nares Strait, along with a diluted upper halocline layer of partial Pacific-origin. These mixed-origin water masses appear to play an important role in driving a modest phytoplankton bloom in Kane Basin, leading to decreased surface pCO2 concentrations in Nares Strait. Although inorganic nitrogen was already limited in the surface mixed layer in northern Nares Strait, the unique properties of mixed Atlantic-Pacific water facilitated upwelling and nutrient supply to the surface. These observations suggest that the positioning of the Transpolar Drift, and hence the balance of Atlantic and Pacific water delivered to Nares Strait, is likely to play an important role in regional biological productivity and carbon uptake from the atmosphere. We also observed water masses from the WGC transported as far north as Kane Basin, contributing to relatively high pCO2 and low pH in the intermediate and deep water column of southern Nares Strait and northern Baffin Bay.

The Late Miocene Biogenic Bloom (LMBB) is a late Miocene to early Pliocene oceanographic event characterized by high accumulation rates of opal from diatoms and calcite from calcareous nannofossils and planktic foraminifera. This multi-million year event has been recognized in sediment cores from the Pacific, Atlantic and Indian Oceans. The numerous studies discussing the LMBB lead us to believe that this event is omnipresent in all oceans, although this hypothesis need to be tested. Moreover, the origin of this event is still widely discussed. In this study we aim to provide a comprehensive overview of the geographical and temporal aspects of the LMBB by compiling published ocean drilling (DSDP, ODP and IODP) records of sedimentation rates, CaCO\textsubscript{3} and opal and terrigenous accumulation rates that cover the late Miocene and early Pliocene interval. Our data compilation shows that traces of the LMBB are present in many different locations but in a very heterogeneous way, highlighting that the LMBB is not a pervasive event. The compilation in addition shows that the sites where the LMBB is recorded are mainly located in areas with a high productivity regime (i.e. upwelling systems). We suggest that the most likely hypothesis to explain the LMBB is a global increase in upwelling intensity due to an increase in wind strength or an increase in deep water formation, ramping up global thermohaline circulation.

Protists represent the majority of the eukaryotic diversity in the oceans. They have different functions in the marine food web, playing essential roles in the biogeochemical cycles. Meanwhile the available data is rich in horizontal and temporal coverage, little is known on their vertical structuring, particularly below the photic zone. The present study applies DNA metabarcoding to samples collected over three years in conjunction with the BATS time-series to assess marine protist communities in the epipelagic and mesopelagic zones. The protist community showed a dynamic seasonality in the epipelagic, responding to hydrographic yearly cycles. Mixotrophic lineages dominated throughout the year; however, autotrophs bloomed during the rapid transition between the winter mixing and the stratified summer, and heterotrophs had their peak at the end of summer, when the base of the thermocline reaches its deepest depth. Below the photic zone, the community, dominated by Rhizaria, is depth-stratified and relatively constant throughout the year, mirroring local hydrographic and biological features such as the oxygen minimum zone. The results suggest a dynamic partitioning of the water column, where the niche vertical position for each community changes throughout the year, likely depending on nutrient availability, the mixed layer depth, and other hydrographic features. Finally, the protist community closely followed mesoscale events (eddies), where the communities mirrored the hydrographic uplift, raising the deeper communities for hundreds of meters, and compressing the communities above.