Seasonal Variability of The Mixed Layer Depth in The Sulawesi Sea

The Sulawesi Sea is a semi-enclosed basin located in the Indonesian Seas and considered as the one of location in the west route of Indonesian Throughflow (ITF). There is less attention on the mixed layer depth investigation in the Sulawesi Sea. Concerning that the mixed layer plays an important role in influencing the ocean in air-sea interaction and affects biological activity, the estimation of mixed layer depth (MLD) in the Sulawesi Sea is important. Seasonal variation of the mixed layer in the Sulawesi Sea between 115°-125°E and 0°-8°N is estimated by using World Ocean Atlas 2013. Forcing elements on the mixed layer in terms of surface-forced turbulent mixing from mechanical forcing of wind stress and buoyancy forcing (from heat flux as well as freshwater flux) in the Sulawesi Sea is provided by using a reanalysis dataset. The MLD is estimated directly on grid profiles with interpolated levels based on chosen density fixed criterion of 0.03 kg.m^-3 and temperature criterion of 0.5°C difference from the surface. The results show that mixed layer depth in the Sulawesi Sea varies both spatially and temporally. Generally, the deepest MLD was occurred during the southwest monsoon (JJA), and the lowest MLD was occurred during the first transition (MAM) and second transition monsoon (SON). Strengthening and weakening MLD are influenced by mechanical forcing from wind stress and buoyancy flux. In the Sulawesi Sea, the mixed layer deepening coincides with the occurrence of a maximum in wind stress, and low buoyancy flux at the surface. This condition is the opposite when mixed layer shallowing occurs.

Under pressure: turbulent plumes in a uniform crossflow

Direct numerical simulation is used to investigate the integral behaviour of buoyant plumes subjected to a uniform crossflow that are infinitely lazy at the source. Neither a plume trajectory defined by the centre of mass of the plume $z_c$ nor a trajectory defined by the central streamline $z_{U}$ is aligned with the average streamlines inside the plume. Both $z_c$ and $z_{U}$ are shown to correlate with field lines of the total buoyancy flux, which implies that a model for the vertical turbulent buoyancy flux is required to faithfully predict the plume angle. A study of the volume conservation equation shows that entrainment due to incorporation of ambient fluid with non-zero velocity due to the increase in the surface area (the Leibniz term) is the dominant entrainment mechanism in strong crossflows. The data indicate that pressure differences between the top and bottom of the plume play a leading role in the evolution of the horizontal and vertical momentum balances and are crucial for appropriately modelling plume rise. By direct parameterisation of the vertical buoyancy flux, the entrainment and the pressure, an integral plume model is developed which is in good agreement with the simulations for sufficiently strong crossflow. A perturbation expansion shows that the current model is an intermediate-range model valid for downstream distances up to $100\ell _b$ – $1000 \ell _b$ , where $\ell _b$ is the buoyancy length scale based on the flow speed and plume buoyancy flux.

The diurnal cycle of entrainment and detrainment in LES of the Southern Ocean driven by observed surface fluxes and waves

Empirical rules for both entrainment and detrainment are developed from LES of the Southern Ocean boundary layer when the turbulence, stratification and shear cannot be assumed to be in equilibrium with diurnal variability in surface flux and wave (Stokes drift) forcing. A major consequence is the failure of down-gradient eddy viscosity, which becomes more serious with Stokes drift and is overcome by relating the angle between the stress and shear vectors to the orientations of Lagrangian shear to the surface and of local Eulerian shear over five meters. Thus, the momentum flux can be parameterized as a stress magnitude and this empirical direction. In addition, the response of a deep boundary layer to sufficiently strong diurnal heating includes boundary layer collapse and the subsequent growth of a morning boundary layer, whose depth is empirically related to the time history of the forcing, as are both morning detrainment and afternoon entrainment into weak diurnal stratification. Below the boundary layer, detrainment rules give the maximum buoyancy flux and its depth, as well a specific stress direction. Another rule relates both afternoon and night-time entrainment depth and buoyancy flux to surface layer turbulent kinetic energy production integrals. These empirical relationships are combined with rules for boundary layer transport to formulate two parameterizations; one based on eddy diffusivity and viscosity profiles and another on flux profiles of buoyancy and of stress magnitude. Evaluations against LES fluxes show the flux profiles to be more representative of the diurnal cycle, especially with Stokes drift.

Wind- and Wave-driven Ocean Surface Boundary Layer in a Frontal Zone: Roles of Submesoscale Eddies and Ekman-Stokes Transport

Large-eddy simulations are used to investigate the influence of a horizontal frontal zone, represented by a stationary uniform background horizontal temperature gradient, on the wind- and wave-driven ocean surface boundary layers. In a frontal zone, the temperature structure, the ageostrophic mean horizontal current, and the turbulence in the ocean surface boundary layer all change with the relative angle among the wind and the front. The net heating and cooling of the boundary layer could be explained by the depth-integrated horizontal advective buoyancy flux, called the Ekman Buoyancy Flux (or the Ekman-Stokes Buoyancy Flux if wave effects are included). However, the detailed temperature profiles are also modulated by the depth-dependent advective buoyancy flux and submesoscale eddies. The surface current is deflected less (more) to the right of the wind and wave when the depth-integrated advective buoyancy flux cools (warms) the ocean surface boundary layer. Horizontal mixing is greatly enhanced by submesoscale eddies. The eddy-induced horizontal mixing is anisotropic and is stronger to the right of the wind direction. Vertical turbulent mixing depends on the superposition of the geostrophic and ageostrophic current, the depth-dependent advective buoyancy flux, and submesoscale eddies.

The Effect of Bottom – Generated Tidal Mixing on Tidally Pulsed River Plumes

The mixing of river plumes into the coastal ocean influences the fate of river-borne tracers over the inner-shelf, though the relative importance of mixing mechanisms under different environmental conditions is not fully understood. In particular, the contribution to plume mixing from bottom generated shear stresses, referred to as tidal mixing, is rarely considered important relative to frontal and stratified shear (interfacial) mixing in surface advected plumes. The effect of different mixing mechanisms is investigated numerically on an idealized, tidally pulsed river plume with varying river discharge and tidal amplitudes. Frontal, interfacial, and tidal mixing are quantified via a mixing energy budget to compare the relative importance of each to the overall buoyancy flux over one tide. Results indicate that tidal mixing can dominate the energy budget when the tidal mixing power exceeds that of the input buoyancy flux. This occurs when the non-dimensional number, RiE (the estuarine Richardson number divided by the mouth Rossby number), is generally less than 1. Tidal mixing accounts for between 60% and 90% of the net mixing when RiE < 1, with the largest contributions during large tides and low discharge. Interfacial mixing varies from 10% to 90% of total mixing and dominates the budget for high discharge events with relatively weaker tides (RiE > 1). Frontal mixing is always less than 10% of total mixing and never dominates the budget. This work is the first to show tidal mixing as an important mixing mechanism in surface advected river plumes.

Analysis of Energy Sources along the Kuroshio in the East of Taiwan Island and East China Sea

&lt;p&gt;The locations and generation mechanisms of energy sources in the Kuroshio were analyzed. The slope of the one-dimensional spectral energy density varies between -5/3 and -3 in the wavenumber range of 0.03-0.1 cpkm (wavelengths of approximately 209 to 63 km, respectively), indicating an inverse energy cascade in the Kuroshio; according to the steady-state energy evolution, an energy source which occurs at scale smaller than Rhines scale must be present. By analyzing the wavenumber-frequency spectrum, the period of higher kinetic energy (KE) is about 89-209 days and spatial scale is less than 0.03 cpkm. The locations of energy sources were identified with using the spectral energy transfer calculated by altimetry and model data. At the sea surface, the KE sources are mainly within 23.2&amp;#176;-25.2&amp;#176;N and 28&amp;#176;-30&amp;#176;N at less than 0.03 cpkm and 23.2&amp;#176;-23.6&amp;#176;N and 26&amp;#176;-30&amp;#176;N at 0.03-0.1 cpkm. The available potential energy (APE) sources are mainly within 22.2&amp;#176;-28&amp;#176;N and 28.6&amp;#176;-30&amp;#176;N at less than 0.03 cpkm and 29.2&amp;#176;-30&amp;#176;N at 0.03-0.1 cpkm. Wind stress and density differences (including buoyancy flux, temperature flux and salinity flux) are primarily responsible for the KE and APE sources, respectively. Beneath the sea surface, the energy sources are mainly above 400 m depth, and buoyancy flux plays a major role in the generation of energy sources. The energy cycle process can be summarized as follows: once an energy source is formed, to maintain a steady state, energy cascades (mainly inverse cascades) will be engendered.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;

Correlation between Buoyancy Flux, Dissipation and Potential Vorticity in Rotating Stratified Turbulence

We study in this paper the correlation between the buoyancy flux, the efficiency of energy dissipation and the linear and nonlinear components of potential vorticity, PV, a point-wise invariant of the Boussinesq equations, contrasting the three identified regimes of rotating stratified turbulence, namely wave-dominated, wave–eddy interactions and eddy-dominated. After recalling some of the main novel features of these flows compared to homogeneous isotropic turbulence, we specifically analyze three direct numerical simulations in the absence of forcing and performed on grids of 10243 points, one in each of these physical regimes. We focus in particular on the link between the point-wise buoyancy flux and the amount of kinetic energy dissipation and of linear and nonlinear PV. For flows dominated by waves, we find that the highest joint probability is for minimal kinetic energy dissipation (compared to the buoyancy flux), low dissipation efficiency and low nonlinear PV, whereas for flows dominated by nonlinear eddies, the highest correlation between dissipation and buoyancy flux occurs for weak flux and high localized nonlinear PV. We also show that the nonlinear potential vorticity is strongly correlated with high dissipation efficiency in the turbulent regime, corresponding to intermittent events, as observed in the atmosphere and oceans.

Diapycnal Transport near a Sloping Bottom Boundary

The diapycnal motion in the stratified ocean near a sloping bottom boundary is studied using analytical solutions from one-dimensional boundary layer theory. Bottom-intensification of the diapycnal mixing intensity ensures that in the stratified mixing layer (SML), where isopycnals are relatively flat, the diapycnal motion is downward toward denser fluid. In contrast, convergence of the diffusive buoyancy flux near the seafloor drives diapycnal upwelling in what we define as the bottom boundary layer (BBL). Much of the one-dimensional BBL is characterized by a stratification only slightly reduced from that in the SML because the maximum in the buoyancy flux at the top of the BBL, where the diapycnal velocity changes sign, must occur in well-stratified fluid. The diapycnal upwelling in the BBL is determined by variations not only in the magnitude of the buoyancy gradient but also in the curvature of isopycnals. The net diapycnal upwelling is concentrated in the bottom half of the BBL where the magnitude of the buoyancy gradient changes most rapidly. The curvature effect drives upwelling near the seafloor that only makes a significant contribution to the net upwelling for steep slopes. The structure of the diapycnal velocity in this stratified BBL differs from the case of a turbulent well-mixed BBL that has been assumed in some recent theoretical studies on bottom-intensified mixing. This work therefore extends recent theories in a way that should be more applicable to abyssal ocean observations where well-mixed BBLs are not common.

Upper mantle mush zones beneath low melt flux ocean island volcanoes: insights from Isla Floreana, Galápagos

The physicochemical characteristics of sub-volcanic magma storage regions have important implications for magma system dynamics and pre-eruptive behaviour. The architecture of magma storage regions located directly above high buoyancy flux mantle plumes (such as Kīlauea, Hawai’i and Fernandina, Galápagos) are relatively well understood. However, far fewer constraints exist on the nature of magma storage beneath ocean island volcanoes that are distal to the main zone of mantle upwelling or above low buoyancy flux plumes, despite these systems representing a substantial proportion of ocean island volcanism globally. To address this, we present a detailed petrological study of Isla Floreana in the Galápagos Archipelago, which lies at the periphery of the upwelling mantle plume and is thus characterised by an extremely low flux of magma into the lithosphere. Detailed in situ major and trace element analyses of crystal phases within exhumed cumulate xenoliths, lavas and scoria deposits, indicate that the erupted crystal cargo is dominated by disaggregated crystal-rich material (i.e., mush or wall rock). Trace element disequilibria between cumulus phases and erupted melts, as well as trace element zoning within the xenolithic clinopyroxenes, reveals that reactive porous flow (previously identified beneath mid-ocean ridges) is an important process of melt transport within crystal-rich magma storage regions. In addition, application of three petrological barometers reveal that the Floreana mush zones are located in the upper mantle, at a depth of 23.7 ± 5.1 km. Our barometric results are compared to recent studies of high melt flux volcanoes in the western Galápagos, and other ocean island volcanoes worldwide, and demonstrate that the flux of magma from the underlying mantle source represents a first-order control on the depth and physical characteristics of magma storage.