Search Results
You are looking at 1 - 10 of 17 items for
- Author or Editor: Amala Mahadevan x
- Refine by Access: All Content x
Abstract
Within the pycnocline, where diapycnal mixing is suppressed, both the vertical movement (uplift) of isopycnal surfaces and upward motion along sloping isopycnals supply nutrients to the euphotic layer, but the relative importance of each of these mechanisms is unknown. We present a method for decomposing vertical velocity w into two components in a Lagrangian frame: vertical velocity along sloping isopycnal surfaces 
Abstract
Within the pycnocline, where diapycnal mixing is suppressed, both the vertical movement (uplift) of isopycnal surfaces and upward motion along sloping isopycnals supply nutrients to the euphotic layer, but the relative importance of each of these mechanisms is unknown. We present a method for decomposing vertical velocity w into two components in a Lagrangian frame: vertical velocity along sloping isopycnal surfaces
Abstract
Using a process study model, the effect of mixed layer submesoscale instabilities on the lateral mixing of passive tracers in the pycnocline is explored. Mixed layer eddies that are generated from the baroclinic instability of a front within the mixed layer are found to penetrate into the pycnocline leading to an eddying flow field that acts to mix properties laterally along isopycnal surfaces. The mixing of passive tracers released on such isopycnal surfaces is quantified by estimating the variance of the tracer distribution over time. The evolution of the tracer variance reveals that the flow undergoes three different turbulent regimes. The first regime, lasting about 3–4 days (about 5 inertial periods) exhibits near-diffusive behavior; dispersion of the tracer grows nearly linearly with time. In the second regime, which lasts for about 10 days (about 14 inertial periods), tracer dispersion exhibits exponential growth because of the integrated action of high strain rates created by the instabilities. In the third regime, tracer dispersion follows Richardson’s power law. The Nakamura effective diffusivity is used to study the role of individual dynamical filaments in lateral mixing. The filaments, which carry a high concentration of tracer, are characterized by the coincidence of large horizontal strain rate with large vertical vorticity. Within filaments, tracer is sheared without being dispersed, and consequently the effective diffusivity is small in filaments. While the filament centers act as barriers to transport, eddy fluxes are enhanced at the filament edges where gradients are large.
Abstract
Using a process study model, the effect of mixed layer submesoscale instabilities on the lateral mixing of passive tracers in the pycnocline is explored. Mixed layer eddies that are generated from the baroclinic instability of a front within the mixed layer are found to penetrate into the pycnocline leading to an eddying flow field that acts to mix properties laterally along isopycnal surfaces. The mixing of passive tracers released on such isopycnal surfaces is quantified by estimating the variance of the tracer distribution over time. The evolution of the tracer variance reveals that the flow undergoes three different turbulent regimes. The first regime, lasting about 3–4 days (about 5 inertial periods) exhibits near-diffusive behavior; dispersion of the tracer grows nearly linearly with time. In the second regime, which lasts for about 10 days (about 14 inertial periods), tracer dispersion exhibits exponential growth because of the integrated action of high strain rates created by the instabilities. In the third regime, tracer dispersion follows Richardson’s power law. The Nakamura effective diffusivity is used to study the role of individual dynamical filaments in lateral mixing. The filaments, which carry a high concentration of tracer, are characterized by the coincidence of large horizontal strain rate with large vertical vorticity. Within filaments, tracer is sheared without being dispersed, and consequently the effective diffusivity is small in filaments. While the filament centers act as barriers to transport, eddy fluxes are enhanced at the filament edges where gradients are large.
Abstract
The nonhydrostatic model with a free surface is numerically implemented in boundary-fitted curvilinear coordinates to model the mesoscale circulation in an ocean basin with natural topography. A semi-implicit numerical scheme is used, and the directional inhomogeneity in the elliptic equation for pressure is exploited to speed up the computation of its solution while using the multigrid method.
The model is used to simulate the circulation in the Gulf of Mexico. We observe the formation of the Loop Current and several eddies. The flow is very strongly controlled by the topography and our numerical experiments reveal that in the bottom layers, the flow along topographic contours is in the opposite direction of the anticyclonic circulation in the top layers.
Abstract
The nonhydrostatic model with a free surface is numerically implemented in boundary-fitted curvilinear coordinates to model the mesoscale circulation in an ocean basin with natural topography. A semi-implicit numerical scheme is used, and the directional inhomogeneity in the elliptic equation for pressure is exploited to speed up the computation of its solution while using the multigrid method.
The model is used to simulate the circulation in the Gulf of Mexico. We observe the formation of the Loop Current and several eddies. The flow is very strongly controlled by the topography and our numerical experiments reveal that in the bottom layers, the flow along topographic contours is in the opposite direction of the anticyclonic circulation in the top layers.
Abstract
The incompressibility and hydrostatic approximations that are traditionally used in large-scale oceanography to make the hydrodynamic equations more amenable to numerical integration result in the primitive equations. These are ill-posed in domains with open boundaries and hence not well-suited to mesoscale or regional modeling. Instead of using the hydrostatic approximation, the authors permit a greater deviation from hydrostatic balance than what exists in the ocean to obtain a system of equations that is well-posed with the specification of pointwise boundary conditions at open or solid boundaries. These equations, formulated with a free-surface boundary, model the mesoscale flow field accurately in all three-dimensions, even the vertical. It is essential to retain the vertical component of the Coriolis acceleration in the model since it is nonhydrostatic.
Abstract
The incompressibility and hydrostatic approximations that are traditionally used in large-scale oceanography to make the hydrodynamic equations more amenable to numerical integration result in the primitive equations. These are ill-posed in domains with open boundaries and hence not well-suited to mesoscale or regional modeling. Instead of using the hydrostatic approximation, the authors permit a greater deviation from hydrostatic balance than what exists in the ocean to obtain a system of equations that is well-posed with the specification of pointwise boundary conditions at open or solid boundaries. These equations, formulated with a free-surface boundary, model the mesoscale flow field accurately in all three-dimensions, even the vertical. It is essential to retain the vertical component of the Coriolis acceleration in the model since it is nonhydrostatic.
Abstract
Changes in axial atmospheric angular momentum M are related to zonal torques on the atmosphere, but studies reveal large imbalances between the estimated torques and M variations on seasonal timescales. The observed imbalances are commonly attributed to uncertainties in the torque estimates. One particularly important torque component at the seasonal period is that due to zonal wind stresses over the ocean T O . The uncertainties in T O are explored by calculating different multiyear time series based on surface wind products derived from passive and active microwave satellite data. The satellite-based T O are compared to available reanalysis products. Results indicate that there are indeed substantial uncertainties in the seasonal T O , and that these uncertainties are related mostly to the wind fields rather than to the particular parameterizations of the surface stress in the boundary layer. Regional analyses point to the need to improve knowledge of the wind fields over extensive areas of the ocean, particularly in many tropical and southern latitude regions. Resolving subweekly variability in surface winds is also found to be important when determining the seasonal cycle in T O . The current satellite-based T O estimates can lead to a better seasonal momentum budget, but results are tempered by the uncertain effects of gravity wave torque in that budget.
Abstract
Changes in axial atmospheric angular momentum M are related to zonal torques on the atmosphere, but studies reveal large imbalances between the estimated torques and M variations on seasonal timescales. The observed imbalances are commonly attributed to uncertainties in the torque estimates. One particularly important torque component at the seasonal period is that due to zonal wind stresses over the ocean T O . The uncertainties in T O are explored by calculating different multiyear time series based on surface wind products derived from passive and active microwave satellite data. The satellite-based T O are compared to available reanalysis products. Results indicate that there are indeed substantial uncertainties in the seasonal T O , and that these uncertainties are related mostly to the wind fields rather than to the particular parameterizations of the surface stress in the boundary layer. Regional analyses point to the need to improve knowledge of the wind fields over extensive areas of the ocean, particularly in many tropical and southern latitude regions. Resolving subweekly variability in surface winds is also found to be important when determining the seasonal cycle in T O . The current satellite-based T O estimates can lead to a better seasonal momentum budget, but results are tempered by the uncertain effects of gravity wave torque in that budget.
Abstract
A cluster of 45 drifters deployed in the Bay of Bengal is tracked for a period of four months. Pair dispersion statistics, from observed drifter trajectories and simulated trajectories based on surface geostrophic velocity, are analyzed as a function of drifter separation and time. Pair dispersion suggests nonlocal dynamics at submesoscales of 1–20 km, likely controlled by the energetic mesoscale eddies present during the observations. Second-order velocity structure functions and their Helmholtz decomposition, however, suggest local dispersion and divergent horizontal flow at scales below 20 km. This inconsistency cannot be explained by inertial oscillations alone, as has been reported in recent studies, and is likely related to other nondispersive processes that impact structure functions but do not enter pair dispersion statistics. At scales comparable to the deformation radius L D , which is approximately 60 km, we find dynamics in agreement with Richardson’s law and observe local dispersion in both pair dispersion statistics and second-order velocity structure functions.
Abstract
A cluster of 45 drifters deployed in the Bay of Bengal is tracked for a period of four months. Pair dispersion statistics, from observed drifter trajectories and simulated trajectories based on surface geostrophic velocity, are analyzed as a function of drifter separation and time. Pair dispersion suggests nonlocal dynamics at submesoscales of 1–20 km, likely controlled by the energetic mesoscale eddies present during the observations. Second-order velocity structure functions and their Helmholtz decomposition, however, suggest local dispersion and divergent horizontal flow at scales below 20 km. This inconsistency cannot be explained by inertial oscillations alone, as has been reported in recent studies, and is likely related to other nondispersive processes that impact structure functions but do not enter pair dispersion statistics. At scales comparable to the deformation radius L D , which is approximately 60 km, we find dynamics in agreement with Richardson’s law and observe local dispersion in both pair dispersion statistics and second-order velocity structure functions.
Abstract
While near-inertial waves are known to be generated by atmospheric storms, recent observations in the Kuroshio Front find intense near-inertial internal-wave shear along sloping isopycnals, even during calm weather. Recent literature suggests that spontaneous generation of near-inertial waves by frontal instabilities could represent a major sink for the subinertial quasigeostrophic circulation. An unforced three-dimensional 1-km-resolution model, initialized with the observed cross-Kuroshio structure, is used to explore this mechanism. After several weeks, the model exhibits growth of 10–100-km-scale frontal meanders, accompanied by O(10) mW m−2 spontaneous generation of near-inertial waves associated with readjustment of submesoscale fronts forced out of balance by mesoscale confluent flows. These waves have properties resembling those in the observations. However, they are reabsorbed into the model Kuroshio Front with no more than 15% dissipating or radiating away. Thus, spontaneous generation of near-inertial waves represents a redistribution of quasigeostrophic energy rather than a significant sink.
Abstract
While near-inertial waves are known to be generated by atmospheric storms, recent observations in the Kuroshio Front find intense near-inertial internal-wave shear along sloping isopycnals, even during calm weather. Recent literature suggests that spontaneous generation of near-inertial waves by frontal instabilities could represent a major sink for the subinertial quasigeostrophic circulation. An unforced three-dimensional 1-km-resolution model, initialized with the observed cross-Kuroshio structure, is used to explore this mechanism. After several weeks, the model exhibits growth of 10–100-km-scale frontal meanders, accompanied by O(10) mW m−2 spontaneous generation of near-inertial waves associated with readjustment of submesoscale fronts forced out of balance by mesoscale confluent flows. These waves have properties resembling those in the observations. However, they are reabsorbed into the model Kuroshio Front with no more than 15% dissipating or radiating away. Thus, spontaneous generation of near-inertial waves represents a redistribution of quasigeostrophic energy rather than a significant sink.
Abstract
While lee-wave generation has been argued to be a major sink for the 1-TW wind work on the ocean’s circulation, microstructure measurements in the Antarctic Circumpolar Currents find dissipation rates as much as an order of magnitude weaker than linear lee-wave generation predictions in bottom-intensified currents. Wave action conservation suggests that a substantial fraction of lee-wave radiation can be reabsorbed into bottom-intensified flows. Numerical simulations are conducted here to investigate generation, reabsorption, and dissipation of internal lee waves in a bottom-intensified, laterally confined jet that resembles a localized abyssal current over bottom topography. For the case of monochromatic topography with |kU 0| ≈ 0.9N, where k is the along-stream topographic wavenumber, |U 0| is the near-bottom flow speed, and N is the buoyancy frequency; Reynolds-decomposed energy conservation is consistent with linear wave action conservation predictions that only 14% of lee-wave generation is dissipated, with the bulk of lee-wave energy flux reabsorbed by the bottom-intensified flow. Thus, water column reabsorption needs to be taken into account as a possible mechanism for reducing the lee-wave dissipative sink for balanced circulation.
Abstract
While lee-wave generation has been argued to be a major sink for the 1-TW wind work on the ocean’s circulation, microstructure measurements in the Antarctic Circumpolar Currents find dissipation rates as much as an order of magnitude weaker than linear lee-wave generation predictions in bottom-intensified currents. Wave action conservation suggests that a substantial fraction of lee-wave radiation can be reabsorbed into bottom-intensified flows. Numerical simulations are conducted here to investigate generation, reabsorption, and dissipation of internal lee waves in a bottom-intensified, laterally confined jet that resembles a localized abyssal current over bottom topography. For the case of monochromatic topography with |kU 0| ≈ 0.9N, where k is the along-stream topographic wavenumber, |U 0| is the near-bottom flow speed, and N is the buoyancy frequency; Reynolds-decomposed energy conservation is consistent with linear wave action conservation predictions that only 14% of lee-wave generation is dissipated, with the bulk of lee-wave energy flux reabsorbed by the bottom-intensified flow. Thus, water column reabsorption needs to be taken into account as a possible mechanism for reducing the lee-wave dissipative sink for balanced circulation.
On Characterizing Ocean Kinematics from Surface DriftersAbstract
Horizontal kinematic properties, such as vorticity, divergence, and lateral strain rate, are estimated from drifter clusters using three approaches. At submesoscale horizontal length scales
Significance Statement
The purpose of this study is to provide insights and guidance for computing horizontal velocity gradients from clusters (i.e., three or more) of Lagrangian surface ocean drifters. The uncertainty in velocity gradient estimates depends strongly on the shape deformation of drifter clusters by the ocean currents. We propose criteria for drifter cluster length scales and aspect ratios to reduce uncertainties and develop ways of estimating the magnitude of the resulting errors. The findings are applied to a real ocean dataset from the Bay of Bengal.
Abstract
Horizontal kinematic properties, such as vorticity, divergence, and lateral strain rate, are estimated from drifter clusters using three approaches. At submesoscale horizontal length scales
Significance Statement
The purpose of this study is to provide insights and guidance for computing horizontal velocity gradients from clusters (i.e., three or more) of Lagrangian surface ocean drifters. The uncertainty in velocity gradient estimates depends strongly on the shape deformation of drifter clusters by the ocean currents. We propose criteria for drifter cluster length scales and aspect ratios to reduce uncertainties and develop ways of estimating the magnitude of the resulting errors. The findings are applied to a real ocean dataset from the Bay of Bengal.
Abstract
A new method is proposed for extrapolating subsurface velocity and density fields from sea surface density and sea surface height (SSH). In this, the surface density is linked to the subsurface fields via the surface quasigeostrophic (SQG) formalism, as proposed in several recent papers. The subsurface field is augmented by the addition of the barotropic and first baroclinic modes, whose amplitudes are determined by matching to the sea surface height (pressure), after subtracting the SQG contribution. An additional constraint is that the bottom pressure anomaly vanishes. The method is tested for three regions in the North Atlantic using data from a high-resolution numerical simulation. The decomposition yields strikingly realistic subsurface fields. It is particularly successful in energetic regions like the Gulf Stream extension and at high latitudes where the mixed layer is deep, but it also works in less energetic eastern subtropics. The demonstration highlights the possibility of reconstructing three-dimensional oceanic flows using a combination of satellite fields, for example, sea surface temperature (SST) and SSH, and sparse (or climatological) estimates of the regional depth-resolved density. The method could be further elaborated to integrate additional subsurface information, such as mooring measurements.
Abstract
A new method is proposed for extrapolating subsurface velocity and density fields from sea surface density and sea surface height (SSH). In this, the surface density is linked to the subsurface fields via the surface quasigeostrophic (SQG) formalism, as proposed in several recent papers. The subsurface field is augmented by the addition of the barotropic and first baroclinic modes, whose amplitudes are determined by matching to the sea surface height (pressure), after subtracting the SQG contribution. An additional constraint is that the bottom pressure anomaly vanishes. The method is tested for three regions in the North Atlantic using data from a high-resolution numerical simulation. The decomposition yields strikingly realistic subsurface fields. It is particularly successful in energetic regions like the Gulf Stream extension and at high latitudes where the mixed layer is deep, but it also works in less energetic eastern subtropics. The demonstration highlights the possibility of reconstructing three-dimensional oceanic flows using a combination of satellite fields, for example, sea surface temperature (SST) and SSH, and sparse (or climatological) estimates of the regional depth-resolved density. The method could be further elaborated to integrate additional subsurface information, such as mooring measurements.
