skip to main content


Title: Revised Estimates of Ocean Surface Drag in Strong Winds
Abstract

Air‐sea drag governs the momentum transfer between the atmosphere and the ocean and remains largely unknown in hurricane winds. We revisit the momentum budget and eddy covariance methods to estimate the surface drag coefficient in the laboratory. Our drag estimates agree with field measurements in low‐to‐moderate winds and previous laboratory measurements in hurricane‐force winds. The drag coefficient saturates at 2.6×10−3andU10≈25 m s−1, in agreement with previous laboratory results by Takagaki et al. (2012,). During our analysis, we discovered an error in the original source code used by Donelan et al. (2004,). We present the corrected data and describe the correction procedure. Although the correction to the data does not change the key finding of drag saturation in strong winds, its magnitude and wind speed threshold are significantly changed. Our findings emphasize the need for an updated and unified drag parameterization based on field and laboratory data.

 
more » « less
Award ID(s):
1745384
NSF-PAR ID:
10447472
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
47
Issue:
10
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    In strong winds, air flow detaches from the ocean surface in the lee of wave crests and creates a low‐pressure zone on the wave’s leeward face. The pressure difference between the wave’s rear and front face modulates the momentum input from wind to waves. Numerical wave models parameterize this effect using a so‐called sheltering coefficient. However, its value and dependence on wind speed are not well understood, particularly with background swell waves. To bridge this gap, we conducted laboratory experiments with winds up to Category 4 hurricane force blown over various mechanically generated wave conditions (pure wind sea, mixed waves with directional spreading, and monochromatic unidirectional waves) and measured the wind, waves, and stress at a sufficient frequency to resolve wind‐wave variability over the long‐wave phase. We analyze the results in the context of Jeffreys’s sheltering theory and find two regimes: (a) from low‐to‐moderate wind forcing (10 m s−1 < U10 < 33 m s−1), the aerodynamic sheltering increases with wind speed, consistent with previous studies; (b) in hurricane conditions (U10 > 33 m s−1), the aerodynamic sheltering decreases with wind at a rate depending on wave state. Further, we isolate the short wind waves from the longer paddle waves and find that the aerodynamic sheltering by longer waves leads to a phase‐dependent variability of the short wind‐waves’ local steepness, which is evidenced by the sheltering coefficient’s value derived from wind and wave measurements. Our results emphasize the need for further measurements of aerodynamic sheltering and improving its representation in models.

     
    more » « less
  2. This dataset contains three netcdf files that pertain to monthly, seasonal, and annual fields of surface wind stress, wind stress curl, and curl-derived upwelling velocities over the Northwest Atlantic (80-45W, 30-45N) covering a forty year period from 1980 to 2019. Six-hourly surface (10 m) wind speed components from the Japanese 55-year reanalysis (JRA-55; Kobayashi et al., 2015) were processed from 1980 to 2019 over a larger North Atlantic domain of 100W to 10E and 10N to 80N. Wind stress was computed using a modified step-wise formulation, originally based on (Gill, 1982) and a non-linear drag coefficient (Large and Pond, 1981), and later modified for low speeds (Trenberth et al., 1989). See Gifford (2023) for more details.   

    After the six-hourly zonal and meridional wind stresses were calculated, the zonal change in meridional stress (curlx) and the negative meridional change in zonal stress (curly) were found using NumPy’s gradient function in Python (Harris et al., 2020) over the larger North Atlantic domain (100W-10E, 10-80N). The curl (curlx + curly) over the study domain (80-45W, 10-80N) is then extracted, which maintain a constant order of computational accuracy in the interior and along the boundaries for the smaller domain in a centered-difference gradient calculation. 

    The monthly averages of the 6-hour daily stresses and curls were then computed using the command line suite climate data operators (CDO, Schulzweida, 2022) monmean function. The seasonal (3-month average) and annual averages (12-month average) were calculated in Python using the monthly fields with NumPy (NumPy, Harris et al., 2020). 

    Corresponding upwelling velocities at different time-scales were obtained from the respective curl fields and zonal wind stress by using the Ekman pumping equation of the study by Risien and Chelton (2008; page 2393). Please see Gifford (2023) for more details.   

    The files each contain nine variables that include longitude, latitude, time, zonal wind stress, meridional wind stress, zonal change in meridional wind stress (curlx), the negative meridional change in zonal wind stress (curly), total curl, and upwelling. Units of time begin in 1980 and are months, seasons (JFM etc.), and years to 2019. The longitude variable extends from 80W to 45W and latitude is 30N to 45N with uniform 1.25 degree resolution.  

    Units of stress are in Pascals, units of curl are in Pascals per meter, and upwelling velocity is described by centimeters per day. The spatial grid is a 29 x 13 longitude x latitude array. 

    Filenames: 

    monthly_windstress_wsc_upwelling.nc: 480 time steps from 80W to 45W and 30N to 45N.

    seasonal_windstress_wsc_upwelling.nc: 160 time steps from 80W to 45W and 30N to 45N.

    annual_windstress_wsc_upwelling.nc: 40 time steps from 80W to 45W and 30N to 45N.

    Please contact igifford@earth.miami.edu for any queries. {"references": ["Gifford, I.H., 2023. The Synchronicity of the Gulf Stream Free Jet and the Wind Induced Cyclonic Vorticity Pool. MS Thesis, University of Massachusetts Dartmouth. 75pp.", "Gill, A. E. (1982). Atmosphere-ocean dynamics (Vol. 30). Academic Press.", "Harris, C.R., Millman, K.J., van der Walt, S.J. et al. Array programming with NumPy. Nature 585, 357\u2013362 (2020). DOI: 10.1038/s41586-020-2649-2.", "Japan Meteorological Agency/Japan (2013), JRA-55: Japanese 55-year Reanalysis, Daily 3-Hourly and 6-Hourly Data, https://doi.org/10.5065/D6HH6H41, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory, Boulder, Colo. (Updated monthly.)", "Kobayashi, S., Ota, Y., Harada, Y., Ebita, A., Moriya, M., Onoda, H., Onogi, K., Kamahori, H., Kobayashi, C., Endo, H. and Miyaoka, K., 2015. The JRA-55 reanalysis: General specifications and basic characteristics.\u202fJournal of the Meteorological Society of Japan. Ser. II,\u202f93(1), pp.5-48.", "Large, W.G. and Pond, S., 1981. Open ocean momentum flux measurements in moderate to strong winds.\u202fJournal of physical oceanography,\u202f11(3), pp.324-336.", "Risien, C.M. and Chelton, D.B., 2008. A global climatology of surface wind and wind stress fields from eight years of QuikSCAT scatterometer data.\u202fJournal of Physical Oceanography,\u202f38(11), pp.2379-2413.", "Schulzweida, Uwe. (2022). CDO User Guide (2.1.0). Zenodo. https://doi.org/10.5281/zenodo.7112925.", "Trenberth, K.E., Large, W.G. and Olson, J.G., 1989. The effective drag coefficient for evaluating wind stress over the oceans.\u202fJournal of Climate,\u202f2(12), pp.1507-1516."]} 
    more » « less
  3. Abstract

    The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,Ljet∼ 1050erg s−1. We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,Ljet∼ 1052erg s−1, but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2.

     
    more » « less
  4. Abstract

    We present measurements of the rest-frame UV spectral slope,β, for a sample of 36 faint star-forming galaxies atz∼ 9–16 discovered in one of the deepest JWST NIRCam surveys to date, the Next Generation Deep Extragalactic Exploratory Public Survey. We use robust photometric measurements for UV-faint galaxies (down toMUV∼ −16), originally published in Leung et al., and measure values of the UV spectral slope via photometric power-law fitting to both the observed photometry and stellar population models obtained through spectral energy distribution (SED) fitting withBagpipes. We obtain a median and 68% confidence interval forβfrom photometric power-law fitting ofβPL=2.70.5+0.5and from SED fitting,βSED=2.30.1+0.2for the full sample. We show that when only two to three photometric detections are available, SED fitting has a lower scatter and reduced biases than photometric power-law fitting. We quantify this bias and find that after correction the medianβSED,corr=2.50.2+0.2. We measure physical properties for our galaxies withBagpipesand find that our faint (MUV=18.10.9+0.7) sample is low in mass (log[M*/M]=7.70.5+0.5), fairly dust-poor (Av=0.10.1+0.2mag), and modestly young (log[age]=7.80.8+0.2yr) with a median star formation rate oflog(SFR)=0.30.4+0.4Myr1. We find no strong evidence for ultrablue UV spectral slopes (β∼ −3) within our sample, as would be expected for exotically metal-poor (Z/Z< 10−3) stellar populations with very high Lyman continuum escape fractions. Our observations are consistent with model predictions that galaxies of these stellar masses atz∼ 9–16 should have only modestly low metallicities (Z/Z∼ 0.1–0.2).

     
    more » « less
  5. This study analyzes high-resolution ship data collected in the Gulf of Mexico during the Lagrangian Submesoscale Experiment (LASER) from January to February 2016 to produce the first reported measurements of dissipative heating in the explicitly nonhurricane atmospheric surface layer. Although typically computed from theory as a function of wind speed cubed, the dissipative heating directly estimated via the turbulent kinetic energy (TKE) dissipation rate is also presented. The dissipative heating magnitude agreed with a previous study that estimated the dissipative heating in the hurricane boundary layer using in situ aircraft data. Our observations that the 10-m neutral drag coefficient parameterized using TKE dissipation rate approaches zero slope as wind increases suggests that TKE dissipation and dissipative heating are constrained to a physical limit. Both surface-layer stability and sea state were observed to be important conditions influencing dissipative heating, with the stability determined via TKE budget terms and the sea state determined via wave steepness and age using direct shipboard measurements. Momentum and enthalpy fluxes used in the TKE budget are determined using the eddy-correlation method. It is found that the TKE dissipation rate and the dissipative heating are largest in a nonneutral atmospheric surface layer with a sea surface comprising steep wind sea and slow swell waves at a given surface wind speed, whereas the ratio of dissipative heating to enthalpy fluxes is largest in near-neutral stability where the turbulent vertical velocities are near zero.

     
    more » « less