An Overall Uniformity Optimization Method of the Spherical Icosahedral Grid Based on the Optimal Transformation Theory

The improvement of overall uniformity and smoothness of spherical icosahedral grids, the basic framework of atmospheric models, is a key to reducing simulation errors. However, most of the existing grid optimization methods have optimized grid from different aspects and not improved overall uniformity and smoothness of grid at the same time, directly affecting the accuracy and stability of numerical simulation. Although a well-defined grid with more than 12 points cannot be constructed on a sphere, the area uniformity and the interval uniformity of the spherical grid can be traded off to enhance extremely the overall grid uniformity and smoothness. To solve this problem, an overall uniformity and smoothness optimization method of the spherical icosahedral grid is proposed based on the optimal transformation theory. The spherical cell decomposition method has been introduced to iteratively update the grid to minimize the spherical transportation cost, achieving an overall optimization of the spherical icosahedral grid. Experiments on the four optimized grids (the spring dynamics optimized grid, the Heikes and Randall optimized grid, the spherical centroidal Voronoi tessellations optimized grid and XU optimized grid) demonstrate that the grid area uniformity of our method has been raised by 22.60% of SPRG grid, −1.30% of HR grid, 38.30% of SCVT grid and 38.20% of XU grid, and the grid interval uniformity has been improved by 2.50% of SPRG grid, 2.80% of HR grid, 11.10% of SCVT grid and 11.00% of XU grid. Although the grid uniformity of the proposed method is similar with the HR grid, the smoothness of grid deformation has been enhanced by 79.32% of grid area and 24.07% of grid length. To some extent, the proposed method may be viewed as a novel optimization approach of the spherical icosahedral grid which can improve grid overall uniformity and smoothness of grid deformation.

[1]  Shin-ichi Iga.  An equatorially enhanced grid with smooth resolution distribution generated by a spring dynamics method , 2017, J. Comput. Phys..

[2]  T. Nakajima,et al.  Comparison and evaluation of the simulated annual aerosol characteristics over China with two global aerosol models. , 2020, The Science of the total environment.

[3]  H. Cheong,et al.  Eigensolutions of the spherical Laplacian for the cubed‐sphere and icosahedral‐hexagonal grids , 2015 .

[4]  Peter Korn,et al.  A conservative discretization of the shallow-water equations on triangular grids , 2018, J. Comput. Phys..

[5]  Teruyuki Nakajima,et al.  Aerosol model evaluation using two geostationary satellites over East Asia in May 2016 , 2019, Atmospheric Research.

[6]  Thomas Dubos,et al.  DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility , 2015 .

[7]  Wei Zeng,et al.  Shape Classification Using Wasserstein Distance for Brain Morphometry Analysis , 2015, IPMI.

[8]  Peter Korn,et al.  Formulation of an unstructured grid model for global ocean dynamics , 2017, J. Comput. Phys..

[9]  M. Gunzburger,et al.  Voronoi-based finite volume methods, optimal Voronoi meshes, and PDEs on the sphere ☆ , 2003 .

[10]  Hirofumi Tomita,et al.  An optimization of the Icosahedral grid modified by spring dynamics , 2002 .

[11]  S. Yau,et al.  Variational Principles for Minkowski Type Problems, Discrete Optimal Transport, and Discrete Monge-Ampere Equations , 2013, 1302.5472.

[12]  Saulo R. M. Barros,et al.  Analysis of grid imprinting on geodesic spherical icosahedral grids , 2013, J. Comput. Phys..

[13]  David A. Randall,et al.  Optimized Icosahedral Grids: Performance of Finite-Difference Operators and Multigrid Solver , 2013 .

[14]  G. Zängl,et al.  The ICON (ICOsahedral Non‐hydrostatic) modelling framework of DWD and MPI‐M: Description of the non‐hydrostatic dynamical core , 2015 .

[15]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[16]  Ning Wang,et al.  Geometric Properties of the Icosahedral-Hexagonal Grid on the Two-Sphere , 2011, SIAM J. Sci. Comput..

[17]  Todd D. Ringler,et al.  A Multiscale Nonhydrostatic Atmospheric Model Using Centroidal Voronoi Tesselations and C-Grid Staggering , 2012 .

[18]  Shing-Tung Yau,et al.  A Geometric View of Optimal Transportation and Generative Model , 2017, Comput. Aided Geom. Des..

[19]  Gary P. T. Choi,et al.  Open and closed anatomical surface description via hemispherical area-preserving map , 2021, Signal Process..

[20]  Ali Mahdavi-Amiri,et al.  Icosahedral Maps for a Multiresolution Representation of Earth Data , 2016, VMV.

[21]  H. Miura.  Application of the Synchronized B-Grid Staggering for Solution of the Shallow-Water Equations on the Spherical Icosahedral Grid , 2019, Monthly Weather Review.

[22]  Dan Chen,et al.  Simulations of the East Asian Winter Monsoon on Subseasonal to Seasonal Time Scales Using the Model for Prediction Across Scales , 2021, Atmosphere.

[23]  Christopher J. Subich.  Higher-order finite volume differential operators with selective upwinding on the icosahedral spherical grid , 2018, J. Comput. Phys..

[24]  Roni Avissar,et al.  The Ocean-Land-Atmosphere Model (OLAM). Part II: Formulation and Tests of the Nonhydrostatic Dynamic Core , 2008 .

[25]  Qiang Du,et al.  Centroidal Voronoi Tessellations: Applications and Algorithms , 1999, SIAM Rev..

[26]  A. E. Nikitin,et al.  High-Resolution Simulation of Polar Lows over Norwegian and Barents Seas Using the COSMO-CLM and ICON Models for the 2019–2020 Cold Season , 2021, Atmosphere.

[27]  T. Takemura,et al.  Global cloud‐system‐resolving simulation of aerosol effect on warm clouds , 2008 .

[28]  Pedro S. Peixoto,et al.  Accuracy analysis of mimetic finite volume operators on geodesic grids and a consistent alternative , 2016, J. Comput. Phys..

[29]  Almut Gassmann,et al.  A global hexagonal C‐grid non‐hydrostatic dynamical core (ICON‐IAP) designed for energetic consistency , 2013 .

[30]  Hiroaki Miura,et al.  A Comparison of Grid Quality of Optimized Spherical Hexagonal–Pentagonal Geodesic Grids , 2005 .

[31]  Jie Zhang,et al.  Isometry Invariant Shape Descriptors for Abnormality Detection on Brain Surfaces Affected by Alzheimer’s Disease , 2018, 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[32]  佐藤 正樹.  Atmospheric circulation dynamics and general circulation models , 2014 .

[33]  Hirofumi Tomita,et al.  Shallow water model on a modified icosahedral geodesic grid by using spring dynamics , 2001 .

[34]  Roni Avissar,et al.  The Ocean-Land-Atmosphere Model (OLAM). Part I: Shallow-Water Tests , 2008 .

[35]  Wei Zeng,et al.  Optimal Mass Transport for Shape Matching and Comparison , 2015, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[36]  Guoliang Xu,et al.  Discrete Laplace-Beltrami Operator on Sphere and Optimal Spherical Triangulations , 2006, Int. J. Comput. Geom. Appl..

[37]  Qiang Du,et al.  Constrained Centroidal Voronoi Tessellations for Surfaces , 2002, SIAM J. Sci. Comput..

[38]  Colin J. Cotter,et al.  Computational Modes and Grid Imprinting on Five Quasi-Uniform Spherical C Grids , 2012 .

[39]  B. Stevens,et al.  Madden‐Julian oscillation as simulated by the MPI Earth System Model: Over the last and into the next millennium , 2013 .

[40]  I. Matsuyama,et al.  Nonlinear tidal dissipation in the subsurface oceans of Enceladus and other icy satellites , 2019, Icarus.

[41]  Ning Wang,et al.  Wavelet Compression Technique for High-Resolution Global Model Data on an Icosahedral Grid , 2015 .

[42]  Min Zhang,et al.  Spherical optimal transportation , 2019, Comput. Aided Des..

[43]  T. Nakajima,et al.  Investigating the assimilation of CALIPSO global aerosol vertical observations using a four-dimensional ensemble Kalman filter , 2019, Atmospheric Chemistry and Physics.

[44]  Hirofumi Tomita,et al.  Improved smoothness and homogeneity of icosahedral grids using the spring dynamics method , 2014, J. Comput. Phys..

[45]  Todd D. Ringler,et al.  Voronoi Tessellations and Their Application to Climate and Global Modeling , 2011 .

[46]  David A. Randall,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid. Part II. A Detailed Description of the Grid and an Analysis of Numerical Accuracy , 1995 .