Browse
Abstract
The formation of a sharp oceanic front located south-southeast of Sri Lanka during the southwest monsoon is examined through in situ and remote observations and high-resolution model output. Remote sensing and model output reveal that the front extends approximately 200 km eastward from the southeast coast of Sri Lanka toward the southern Bay of Bengal (BoB). This annually occurring front is associated with the boundary between the southwest monsoon current with high-salinity water to the south, and a weak flow field comprised of relatively fresh BoB water to the north. The front contains a line of high chlorophyll extending from the coastal upwelling zone, often for several hundred kilometers. Elevated turbulent diffusivities ∼10−2 m2 s−1 along with large diapycnal fluxes of heat and salt were found within the front. The formation of the front and vertical transports are linked to local wind stress curl. Large vertical velocities (∼50 m day−1) indicate the importance of ageostrophic, submesoscale processes. To examine these processes, the Ertel potential vorticity (PV) was computed using the observations and numerical model output. The model output shows a ribbon of negative PV along the front between the coastal upwelling zone and two eddies (Sri Lanka Dome and an anticyclonic eddy) typically found in the southern BoB. PV estimates support the view that the flow is susceptible to submesoscale instabilities, which in turn generate high vertical velocities within the front. Frontal upwelling and heightened mixing show that the seasonal front is regionally important to linking the fresh surface water of the BoB with the Arabian Sea.
Significance Statement
Within the ocean, motions span extraordinarily wide ranges of sizes and time scales. In this study we focus on a narrow, intensified feature called a front. This front occurs in the southern Bay of Bengal during the summer monsoon and forms a boundary between fresher water to the north and saltier water to the south. Features such as this are difficult to study, however, by combining observations made from ships and satellites with output from numerical models of the ocean, we are able to better understand the front. This is important because fronts like the one studied here play a role in determining the pathways of heat within the ocean, which, in turn, may feedback into the atmosphere and weather patterns.
Abstract
The formation of a sharp oceanic front located south-southeast of Sri Lanka during the southwest monsoon is examined through in situ and remote observations and high-resolution model output. Remote sensing and model output reveal that the front extends approximately 200 km eastward from the southeast coast of Sri Lanka toward the southern Bay of Bengal (BoB). This annually occurring front is associated with the boundary between the southwest monsoon current with high-salinity water to the south, and a weak flow field comprised of relatively fresh BoB water to the north. The front contains a line of high chlorophyll extending from the coastal upwelling zone, often for several hundred kilometers. Elevated turbulent diffusivities ∼10−2 m2 s−1 along with large diapycnal fluxes of heat and salt were found within the front. The formation of the front and vertical transports are linked to local wind stress curl. Large vertical velocities (∼50 m day−1) indicate the importance of ageostrophic, submesoscale processes. To examine these processes, the Ertel potential vorticity (PV) was computed using the observations and numerical model output. The model output shows a ribbon of negative PV along the front between the coastal upwelling zone and two eddies (Sri Lanka Dome and an anticyclonic eddy) typically found in the southern BoB. PV estimates support the view that the flow is susceptible to submesoscale instabilities, which in turn generate high vertical velocities within the front. Frontal upwelling and heightened mixing show that the seasonal front is regionally important to linking the fresh surface water of the BoB with the Arabian Sea.
Significance Statement
Within the ocean, motions span extraordinarily wide ranges of sizes and time scales. In this study we focus on a narrow, intensified feature called a front. This front occurs in the southern Bay of Bengal during the summer monsoon and forms a boundary between fresher water to the north and saltier water to the south. Features such as this are difficult to study, however, by combining observations made from ships and satellites with output from numerical models of the ocean, we are able to better understand the front. This is important because fronts like the one studied here play a role in determining the pathways of heat within the ocean, which, in turn, may feedback into the atmosphere and weather patterns.
Abstract
Around Hopen Island, the satellite images and experiments with drifting buoys describe the movement of the drifting ice and depict tidally generated trapped motion. An analytical solution is applied to investigate the trapping phenomenon. A general solution is achieved by the superposition of the incident and reflected (scattered) waves for an elliptically shaped island above the critical latitude. The incident wave simulates the tidal wave propagation toward the island and its prominent feature, an amphidromic point located to the southeast from Hopen Island. The analytical solution for the reflected wave is constructed in elliptic coordinates. Tide amplitudes and cophase lines are analyzed in the island’s vicinity and compared to observations and numerical model results. A simulated drift of Lagrangian water particles constructed with the help of analytical solutions reproduces well the observed clockwise trapped motion of the drifting buoy near Hopen Island. Since the resonance may amplify the semidiurnal incident tide, we have also investigated the natural modes of water oscillations near the island. While this paper focuses on the details of the model used at the specific site of Hopen Island, a similar trapping analysis can be applied to circular or elliptic islands that have a small scale relative to the barotropic Rossby deformation radius.
Significance Statement
This study aims to understand how the semidiurnal tide propagates and generates strong currents near Hopen Island in the Barents Sea. The trapping of the semidiurnal (M2) tide around Hopen Island leads to an organized dipole structure in sea level, which rotates clockwise. The dipole generates maximum amplitudes of water surface elevation and the strong current near the south and north tips of the island. The abrupt sea level change induced by the dipole sets up often violent currents, which, together with drifting ice, can be dangerous for navigation. The strong tidal currents generate permanent clockwise circulation around the islands, which is essential for biological life and waste disposal as material disposed of near the islands will be trapped for an extended time. Our investigation elucidates the role of dipoles in the local enhancement of tides around the islands.
Abstract
Around Hopen Island, the satellite images and experiments with drifting buoys describe the movement of the drifting ice and depict tidally generated trapped motion. An analytical solution is applied to investigate the trapping phenomenon. A general solution is achieved by the superposition of the incident and reflected (scattered) waves for an elliptically shaped island above the critical latitude. The incident wave simulates the tidal wave propagation toward the island and its prominent feature, an amphidromic point located to the southeast from Hopen Island. The analytical solution for the reflected wave is constructed in elliptic coordinates. Tide amplitudes and cophase lines are analyzed in the island’s vicinity and compared to observations and numerical model results. A simulated drift of Lagrangian water particles constructed with the help of analytical solutions reproduces well the observed clockwise trapped motion of the drifting buoy near Hopen Island. Since the resonance may amplify the semidiurnal incident tide, we have also investigated the natural modes of water oscillations near the island. While this paper focuses on the details of the model used at the specific site of Hopen Island, a similar trapping analysis can be applied to circular or elliptic islands that have a small scale relative to the barotropic Rossby deformation radius.
Significance Statement
This study aims to understand how the semidiurnal tide propagates and generates strong currents near Hopen Island in the Barents Sea. The trapping of the semidiurnal (M2) tide around Hopen Island leads to an organized dipole structure in sea level, which rotates clockwise. The dipole generates maximum amplitudes of water surface elevation and the strong current near the south and north tips of the island. The abrupt sea level change induced by the dipole sets up often violent currents, which, together with drifting ice, can be dangerous for navigation. The strong tidal currents generate permanent clockwise circulation around the islands, which is essential for biological life and waste disposal as material disposed of near the islands will be trapped for an extended time. Our investigation elucidates the role of dipoles in the local enhancement of tides around the islands.
Abstract
Recent work has shown that, when nontraditional (NT) effects associated with the horizontal component of the Coriolis parameter are taken into account, equatorial waves (EWs) experience critical reflection when they reflect off the seafloor at the latitude where their frequency is equal to the inertial frequency. As a result, the vertical shear associated with the wave is strongly enhanced locally and results in bottom-intensified mixing. Using an off-the-shelf parameterization for mixing, these studies have shown that this process could play an important role in driving diapycnal upwelling in the abyssal ocean, but the specific mechanisms generating the mixing have not been studied yet. In this work, we address this limitation by running two-dimensional, high-resolution, nonhydrostatic simulations of the critical reflection of internal waves modified by NT effects. These simulations can resolve the instabilities triggered when the wave reflects off the bottom, allowing us to characterize the energy cascade to smaller scales and to estimate the mixing it generates. We find that shear instabilities drive elevated turbulent diffusivities between 10−1 and −10−3 m2 s−1 over a critical layer of 100–300 m thick. The shear instabilities result directly from the enhancement of kinetic energy in the reflected wave that is confined against the seafloor during the critical reflection process. Simultaneously, higher harmonics are generated and flux energy upward in the water column. These higher harmonics are unstable to parametric subharmonic instability, which absorbs their energy and drive enhanced dissipation above the critical layer, to a height of O(1000) m off the bottom. We show how these results depend on key elements of the EWs and of the medium and discuss the implementation of a parameterization of these effects in global ocean models.
Abstract
Recent work has shown that, when nontraditional (NT) effects associated with the horizontal component of the Coriolis parameter are taken into account, equatorial waves (EWs) experience critical reflection when they reflect off the seafloor at the latitude where their frequency is equal to the inertial frequency. As a result, the vertical shear associated with the wave is strongly enhanced locally and results in bottom-intensified mixing. Using an off-the-shelf parameterization for mixing, these studies have shown that this process could play an important role in driving diapycnal upwelling in the abyssal ocean, but the specific mechanisms generating the mixing have not been studied yet. In this work, we address this limitation by running two-dimensional, high-resolution, nonhydrostatic simulations of the critical reflection of internal waves modified by NT effects. These simulations can resolve the instabilities triggered when the wave reflects off the bottom, allowing us to characterize the energy cascade to smaller scales and to estimate the mixing it generates. We find that shear instabilities drive elevated turbulent diffusivities between 10−1 and −10−3 m2 s−1 over a critical layer of 100–300 m thick. The shear instabilities result directly from the enhancement of kinetic energy in the reflected wave that is confined against the seafloor during the critical reflection process. Simultaneously, higher harmonics are generated and flux energy upward in the water column. These higher harmonics are unstable to parametric subharmonic instability, which absorbs their energy and drive enhanced dissipation above the critical layer, to a height of O(1000) m off the bottom. We show how these results depend on key elements of the EWs and of the medium and discuss the implementation of a parameterization of these effects in global ocean models.
Abstract
The Gulf Stream (GS) is one of the strongest ocean currents on the planet. Eddy-rich resolution models are needed to properly represent the dynamics of the GS; however, kinetic energy (KE) can be in excess in these models if not dissipated efficiently. The question of how and how much energy is dissipated and in particular how it flows through ocean scales thus remains an important and largely unanswered question. Using a high-resolution (∼2 km) ocean model [Coastal and Regional Ocean Community (CROCO)], we characterize the spatial and temporal distribution of turbulent cascades in the GS based on a coarse-grained method. We show that the balanced flow is associated with an inverse cascade while the forward cascade is explained by ageostrophic advection associated with frontogenesis. Downscale fluxes are dominant at scales smaller than about 20 km near the surface and most intense at the GS North Wall. There is also strong seasonal variability in KE flux, with the forward cascade intensifying in winter and the inverse cascade later in spring. The forward cascade, which represents an interior route to dissipation, is compared with both numerical and boundary dissipation processes. The contribution of interior dissipation is an order of magnitude smaller than that of the other energy sinks. We thus evaluate the sensitivity of horizontal momentum advection schemes on energy dissipation and show that the decrease in numerical dissipation in a high-order scheme leads to an increase in dissipation at the boundaries, not in the downscale flux.
Abstract
The Gulf Stream (GS) is one of the strongest ocean currents on the planet. Eddy-rich resolution models are needed to properly represent the dynamics of the GS; however, kinetic energy (KE) can be in excess in these models if not dissipated efficiently. The question of how and how much energy is dissipated and in particular how it flows through ocean scales thus remains an important and largely unanswered question. Using a high-resolution (∼2 km) ocean model [Coastal and Regional Ocean Community (CROCO)], we characterize the spatial and temporal distribution of turbulent cascades in the GS based on a coarse-grained method. We show that the balanced flow is associated with an inverse cascade while the forward cascade is explained by ageostrophic advection associated with frontogenesis. Downscale fluxes are dominant at scales smaller than about 20 km near the surface and most intense at the GS North Wall. There is also strong seasonal variability in KE flux, with the forward cascade intensifying in winter and the inverse cascade later in spring. The forward cascade, which represents an interior route to dissipation, is compared with both numerical and boundary dissipation processes. The contribution of interior dissipation is an order of magnitude smaller than that of the other energy sinks. We thus evaluate the sensitivity of horizontal momentum advection schemes on energy dissipation and show that the decrease in numerical dissipation in a high-order scheme leads to an increase in dissipation at the boundaries, not in the downscale flux.
Abstract
Long waves play an important role in coastal inundation and shoreline and dune erosion, requiring a detailed understanding of their evolution in nearshore regions and interaction with shorelines. While their generation and dissipation mechanisms are relatively well understood, there are fewer studies describing how reflection processes govern their propagation in the nearshore. We propose a new approach, accounting for partial reflections, which leads to an analytical solution to the free wave linear shallow-water equations at the wave-group scale over general varying bathymetry. The approach, supported by numerical modeling, agrees with the classic Bessel standing solution for a plane sloping beach but extends the solution to arbitrary alongshore uniform bathymetry profiles and decomposes it into incoming and outgoing wave components, which are a combination of successively partially reflected waves lagging each other. The phase lags introduced by partial reflections modify the wave amplitude and explain why Green’s law, which describes the wave growth of free waves with decreasing depth, breaks down in very shallow water. This reveals that the wave amplitude at the shoreline is highly dependent on partial reflections. Consistent with laboratory and field observations, our analytical model predicts a reflection coefficient that increases and is highly correlated with the normalized bed slope (bed slope relative to wave frequency). Our approach shows that partial reflections occurring due to depth variations in the nearshore are responsible for the relationship between the normalized bed slope and the amplitude of long waves in the nearshore, with direct implications for determining long-wave amplitudes at the shoreline and wave runup.
Abstract
Long waves play an important role in coastal inundation and shoreline and dune erosion, requiring a detailed understanding of their evolution in nearshore regions and interaction with shorelines. While their generation and dissipation mechanisms are relatively well understood, there are fewer studies describing how reflection processes govern their propagation in the nearshore. We propose a new approach, accounting for partial reflections, which leads to an analytical solution to the free wave linear shallow-water equations at the wave-group scale over general varying bathymetry. The approach, supported by numerical modeling, agrees with the classic Bessel standing solution for a plane sloping beach but extends the solution to arbitrary alongshore uniform bathymetry profiles and decomposes it into incoming and outgoing wave components, which are a combination of successively partially reflected waves lagging each other. The phase lags introduced by partial reflections modify the wave amplitude and explain why Green’s law, which describes the wave growth of free waves with decreasing depth, breaks down in very shallow water. This reveals that the wave amplitude at the shoreline is highly dependent on partial reflections. Consistent with laboratory and field observations, our analytical model predicts a reflection coefficient that increases and is highly correlated with the normalized bed slope (bed slope relative to wave frequency). Our approach shows that partial reflections occurring due to depth variations in the nearshore are responsible for the relationship between the normalized bed slope and the amplitude of long waves in the nearshore, with direct implications for determining long-wave amplitudes at the shoreline and wave runup.
Abstract
Mixing along sloping isopycnals plays a key role in the transport and uptake of heat and carbon by the ocean. This mixing is quantified by a lateral diffusivity, which can be measured by tracking the lateral spreading of point release tracer patches. We present a definition for the area of a tracer patch, the time derivative of which provides the lateral diffusivity. To accurately estimate the diffusivity, an ensemble mean concentration field of many tracer release experiments is required. We use numerical experiments to quantify how accurately the “true” lateral diffusivity (obtained from the ensemble mean concentration field) can be estimated from a single tracer release experiment (one ensemble member). To simulate observational campaigns, we also estimate the diffusivity from a single tracer release that is spatially and/or temporally subsampled, quantifying how the error between the estimated diffusivity and the true diffusivity grows as this sampling resolution worsens. We perform these numerical experiments in a two-layer quasigeostrophic model of turbulent flow on a β plane, using an ensemble of 50 passive tracer release experiments, each initialized as a 2D Gaussian but with differing realizations of the turbulent flow. We find that the diffusivity estimates from the single tracer releases have a relative root-mean-square error (RMSE) of 1.43% from the true diffusivity. Subsampling a single tracer release experiment every 956 km increases the relative RMSE from the true diffusivity to 3.1%; also subsampling every 277 days raises this figure to 6.5%.
Abstract
Mixing along sloping isopycnals plays a key role in the transport and uptake of heat and carbon by the ocean. This mixing is quantified by a lateral diffusivity, which can be measured by tracking the lateral spreading of point release tracer patches. We present a definition for the area of a tracer patch, the time derivative of which provides the lateral diffusivity. To accurately estimate the diffusivity, an ensemble mean concentration field of many tracer release experiments is required. We use numerical experiments to quantify how accurately the “true” lateral diffusivity (obtained from the ensemble mean concentration field) can be estimated from a single tracer release experiment (one ensemble member). To simulate observational campaigns, we also estimate the diffusivity from a single tracer release that is spatially and/or temporally subsampled, quantifying how the error between the estimated diffusivity and the true diffusivity grows as this sampling resolution worsens. We perform these numerical experiments in a two-layer quasigeostrophic model of turbulent flow on a β plane, using an ensemble of 50 passive tracer release experiments, each initialized as a 2D Gaussian but with differing realizations of the turbulent flow. We find that the diffusivity estimates from the single tracer releases have a relative root-mean-square error (RMSE) of 1.43% from the true diffusivity. Subsampling a single tracer release experiment every 956 km increases the relative RMSE from the true diffusivity to 3.1%; also subsampling every 277 days raises this figure to 6.5%.
Abstract
The bottom boundary layer (BBL) contributes significantly to the global energy dissipation of low-frequency flows in the abyssal ocean, but how this dissipation occurs remains poorly understood. Using in situ data collected near the BBL at an abyssal seamount in the western Pacific Ocean, we demonstrate that strong bottom-trapped flows over sloping topography can lose their energy to near-inertial waves (NIWs) generated via the adjustment of the bottom Ekman layer. The NIWs with near-resonant frequencies corresponding to internal waves with propagation direction parallel to the topographic slope are observed. These waves are strongest in the BBL and have a correlation with the off-seamount subinertial flows largely attributed to the Ekman transport driven by the bottom-trapped anticyclonic circulation over the seamount. The bottom-intensified NIWs are observed to have dominant upward-propagating energy and hypothesized to be generated via Ekman flow–topography interactions in the BBL. Energy loss from the near-bottom flows to radiating NIWs (∼8 × 10−4 W m−2) is estimated to be substantially larger than that due to bottom drag dissipation (∼2 × 10−4 W m−2), suggesting the important role of internal-wave generation via the Ekman transport adjustment in damping the subinertial flows over the sloping seafloor.
Significance Statement
Dissipation of geostrophic currents and eddies via the oceanic bottom boundary layer (BBL) plays an important role in modulating the global oceanic mechanical energy budget. The bottom drag has long been considered a key process in inducing such dissipation, but it is recently suggested to be less effective at a sloping bottom. This study suggests another potentially important mechanism that can remove energy from the geostrophic flows at the sloping bottom. This mechanism is depicted as the resonant generation of near-inertial internal waves by near-bottom flows. Generation of internal waves due to near-bottom flows over sloping topography should be common in the global ocean and therefore have significant implications for the oceanic mechanical energy budget.
Abstract
The bottom boundary layer (BBL) contributes significantly to the global energy dissipation of low-frequency flows in the abyssal ocean, but how this dissipation occurs remains poorly understood. Using in situ data collected near the BBL at an abyssal seamount in the western Pacific Ocean, we demonstrate that strong bottom-trapped flows over sloping topography can lose their energy to near-inertial waves (NIWs) generated via the adjustment of the bottom Ekman layer. The NIWs with near-resonant frequencies corresponding to internal waves with propagation direction parallel to the topographic slope are observed. These waves are strongest in the BBL and have a correlation with the off-seamount subinertial flows largely attributed to the Ekman transport driven by the bottom-trapped anticyclonic circulation over the seamount. The bottom-intensified NIWs are observed to have dominant upward-propagating energy and hypothesized to be generated via Ekman flow–topography interactions in the BBL. Energy loss from the near-bottom flows to radiating NIWs (∼8 × 10−4 W m−2) is estimated to be substantially larger than that due to bottom drag dissipation (∼2 × 10−4 W m−2), suggesting the important role of internal-wave generation via the Ekman transport adjustment in damping the subinertial flows over the sloping seafloor.
Significance Statement
Dissipation of geostrophic currents and eddies via the oceanic bottom boundary layer (BBL) plays an important role in modulating the global oceanic mechanical energy budget. The bottom drag has long been considered a key process in inducing such dissipation, but it is recently suggested to be less effective at a sloping bottom. This study suggests another potentially important mechanism that can remove energy from the geostrophic flows at the sloping bottom. This mechanism is depicted as the resonant generation of near-inertial internal waves by near-bottom flows. Generation of internal waves due to near-bottom flows over sloping topography should be common in the global ocean and therefore have significant implications for the oceanic mechanical energy budget.
Abstract
The ice shelf–ocean boundary current has an important control on heat delivery to the base of an ice shelf. Climate and regional models that include a representation of ice shelf cavities often use a coarse grid, and results have a strong dependence on resolution near the ice shelf–ocean interface. This study models the ice shelf–ocean boundary current with a nonhydrostatic z-level configuration at turbulence-permitting resolution (1 m). The z-level model performs well when compared against state-of-the-art large-eddy simulations, showing its capability in representing the correct physics. We show that theoretical results from a one-dimensional model with parameterized turbulence reproduce the z-level model results to a good degree, indicating possible utility as a turbulence closure. The one-dimensional model evolves to a state of marginal instability, and we use the z-level model to demonstrate how this is represented in three dimensions. Instabilities emerge that regulate the strength of the pycnocline and coexist with persistent Ekman rolls, which are identified prior to the flow becoming intermittently unstable. When resolution of the z-level model is degraded to understand the gridscale dependencies, the degradation is dominated by the established problem of excessive numerical diffusion. We show that at intermediate resolutions (2–4 m), the boundary layer structure can be partially recovered by tuning diffusivities. Last, we compare replacing prescribed melting with interactive melting that is dependent on the local ocean conditions. Interactive melting results in a feedback such that the system evolves more slowly, which is exaggerated at lower resolution.
Abstract
The ice shelf–ocean boundary current has an important control on heat delivery to the base of an ice shelf. Climate and regional models that include a representation of ice shelf cavities often use a coarse grid, and results have a strong dependence on resolution near the ice shelf–ocean interface. This study models the ice shelf–ocean boundary current with a nonhydrostatic z-level configuration at turbulence-permitting resolution (1 m). The z-level model performs well when compared against state-of-the-art large-eddy simulations, showing its capability in representing the correct physics. We show that theoretical results from a one-dimensional model with parameterized turbulence reproduce the z-level model results to a good degree, indicating possible utility as a turbulence closure. The one-dimensional model evolves to a state of marginal instability, and we use the z-level model to demonstrate how this is represented in three dimensions. Instabilities emerge that regulate the strength of the pycnocline and coexist with persistent Ekman rolls, which are identified prior to the flow becoming intermittently unstable. When resolution of the z-level model is degraded to understand the gridscale dependencies, the degradation is dominated by the established problem of excessive numerical diffusion. We show that at intermediate resolutions (2–4 m), the boundary layer structure can be partially recovered by tuning diffusivities. Last, we compare replacing prescribed melting with interactive melting that is dependent on the local ocean conditions. Interactive melting results in a feedback such that the system evolves more slowly, which is exaggerated at lower resolution.
Abstract
Examining eddy–mean flow interactions in western boundary currents is crucial for understanding the mechanisms of mesoscale eddy generation and the role of eddies in the large-scale circulation. However, this analysis is lacking in the East Australian Current (EAC) system. Here we show the detailed three-dimensional structure of the eddy–mean flow interactions and energy budget in the EAC system. The energy reservoirs and conversions are greatest in the upper 500 m, with complex vertical structures. Strong mean kinetic energy is confined within a narrow band (24.5°–32.5°S) in the EAC jet. Most energy is contained in the eddy fields instead of the mean flow in the EAC typical separation and extension regions (south of 32.5°S). Strong barotropic instability is the primary source of eddy kinetic energy north of 36°S, while baroclinic instability dominates the eddy kinetic energy production in the EAC southern extension, which peaks in the subsurface. The mean flow transfers 5.22 GW of kinetic energy and 3.33 GW of available potential energy to the eddy field in the EAC typical separation region. The largest conversion term is from available potential energy conversion from the mean flow to the eddy field through baroclinic instability, dominating between 29° and 35.5°S. Nonlocal eddy–mean flow interactions also play a role in the energy exchange between the mean flow and the eddy fields. This study provides the mean state of the eddy–mean flow interactions in the EAC system, paving the way for further studies exploring seasonal and interannual variability and provides a baseline for assessing the impact of environmental change.
Abstract
Examining eddy–mean flow interactions in western boundary currents is crucial for understanding the mechanisms of mesoscale eddy generation and the role of eddies in the large-scale circulation. However, this analysis is lacking in the East Australian Current (EAC) system. Here we show the detailed three-dimensional structure of the eddy–mean flow interactions and energy budget in the EAC system. The energy reservoirs and conversions are greatest in the upper 500 m, with complex vertical structures. Strong mean kinetic energy is confined within a narrow band (24.5°–32.5°S) in the EAC jet. Most energy is contained in the eddy fields instead of the mean flow in the EAC typical separation and extension regions (south of 32.5°S). Strong barotropic instability is the primary source of eddy kinetic energy north of 36°S, while baroclinic instability dominates the eddy kinetic energy production in the EAC southern extension, which peaks in the subsurface. The mean flow transfers 5.22 GW of kinetic energy and 3.33 GW of available potential energy to the eddy field in the EAC typical separation region. The largest conversion term is from available potential energy conversion from the mean flow to the eddy field through baroclinic instability, dominating between 29° and 35.5°S. Nonlocal eddy–mean flow interactions also play a role in the energy exchange between the mean flow and the eddy fields. This study provides the mean state of the eddy–mean flow interactions in the EAC system, paving the way for further studies exploring seasonal and interannual variability and provides a baseline for assessing the impact of environmental change.
Abstract
This study reveals the role of the tropical Atlantic variability in modulating barrier layer thickness (BLT) in peak seasons. Based on reanalysis data during 1980–2016, statistical and dynamical analyses are performed to investigate the mechanism of BLT variability associated with the tropical Atlantic modes. The regions with significant correlation between BLT and tropical Atlantic modes are located on the northwest and southeast coasts of the tropical Atlantic, which are consistent with BLT maximum variability regions. In boreal spring, BLT decreases in the northwest because less latent heat release affected by weak trade wind related to the Atlantic meridional mode (AMM) shoals the isothermal layer depth (ITLD). In the south equatorial Atlantic, deepened mixed layer depth (MLD) is controlled by the decreasing freshwater input brought by a northward shift of the intertropical convergence zone (ITCZ) and further leads to a thinner barrier layer (BL). However, a shoaling MLD appears in the north equatorial Atlantic, which results from excessive freshwater input, causing a thick BL there. In boreal summer, positive runoff anomaly caused by the Atlantic equatorial mode (AEM) leads to upper warming of the tropical northwest Atlantic and a shallowing ITLD, favoring a thinner BL there. However, a southward shift of ITCZ brings more freshwater into the south equatorial Atlantic, inducing a shallowing MLD as well as a thicker BL. AEM-driven horizontal heat advection of the south equatorial current contributes to a thick ITLD in the central southern tropical Atlantic and thus increases BLT.
Significance Statement
This research aims to reveal how the tropical Atlantic meridional and equatorial interannual climatic modes affect barrier layer thickness (BLT). These two climate modes can affect the wind field, ocean current, and precipitation through air–sea interaction processes, and further affect mixing, heat–salt transport, and stratification in the upper ocean and thus BLT. This finding is important because the barrier layer restricts the exchange of heat, momentum, mass, and nutrients between the mixed layer and the thermocline, thereby impacting local and remote weather events, the ecological environment, and the climate. Our results provide guidance for interpreting the interannual variability of BLT in the tropical Atlantic.
Abstract
This study reveals the role of the tropical Atlantic variability in modulating barrier layer thickness (BLT) in peak seasons. Based on reanalysis data during 1980–2016, statistical and dynamical analyses are performed to investigate the mechanism of BLT variability associated with the tropical Atlantic modes. The regions with significant correlation between BLT and tropical Atlantic modes are located on the northwest and southeast coasts of the tropical Atlantic, which are consistent with BLT maximum variability regions. In boreal spring, BLT decreases in the northwest because less latent heat release affected by weak trade wind related to the Atlantic meridional mode (AMM) shoals the isothermal layer depth (ITLD). In the south equatorial Atlantic, deepened mixed layer depth (MLD) is controlled by the decreasing freshwater input brought by a northward shift of the intertropical convergence zone (ITCZ) and further leads to a thinner barrier layer (BL). However, a shoaling MLD appears in the north equatorial Atlantic, which results from excessive freshwater input, causing a thick BL there. In boreal summer, positive runoff anomaly caused by the Atlantic equatorial mode (AEM) leads to upper warming of the tropical northwest Atlantic and a shallowing ITLD, favoring a thinner BL there. However, a southward shift of ITCZ brings more freshwater into the south equatorial Atlantic, inducing a shallowing MLD as well as a thicker BL. AEM-driven horizontal heat advection of the south equatorial current contributes to a thick ITLD in the central southern tropical Atlantic and thus increases BLT.
Significance Statement
This research aims to reveal how the tropical Atlantic meridional and equatorial interannual climatic modes affect barrier layer thickness (BLT). These two climate modes can affect the wind field, ocean current, and precipitation through air–sea interaction processes, and further affect mixing, heat–salt transport, and stratification in the upper ocean and thus BLT. This finding is important because the barrier layer restricts the exchange of heat, momentum, mass, and nutrients between the mixed layer and the thermocline, thereby impacting local and remote weather events, the ecological environment, and the climate. Our results provide guidance for interpreting the interannual variability of BLT in the tropical Atlantic.
