Search Results
You are looking at 1 - 6 of 6 items for
- Author or Editor: Michael D. Toy x
- Refine by Access: All Content x
Abstract
Using isentropic coordinates in atmospheric models has the advantage of eliminating the cross-coordinate vertical mass flux for adiabatic flow, and virtually eliminating the associated numerical error in the vertical transport. This is a significant benefit since much of the flow in the atmosphere is approximately adiabatic. Nonadiabatic processes, such as condensational heating, result in a nonzero vertical velocity 
Abstract
Using isentropic coordinates in atmospheric models has the advantage of eliminating the cross-coordinate vertical mass flux for adiabatic flow, and virtually eliminating the associated numerical error in the vertical transport. This is a significant benefit since much of the flow in the atmosphere is approximately adiabatic. Nonadiabatic processes, such as condensational heating, result in a nonzero vertical velocity
Abstract
A three-dimensional simulation of a supercell storm is performed with a nonhydrostatic model based on a hybrid isentropic-sigma vertical coordinate. The coordinate is a terrain-following, height-based coordinate near the surface that smoothly transitions to potential temperature with height. Using isentropic coordinates provides the advantage of having zero cross-coordinate vertical mass flux for adiabatic flow, which virtually eliminates the numerical error in the vertical transport. The model uses an adaptive grid algorithm by which the coordinate surfaces may deviate from their target isentropes to maintain a sufficiently smooth mesh, while allowing the turbulence and vertical motion associated with convection to develop. The storm simulated by the hybrid-coordinate model compares well with simulations by Eulerian-coordinate models, but with the key difference being that the cross-coordinate mass flux is significantly smaller in much of the domain with the hybrid-coordinate model. A semi-implicit time-differencing scheme for numerically stabilizing vertically propagating acoustic modes in isentropic coordinates is also presented in the paper.
Abstract
A three-dimensional simulation of a supercell storm is performed with a nonhydrostatic model based on a hybrid isentropic-sigma vertical coordinate. The coordinate is a terrain-following, height-based coordinate near the surface that smoothly transitions to potential temperature with height. Using isentropic coordinates provides the advantage of having zero cross-coordinate vertical mass flux for adiabatic flow, which virtually eliminates the numerical error in the vertical transport. The model uses an adaptive grid algorithm by which the coordinate surfaces may deviate from their target isentropes to maintain a sufficiently smooth mesh, while allowing the turbulence and vertical motion associated with convection to develop. The storm simulated by the hybrid-coordinate model compares well with simulations by Eulerian-coordinate models, but with the key difference being that the cross-coordinate mass flux is significantly smaller in much of the domain with the hybrid-coordinate model. A semi-implicit time-differencing scheme for numerically stabilizing vertically propagating acoustic modes in isentropic coordinates is also presented in the paper.
Abstract
An energy and potential enstrophy conserving finite-difference scheme for the shallow-water equations is derived in generalized curvilinear coordinates. This is an extension of a scheme formulated by Arakawa and Lamb for orthogonal coordinate systems. The starting point for the present scheme is the shallow-water equations cast in generalized curvilinear coordinates, and tensor analysis is used to derive the invariant conservation properties. Preliminary tests on a flat plane with doubly periodic boundary conditions are presented. The scheme is shown to possess similar order-of-convergence error characteristics using a nonorthogonal coordinate compared to Cartesian coordinates for a nonlinear test of flow over an isolated mountain. A linear normal mode analysis shows that the discrete form of the Coriolis term provides stationary geostrophically balanced modes for the nonorthogonal coordinate and no unphysical computational modes are introduced. The scheme uses centered differences and averages, which are formally second-order accurate. An empirical test with a steady geostrophically balanced flow shows that the convergence rate of the truncation errors of the discrete operators is second order. The next step will be to adapt the scheme for use on the cubed sphere, which will involve modification at the lateral boundaries of the cube faces.
Abstract
An energy and potential enstrophy conserving finite-difference scheme for the shallow-water equations is derived in generalized curvilinear coordinates. This is an extension of a scheme formulated by Arakawa and Lamb for orthogonal coordinate systems. The starting point for the present scheme is the shallow-water equations cast in generalized curvilinear coordinates, and tensor analysis is used to derive the invariant conservation properties. Preliminary tests on a flat plane with doubly periodic boundary conditions are presented. The scheme is shown to possess similar order-of-convergence error characteristics using a nonorthogonal coordinate compared to Cartesian coordinates for a nonlinear test of flow over an isolated mountain. A linear normal mode analysis shows that the discrete form of the Coriolis term provides stationary geostrophically balanced modes for the nonorthogonal coordinate and no unphysical computational modes are introduced. The scheme uses centered differences and averages, which are formally second-order accurate. An empirical test with a steady geostrophically balanced flow shows that the convergence rate of the truncation errors of the discrete operators is second order. The next step will be to adapt the scheme for use on the cubed sphere, which will involve modification at the lateral boundaries of the cube faces.
Abstract
The isentropic system of equations has particular advantages in the numerical modeling of weather and climate. These include the elimination of the vertical velocity in adiabatic flow, which simplifies the motion to a two-dimensional problem and greatly reduces the numerical errors associated with vertical advection. The mechanism for the vertical transfer of horizontal momentum is simply the pressure drag acting on isentropic coordinate surfaces under frictionless, adiabatic conditions. In addition, vertical resolution is enhanced in regions of high static stability, which leads to better resolution of features such as the tropopause. Negative static stability and isentropic overturning frequently occur in finescale atmospheric motion. This presents a challenge to nonhydrostatic modeling with the isentropic vertical coordinate. This paper presents a new nonhydrostatic atmospheric model based on a generalized vertical coordinate. The coordinate is specified in a manner similar to that of Konor and Arakawa, but “arbitrary Eulerian–Lagrangian” (ALE) methods are used to maintain coordinate monotonicity in regions of negative static stability and to return the coordinate surfaces to their isentropic “targets” in statically stable regions. The model is mass conserving and implements a vertical differencing scheme that satisfies two additional integral constraints for the limiting case of z coordinates. The hybrid vertical coordinate model is tested with mountain-wave experiments including a downslope windstorm with breaking gravity waves. The results show that the advantages of the isentropic coordinate are realized in the model with regard to vertical tracer and momentum transport. Also, the isentropic overturning associated with the wave breaking is successfully handled by the coordinate formulation.
Abstract
The isentropic system of equations has particular advantages in the numerical modeling of weather and climate. These include the elimination of the vertical velocity in adiabatic flow, which simplifies the motion to a two-dimensional problem and greatly reduces the numerical errors associated with vertical advection. The mechanism for the vertical transfer of horizontal momentum is simply the pressure drag acting on isentropic coordinate surfaces under frictionless, adiabatic conditions. In addition, vertical resolution is enhanced in regions of high static stability, which leads to better resolution of features such as the tropopause. Negative static stability and isentropic overturning frequently occur in finescale atmospheric motion. This presents a challenge to nonhydrostatic modeling with the isentropic vertical coordinate. This paper presents a new nonhydrostatic atmospheric model based on a generalized vertical coordinate. The coordinate is specified in a manner similar to that of Konor and Arakawa, but “arbitrary Eulerian–Lagrangian” (ALE) methods are used to maintain coordinate monotonicity in regions of negative static stability and to return the coordinate surfaces to their isentropic “targets” in statically stable regions. The model is mass conserving and implements a vertical differencing scheme that satisfies two additional integral constraints for the limiting case of z coordinates. The hybrid vertical coordinate model is tested with mountain-wave experiments including a downslope windstorm with breaking gravity waves. The results show that the advantages of the isentropic coordinate are realized in the model with regard to vertical tracer and momentum transport. Also, the isentropic overturning associated with the wave breaking is successfully handled by the coordinate formulation.
Abstract
A long-lived heavy precipitation area was observed along the southwest coast of Taiwan from 13 to 18 June 2008 during the Terrain-Influenced Monsoon Rainfall Experiment (TiMREX). Rainfall amounts exceeded 500 mm along portions of the coast, and the coastal plains experienced severe flooding. The precipitation systems were influenced by blocking effects, as the southerly moist monsoon flow impinged on the island. A relatively strong gradient in the sea surface temperature (SST) off the southwest coast of Taiwan existed during the rainfall event. Mesoscale SST fronts are known to influence the planetary boundary layer (PBL) such that low-level convergence and precipitation are enhanced under certain circumstances. In this study, the authors investigate the role of the SST front in enhancing the 13–18 June 2008 precipitation event over Taiwan using the Weather Research and Forecasting (WRF) Model. In control simulations with the observed SST, there is a transition from a well-mixed to a stable PBL across the front, causing the low-level flow to decelerate, resulting in an enhancement of horizontal convergence. Such a transition in the PBL and the associated convergence is greatly reduced in smoothed SST gradient model simulations, which produce over 20% less precipitation over southwest Taiwan. Sensitivity tests show that, qualitatively, the results are independent of the existence of the island of Taiwan. These findings indicate that the SST gradient over the northern South China Sea during the early summer monsoon can have a significant impact on the intensity of rainfall over Taiwan.
Abstract
A long-lived heavy precipitation area was observed along the southwest coast of Taiwan from 13 to 18 June 2008 during the Terrain-Influenced Monsoon Rainfall Experiment (TiMREX). Rainfall amounts exceeded 500 mm along portions of the coast, and the coastal plains experienced severe flooding. The precipitation systems were influenced by blocking effects, as the southerly moist monsoon flow impinged on the island. A relatively strong gradient in the sea surface temperature (SST) off the southwest coast of Taiwan existed during the rainfall event. Mesoscale SST fronts are known to influence the planetary boundary layer (PBL) such that low-level convergence and precipitation are enhanced under certain circumstances. In this study, the authors investigate the role of the SST front in enhancing the 13–18 June 2008 precipitation event over Taiwan using the Weather Research and Forecasting (WRF) Model. In control simulations with the observed SST, there is a transition from a well-mixed to a stable PBL across the front, causing the low-level flow to decelerate, resulting in an enhancement of horizontal convergence. Such a transition in the PBL and the associated convergence is greatly reduced in smoothed SST gradient model simulations, which produce over 20% less precipitation over southwest Taiwan. Sensitivity tests show that, qualitatively, the results are independent of the existence of the island of Taiwan. These findings indicate that the SST gradient over the northern South China Sea during the early summer monsoon can have a significant impact on the intensity of rainfall over Taiwan.
Abstract
The primary goal of the Second Wind Forecast Improvement Project (WFIP2) is to advance the state-of-the-art of wind energy forecasting in complex terrain. To achieve this goal, a comprehensive 18-month field measurement campaign was conducted in the region of the Columbia River basin. The observations were used to diagnose and quantify systematic forecast errors in the operational High-Resolution Rapid Refresh (HRRR) model during weather events of particular concern to wind energy forecasting. Examples of such events are cold pools, gap flows, thermal troughs/marine pushes, mountain waves, and topographic wakes. WFIP2 model development has focused on the boundary layer and surface-layer schemes, cloud–radiation interaction, the representation of drag associated with subgrid-scale topography, and the representation of wind farms in the HRRR. Additionally, refinements to numerical methods have helped to improve some of the common forecast error modes, especially the high wind speed biases associated with early erosion of mountain–valley cold pools. This study describes the model development and testing undertaken during WFIP2 and demonstrates forecast improvements. Specifically, WFIP2 found that mean absolute errors in rotor-layer wind speed forecasts could be reduced by 5%–20% in winter by improving the turbulent mixing lengths, horizontal diffusion, and gravity wave drag. The model improvements made in WFIP2 are also shown to be applicable to regions outside of complex terrain. Ongoing and future challenges in model development will also be discussed.
Abstract
The primary goal of the Second Wind Forecast Improvement Project (WFIP2) is to advance the state-of-the-art of wind energy forecasting in complex terrain. To achieve this goal, a comprehensive 18-month field measurement campaign was conducted in the region of the Columbia River basin. The observations were used to diagnose and quantify systematic forecast errors in the operational High-Resolution Rapid Refresh (HRRR) model during weather events of particular concern to wind energy forecasting. Examples of such events are cold pools, gap flows, thermal troughs/marine pushes, mountain waves, and topographic wakes. WFIP2 model development has focused on the boundary layer and surface-layer schemes, cloud–radiation interaction, the representation of drag associated with subgrid-scale topography, and the representation of wind farms in the HRRR. Additionally, refinements to numerical methods have helped to improve some of the common forecast error modes, especially the high wind speed biases associated with early erosion of mountain–valley cold pools. This study describes the model development and testing undertaken during WFIP2 and demonstrates forecast improvements. Specifically, WFIP2 found that mean absolute errors in rotor-layer wind speed forecasts could be reduced by 5%–20% in winter by improving the turbulent mixing lengths, horizontal diffusion, and gravity wave drag. The model improvements made in WFIP2 are also shown to be applicable to regions outside of complex terrain. Ongoing and future challenges in model development will also be discussed.
