scholarly journals The Matsuno baroclinic wave test case

2018 ◽  
Author(s):  
Ofer Shamir ◽  
Itamar Yacoby ◽  
Nathan Paldor
Keyword(s):  
Shallow Water ◽  
Gravity Waves ◽  
Analytic Solutions ◽  
Wave Mode ◽  
Wave Number ◽  
Water Model ◽  
Test Case ◽  
Wave Solutions ◽  
Barotropic Vorticity ◽  
Spatial Spectra

Abstract. The analytic wave-solutions obtained by Matsuno (1966) in his seminal work on equatorial waves provide a simple and informative way of assessing atmospheric and oceanic models by measuring the accuracy with which they simulate these waves. These solutions approximate the solutions of the shallow water equations on the sphere for small speeds of gravity waves such as those of the baroclinic modes in the atmosphere and ocean. This is in contrast to the solutions of the non-divergent barotropic vorticity equation, used in the Rossby-Haurwitz test case, which are only accurate for large speeds of gravity waves such as those of the barotropic mode. The proposed test case assigns specific values to the wave-parameters (gravity wave speed, zonal wave-number, meridional wave-mode and amplitude) for both planetary and inertia gravity waves, and confirms the accuracy of the simulation by employing Hovmöller diagrams and temporal and spatial spectra. The proposed test case is successfully applied to a standard finite-difference, equatorial, non-linear, shallow water model in spherical coordinates, which demonstrates that Matsuno’s wave-solutions can be accurately simulated for at least 10 wave-periods, which for oceanic planetary waves is nearly 1300 days. In order to facilitate the use of the proposed test case, we provide Matlab, Python and Fortran codes for computing the analytic solutions at any time on arbitrary latitude-longitude grids.

scholarly journals The Matsuno baroclinic wave test case

2019 ◽  
Vol 12 (6) ◽  
pp. 2181-2193
Author(s):  
Ofer Shamir ◽  
Itamar Yacoby ◽  
Shlomi Ziskin Ziv ◽  
Nathan Paldor
Keyword(s):  
Shallow Water ◽  
Gravity Wave ◽  
Global Scale ◽  
Wave Mode ◽  
Wave Number ◽  
Scale Model ◽  
Water Model ◽  
Test Case ◽  
Wave Speeds ◽  
Wave Solutions

Abstract. The analytic wave solutions obtained by Matsuno (1966) in his seminal work on equatorial waves provide a simple and informative way of assessing the performance of atmospheric models by measuring the accuracy with which they simulate these waves. These solutions approximate the solutions of the shallow-water equations on the sphere for low gravity-wave speeds such as those of the baroclinic modes in the atmosphere. This is in contrast to the solutions of the non-divergent barotropic vorticity equation, used in the Rossby–Haurwitz test case, which are only accurate for high gravity-wave speeds such as those of the barotropic mode. The proposed test case assigns specific values to the wave parameters (gravity-wave speed, zonal wave number, meridional wave mode and wave amplitude) for both planetary and inertia-gravity waves, and suggests simple assessment criteria suitable for zonally propagating wave solutions. The test is successfully applied to a spherical shallow-water model in an equatorial channel and to a global-scale model. By adding a small perturbation to the initial fields it is demonstrated that the chosen initial waves remain stable for at least 100 wave periods. The proposed test case can also be used as a resolution convergence test.

Asymptotic stability of a polar vortex perturbed by harmonic waves describing atmospheric gravity waves circulating in an equatorial plane of a spherical planet

2018 ◽  
Vol 13 (4) ◽  
pp. 36
Author(s):  
Ranis Ibragimov ◽  
Pirooz Mohazzabi ◽  
Rebecca Roembke ◽  
Justin Van Ee
Keyword(s):  
Shallow Water ◽  
Gravity Waves ◽  
Equatorial Plane ◽  
Polar Vortex ◽  
Water Model ◽  
Large Radius ◽  
Atmospheric Gravity Waves ◽  
Nonstationary Motion ◽  
Two Parameters ◽  
Atmospheric Gravity

We examine stability of the vortex that represents one particular class of exact solution of a a nonlinear shallow water model describing atmospheric gravity waves circulating in an equatorial plane of a spherical planet. The mathematical model is represented by a two-dimensional free boundary Cauchy–Poisson problem on the nonstationary motion of a perfect uid around a solid circle with a sufficiently large radius so that the gravity is directed to the center of the circle. It is shown that the model admits two functionally independent nonlinear systems of shallow water equations. Two essential parameters that control stability of the vortex for both systems are identified. The order of their importance is analyzed and it is shown that one of the systems is more resistant to small perturbations and remains stable for larger range of these two parameters.

scholarly journals A Shallow Water Model for Convective Self-Aggregation

2020 ◽  
Author(s):  
Da Yang
Keyword(s):  
Shallow Water ◽  
Gravity Waves ◽  
Large Scale ◽  
Wave Pattern ◽  
Water Model ◽  
Shallow Water Model ◽  
Modeling Framework ◽  
Convective Storms ◽  
Wide Range ◽  
Self Aggregation

AbstractRandomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. We show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping timescale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

Existence and properties of ageostrophic modons and coherent tripoles in the two-layer rotating shallow water model on the -plane

2012 ◽  
Vol 706 ◽  
pp. 71-107 ◽  
Author(s):  
Noé Lahaye ◽  
Vladimir Zeitlin
Keyword(s):  
Shallow Water ◽  
Gravity Waves ◽  
Coherent Structures ◽  
Water Model ◽  
Shallow Water Model ◽  
Scatter Plot ◽  
Weak Stratification ◽  
Rotating Shallow Water ◽  
Frontal Collisions ◽  
Inertia Gravity Waves

AbstractWe study formation and properties of new coherent structures: ageostrophic modons in the two-layer rotating shallow water model. The ageostrophic modons are obtained by ‘ageostrophic adjustment’ of the exact modon solutions of the two-layer quasi-geostrophic equations with the free surface, which are used to initialize the full two-layer shallow water model. Numerical simulations are performed using a well-balanced high-resolution finite volume numerical scheme. For large enough Rossby numbers, the initial configurations undergo ageostrophic adjustment towards asymmetric ageostrophic quasi-stationary coherent dipoles. This process is accompanied by substantial emission of inertia–gravity waves. The resulting dipole is shown to be robust and survives frontal collisions. It contains captured inertia–gravity waves and, for higher Rossby numbers and weak stratification, carries a (baroclinic) hydraulic jump at its axis. For stronger stratifications and high enough Rossby numbers, ‘rider’ coherent structures appear as a result of adjustment, with a monopole in one layer and a dipole in another. Other ageostrophic coherent structures, such as two-layer tripoles and two-layer modons with nonlinear scatter plot, result from the collisions of ageostrophic modons. They are shown to be long-living and robust, and to capture waves.

scholarly journals Towards fast large-scale flood simulations using 2D Shallow water modelling with depth-dependant porosity

2020 ◽  
Author(s):  
Vita Ayoub ◽  
Carole Delenne ◽  
Patrick Matgen ◽  
Pascal Finaud-Guyot ◽  
Renaud Hostache
Keyword(s):  
Shallow Water ◽  
Spatial Data ◽  
Large Scale ◽  
Model Performance ◽  
Water Model ◽  
Computational Time ◽  
Small Scale ◽  
Test Case ◽  
Strong Impact ◽  
Flood Simulations

<p><span>In hydrodynamic modelling, the mesh resolution has a strong impact on run time and result accuracy. Coarser meshes allow faster simulations but often at the cost of accuracy. Conversely, finer meshes offer a better description of complex geometries but require much longer computational time, which makes their use at a large scale challenging. In this context, we aim to assess the potential of a two-dimensional shallow water model with depth-dependant porosity (SW2D-DDP) for flood simulations at a large scale. This modelling approach relies on nesting a sub-grid mesh containing high-resolution topographic and bathymetric data within each computational cell via a so-called depth-dependant storage porosity. It enables therefore faster simulations on rather coarse grids while preserving small-scale topography information. The July 2007 flood event in the Severn River basin (UK) is used as a test case, for which hydrometric measurements and spatial data are available for evaluation. A sensitivity analysis is carried out to investigate the porosity influence on the model performance in comparison with other classical parameters such as boundary conditions.</span></p>

scholarly journals A Discontinuous Galerkin Global Shallow Water Model

2005 ◽  
Vol 133 (4) ◽  
pp. 876-888 ◽  
Author(s):  
Ramachandran D. Nair ◽  
Stephen J. Thomas ◽  
Richard D. Loft
Keyword(s):  
Shallow Water ◽  
Discontinuous Galerkin ◽  
Numerical Solutions ◽  
Time Integration ◽  
Curvilinear Coordinates ◽  
Basis Set ◽  
Water Model ◽  
Test Case ◽  
Shallow Water Model ◽  
Cubed Sphere

A discontinuous Galerkin shallow water model on the cubed sphere is developed, thereby extending the transport scheme developed by Nair et al. The continuous flux form nonlinear shallow water equations in curvilinear coordinates are employed. The spatial discretization employs a modal basis set consisting of Legendre polynomials. Fluxes along the element boundaries (internal interfaces) are approximated by a Lax–Friedrichs scheme. A third-order total variation diminishing Runge–Kutta scheme is applied for time integration, without any filter or limiter. Numerical results are reported for the standard shallow water test suite. The numerical solutions are very accurate, there are no spurious oscillations in test case 5, and the model conserves mass to machine precision. Although the scheme does not formally conserve global invariants such as total energy and potential enstrophy, conservation of these quantities is better preserved than in existing finite-volume models.

Equilibrium dynamics in a forced-dissipative f-plane shallow-water system

1994 ◽  
Vol 280 ◽  
pp. 369-394 ◽  
Author(s):  
Li Yuan ◽  
Kevin Hamilton
Keyword(s):  
Energy Spectrum ◽  
Shallow Water ◽  
Gravity Waves ◽  
Water System ◽  
Water Model ◽  
Dimensional System ◽  
Linear Potential ◽  
Two Dimensional ◽  
Equilibrium Dynamics ◽  
Shallow Water System

The equilibrium dynamics in a homogeneous forced-dissipative f-plane shallow-water system is investigated through numerical simulations. In addition to classical two-dimensional turbulence, inertio-gravity waves also exist in this system. The dynamics is examined by decomposing the full flow field into a dynamically balanced potential-vortical component and a residual ‘free’ component. Here the potential-vortical component is defined as part of the flow that satisfies the gradient-wind balance equation and that contains all the linear potential vorticity of the system. The residual component is found to behave very nearly as linear inertio-gravity waves. The forcing employed is a mass and momentum source balanced so that only the large-scale potential-vortical component modes are directly excited. The dissipation is provided by a linear relaxation applied to the large scales and by an eighth-order linear hyperdiffusion. The statistical properties of the potential-vortical component in the fully developed flow were found to be very similar to those of classical two-dimensional turbulence. In particular, the energy spectrum of the potential-vortical component at scales smaller than the forcing is close to the ∼ k−3 expected for a purely two-dimensional system. Detailed analysis shows that the downscale enstrophy cascade into any wavenumber is dominated by very elongated triads involving interactions with large scales. Although not directly forced, a substantial amount of energy is found in the inertio-gravity modes and interactions among inertio-gravity modes are principally responsible for transferring energy to the small scales. The contribution of the inertio-gravity modes to the flow leads to a shallow tail at the high-wavenumber end of the total energy spectrum. For parameters roughly appropriate for the midlatitude atmosphere (notably Rossby number ∼ 0.5), the break between the roughly ∼ k−3 regime and this shallower regime occurs at scales of a few hundred km. This is similar to the observed mesoscale regime in the atmosphere. The nonlinear interactions among the inertio-gravity modes are extremely broadband in spectral space. The implications of this result for the subgrid-scale closure in the shallow-water model are discussed.

scholarly journals A novel method for the extraction of local gravity wave parameters from gridded three-dimensional data: description, validation, and application

2018 ◽  
Vol 18 (9) ◽  
pp. 6971-6983 ◽  
Author(s):  
Lena Schoon ◽  
Christoph Zülicke
Keyword(s):  
Gravity Wave ◽  
Gravity Waves ◽  
Middle Atmosphere ◽  
Grid Point ◽  
Three Dimensional ◽  
Wave Number ◽  
Test Case ◽  
Synthetic Test ◽  
Novel Method ◽  
Wave Properties

Abstract. For the local diagnosis of wave properties, we develop, validate, and apply a novel method which is based on the Hilbert transform. It is called Unified Wave Diagnostics (UWaDi). It provides the wave amplitude and three-dimensional wave number at any grid point for gridded three-dimensional data. UWaDi is validated for a synthetic test case comprising two different wave packets. In comparison with other methods, the performance of UWaDi is very good with respect to wave properties and their location. For a first practical application of UWaDi, a minor sudden stratospheric warming on 30 January 2016 is chosen. Specifying the diagnostics for hydrostatic inertia–gravity waves in analyses from the European Centre for Medium-Range Weather Forecasts, we detect the local occurrence of gravity waves throughout the middle atmosphere. The local wave characteristics are discussed in terms of vertical propagation using the diagnosed local amplitudes and wave numbers. We also note some hints on local inertia–gravity wave generation by the stratospheric jet from the detection of shallow slow waves in the vicinity of its exit region.

scholarly journals Data assimilation of two-dimensional geophysical flows with a Variational Ensemble Kalman Filter

2014 ◽  
Vol 1 (1) ◽  
pp. 403-446 ◽  
Author(s):  
Z. Mussa ◽  
I. Amour ◽  
A. Bibov ◽  
T. Kauranne
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Shallow Water ◽  
Ensemble Kalman Filter ◽  
Geophysical Flows ◽  
Water Model ◽  
Test Case ◽  
Two Dimensional ◽  
Shallow Water Model ◽  
Model Bias

Abstract. The Variational Ensemble Kalman Filter (VEnKF), a recent data assimilation method that combines a variational assimilation of the Bayesian estimation problem with an ensemble of forecasts, is demonstrated in two-dimensional geophysical flows using a Quasi-Geostrophic (QG) model and a shallow water model. Using a synthetic experiment, a two layer QG model with model bias is solved on a cylindrical 40 x 20 domain. The performance of VEnKF on the QG model with increasing ensemble size is compared with the classical Extended Kalman Filter (EKF). It is shown that although convergence can be achieved with just 20 ensemble members, increasing the number of members results in a better estimate that approaches the one produced by EKF. In the second test case, a 2-D shallow water model is described using a real dam-break experiment. The VEnKF algorithm was used to assimilate observations obtained from a modified laboratory dam-break experiment with a two-dimensional setup of sensors at the downstream end. The wave meters are placed parallel to the direction of the flow alongside the flume walls to capture both cross flow and stream flow. In both test cases, VEnKF was able to predict genuinely two-dimensional flow patterns when the sensors had a two-dimensional geometry and was stable against model bias in the first test case. In the second test case, the experiments are complemented with an empirical study of the impact of observation interpolation on the stability of the VEnKF filter. In this study, a novel Courant–Friedrichs–Lewy type filter stability condition is observed that relates ensemble variance to the time interpolation distance between observations. The results of the two experiments shows that VEnKF is a good candidate for data assimilation problems and can be implemented in higher dimensional nonlinear models.

Sign in / Sign up

Export Citation Format

Share Document