Search NASA⌕ Search

SEARCH · Search NASA

Results for “air-sea”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

North Pacific meridional mode has larger impacts on El Niño evolution than the March Madden-Julian Oscillation

The El Niño–Southern Oscillation (ENSO) is a key driver of global climate variability. Early-season westerly wind bursts (WWBs) have long been suggested to be important for ENSO evolution and diversity, with the Madden-Julian Oscillation (MJO) among the main sources of WWBs. However, MJO’s contribution to ENSO evolution has been difficult to quantify. Here, using an ensemble hindcast approach specifically designed to isolate internal atmospheric variability, we evaluate the influences of March MJO on subsequent ENSO development. Our results show that the March MJO, under favorable background conditions, by itself has limited impacts on ENSO due to weak equatorial air-sea coupling in spring. In comparison, the North Pacific Oscillation–induced meridional mode exerts a more sustained influence on ENSO evolution. A cyclonic circulation anomaly over Hawaii, associated with the Pacific-North American pattern, also plays a role. These findings suggest that March MJO activity alone may not be a reliable predictor for ENSO evolution, but underscore the importance of North Pacific atmospheric variability.

Science & Technology - Other Topics↗

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗

Mantaray: A Rust Package for Ray Tracing Ocean Surface Gravity Waves

Ocean surface gravity waves are an important component of air-sea interaction, influencing energy, momentum, and gas exchanges across the ocean-atmosphere interface. In specific applications such as refraction by ocean currents or bathymetry, ray tracing provides a computationally efficient way to gain insight into wave propagation. In this paper, we introduce Mantaray, an open-source software package implemented in Rust, with a Python interface, that solves the ray equations for ocean surface gravity waves. Mantaray is designed for performance, robustness, and ease of use. The package is modular to facilitate further development and can currently be applied to both idealized and realistic wave propagation problems (Fig. 1).

16 TIDAL AND WAVE POWER↗

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗

Mid-Atlantic Turbulence Intensity Model Data

The dataset is a subset of simulations carried out for the mid-Atlantic region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the Air-Sea Interaction Tower (ASIT) of Woods Hole Oceanographic Institution’s Martha's Vineyard Coastal Observatory are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing. Simulations are performed from February through June in 2020 and output every 10 minutes, which permits comparison with corresponding lidar observations.

17 WIND ENERGY↗

Mid-Atlantic Turbulence Intensity Data in Observational Space

The dataset archives observed and model-simulated turbulence intensity and meteorological profiles and timeseries at the Air-Sea Interaction Tower (ASIT) of Woods Hole Oceanographic Institution’s Martha’s Vineyard Coastal Observatory (MVCO). The observational data were measured by lidar and buoy deployed at ASIT. The simulated profiles and timeseries data are interpolated in time and/or space according to observations. Simulations were carried out for the mid-Atlantic region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles at the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing.

17 WIND ENERGY↗

The Impact of Model Resolution on the Southern Hemisphere in Community Climate System Model Version 4 Idealized Climate Simulations

Model resolution plays a large role in accurately simulating the Southern Hemisphere circulation in both the ocean and atmosphere. Resolving the mesoscale field is important as it has been shown to have a significant impact on the large-scale climate in eddy-rich regions, which are regions of large CO2 absorption. The presence of ocean and atmospheric mesoscale features can affect sea surface temperatures, the strength and location of storm tracks, and many other air-sea processes. Additionally, with an improvement in resolution, the eddy kinetic energy in the ocean can be expected to change considerably. The significance model resolution has on the Southern Hemisphere is examined using the Community Climate System Model, Version 4, eddy-parameterizing and eddy-resolving simulations. The CO2 concentrations and ozone levels are specified independently to better understand how the mesoscale field responds to extreme changes in external forcing and the resulting climate impacts. Overall, in the eddy-parameterizing simulations, the ozone forcing is found to be more important than the changes in CO2 concentrations. However, in the case of the eddy-resolving simulations, the CO2 concentrations are found to be more dominant, especially in eddy-rich regions. These results demonstrate the need for an increase in model resolution for climate prediction.

Daher, Houraa (ORCID:0000000200177346)↗

Annual carbon dioxide flux over seasonal sea ice in the Canadian Arctic

Continuous measurements of carbon dioxide (CO 2 ) flux were collected from a 10 m eddy covariance tower in a coastal-marine environment in the Canadian Arctic Archipelago over the course of a 17-month period. The extended length of data collection resulted in a unique dataset that includes measurements from two spring melt and summer seasons and one autumn freeze-up. These field observations were used to verify findings from previous theoretical and laboratory experiments investigating air-sea gas exchange in connection with sea ice. The results corroborated previous findings showing that thick ice cover under winter conditions acts as a barrier to gas exchange. In the spring, CO 2 fluxes were downward (uptake) in both the presence of melt ponds and during ice break-up. However, diurnal cycles were present throughout the early spring melt period, corresponding to the opposing influences of freezing and melting at the ice surface. Fluxes measured during melt periods confirmed previous laboratory tank measurements that showed a gas transfer coefficient of melting ice of 0.4 mol m −2 d −1 atm −1 . Open water CO 2 fluxes showed outgassing in early summer and uptake in mid-to-late summer, tied closely to trends in surface water temperature and its effect on the partial pressure of CO 2 in the water. The autumn period of the field campaign represents the first eddy covariance CO 2 fluxes measured over naturally forming sea ice. Our measurements showed mean upward fluxes (outgassing) of 1.1±1.5 mmol m −2 d −1 associated with the freezing of ice – the same order of magnitude found by previous laboratory tank experiments. However, peak flux periods during ice formation had measured fluxes that were a factor of 3 higher than the tank experiments, suggesting the importance of natural conditions (e.g., wind) on air-ice gas exchange. Conducting an Arctic-wide extrapolation we estimate CO 2 outgassing from the freezing period to be a counterbalance equivalent to 5 to 15 % of the magnitude of the estimated Arctic CO 2 sink. Overall, there was no evidence of dramatically enhanced gas exchange in marginal ice conditions as proposed by previous studies. Although the different seasons showed active CO 2 exchange, there was a balance between upward and downward fluxes at this specific location, resulting in a small net CO 2 uptake over the annual cycle of −0.3 g C m −2 .

54 ENVIRONMENTAL SCIENCES↗

CAPE-K_ISO: Water vapor isotopic measurements from Western Tasmania

This dataset provides processed 1-minute water vapor isotope observations from the Water Vapor Isotopic Measurements During CAPE-K project at the Kennaook/Cape Grim Baseline Air Pollution Station in northwestern Tasmania. Measurements were made with a Picarro L-2130 cavity ring-down spectroscopy analyzer during the CAPE-K deployment. The dataset includes water vapor mixing ratio, d18O, dD, deuterium excess, and propagated 1-sigma uncertainties. Isotope ratios are calibrated to the VSMOW/SLAP scale. Calibration intervals, transition periods, and obvious isotope outliers have been removed. The processed record spans 11 June 2024 to 15 October 2025 and is intended for studies of Southern Ocean marine boundary-layer moisture, isotope variability, air-sea exchange, and cloud and precipitation processes.

Water Vapor Isotopic Composition↗

Future North Atlantic tropical cyclone intensities in thermodynamically modified historical environments

Tropical cyclones (TCs) rank as the deadliest and most financially crippling natural disasters in the United States for the last half-century. It is imperative to assess potential shifts in TC intensity within the paradigm of an evolving climate. In this study, we have modeled the intensities of 620 historical TC events in the North Atlantic Basin using the Risk Analysis Framework for Tropical Cyclones (RAFT)'s deep learning intensity model. By applying a thermodynamic warming signal extrapolated from Global Climate Models, we rerun historical events under eight different future climate scenarios, providing a spectrum of potential TC intensity outcomes. One of the future simulations indicates a staggering 43% increase in the number of major hurricanes, underscoring the critical impact of climate change on TC intensity. Additionally, an interactive dashboard has been created to enable users to explore individual storm simulations and understand the influence of future climate signals on environmental conditions of TC development and resulting TC intensities. This dataset and the user-friendly tool offer invaluable resources for systematic exploration of the discrete effects that changes in the air-sea thermodynamic state have on the intensities of TCs.

Climate Change↗

Future North Atlantic tropical cyclone intensities in thermodynamically modified historical environments

Tropical cyclones (TCs) have ranked as the deadliest and most financially crippling natural disasters in the United States. It is imperative to assess potential shifts in TC intensity within the paradigm of an evolving climate. In this study we apply a fixed-constraint storyline approach that holds storm tracks and initial conditions constant to probe future TC intensity in the North Atlantic Basin. First, we simulate 618 historical TC events using the Risk Analysis Framework for Tropical Cyclones (RAFT)'s deep-learning intensity model. Next, we apply warming signals derived from eight CMIP6 climate scenarios and rerun each event to explore how intensities respond across scenarios. Finally, we develop an interactive dashboard that allows users to explore individual storm simulations and the scenario-modified environmental drivers. Together, this dataset and tool provide a clear, illustrative way to investigate how TC intensity responds to changes in air-sea state.

Climate Change↗

E3SM‐Arctic: Regionally Refined Coupled Model for Advanced Understanding of Arctic Systems Interactions

Earth system models are essential tools for climate projections, but coarse resolutions limit regional accuracy, especially in the Arctic. Regionally refined meshes (RRMs) enhance resolution in key areas while maintaining computational efficiency. This paper provides an overview of the United States (U.S.) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2.1 with an Arctic RRM, hereafter referred to as E3SMv2.1-Arctic, for the atmosphere (25 km), land (25 km), and ocean/ice (10 km) components. We evaluate the atmospheric component and its interactions with land, ocean, and cryosphere by comparing the RRM (E3SM2.1-Arctic) historical simulations (1950–2014) with the uniform low-resolution (LR) counterpart, reanalysis products, and observational data sets. The RRM generally reduces biases in the LR model, improving simulations of Arctic large-scale mean fields, such as precipitation, atmospheric circulation, clouds, atmospheric river frequency, and sea ice thickness. However, it introduces a seasonally dependent surface air temperature bias, reducing the LR cold bias in summer but enhancing the LR warm bias in winter, which contributes to the underestimated winter sea ice area and volume. Radiative feedback analysis shows similar climate feedback strengths in both model configurations, with the RRM exhibiting a more positive surface albedo feedback and contributing to a stronger surface warming than LR. These findings underscore the importance of high-resolution modeling for advancing our understanding of Arctic climate changes and their broader global impacts, although some persistent biases appear to be independent of model resolution at 10–100 km scales.

Energy Exascale Earth System Model (E3SM)↗

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗

Stronger Tropical Cyclone–Induced Ocean Cooling in Near-Coastal Regions Compared to the Open Ocean

Abstract Tropical cyclones (TC) often induce strong mixing in the upper ocean that generates a trail of cooler sea surface temperature (Twake) in their wakes. The Twake can affect TC intensity, so its prediction is important, especially in coastal regions where TCs can make landfall. Coastal Twakes are often more complex than those in the open ocean due to the influences of coastline geometry, highly variable water depth, continental runoff, and shelf processes. Using observational data since 2002, here we show a significantly stronger global mean Twake in coastal regions compared to offshore regions. Temperature stratification is the main driver of stronger coastal Twakes in the North Atlantic and east Pacific. In the northwest Pacific and north Indian Ocean, the differences between coastal and offshore Twakes are smaller due to compensation between TC forcings and ocean stratification. The north Indian Ocean is unique in the Northern Hemisphere because salinity stratification plays a major role on the spatial distribution of Twake. In the South Pacific Ocean, TC intensity and translation speed are crucial for explaining coastal–offshore Twake differences, while ocean stratification and mixed layer depth are more important for the coastal–offshore Twake differences in the south Indian Ocean. These findings suggest that coastal–offshore differences in ocean stratification need to be properly represented in models in order to capture changes in TC-induced ocean cooling as storms approach landfall. Significance Statement Landfalling tropical cyclones (TCs) often cause considerable damage in coastal regions with dense human populations. Understanding TC–ocean interaction and how it differs between coastal and offshore regions can help predict TC intensity prior to landfall. Sea surface cooling after TC passage is an important proxy for TC–ocean interaction. A global evaluation of coastal TC-induced cooling has not been conducted. Using data covering two decades, we show significantly stronger TC-induced surface cooling in coastal regions compared to offshore regions at the global scale and in all basins except the northwest Pacific and north Indian Ocean. The difference is driven mainly by upper-ocean conditions in the North Atlantic, east Pacific, and south Indian Ocean, and by TC characteristics in the South Pacific.

54 ENVIRONMENTAL SCIENCES↗

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

54 ENVIRONMENTAL SCIENCES↗

Observations of Coastal Wind Momentum Flux: Dependence on Fetch and Waves with Comparisons to COARE

Here, using observations from the Martha’s Vineyard Coastal Observatory, this paper investigates how momentum flux in the marine atmospheric surface layer over the coastal ocean varies with fetch, wave age, and wave slope and assesses the performance of the COARE 3.5 bulk parameterization. Long-fetch (at least 300 km) and short-fetch (3–6 km away from land) conditions have very similar momentum flux, with the latter being just 15% higher. The COARE 3.5 wind speed–dependent formulation closely matches the observations. The sea state dependence of wave age and wave slope is analyzed by considering both peak frequency and mean frequency in a wave spectrum. The slope of regression lines between normalized roughness and wave age is sensitive to wind speed ranges and the scatter of momentum flux, potentially explaining why earlier studies did not find a universal formula to characterize the wave dependence. Although the observed momentum flux exhibits an obvious dependence on wave age, a robust quadratic fit between momentum flux and the 10-m neutral wind speed exists only for young waves. For the short-fetch conditions, the momentum flux does not increase as wave slope increases. This may be because waves generated by local wind are still weak and swell that propagates from other areas dominates the wave spectra. In other words, the wave slope computed by considering the spectra does not properly reflect the local wind–wave interaction.

Air-sea interaction↗

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗

A Large-Eddy Simulation Study in WRF on Wind over Broadband Waves of Different Directions and Spreading Widths

We use the Weather Research and Forecasting (WRF) Model coupled with moving waves to conduct large-eddy simulations (LESs) of wind over broadband waves with different propagating directions and spreading widths. Our results show that wind-opposing waves can double the form drag, and the wave propagating direction affects the mean wind velocity, velocity variances, and pressure stress. On the other hand, waves with wider spread tend to reduce the form drag in the streamwise direction. Results further indicate that the wave direction can impact the bulk drag coefficient by as much as 25%, while the wave directional spreading width can change it by 5%. Based on our wave-phase-resolved simulations, we demonstrate that the parameterization of sea surface roughness is significantly influenced by the direction of wave propagation relative to the wind, a factor that the commonly used Charnock relation does not account for.

54 ENVIRONMENTAL SCIENCES↗