The Euler equations are solved for non-hydrostatic atmospheric flow problems in two dimensions using a high-resolution Godunov-type scheme. The Riemann problem is solved using a flux-based wave decomposition suggested by LeVeque. This paper describes in detail, the design and implementation of the Riemann solver used for computing the Godunov fluxes. The methodology is then validated against benchmark cases for non-hydrostatic atmospheric flows. Comparisons are made with solutions obtained from the National Center for Atmospheric Research's state-of-the-art numerical model. The method shows promise in simulating non-hydrostatic flows, which are characterized by steep gradients on the meso-, micro-and urban-scales.

Fully nonlinear mesoscale model simulations are used to investigate the momentum fluxes of gravity waves that emerge at a “far-field” height of 6 km from steady unsheared flow over both an axisymmetric and elliptical obstacle for nondimensional mountain heights ĥm = Fr−1 in the range 0.1–5, where Fr is the surface Froude number. Fourier- and Hilbert-transform diagnostics of model output yield local estimates of phase-averaged momentum flux, while area integrals of momentum flux quantify the amount of surface pressure drag that translates into far-field gravity waves, referred to here as the “wave drag” component. Estimates of surface and wave drag are compared to parameterization predictions and theory. Surface dynamics transition from linear to high-drag (wave breaking) states at critical inverse Froude numbers Frc−1 predicted to within 10% by transform relations. Wave drag peaks at Frc−1 < ĥm ≲ 2, where for the elliptical obstacle both surface and wave drag vacillate owing to cyclical buildup and breakdown of waves. For the axisymmetric obstacle, this occurs only at ĥm = 1.2. At ĥm ≳ 2–3 vacillation abates and normalized pressure drag assumes a common normalized form for both obstacles that varies approximately as ĥm−1.3. Wave drag in this range asymptotes to a constant absolute value that, despite its theoretical shortcomings, is predicted to within 10%–40% by an analytical relation based on linear clipped-obstacle drag for a Sheppard-based prediction of dividing streamline height. Constant wave drag at ĥm ∼ 2–5 arises despite large variations with ĥm in the three-dimensional morphology of the local wave momentum fluxes. Specific implications of these results for the parameterization of subgrid-scale orographic drag in global climate and weather models are discussed.

A time-dependent generalization of a Fourier-ray method is presented and tested for fast numerical computation of high-resolution nonhydrostatic mountain-wave fields. The method is used to model mountain waves from Jan Mayen on 25 January 2000, a period when wavelike cloud banding was observed long distances downstream of the island by the Advanced Very High Resolution Radiometer Version 3 (AVHRR-3). Surface weather patterns show intensifying surface geostrophic winds over the island at 1200 UTC caused by rapid eastward passage of a compact low pressure system. The 1200 UTC wind profiles over the island increase with height to a jet maximum of ∼60–70 m s−1, yielding Scorer parameters that indicate vertical trapping of any short wavelength mountain waves. Separate Fourier-ray solutions were computed using high-resolution Jan Mayen orography and 1200 UTC vertical profiles of winds and temperatures over the island from a radiosonde sounding and an analysis system. The radiosonde-based simulations produce a purely diverging trapped wave solution that reproduces the salient features in the AVHRR-3 imagery. Differences in simulated wave patterns governed by the radiosonde and analysis profiles are explained in terms of resonant modes and are corroborated by spatial ray-group trajectories computed for wavenumbers along the resonant mode curves. Output from a nonlinear Lipps–Hemler orographic flow model also compares well with the Fourier-ray solution horizontally. Differences in vertical cross sections are ascribed to the Fourier-ray model’s current omission of tunneling of trapped wave energy through evanescent layers.

A Fourier method is combined with a mesoscale model to simulate mountain waves. The mesoscale model describes the nonlinear low-level flow and predicts the emerging wave field above the mountain. This solution serves as the lower boundary condition for the Fourier method, which follows the waves upward to much higher altitudes and downward to the ground to examine parameterizations for the orography and the lower boundary condition. A high-drag case with a Froude number of 2 ⁄3 is presented.