Search Results
You are looking at 11 - 14 of 14 items for
- Author or Editor: Leonel Romero x
- Refine by Access: All Content x
Abstract
A large-eddy simulation (LES) model, which adopts wave-averaged equations with vortex force, is used to investigate Langmuir turbulence and ocean boundary layer (OBL) dynamics in high-wind hurricane conditions. The temporally evolving spatially asymmetric wind and wave Stokes drift velocity imposed in the LES are generated by a spectral wave prediction model adapted to Hurricane Frances traveling at a speed of 5.5 m s−1. The potency of Langmuir turbulence depends on the turbulent Langmuir number, the wind–Stokes drift alignment, and the depth scale of the Stokes profile Ds relative to the OBL depth h. At the time of maximum winds, large-scale vigorous coherent cells develop on the right-hand side of the storm under the inertially rotating winds; the Stokes drift velocity is well tuned to the surface winds. Much weaker cells develop on the left-hand side of the storm, partly because of reduced Stokes production. With misaligned winds and waves the vertical momentum fluxes can be counter to the gradient of Stokes drift, and the cell orientation tracks the direction of the mean Lagrangian shear. The entrainment flux is increased by 20% and the sea surface temperature is 0.25 K cooler on the right-hand side of the storm in the presence of Langmuir turbulence. Wave effects impact entrainment when the ratio Ds /|h| > 0.75. Because of wind–wave asymmetry Langmuir cells add quantitatively to the left–right asymmetry already understood for hurricanes due to resonance. And the transient evolution of the OBL cannot be understood simply in terms of equilibrium snapshots.
Abstract
A large-eddy simulation (LES) model, which adopts wave-averaged equations with vortex force, is used to investigate Langmuir turbulence and ocean boundary layer (OBL) dynamics in high-wind hurricane conditions. The temporally evolving spatially asymmetric wind and wave Stokes drift velocity imposed in the LES are generated by a spectral wave prediction model adapted to Hurricane Frances traveling at a speed of 5.5 m s−1. The potency of Langmuir turbulence depends on the turbulent Langmuir number, the wind–Stokes drift alignment, and the depth scale of the Stokes profile Ds relative to the OBL depth h. At the time of maximum winds, large-scale vigorous coherent cells develop on the right-hand side of the storm under the inertially rotating winds; the Stokes drift velocity is well tuned to the surface winds. Much weaker cells develop on the left-hand side of the storm, partly because of reduced Stokes production. With misaligned winds and waves the vertical momentum fluxes can be counter to the gradient of Stokes drift, and the cell orientation tracks the direction of the mean Lagrangian shear. The entrainment flux is increased by 20% and the sea surface temperature is 0.25 K cooler on the right-hand side of the storm in the presence of Langmuir turbulence. Wave effects impact entrainment when the ratio Ds /|h| > 0.75. Because of wind–wave asymmetry Langmuir cells add quantitatively to the left–right asymmetry already understood for hurricanes due to resonance. And the transient evolution of the OBL cannot be understood simply in terms of equilibrium snapshots.
Abstract
Monin–Obukhov similarity theory is applied to the surface layer of large-eddy simulations (LES) of deep Southern Ocean boundary layers. Observations from the Southern Ocean Flux Station provide a wide range of wind, buoyancy, and wave (Stokes drift) forcing. Two No-Stokes LES are used to determine the extent of the ocean surface layer and to adapt the nondimensional momentum and buoyancy gradients, as functions of the stability parameter. Stokes-forced LES are used to modify this parameter for wave effects, then to formulate dependencies of Stokes similarity functions on a Stokes parameter ξ. To account for wind-wave misalignment, the dimensional analysis is extended with two independent variables, namely, the production of turbulent kinetic energy in the surface layer due to Stokes shear and the total production, so that their ratio gives ξ. Stokes forcing is shown to reduce vertical shear more than stratification, and to enhance viscosity and diffusivity by factors up to 5.8 and 4.0, respectively, such that the Prandtl number can exceed unity. A practical parameterization is developed for ξ in terms of the meteorological forcing plus a Stokes drift profile, so that the Stokes and stability similarity functions can be combined to give turbulent velocity scales. These scales for both viscosity and diffusivity are evaluated against the LES, and the correlations are nearly 0.97. The benefit of calculating Stokes drift profiles from directional wave spectra is demonstrated by similarly evaluating three alternatives.
Abstract
Monin–Obukhov similarity theory is applied to the surface layer of large-eddy simulations (LES) of deep Southern Ocean boundary layers. Observations from the Southern Ocean Flux Station provide a wide range of wind, buoyancy, and wave (Stokes drift) forcing. Two No-Stokes LES are used to determine the extent of the ocean surface layer and to adapt the nondimensional momentum and buoyancy gradients, as functions of the stability parameter. Stokes-forced LES are used to modify this parameter for wave effects, then to formulate dependencies of Stokes similarity functions on a Stokes parameter ξ. To account for wind-wave misalignment, the dimensional analysis is extended with two independent variables, namely, the production of turbulent kinetic energy in the surface layer due to Stokes shear and the total production, so that their ratio gives ξ. Stokes forcing is shown to reduce vertical shear more than stratification, and to enhance viscosity and diffusivity by factors up to 5.8 and 4.0, respectively, such that the Prandtl number can exceed unity. A practical parameterization is developed for ξ in terms of the meteorological forcing plus a Stokes drift profile, so that the Stokes and stability similarity functions can be combined to give turbulent velocity scales. These scales for both viscosity and diffusivity are evaluated against the LES, and the correlations are nearly 0.97. The benefit of calculating Stokes drift profiles from directional wave spectra is demonstrated by similarly evaluating three alternatives.
Abstract
Knowledge of horizontal relative dispersion in nearshore oceans is important for many applications including the transport and fate of pollutants and the dynamics of nearshore ecosystems. Two-particle dispersion statistics are calculated from millions of synthetic particle trajectories from high-resolution numerical simulations of the Southern California Bight. The model horizontal resolution of 250 m allows the investigation of the two-particle dispersion, with an initial pair separation of 500 m. The relative dispersion is characterized with respect to the coastal geometry, bathymetry, eddy kinetic energy, and the relative magnitudes of strain and vorticity. Dispersion is dominated by the submesoscale, not by tides. In general, headlands are more energetic and dispersive than bays. Relative diffusivity estimates are smaller and more anisotropic close to shore. Farther from shore, the relative diffusivity increases and becomes less anisotropic, approaching isotropy ~10 km from the coast. The degree of anisotropy of the relative diffusivity is qualitatively consistent with that for eddy kinetic energy. The total relative diffusivity as a function of pair separation distance R is on average proportional to R 5/4. Additional Lagrangian experiments at higher horizontal numerical resolution confirmed the robustness of these results. Structures of large vorticity are preferably elongated and aligned with the coastline nearshore, which may limit cross-shelf dispersion. The results provide useful information for the design of subgrid-scale mixing parameterizations as well as quantifying the transport and dispersal of dissolved pollutants and biological propagules.
Abstract
Knowledge of horizontal relative dispersion in nearshore oceans is important for many applications including the transport and fate of pollutants and the dynamics of nearshore ecosystems. Two-particle dispersion statistics are calculated from millions of synthetic particle trajectories from high-resolution numerical simulations of the Southern California Bight. The model horizontal resolution of 250 m allows the investigation of the two-particle dispersion, with an initial pair separation of 500 m. The relative dispersion is characterized with respect to the coastal geometry, bathymetry, eddy kinetic energy, and the relative magnitudes of strain and vorticity. Dispersion is dominated by the submesoscale, not by tides. In general, headlands are more energetic and dispersive than bays. Relative diffusivity estimates are smaller and more anisotropic close to shore. Farther from shore, the relative diffusivity increases and becomes less anisotropic, approaching isotropy ~10 km from the coast. The degree of anisotropy of the relative diffusivity is qualitatively consistent with that for eddy kinetic energy. The total relative diffusivity as a function of pair separation distance R is on average proportional to R 5/4. Additional Lagrangian experiments at higher horizontal numerical resolution confirmed the robustness of these results. Structures of large vorticity are preferably elongated and aligned with the coastline nearshore, which may limit cross-shelf dispersion. The results provide useful information for the design of subgrid-scale mixing parameterizations as well as quantifying the transport and dispersal of dissolved pollutants and biological propagules.
Abstract
The observation-based source terms available in the third-generation wave model WAVEWATCH III (i.e., the ST6 package for parameterizations of wind input, wave breaking, and swell dissipation terms) are recalibrated and verified against a series of academic and realistic simulations, including the fetch/duration-limited test, a Lake Michigan hindcast, and a 1-yr global hindcast. The updated ST6 not only performs well in predicting commonly used bulk wave parameters (e.g., significant wave height and wave period) but also yields a clearly improved estimation of high-frequency energy level (in terms of saturation spectrum and mean square slope). In the duration-limited test, we investigate the modeled wave spectrum in a detailed way by introducing spectral metrics for the tail and the peak of the omnidirectional wave spectrum and for the directionality of the two-dimensional frequency–direction spectrum. The omnidirectional frequency spectrum E(f) from the recalibrated ST6 shows a clear transition behavior from a power law of approximately f −4 to a power law of about f −5, comparable to previous field studies. Different solvers for nonlinear wave interactions are applied with ST6, including the Discrete Interaction Approximation (DIA), the more expensive Generalized Multiple DIA (GMD), and the very expensive exact solutions [using the Webb–Resio–Tracy method (WRT)]. The GMD-simulated E(f) is in excellent agreement with that from WRT. Nonetheless, we find the peak of E(f) modeled by the GMD and WRT appears too narrow. It is also shown that in the 1-yr global hindcast, the DIA-based model overestimates the low-frequency wave energy (wave period T > 16 s) by 90%. Such model errors are reduced significantly by the GMD to ~20%.
Abstract
The observation-based source terms available in the third-generation wave model WAVEWATCH III (i.e., the ST6 package for parameterizations of wind input, wave breaking, and swell dissipation terms) are recalibrated and verified against a series of academic and realistic simulations, including the fetch/duration-limited test, a Lake Michigan hindcast, and a 1-yr global hindcast. The updated ST6 not only performs well in predicting commonly used bulk wave parameters (e.g., significant wave height and wave period) but also yields a clearly improved estimation of high-frequency energy level (in terms of saturation spectrum and mean square slope). In the duration-limited test, we investigate the modeled wave spectrum in a detailed way by introducing spectral metrics for the tail and the peak of the omnidirectional wave spectrum and for the directionality of the two-dimensional frequency–direction spectrum. The omnidirectional frequency spectrum E(f) from the recalibrated ST6 shows a clear transition behavior from a power law of approximately f −4 to a power law of about f −5, comparable to previous field studies. Different solvers for nonlinear wave interactions are applied with ST6, including the Discrete Interaction Approximation (DIA), the more expensive Generalized Multiple DIA (GMD), and the very expensive exact solutions [using the Webb–Resio–Tracy method (WRT)]. The GMD-simulated E(f) is in excellent agreement with that from WRT. Nonetheless, we find the peak of E(f) modeled by the GMD and WRT appears too narrow. It is also shown that in the 1-yr global hindcast, the DIA-based model overestimates the low-frequency wave energy (wave period T > 16 s) by 90%. Such model errors are reduced significantly by the GMD to ~20%.
