A severe derecho impacted the Midwestern United States on 10 August 2020, causing over $12 billion (U.S. dollars) in damage, and producing peak winds estimated at 63 m s−1, with the worst impacts in Iowa. The event was not forecast well by operational forecasters, nor even by operational and quasi-operational convection-allowing models. In the present study, nine simulations are performed using the Limited Area Model version of the Finite-Volume-Cubed-Sphere model (FV3-LAM) with three horizontal grid spacings and two physics suites. In addition, when a prototype of the Rapid Refresh Forecast System (RRFS) physics is used, sensitivity tests are performed to examine the impact of using the Grell–Freitas (GF) convective scheme. Several unusual results are obtained. With both the RRFS (not using GF) and Global Forecast System (GFS) physics suites, simulations using relatively coarse 13- and 25-km horizontal grid spacing do a much better job of showing an organized convective system in Iowa during the daylight hours of 10 August than the 3-km grid spacing runs. In addition, the RRFS run with 25-km grid spacing becomes much worse when the GF convective scheme is used. The 3-km RRFS run that does not use the GF scheme develops spurious nocturnal convection the night before the derecho, removing instability and preventing the derecho from being simulated at all. When GF is used, the spurious storms are removed and an excellent forecast is obtained with an intense bowing echo, exceptionally strong cold pool, and roughly 50 m s−1surface wind gusts.
The dynamical core that predicts the three‐dimensional vorticity rather than the momentum, which is called Vector‐Vorticity Model (VVM), is implemented on a cubed sphere. Its horizontal coordinate system is not restricted to orthogonal, while the vertical coordinate is orthogonal to the horizontal surface. Accordingly, all the governing equations of the VVM, which are originally developed with Cartesian coordinates, are rewritten in terms of general curvilinear coordinates. The local coordinates on each cube surface are constructed with the gnomonic equiangular projection. Using global channel domains, the VVM on the cubed sphere has been evaluated by (1) advecting a passive tracer with a bell‐shaped initial perturbation along an east‐west latitude circle and along a north‐south meridional circle and (2) simulating the evolution of barotropic and baroclinic instabilities. The simulated results with the cubed‐sphere grids are compared to analytic solutions or those with the regular longitude‐latitude grids. The convergence with increasing spatial resolution is also quantified using standard error norms. The comparison shows that the solutions with the cubed‐sphere grids are quite reasonable for both linear and nonlinear problems when high resolutions are used. With coarse resolution, degeneracy appears in the solutions of the nonlinear problems such as spurious wave growth; however, it is effectively reduced with increased resolution. Based on the encouraging results in this study, we intend to use this model as the cloud‐resolving component in a global Quasi‐Three‐Dimensional Multiscale Modeling Framework.more » « less
- NSF-PAR ID:
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Page Range / eLocation ID:
- p. 560-577
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
The mixing of tracers by mesoscale eddies, parameterized in many ocean general circulation models (OGCMs) as a diffusive‐advective process, contributes significantly to the distribution of tracers in the ocean. In the ocean interior, diffusive contribution occurs mostly along the direction parallel to local neutral density surfaces. However, near the surface of the ocean, small‐scale turbulence and the presence of the boundary itself break this constraint and the mesoscale transport occurs mostly along a plane parallel to the ocean surface (horizontal). Although this process is easily represented in OGCMs with geopotential vertical coordinates, the representation is more challenging in OGCMs that use a general vertical coordinate, where surfaces can be tilted with respect to the horizontal. We propose a method for representing the diffusive horizontal mesoscale fluxes within the surface boundary layer of general vertical coordinate OGCMs. The method relies on regridding/remapping techniques to represent tracers in a geopotential grid. Horizontal fluxes are calculated on this grid and then remapped back to the native grid, where fluxes are applied. The algorithm is implemented in an ocean model and tested in idealized and realistic settings. Horizontal diffusion can account for up to 10% of the total northward heat transport in the Southern Ocean and Western boundary current regions of the Northern Hemisphere. It also reduces the vertical stratification of the upper ocean, which results in an overall deepening of the surface boundary layer depth. Finally, enabling horizontal diffusion leads to meaningful reductions in the near‐surface global bias of potential temperature and salinity.
We present three new semi-Lagrangian methods based on radial basis function (RBF) interpolation for numerically simulating transport on a sphere. The methods are mesh-free and are formulated entirely in Cartesian coordinates, thus avoiding any irregular clustering of nodes at artificial boundaries on the sphere and naturally bypassing any apparent artificial singularities associated with surface-based coordinate systems. For problems involving tracer transport in a given velocity field, the semi-Lagrangian framework allows these new methods to avoid the use of any stabilization terms (such as hyperviscosity) during time-integration, thus reducing the number of parameters that have to be tuned. The three new methods are based on interpolation using 1) global RBFs, 2) local RBF stencils, and 3) RBF partition of unity. For the latter two of these methods, we find that it is crucial to include some low degree spherical harmonics in the interpolants. Standard test cases consisting of solid body rotation and deformational flow are used to compare and contrast the methods in terms of their accuracy, efficiency, conservation properties, and dissipation/dispersion errors. For global RBFs, spectral spatial convergence is observed for smooth solutions on quasi-uniform nodes, while high-order accuracy is observed for the local RBF stencil and partition of unity approaches.more » « less
Solving the shallow water equations efficiently is critical to the study of natural hazards induced by tsunami and storm surge, since it provides more response time in an early warning system and allows more runs to be done for probabilistic assessment where thousands of runs may be required. Using adaptive mesh refinement speeds up the process by greatly reducing computational demands while accelerating the code using the graphics processing unit (GPU) does so through using faster hardware. Combining both, we present an efficient CUDA implementation of GeoClaw, an open source Godunov‐type high‐resolution finite volume numerical scheme on adaptive grids for shallow water system with varying topography. The use of adaptive mesh refinement and spherical coordinates allows modeling transoceanic tsunami simulation. Numerical experiments on the 2011 Japan tsunami and a local tsunami triggered by a hypothetical
M w7.3 earthquake on the Seattle Fault illustrate the correctness and efficiency of the code, which implements a simplified dimensionally split version of the algorithms. Both numerical simulations are conducted on subregions on a sphere with adaptive grids that adequately resolve the propagating waves. The implementation is shown to be accurate and faster than the original when using Central Processing Units (CPUs) alone. The GPU implementation, when running on a single GPU, is observed to be 3.6 to 6.4 times faster than the original model running in parallel on a 16‐core CPU. Three metrics are proposed to evaluate relative performance of the model, which shows efficient usage of hardware resources.
AeolusLevel‐2B (L2B) horizontal line‐of‐sight (HLOS) retrieved Rayleigh winds, produced by the European Space Agency (ESA), utilize European Centre for Medium‐Range Weather Forecasts (ECMWF) short‐term forecasts of temperature, pressure, and horizontal winds in the Rayleigh–Brillouin and M1 correction procedures. These model fields or backgrounds can contain ECMWF model‐specific errors, which may propagate to the retrieved Rayleigh winds. This study examines the sensitivity of the retrieved Rayleigh winds to the changes in the model backgrounds, and the potential benefit of using the same system, in this case the National Oceanic and Atmospheric Administration's Finite‐Volume Cubed Sphere Global Forecast System (FV3GFS), for both the corrections and the data assimilation and forecast procedures. It is shown that the differences in the model backgrounds (FV3GFS minus ECMWF) can propagate through the Level‐2B horizontal line‐of‐sight Rayleigh wind retrieval process, mainly the M1 correction, resulting in differences in the retrieved Rayleigh winds with mean and standard deviation of magnitude as large as 0.2 m·s−1. The differences reach up to 0.4, 0.6, and 0.7 m·s−1for the 95th, 99th, and 99.5th percentiles of the sample distribution with maxima of ∼1.4 m·s−1. The numbers of the large differences for the combined lower and upper 5th, 1st, and 0.5th percentile pairs are ∼6,100, 1,220, and 610 between 2.5 and 25 km height globally per day respectively. The ESA‐disseminated Rayleigh wind product (based on the ECMWF corrections) already shows a significant positive impact on the FV3GFS global forecasts. In the observing system experiments performed, compared with the ESA Rayleigh winds, the use of the FV3GFS‐corrected Rayleigh winds lead to ∼0.5% more Rayleigh winds assimilated in the lower troposphere and show enhanced positive impact on FV3GFS forecasts at the day 1–10 range but limited to the Southern Hemisphere.