Mixed Layer Research Articles

Changes in the SST Seasonal Cycle in a Warmer North Pacific without Ocean Dynamical Feedbacks

Abstract Climate models project a significant intensification of the sea surface temperature (SST) seasonal cycle over the subpolar North Pacific due to global warming, with the shallower mixed layer widely recognized as the dominant factor. However, employing slab ocean experiments with only ocean–atmosphere thermal coupling, we find a substantial contribution from changes in surface heat flux to this seasonal cycle intensification. In particular, the stronger Newtonian cooling effect in winter acts as a more potent damping than in summer. This differential damping inhibits the warming in colder seasons, significantly contributing to the intensified SST seasonal cycle in the subpolar North Pacific. In addition, consistent phase shifts in the North Pacific are identified across CMIP6 models. In the northwest North Pacific, a phase advance is associated with anomalous heating in early spring, driven by enhanced warm atmospheric advection from lower latitudes and sea ice melting in marginal seas. In contrast, the southeast North Pacific exhibits a phase delay attributed to the anomalous cooling in spring relative to autumn. This cooling is due to weakened trade winds and increased presence of high clouds. The former leads to stronger evaporative cooling in spring, while the latter impedes shortwave radiation from reaching the ocean.

The effect of Atmospheric Boundary Layer Meteorology in climate-related gases and VOCs concentrations over 3 European cities using an aerial platform

In the frame of the EU Horizon2020 ICARUS project, a N.A.S.A Awarded Light Manned Aircraft equipped with high-tech scientific instrumentation was used to perform an aerial mapping over Athens, Thessaloniki and Ljubljana greater areas. This study aimed to evaluate the effect of Atmospheric Boundary Layer (ABL) on Green House Gases (GHGs) and Volatile Organic Compounds (VOCs) concentrations over urban and rural areas of the above cities. Simultaneous ground-based measurements were performed in the respective regions. Air samples were pressurized with a stainless-steel bellows compressor into electropolished stainless steel canisters and analyzed by the use of a novel gas chromatographic system. The estimation of the mixing height (ΜΗ) of the ABL was based on the synoptic scale atmospheric circulation and the prevailing background wind. It was found that the MH variation, following the prevailing meteorological conditions, results in different concentration profiles in the lower troposphere over the examined regions. The pollutants concentrations were generally decreasing with altitude in the ABL. Under certain meteorological conditions, vertical mixing plus horizontal transport can cause a high pollution level at the top of ABL. The pollutants concentrations were low over less industrialized and upwind regions, suggesting that local emission sources play significant role on the GHGs and VOCs levels over the regions. However, due to the large scale of sampling area that the aircraft covers the above gradients in concentrations are relative low. AQ modelling activities for simulating the cases studied in the current or any relative future work could reduce operating costs and allow projections of potential impacts.

The Feedback of Ocean Cold Wakes on Three Successive Northward-Moving Tropical Cyclones in 2022

Abstract In 2022, three tropical cyclones (TCs) successively moved northward in the Northwest Pacific. Both the preceding and subsequent TCs can cause obvious sea surface temperature cooling (SSTC), which may last from a few days to a month. In this study, an uncoupled atmosphere-only model and a coupled ocean-atmosphere model were employed to investigate the effects of the SSTC caused by the preceding TCs on the track and intensity of the subsequent TCs. The results indicated that in the atmosphere-only simulations, the subsequent TC had a northeast shift in the track and a reduced intensity. The ocean-atmosphere coupled simulations showed similar TC track adjustments, but different TC intensity change results. The analysis suggested that the deviation of the TC track is the combined effect of the changes in the steering flow, beta drift, and diabatic heating. For the change of TC intensity, not only the SSTC but also the thickness of the oceanic mixed layer needs to be considered. The related physical mechanisms are analyzed in this paper.

Impact of internal wave drag on Arctic sea ice

Abstract A parameterization of the impact of internal waves on momentum transfer at the sea-ice–ocean interface based on previous work by McPhee has been implemented in a sea-ice model for the first time. The ice–ocean drag from internal waves is relevant for shallow mixed layer depth and the presence of a density jump at the pycnocline and is also a function of the strength of the stratification beneath the ocean mixed layer and geometry of the ice interface. We present results from a coupled sea-ice–ocean model where the parameterization of internal wave drag has been implemented. We conducted simulations spanning the years from 2000 to 2017. We find a deceleration of ice drift by 5–8% in both winter and summer, but with significant spatial and temporal variation reaching seasonal average values of ~10%. The spatial variation of ice transport leads to local impacts on deformed ice of magnitude ~0.05 m (2–5%), and reductions in ocean-to-ice heat fluxes of ~1 W m−2, and a decrease in bottom melt of ~0.02–0.04 cm d−1. There is an increase of up to 15% in thickness and ice concentration in the Canadian Arctic and a 10% overall impact on the total sea-ice volume.

The Control of Isolated Kerogen on Pore Structure and Heterogeneity in Marine-Continental Transitional Shale: A Case Study on the Taiyuan Formation, Northern Ordos Basin

Accurately determining the pore structure and heterogeneity characteristics of marine-continental transitional shale in the Taiyuan Formation is crucial for evaluating the shale gas resources in the northern Ordos Basin. However, the studies on pore characteristics and heterogeneity of marine-continental transitional shales and isolated kerogen are limited. This study collected Taiyuan Formation shale in the northern Ordos Basin, and corresponding kerogen isolated from shale and used N2 and CO2 adsorption experiment and Frenkel–Halsey–Hill and Volume-Specific Surface Area model to investigate the pore structure and heterogeneity of both. The results show that the isolated kerogen is dominated by micropores, and the micropore’s specific surface area and volume are 4.7 and 3.5 times the corresponding shale, respectively. In addition, the microporous heterogeneity of the isolated kerogen is stronger than that of shale, while the mesoporous heterogeneity is exactly the opposite. Meanwhile, the micropores fractal dimension Dm is positively correlated with organic matter (OM) content, while mesopores fractal dimension D1 and D2 are negatively linearly correlated with TOC content and have no significant relationship with clay mineral and quartz content (but show a significant positive correlation with illite and illite/smectite mixed layer). Isolated kerogen plays an important role in the pore (especially micropores) heterogeneity of shale, while other minerals (such as clay minerals) have a controlling effect on the mesopores heterogeneity of shale. Compared with marine shale, the marine-continental transitional shale of the Taiyuan Formation has a lower fractal dimension and better connectivity, which is conducive to shale gas seepage and migration. The final result can provide a significant basis for the reserve evaluation and the optimization of desert areas in the marine-continental transitional shale gas in the northern Ordos Basin.

Large Mixed Layer Salinity Variation in the Southern Tropical Indian Ocean Due To the Blending of Water Masses

AbstractThe southern tropical Indian Ocean (TIO) displays large mixed layer salinity (MLS) variation. Circulation in this region is governed by the Indian Ocean tropical gyre (IOTG), where the source water proportion and associated mixing remain unclear. Particles integrating into the IOTG and entering the central southern TIO originate from the Bay of Bengal, Malacca Strait, western Indian Ocean, and Indonesian Throughflow. Surprisingly, cross‐equatorial advection is particularly important, implying a significant connection between both the Bay of Bengal and the South China Sea via Malacca Strait into the southern TIO. The anomalous anticlockwise circulation weakens the IOTG during positive Indian Ocean Dipole (IOD). An opposite pattern is observed in the negative IOD. A particle experiment reveals that water masses are modulated by the anomalous circulation that drives the redistribution of MLS by changing the proportion of the different source waters. This represents a potential predictability for the southern TIO MLS variability.

Prediction of the Relative Resource Abundance of the Argentine Shortfin Squid Illex argentinus in the High Sea in the Southwest Atlantic Based on a Deep Learning Model.

To analyze the impact of the marine environment on the relative abundance of Illex argentinus (high and low categories) in the southwest Atlantic, this study collected logbook data from Chinese pelagic trawlers from December 2014 to June 2024, including vessel position data and oceanographic variables such as sea surface temperature, 50 m and 100 m water temperature, sea surface salinity, sea surface height, chlorophyll-a concentration, and mixed layer depth. Vessel positions were used to enhance the logbook data quality, allowing an analysis of the annual trends in the resource center of this squid at a spatial resolution of 0.1° × 0.1° and a temporal resolution of ten days. The findings showed that the resource center is primarily located around 42° S in the north and between 45° S and 47° S in the south, with a trend of northward movement during the study period. Additionally, we constructed two ensemble learning models based on decision trees-AdaBoost and PSO-RF-aiming to identify the most critical environmental factors that affect its resource abundance; we found that the optimal model was the PSO-RF model with max_depth of 5 and n_estimators of 46. The importance analysis revealed that sea surface temperature, mixed layer depth, sea surface height, sea surface salinity, and 50 m water temperature are critical environmental factors affecting this species' resources. Given that deep learning models generally have shorter running times and higher accuracy than other models, we developed a CNN-Attention model based on the five most important input factors. This model achieved an accuracy of 73.6% in forecasting this squid for 2024, predicting that the population would first appear near the Argentine exclusive economic zone around mid-December 2023 and gradually move east and south thereafter. The predictions of the model, validated through log data, maintained over 70% accuracy during most periods at a time scale of ten days. The successful construction of the resource abundance forecasting model and its accuracy improvements can help enterprises save fuel and time costs associated with blind searches for target species. Moreover, this research contributes to improving resource utilization efficiency and reducing fishing duration, thereby aiding in lowering carbon emissions from pelagic trawling activities, offering valuable insights for the sustainable development of this species' resources.

Evaluation of the Regulation and Response of Ocean Heat Content to Typhoons in the Northwest Pacific Based on GDCSM_Argo

The Northwest Pacific is a source of strong typhoons accompanied by significantly changing ocean heat content (OHC). Based on the GDCSM_Argo dataset, all typhoons within twenty years are considered to explore the regulation of OHC between 0 and 2000 m with typhoons in this study. The results show a weak correlation between OHC and typhoons but a strong correlation between ocean heat content uptake anomaly (OHUA) and typhoons by wavelet, lag/lead, and EOF analyses. The energy reserve effect by the OHUA on typhoon generation is most significant at a 3-month lag. The OHUA above 300 m primarily influences typhoons’ temporal and spatial characteristics. On the other hand, typhoons enhance vertical transmission of OHC through entrainment. A modified version of the Richardson number (Rohc) is used for the characterization of the hydrostatic stability between ocean layers. The entrainment is particularly strong in the Northwest Pacific north of 20° N, where Rohc is between −0.2 and 0. The OHC below the mixed layer also responds to typhoon disturbances but transports less heat vertically than the near-surface layers. Our results suggest that GDCSM_Argo provides reliable data support for the study of the relationship between OHC and typhoons.

Impact of ocean vertical-mixing parameterization on Arctic sea ice and upper-ocean properties using the NEMO-SI3 model

Abstract. We evaluate the vertical turbulent-kinetic-energy (TKE) mixing scheme of the NEMO-SI3 ocean–sea-ice model in sea-ice-covered regions of the Arctic Ocean. Specifically, we assess the parameters involved in TKE mixed-layer-penetration (MLP) parameterization. This ad hoc parameterization aims to capture processes that impact the ocean surface boundary layer, such as near-inertial oscillations, ocean swells, and waves, which are often not well represented in the default TKE scheme. We evaluate this parameterization for the first time in three regions of the Arctic Ocean: the Makarov, Eurasian, and Canada basins. We demonstrate the strong effect of the scaling parameter that accounts for the presence of sea ice. Our results confirm that TKE MLP must be scaled down below sea ice to avoid unrealistically deep mixed layers. The other parameters evaluated are the percentage of energy penetrating below the mixed layer and the length scale of its decay with depth. All these parameters affect mixed-layer depth and its seasonal cycle, surface temperature, and salinity, as well as underlying stratification. Shallow mixed layers are associated with stronger stratification and fresh surface anomalies, and deeper mixed layers correspond to weaker stratification and salty surface anomalies. Notably, we observe significant impacts on sea-ice thickness across the Arctic Ocean in two scenarios: when the scaling parameter due to sea ice is absent and when the TKE mixed-layer-penetration process vanishes. In the former case, we observe an increase of several meters in mixed-layer depth, along with a reduction in sea-ice thickness ranging from 30 to 40 cm, reflecting the impact of stronger mixing. Conversely, in the latter case, we notice that a shallower mixed layer is accompanied by a moderate increase in sea-ice thickness, ranging from 10 to 20 cm, as expected from weaker mixing. Additionally, interannual variability suggests that experiments incorporating a scaling parameter based on sea-ice concentration display an increased mixed-layer depth during periods of reduced sea ice, which is consistent with observed trends. These findings underscore the influence of enhanced ocean mixing, through specific parameterizations, on the physical properties of the upper ocean and sea ice.

Effects of Anisotropy, Convection, and Relaxation on Nonlinear Reaction-Diffusion Systems

The effects of relaxation, convection, and anisotropy on a two-dimensional, two-equation system of nonlinearly coupled, second-order hyperbolic, advection–reaction–diffusion equations are studied numerically by means of a three-time-level linearized finite difference method. The formulation utilizes a frame-indifferent constitutive equation for the heat and mass diffusion fluxes, taking into account the tensorial character of the thermal diffusivity of heat and mass diffusion. This approach results in a large system of linear algebraic equations at each time level. It is shown that the effects of relaxation are small although they may be noticeable initially if the relaxation times are smaller than the characteristic residence, diffusion, and reaction times. It is also shown that the anisotropy associated with one of the dependent variables does not have an important role in the reaction wave dynamics, whereas the anisotropy of the other dependent variable results in transitions from spiral waves to either large or small curvature reaction fronts. Convection is found to play an important role in the reaction front dynamics depending on the vortex circulation and radius and the anisotropy of the two dependent variables. For clockwise-rotating vortices of large diameter, patterns similar to those observed in planar mixing layers have been found for anisotropic diffusion tensors.

Aerosol Vertical Turbulent Mass Flux Retrievals Through Novel Remote Sensing Algorithm

AbstractIntegrated measurements of aerosol, radiation, cloud, and turbulent transport in the planetary boundary layer (PBL) are essential for understanding and modeling climate and air quality. Here, we developed a new technique for the identification of convective turbulent regions and deriving the vertical distribution of aerosol turbulent mass fluxes within PBL. The algorithm uses retrievals from coherent Doppler lidars and a high spectral resolution lidar. The technique was applied to study particle mass fluxes over 2 months (November–December 2020) during the campaign conducted at the DOE Atmospheric Radiation Measurement Southern Great Plains (SGP) site in Lamont, Oklahoma. The algorithm developed here is capable of continuously deriving vertically resolved (curtains) aerosol mass fluxes. Our data analysis shows that at the site, the 30‐min averaged fluxes at 135 m above the surface were mainly positive (upward) at ∼1 μg m−2 s−1, suggesting that the surface is the primary source of the particle mass supplied to the boundary layer at the SGP site. Analyses of the individual case studies have revealed that not all the derived fluxes can be linked to surface emissions. Both positive and negative values in a range of ±5 μg m−2 s−1 can be caused by convective thermals interacting between the residual layer and the mixed layer and by rotation of the horizontal wind with the height. Large erroneous negative fluxes can also be caused by drizzling/precipitating clouds. We anticipate that the application of the current technique will lead to a more realistic representation of aerosol mass budgets and bidirectional mixing rates.

Planktonic to sessile: drivers of spatial and temporal variability across barnacle life stages and indirect effects of the Pacific Marine Heatwave

ABSTRACT Barnacles are a foundation species in intertidal habitats. During the Pacific Marine Heatwave (PMH), intertidal barnacle cover increased in the northern Gulf of Alaska (GoA); however, the role of pelagic larval supply in this increase was unknown. Using long-term monitoring data on intertidal benthic (percent cover) and pelagic larval populations (nauplii and cyprid concentrations), we examined potential environmental drivers (temperature, chlorophyll-a, mixed layer depth) of larval concentration and whether including larval concentration at regional and annual scales improved intertidal barnacle percent cover models in two study regions in the GoA. In both regions, larval concentrations were slightly higher following the PMH. Percent cover models were improved by including cyprid concentrations (but not nauplii), and the effect strength varied by site and tidal elevation. This indicates that larval concentration contributes as a bottom–up driver of benthic barnacle abundance. There is little evidence of a direct effect of the PMH on either life stage. Instead, our results may illustrate the positive feedback between life stages, where higher adult benthic abundance increased larval concentrations, which then supplied more new recruits to the benthos. As heatwaves continue to occur, integrating various data types can provide insights into factors influencing both benthic and pelagic communities.

Geochemical characterization of trace fossil assemblages in spotted marls and limestones of the Lower Jurassic of the Western Carpathians: environmental implications

AbstractThe trace fossils of the spotted limestones and marls from upper Pliensbachian of the Skladaná Skala section, Western Carpathians (Slovakia), are characteristic of Zoophycos ichnofacies preserved in outer shelf deposits. The ichnoassociation is dominated by Lamellaeichnus and Chondrites, while Palaeophycus, Planolites, Teichichnus, Thalassinoides, Trichichnus, and Zoophycos, are less abundant. The presence of relatively large sized traces (reaching a decimetre in horizontal dimension for Thalassinoides and Zoophycos) and continuously well bioturbated deposits point to well oxygenated shallow zone (mixed layer) and deeper part of stiffed substrate of the transitional layer of bottom substrate up to the first tens of centimetres to depth. A first tier (tier I) corresponds to the shallowest sediment occupied by Bathysiphon tests and other epifaunal body fossils. The tier II is characterized by shallow infaunal burrows of deposit feeders (Planolites) and permanent dwelling, open burrows (Palaeophycus, Thalassinoides) located in a Ca-rich sediment (corresponding to CaCO3). The tier III is dominated by deposit feeders (Lamellaeichnus, Teichichnus) that burrowed stiff Ca-rich sediment below the redox boundary. Trace fossils of tier III are rich in fine siliciclastics (high content in Si, Al, K mainly associated to clays) and pyrite framboids (enrichment in Fe, S and Co). The deposit-feeders of tier III maintain connection to the sediment-water interface. Tier IV (Zoophycos, Teichichnus and large Chondrites), and deepest infaunal forms of the tier V. (Trichichnus, Pilichnus and tiny forms of Chondrites) comprises trace fossils commonly rich in pyrite framboids with relatively high content of Fe, S, and Co, congruent with chemosymbiotic behaviour commonly inferred for Chondrites, Trichichnus, and the organic-matter storage of Zoophycos.

Response of the upper ocean to northeast Pacific atmospheric rivers under climate change

Atmospheric rivers are important transport vehicles for Earth’s water cycle. Using a high-resolution, eddy-resolving Earth System Model, atmospheric river impacts on the upper ocean are investigated by analyzing historical and climate change simulations. For atmospheric rivers along the North American coastline, strong winds cause significant dynamic and thermodynamic upper ocean responses. They push ocean water towards the coast, measured by sea surface height, a process that is amplified under climate change. Mixed layers are deeper upstream of atmospheric rivers, and shallower downstream, however for climate change, shoaling downstream is subdued. Air-sea heat fluxes tend to promote ocean cooling upstream and warming downstream, although different regions have different climate change heat flux signals. Southern California heat flux changes due to warming are driven by evaporative processes and strengthen the ocean responses seen in historical simulations. The regions north are primarily dominated by sensible heat flux changes and counter the historical patterns.

Effects of internal intervention in tubes on shock wave attenuation, formation and propagation of flame during high-pressure hydrogen release

Experiments on the impacts of the combination of fin structure and expansion chamber on shock wave attenuation, hydrogen ignition, and flame characteristics during high-pressure hydrogen release were carried out. The fin is a four-edged platform hollowed out on all sides. The inflated upper and lower tube walls are located 85, 95, and 105 mm away from the nickel burst disk. With the incident shock wave traveling through the fin, the reflected shock and the superposition of multiple reflected shock waves appear, and the reflected shock wave intensity can equal that of the release pressure. The decrease of the incident shock wave intensity behind the fin leads to the decrease of the non-uniformity of the mixed layer temperature, further causing the decrease of the flame propagation velocity. The presence of fins increases the probability of ignition in the fin region and the critical release pressure for ignition is significantly lower than that in the barrier free tube, mainly because the fin inside intensifies the hydrogen/air mixture, resulting in a lower minimum ignition energy required for self-ignition. Furthermore, the decrease of flame propagation velocity leads to the increase of hydrogen consumption inside the tube and the decrease of shock overpressure outside.

TriMLP : A Foundational MLP-Like Architecture for Sequential Recommendation

In this work, we present TriMLP as a foundational MLP-like architecture for the sequential recommendation, simultaneously achieving computational efficiency and promising performance. First, we empirically study the incompatibility between existing purely MLP-based models and sequential recommendation, that the inherent fully-connective structure endows historical user–item interactions (referred as tokens) with unrestricted communications and overlooks the essential chronological order in sequences. Then, we propose the MLP-based Triangular Mixer to establish ordered contact among tokens and excavate the primary sequential modeling capability under the standard auto-regressive training fashion. It contains (1) a global mixing layer that drops the lower-triangle neurons in MLP to block the anti-chronological connections from future tokens and (2) a local mixing layer that further disables specific upper-triangle neurons to split the sequence as multiple independent sessions. The mixer serially alternates these two layers to support fine-grained preferences modeling, where the global one focuses on the long-range dependency in the whole sequence, and the local one calls for the short-term patterns in sessions. Experimental results on 12 datasets of different scales from 4 benchmarks elucidate that TriMLP consistently attains favorable accuracy/efficiency tradeoff over all validated datasets, where the average performance boost against several state-of-the-art baselines achieves up to 14.88%, and the maximum reduction of inference time reaches 23.73%. The intriguing properties render TriMLP a strong contender to the well-established RNN-, CNN-, and Transformer-based sequential recommenders. Code is available at https://github.com/jiangyiheng1/TriMLP .

Organised large structure in the post-transition mixing layer. Part 3. Dynamics of the spatial growth

Organised large structure in the post-transition mixing layer. Part 3. Dynamics of the spatial growth

Advanced Peak Phase of ENSO under Global Warming

Abstract El Niño–Southern Oscillation (ENSO) is the leading mode of interannual ocean–atmosphere coupling in the tropical Pacific, greatly influencing the global climate system. Seasonal phase locking, which means that ENSO events usually peak in boreal winter, is a distinctive feature of ENSO. In model future projections, the ENSO sea surface temperature (SST) amplitude in winter shows no significant change with a large intermodel spread. However, whether and how ENSO phase locking will respond to global warming are not fully understood. In this study, using Community Earth System Model Large Ensemble (CESM-LE) projections, we found that the seasonality of ENSO events, especially its peak phase, has advanced under global warming. This phenomenon corresponds to the seasonal difference in the changes in the ENSO SST amplitude with an enhanced (weakened) amplitude from boreal summer to autumn (winter). Mixed layer ocean heat budget analysis revealed that the advanced ENSO seasonality is due to intensified positive meridional advective and thermocline feedback during the ENSO developing phase and intensified negative thermal damping during the ENSO peak phase. Furthermore, the seasonal variation in the mean El Niño–like SST warming in the tropical Pacific favors a weakened zonal advective feedback in boreal autumn–winter and earlier decay of ENSO. The advance of the ENSO peak phase is also found in most CMIP5/6 models that simulate the seasonal phase locking of ENSO well in the present climate. Thus, even though the amplitude response in the winter shows no model consensus, ENSO also significantly changes during different stages under global warming.

Phytoplankton Spring Bloom Inhibited by Marine Heatwaves in the North‐Western Mediterranean Sea

AbstractMarine heatwaves (MHWs) represent anomalously warm temperature conditions of seawater that may affect marine life and ocean biogeochemistry. Under such conditions, phytoplankton communities may modify their structure and functions, and their resilience is not assured. This study characterizes the impact of MHWs on the phytoplankton spring bloom in the North‐Western Mediterranean Sea. Here, we synergistically combine autonomous observations from BioGeoChemical‐Argo floats, satellite‐based and marine ecosystem model data, and show that MHW events occurring during winter drastically inhibit phytoplankton carbon biomass in spring by up to 70%. Such reduction is related to the enhanced stratification of the water column under MHWs which hinders the renewal of nutrients from deep‐ocean reservoirs, thus preventing surface phytoplankton from blooming. This process negatively impacts particulate organic carbon stocks within the mixed layer, while severe events cause an earlier shift of phytoplankton phenology that provokes changes in zooplankton biomass distribution.

Drivers of phytoplankton distribution, abundance and community composition off East Antarctica, from 55-80°E (CCAMLR Division 58.4.2 East)

Southern Ocean phytoplankton form the base of the Antarctic food web, influencing higher trophic levels through biomass and community structure. We examined phytoplankton distribution and abundance in the Indian Sector of the Southern Ocean during austral summer as part a multidisciplinary ecosystem survey: Trends in Euphausiids off Mawson, Predators and Oceanography (TEMPO, 2021). Sampling covered six meridional transects from 55-80°E, and from 62°S or 63°S to the ice edge. To determine phytoplankton groups, CHEMTAX analysis was undertaken on pigments measured using HPLC. Diatoms were the dominant component of phytoplankton communities, explaining 56% of variation in chlorophyll a (Chl a), with haptophytes also being a major component. Prior to sampling the sea ice had retreated in a south-westerly direction, leading to shorter ice-free periods in the west (< 44 days, ≤65°E) compared to east (> 44 days, ≥70°E), inducing a strong seasonal effect. The east was nutrient limited, indicated by low-iron forms of haptophytes, and higher silicate:nitrate drawdown ratios (5.1 east vs 4.3 west), pheophytin a (phaeo) concentrations (30.0 vs 18.4 mg m-2) and phaeo:Chl a ratios (1.06 vs 0.53). Biological influences were evident at northern stations between 75-80°E, where krill “super-swarms” and feeding whales were observed. Here, diatoms were depleted from surface waters likely due to krill grazing, as indicated by high phaeo:Chl a ratios (> 0.75), and continued presence of haptophytes, associated with inefficient filtering or selective grazing by krill. Oceanographic influences included deeper mixed layers reducing diatom biomass, and a bloom to the north of the southern Antarctic Circumpolar Current Front in the western survey area thought to be sinking as waters flowed from west to east. Haptophytes were influenced by the Antarctic Slope Front with high-iron forms prevalent to the south only, showing limited iron transfer from coastal waters. Cryptophytes were associated with meltwater, and greens (chlorophytes + prasinophytes) were prevalent below the mixed layer. The interplay of seasonal, biological and oceanographic influences on phytoplankton populations during TEMPO had parallels with processes observed in the BROKE and BROKE-West voyages conducted 25 and 15 years earlier, respectively. Our research consolidates understanding of the krill ecosystem to ensure sustainable management in East Antarctic waters.