**Welcome to the official ADCIRCWiki site! The site is currently under construction, with limited information. ---**

# Difference between revisions of "Wind Stress"

Taylorgasher (talk | contribs) (renamed section) |
Taylorgasher (talk | contribs) (→Roughness Reductions: fixed math typos) |
||

Line 17: | Line 17: | ||

:<math>f=\left ( \frac{z_{0l}}{z_{0m}}\right ) ^{0.0706} \left ( \frac{\ln \frac{10}{z_{0l}}}{\ln \frac{10}{z_{0m}}} \right ) </math> | :<math>f=\left ( \frac{z_{0l}}{z_{0m}}\right ) ^{0.0706} \left ( \frac{\ln \frac{10}{z_{0l}}}{\ln \frac{10}{z_{0m}}} \right ) </math> | ||

for marine roughness length <math>z_{0m}</math> and reduced ("land") roughness length <math>z_{0l}</math>. Wind speed is then reduced as, | for marine roughness length <math>z_{0m}</math> and reduced ("land") roughness length <math>z_{0l}</math>. Wind speed is then reduced as, | ||

− | :<math>\mathbf{w}'=f\mathbf{w}=f | + | :<math>\mathbf{w}'=f\mathbf{w}=f [u,v]</math> |

for x- and y- wind vector components <math>u</math> and <math>v</math>. The marine roughness length is, | for x- and y- wind vector components <math>u</math> and <math>v</math>. The marine roughness length is, | ||

− | :<math> | + | :<math>z_{0m}=\frac{0.018}{g} c_d \left \Vert \mathbf{w} \right \| </math> |

for Charnock parameter <math>0.018</math>, drag coefficient <math>c_d</math> and acceleration due to gravity <math>g</math>. The wind drag coefficient is addressed below in [[#Converting Wind Velocity to Wind Stress|this section]]. As previously noted, the reduced ("land") roughness length <math>z_{0l}</math> is specified by the user via the [[Fort.13_file#Surface_Roughness|surface roughness]] nodal attribute. The fraction <math>f</math> is bounded on <math>[0,1]</math>, meaning the winds cannot be increased, nor change direction. | for Charnock parameter <math>0.018</math>, drag coefficient <math>c_d</math> and acceleration due to gravity <math>g</math>. The wind drag coefficient is addressed below in [[#Converting Wind Velocity to Wind Stress|this section]]. As previously noted, the reduced ("land") roughness length <math>z_{0l}</math> is specified by the user via the [[Fort.13_file#Surface_Roughness|surface roughness]] nodal attribute. The fraction <math>f</math> is bounded on <math>[0,1]</math>, meaning the winds cannot be increased, nor change direction. | ||

## Revision as of 13:59, 8 July 2020

When wind blows over the water, it exerts a shear stress at the water surface that transfers horizontal momentum vertically downward across the air–sea interface, driving circulation. In ADCIRC, wind stress is an input forcing term, with several different formats provided. See the `NWS`

parameter for available formats. In most cases, the exact wind stress to be applied to the model is not provided, therefore ADCIRC must determine how to convert a given wind speed to the actual stress applied at the ocean surface. This page covers the various aspects of this process, as well as the options available to the user.

## Contents

## Definition of Winds

The characteristics of wind forcing are often broken down in three ways:

- Whether the winds are considered to be over-water (termed "marine exposure") or over-land
- The elevation above the sea (or ground) surface of the winds
- The time-averaging (if any) that has been applied

ADCIRC generally expects 10-meter, 10-minute winds at their actual exposure, although the exact expectations vary depending on the input type. For instance, when Holland-type wind inputs are provided (e.g. `NWS=8`

or `NWS=20`

), the wind speed is expected to be the **1-minute maximum sustained wind** at 10 meters elevation. If marine-exposure winds are provided, then surface roughness reductions may be needed
^{[1]}
^{[2]}
. If winds are provided with a different averaging time, then an appropriate correction may be needed, though winds with averaging times of 10 to 60 minutes are generally considered to be quite similar; this is the so-called mesoscale gap. For recommendations on wind time-scale conversions not handled internally by ADCIRC, for tropical cyclones, see the WMO guidelines of Harper et al.
^{[3]}

## Roughness Reductions

Reductions in wind speed to convert to the appropriate exposure come from a logarithmic boundary layer formulation (see, e.g. ^{[1]} ^{[2]}) to determine a fraction to reduce the winds,

for marine roughness length and reduced ("land") roughness length . Wind speed is then reduced as,

for x- and y- wind vector components and . The marine roughness length is,

for Charnock parameter , drag coefficient and acceleration due to gravity . The wind drag coefficient is addressed below in this section. As previously noted, the reduced ("land") roughness length is specified by the user via the surface roughness nodal attribute. The fraction is bounded on , meaning the winds cannot be increased, nor change direction.

### Older Behavior

ADCIRC version: | < 54 |

Note that before version 54, there was a bug in this calculation. The mistake and its effects are addressed in this PDF document: __ ADD A LINK TO A PDF HERE YOOOOO__.

## Converting Wind Velocity to Wind Stress

In ADCIRC, four formulations are available to convert wind speeds to the wind stresses applied in the momentum equations. Although there are several ways to control this, users are generally encouraged to use the metControl namelist in the fort.15 file. The default drag formulation is the Garratt
^{[4]}
linear formula. An alternative for use with tropical cyclones is the Powell formulation,
^{[5]}
which varies drag by the sector of the tropical cyclone. When ice coverage is included in the model, a wind drag formulation that accounts for this effect should be used. By default, if ice coverage input data are supplied, ADCIRC uses a cubic function of ice coverage, termed the "IceCube" drag formulation. Lastly, the "swell" drag law option allows users to utilize SWAN's drag formulation when employing the coupled model.

In all cases, the actual wind drag coefficients determined by ADCIRC can be output to a fort.63-type file named winddrag.173. Output settings (file format, output start/end times, and output interval) match those of either the fort.63 or fort.73/74 files, *I think* fort.73/74. Outputting of this file is enabled by setting `outputWindDrag=.TRUE.`

in the `&metControl`

namelist of the fort.15 file's namelist section.

### Garratt Drag Formulation

This is the default wind drag formulation in ADCIRC. From Garratt (1977)^{[4]}, the formula is,

By default, ADCIRC puts an upper bound on the drag coefficient of . This upper bound `WindDragLimit`

can be changed via the metControl fort.15 namelist.

### Powell Drag Formulation

CONTRIBUTOR NEEDED

### IceCube Drag Formulation

CONTRIBUTOR NEEDED

## References

- ↑
^{1.0}^{1.1}Simiu, E., Scanlan, R.H., 1996. Wind effects on structures: fundamentals and applications to design, 3rd ed. ed. John Wiley, New York. - ↑
^{2.0}^{2.1}Simiu, E., Yeo, D., 2018. Wind effects on structures: modern structural design for wind, Fourth edition. ed. John Wiley & Sons, Hoboken, NJ. - ↑ Harper, B., Kepert, J., Ginger, J., 2010. Guidelines for converting between various wind averaging periods in tropical cyclone conditions (No. WMO/TD-No. 1555). WMO, Geneva, Switzerland.
- ↑
^{4.0}^{4.1}Garratt, J.R., 1977. Review of Drag Coefficients over Oceans and Continents. Mon. Wea. Rev. 105, 915–929. https://doi.org/10.1175/1520-0493(1977)105<0915:RODCOO>2.0.CO;2 - ↑ Powell, M.D., Vickery, P.J., Reinhold, T.A., 2003. Reduced drag coefficient for high wind speeds in tropical cyclones. Nature 422, 279–283. https://doi.org/10.1038/nature01481