-
Notifications
You must be signed in to change notification settings - Fork 10
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Stokes drift for constant momentum fluxes #131
Merged
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
cc @qingli411 |
glwagner
changed the title
Stokes drift for constant winds
Stokes drift for constant momentum fluxes
Mar 5, 2022
Closed
…MA/LESbrary.jl into glw/infinite-impulsive-storm
navidcy
approved these changes
Aug 29, 2022
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR adds a new type of
AbstractStokesDrift
for use withOceananigans
that's appropriate for constant air-sea momentum fluxes over an infinite or isolated patch of ocean. See CliMA/Oceananigans.jl#2290. for preliminary discussion. I've decided to implement this feature in LESbrary.jl because I think it's currently experimental nature means it's not yet appropriate for Oceananigans.jl.The basic idea is that our idealized cases in
IdealizedExperiments/three_layer_constant_fluxes.jl
are meant to represent a scenario in which either infinite or large-scale winds or cooling associated with constant air-sea fluxes act "instantaneously" on a quiescent patch of ocean.The surface wave field that's consistent with this scenario is one that develops in time in response to atmospheric winds that impose constant momentum fluxes into the ocean.
We have to make a choice about the initial sea state because there are many plausible sea states consistent with an initially quiescent ocean interior. This reflects the fact that the surface wave field is, in reality, partly non-local due to distant swell generation and propagation. Here I think we are best served by choosing a "near-quiescent" (not entirely quiescent, since this is a singularity of our model as we will see below) initial sea state as well. This biases our cases towards "calm before storm" physical scenarios, rather than "train of storm" physical scenarios. My initial hunch is that moving towards more complex forcing / sea states is better served with "realistic" cases, rather than attempting to describe all possible ocean conditions with an idealized setup.
Note that we can't really comment about whether atmospheric wind is constant, because air-sea momentum fluxes depend on both the surface wave field and the atmospheric wind, in general (despite the current state-of-the-art in bulk formulae). Yet we don't need to comment on this, because the only thing that matters for ocean dynamics (and boundary layer parameterization) is the momentum flux into the ocean. This is important because it means that our results are independent of how the air-sea momentum flux is parameterized, and thus a bit more robust to future developments and improvements in air-sea momentum flux parameterization.
The
ConstantFluxStokesDrift
modelIt has two components:
The Stokes drift profile
Some preliminary results are here:
which compare to a profile proposed by Breivik et al., 2016.
The main difference between Lenain and Pizzo 2020's parameterization and Breivik et al., 2016's is stronger shear near the surface.
This plot shows the dependence of the surface Stokes drift for different peak wavelengths:
It's being compared to a model proposed by Li et al. 2017 that assumes a constant peak wavenumber for a given momentum flux / atmospheric wind speed.
The peak wavenumber
By manipulating the fetch relationships in Lenain and Melville 2017 we come up with a formula for how the peak wavenumber decreases in time:
However, a few experiments with this formula suggest that we probably can't use it for the constant flux simulations. I think the reason is the non-realism of infinite, constant winds blowing in a constant direction. So either we need to come up with a more realistic model for momentum fluxes, or we need to make approximations with the surface wave model.
We should comment on how this compares to, say, a 0D simulation using wave watch III or something similar.