Next Article in Journal
The Characteristics of Heavy Ozone Pollution Episodes and Identification of the Primary Driving Factors Using a Generalized Additive Model (GAM) in an Industrial Megacity of Northern China
Previous Article in Journal
Utility and Triggers in Uptake of Agricultural Weather and Climate Information Services in Senegal, West Africa
Previous Article in Special Issue
Spatial and Temporal Characteristics of Rainfall Anomalies in 1961–2010 in the Yangtze River Basin, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

School of Geoscience and Surveying Engineering, China University of Mining and Technology (Beijing), D11 Xueyuan Road, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Submission received: 22 September 2021 / Revised: 5 November 2021 / Accepted: 14 November 2021 / Published: 17 November 2021

Abstract

:
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. Introduction

Spherical icosahedral grid is often used as the basic framework of atmospheric models, such as the ICosahedral Nonhydrostatic model (ICON) [1,2,3], Models for Predication Across Scales (MPAS) [4], the Non-hydrostatic ICosahedral Atmospheric Model (NICAM) [5], DYNAMICO [6], Ocean–Land–Atmosphere Model (OLAM) [7,8], etc. These models are wieldy used for polar lows simulation [9], monsoon simulation [10], ocean tides dissipation [11], etc. Recently, the Spectral Radiation Transport Model for Aerosol Species coupled with NICAM (NICAM-SPRINTARS) [12] has been developed for studying and simulating atmosphere-aerosols interactions and their effects on atmospheric pollution and climate change. Related studies include simulating the annual aerosol characteristics over China [13], capturing the horizontal distribution of aerosol optical thickness fields [14], studying the resolution dependency of the transport process of black carbon from Siberia to Japan [15], and assimilating global aerosol vertical observations [16], among others.
The accuracy of numerical simulations is affected by the quality of the spherical icosahedral grid [17,18,19]. Low-quality grids may cause large numerical approximation errors. Some noises are leaded by slight irregularities of the grid, which the long-term simulation will become unstable [20]. This is because nonuniform grids may bring grid imprinting in the numerical solutions and speed up their spread in meteorological models [21,22]. In this way, the nonuniformity of the control volume becomes the main source of numerical simulation error [23]. For example, the non-orthogonality of the grid edge and its dual edge, and the non-coincidence of the midpoints of the two edges increases the truncation error of the spherical differential operators [22,23]. Grid misalignment also leads to errors in some low-order algorithms of the spherical surface [19,24], resulting in discretization of differential operators to only meet first-order convergence [17]. The grid with good quality can reduce the truncation error of these operators on the sphere [25].
At present, some grid optimization methods have been proposed to improve grid quality and to reduce the simulation errors [26]. The spring dynamics method (i.e., the SPRing dynamics Grid optimization (SPRG) method) [27,28] is utilized to improve grid interval uniformity and to reduce geometric variations [19]. Iga [29] found that the interval distributions of grids near the icosahedron vertices, optimized by the SPRG method with a zero natural spring length (SPR0), are inversely proportional to the map factor of the Lambert conformal conic projection. An analytical transformation was proposed to stretch the grids near the icosahedral vertices to decrease grid aggregations [30]. Heikes and Randall [31,32] (HR) have minimized the distances between the midpoints of hexagonal/pentagonal grids edges and their dual (triangles) edges to optimize grid. Then, they have re-optimized the grid (named tweaked grid [33]) based on previous result to avoid the hemispheric twisting, where the grid quality and the convergence of PDE operators have been further improved. Spherical Centroidal Voronoi Tessellations (SCVT) [34,35,36] has been put forward, in which grids’ centers were iteratively moved by Llyod’s method [37] to their mass centroids to minimize the distances between these two points. Xu [25] has updated iteratively grids based on the Laplace-Beltrami solution, and the optimization grid and SCVT grid are very close to each other.
Comparations have been done with the main optimization grids, such as HR grid, SCVT grid, SPRG grid, etc. for grid quality [23,33]. Some indicators, for instance the ratios between minimum and maximum grid area, the ratios between minimum and maximum grid intervals, the grid area ranges, grid interval ranges, area relative deviations, length relative deviations, etc. are used to evaluate grid quality. The grid interval ratio of SPRG grid is bigger than HR grid, that means SPRG grid’s interval uniformity is better than that of HR grid. However, its grid area ratio and grid interval ratio are all less than the non-optimized grid (NOPT grid). Its grid quality is limited by the spring constant. When the spring constant is large, some grids near the icosahedral vertices may collapse at high resolutions. When the constant is small, grids do not collapse, but some grids near the icosahedral vertices may aggregate, such that grid intervals are shortened and the amount of calculation is increased [29]. The grid intervals by Iga [29] only be increased along the stretch paths, decreasing the regularity of the grids near the icosahedral vertices. Although area uniformity of HR gird has been enhanced to the most extent, its grid interval ratio is slightly smaller than SPRG grid and NOPT grid. Among these grids, the quality of SCVT grid is the worst, its two ratios are divergent with the increasing grid resolution, and discretization accuracy of basic operator is also lowest. There is none of them being clearly better than the others [19].
Tomita [27] has implied that simulation error appears in regions with high gradient of grid area and interval deviations, and the smoothness is important for high-accuracy and stable simulation [29]. The node interpolation can follow a second-order convergence on a uniform spherical grid, but if there is any deformation, it has only a first-order convergence [38]. The grid uniformity is also important to determine the maximum time step for numerical integration [20] and to maintain the consistency of physical parameterization in atmospheric simulations [23]. Meanwhile, improvement of the uniformity is helpful for wavelet transform to improve the compression technique for weather and climate data [39,40]. The grid uniformity has been quantified from grid area and grid interval deviations in some reviews [20,23,33]. The grid area uniformity and interval uniformity are not independent of each other. As the spherical surface is non-Euclidean, a well-defined grid with more than 12 points, meting the two attributes simultaneously, cannot be constructed on a sphere [20]. However, the area uniformity and the interval uniformity of the spherical grid can be trade off to enhance extremely the overall grid uniformity and smoothness of the spherical distribution of the grid area and interval deviation.
The present study is devoted to the investigation of a novel grid overall uniformity and smoothness optimization approach rooted in the optimal transportation theory. The spherical cell decomposition method was introduced to iteratively update the grid to minimize the spherical transportation cost, achieving an optimization grid. We discuss the details of the proposed method and its effectiveness from the grid geometry quality and numerical accuracy. Unlike existing optimization approaches, the proposed method is optimal in that it does not only improve the grid overall uniformity but also reduce the grid deformation, enhancing the smoothness of grid deformation. Both the smoothness of deformation and uniformity of grid interval have been greatly improved. The smoothness of grid area deformation has also been heightened even the grid area uniformity is comparable to those obtained by the Heikes and Randall grid.
The rest of this paper is organized as follows: the optimal transportation theory, grid uniformity optimization method and the core algorithm are introduced in the Theory and Methods section. The Results and Discussions section presents the comparative results and discussions from the grid quality and numerical accuracy. The conclusions are shown in the last section.

2. Theory and Methods

Optimal transportation theory is used for the measure-preserving mapping between two probability measure spaces, such that the target probability measure is infinitely close to the source (or true) probability measure. In all measure-preserving mappings, a mapping minimizing the transportation cost is an optimal transportation mapping [41]. Gu [41] proposed a discrete spherical optimal transportation mapping based on a purely geometric method and defined the measures as areas to achieve an area-preserving mapping from the topological sphere to the unit sphere [42,43]. The mapping is global and has been applied in biomedicine [44,45], face recognition [46], generative adversarial networks [47], and other fields.
Given a continuous space Rn, two subspaces with measures, (X, ψ) ∈ Rn and (Y, φ) ∈ Rn, having the same total measures, that is X ψ d x = Y φ d y , there is a measure-preserving mapping T: XY to make φ(B) = ψ(T−1(B)), where BY, T−1(B) ⊂ X. Then the transportation cost can be formalized as in Equation (1):
C ( T ) = X c ( x , T ( x ) ) d ψ ( x )
where C: X × YR is a cost function of T; Topt = argmin{C(T)} is the optimal transportation mapping.
Spherical optimal transportation mapping realizes the area-preserving mapping from the topological sphere to the unit sphere. There is a probability measure ψ in the continuous space defined on the unit sphere and a Dirac measure φ = {φ1, φ2, …, φnv} corresponding to the point set {pi} ⊂ S2, which cannot be covered by any hemisphere, such that ψ ( S 2 ) = i = 1 n v φ i . A set of cells, W = {wi}, decomposing the unit sphere, that is S 2 = i = 1 n v w i , can be found to make φ(wi) = ψi. Then a mapping T: wipi is the optimal transportation mapping minimizing the transportation cost. The existence and uniqueness of the complete solutions in the continuous space have been proved by Gu (2013). Surfaces can be discretized by triangular meshes. The discrete spherical optimal transportation mapping is an approximate preserved area mapping from a topological spherical triangular mesh to a spherical triangular mesh. A topological spherical triangular mesh is M {V, E, F}, with point set V = {v1, v2, …, vnv}, edge set E = {e1, e2, …, ene} and face set F = {f1, f2, …, fnf}. ∀viV, there is a set of first-order adjacent vertices vadj = { v i 1 , v i 2 , …, v i N s }, Ns is the number of first-order adjacent vertices. ∀ψiψ, it can be defined as one-third of total area of triangles formed by vi and vadj, its construction is illustrated in Figure 1 and its formulization is defined in Equation (2),
ψ i = 1 3 j = 1 N s S Δ v i v j v j + 1
where S is the spherical area of ∆vivjvj+1.
There is an optimal transportation mapping T: (M, ψ) → (N, φ), such that φi = ψi and C(T) in Equation (3) is minimal, in which N {V*, E, F} is the image of M, φ is a measure set of N.
C ( T ) = 1 2 i = 1 n v ( φ i ψ i ) 2
Based on the point set V, a cell decomposition set W = {wi} can be constructed, where wi’s area is φi. When φi equals ψi, the transportation cost is minimal equal to zero. Updating the center of cell wi as v i obtains the image of the approximate preserved area mapping.
Cell decomposition on a sphere is a spherical power diagram generation. There is a point set V = {v1, v2, …, vnv} ∈ S2 and its weight set r = {r1, r2, …, rnv} ∈ R, ci(vi, ri) is a circle on a sphere with center vi and radius ri. Spherical power distance between any point pV and ci is defined in Equation (4):
pow ( p , v i ) = cos d ( p , v i ) cos r i
where d(∙) means the spherical great circle, pow(p, vi) is a geodesic distance between p and tangent point, intersection of a line through p tangent to circle and the circle, as shown in Figure 2.
The spherical power diagram of {(pi, ri)} is a cell decomposition of sphere, that is S 2 = i = 1 n v w i , where wi = {viV | pow(p, vi) < pow(p, vj), ∀vjV−{vi}}. A spherical power diagram of a random point set is shown in Figure 3.
Where the green circle is spherical circle of each point with center vi and radius ri, the red spherical polygon is spherical power cell of each point, and the blue spherical triangle is cell’s dual triangle.
It has been proved that the optimal transportation mapping can be achieved by adjusting the weight of the spherical power diagram. When the area of each power cell is equal the predefined weight, defined by the area measure, the optimal transportation mapping can be obtained. According to the spherical power diagram, the transportation cost can be defined as follows:
C ( h ) = 1 2 i = 1 n v ( φ ( ω ( h i ) ) ψ i ) 2
where h is a function of radius, h = −ln(cos(r)), φ is the power cell area.
The power cell area is an analytic function of radius [42]. Let q is a center point of line between vi and vj, there is pow (q, vi) = pow (q, vj) (Figure 4), then:
cos d ( q , v i ) e h i = cos d ( q , v j ) e h j
Where Rl and Rk are triangle power radius of ∆vivjvk and ∆vivjvl.. dl and dk are vertical distances from the triangle center ol and ok to edge [vi, vj], respectively. γi and γj are the distances between vi and q, vj and q, respectively.
Let γi = d(q, vi), γj = d(q, vj), and γi + γj = γij, there is:
cos ( γ i j γ j ) e h i = cos γ j e h j
tan γ j = 1 sin γ i j ( e h j h i cos γ j )
According to the partial derivative of tanγj with respect to hi, there is:
d γ j d h i = cos 2 γ j tan γ j h i = cos 2 γ j 1 sin γ i j e h j h i = 1 sin γ i j e h j h i ( cos 2 γ j e h j ) 2 = 1 sin γ i j cos r i cos r j ( cos R l cos d l ) 2
According to area of infinitesimal spherical quadrilateral, φ i j l = area(wi ∩ ∆vivjvl), φ i j k = area(wi ∩ ∆vivjvk), φ is the power cell area. There is:
φ j k h i = cos 2 R l sin d l sin γ i j cos 2 d l cos r i cos r j
φ i i k h j = cos 2 R k sin d k sin γ i j cos 2 d k cos r i cos r j
φ i h j = φ j h i = cos r i cos r j sin γ i j ( cos 2 R l sin d l cos 2 d l + cos 2 R k sin d k cos 2 d k )
The convexity of the transportation cost has been proved by Gu (2013) and it can be minimized by means of Newton’s method. The gradient of cost can be defined as ∇C = (φ(w1) − ψ1, φ(w2) − ψ2, …, φ(wnv) − ψnv), and the Hessian matrix can be expressed as follows:
H = [ i 1 φ i h 1 φ 2 h 1 φ n v h 1 φ 1 h 2 i 2 φ i h 2 φ n v h 2 φ 1 h n v φ 2 h n v i n v φ i h n v ]
For the spherical icosahedral grid, the measure-preserving mapping becomes a self-mapping. The measure is replaced by a virtual measure ψvirtual, where ψi_virtual = 4πR2/nv for any spherical icosahedral grid point; R is the earth radius. There is a mapping T: (Ico-Grid, ψvirtual) → (AREA-Ico-Grid, φ), where Ico-Grid is composed of a grid point set V = {vi}, a grid edge set E, a grid set F and a grid point set of AREA-Ico-Grid, V* = { v i } is the image of V under the mapping T, its grid edge set and grid set are the same as Ico-Grid’s. When the areas of the grids by cell decomposition are equal to each other, that is φ(v*) = ψvirtual(T−1(v*)), v* ⸦ AREA-Ico-Grid, T−1(v*) ⸦ Ico-Grid, T is the approximate area-preserved mapping. The Algorithm 1 of area quasi-uniformity optimization for the spherical icosahedral grid is described as follows.
Algorithm 1: Area Uniformity Optimization for the Spherical Icosahedral Grid
Input: a spherical icosahedral grid Ico-Grid {V, E, F}, step λ, initial height vector h0, threshold δC
Output: a spherical icosahedral grid with area quasi-uniformity AREA-Ico-Grid {V*, E, F}
(1) Compute the virtual measure ψvirtual = {ψi_virtual};
(2) According to h0 and V, decompose the sphere into a cell set W = {wi}, and calculate each power cell area to get the original measure φ of Ico-Grid by, φ0φ;
(3) Calculate the transportation cost C between φ0 and ψvirtual by Equation (3). If C < δC, V* ← V and move to step 5. If not, move to the next step;
(4) Update cell decomposition of S2
 (4.1) Calculate gradient ∇C = (φ(w1) − ψ1, φ(w2) − ψ2, …, φ(wnv) − ψnv)T;
 (4.2) Compute the transportation cost C via ∇C. If C<δC, move to step 5. If not, move to step 4.3;
 (4.3) According to Equations (10)–(13), calculate a Hessian matrix H;
 (4.4) Establish the relationship between hessian matrix and gradient, that is Hδh = ∇C;
 (4.5) Update height vector, h = h + λδh;
 (4.6) Decompose the sphere into a new cell set W = {wi} based on h; compute cells’ area φW = {φwi};
(5) Compute the centers of cells in W, and obtain CW = {cwi};
(6) V* ← CW, and output result;
(7) End.

3. Results and Discussions

To verify the effectiveness of this algorithm, tests were performed to compare the grid quality and numerical accuracy among the proposed algorithm optimized grid (OURS grid), the non-optimized grid (NOPT grid), the spring dynamic method optimized grid (SPRG grid), tweaked grid optimized by Heikes and Randall grid (HR grid) Spherical Centroidal Voronoi Tessellations grid (SCVT grid) and grid optimized by Xu (XU grid). Here, the construction of the NOPT grid was through recursive division, and the parameter β in the SPRG grid was set to 1.1 according to [33]. All experiments were implemented in C++ and executed on a PC with Intel Core i5-6400 [email protected] GHz, 8 GB RAM.

3.1. The Grid Quality Evaluation

3.1.1. The Grid Area Uniformity

The relative area deviation Darea was used to measure the grid area uniformity, which can be calculated with Equation (14):
D area = A A avg A avg
where A is the area of a grid; Aavg is the average area of the spherical icosahedral grid; earth radius was set to 6, 371.007 km.
Maximum grid area Amax and minimum grid area Amin and average area Aavg of NOPT grid and OURS grid at certain resolutions are listed in Table 1.
The grid area was normalized according to Equation (15) to compare the differences at various resolutions.
A ˜ = A × n v 12 × R 2
The normalized minimum area A ˜ min and maximum area A ˜ max of different spherical icosahedral grids are listed in Table 2, and their curves are illustrated in Figure 5. Ratios between the minimum and maximum area, rA = Amin/Amax, of different spherical icosahedral grids at different resolutions are presented in Figure 6.
The maximum grid area, minimum grid area and their ratios indicate the ranges of grid area. The closer the ratio is to one, the smaller the area range is, which indicates a more uniform grid area. The normalized minimum and maximum areas of all grids except the SCVT grid and XU grid gradually have converged as the resolution increases (Figure 5 and Figure 6 and Table 1 and Table 2). The area range of SCVT grid and XU grid are the same and the biggest among all grids at level 9. The one of NOPT grid is followed and that of HR grid is the smallest. The range of OURS grid is only expanded by 1.20% compared with HR grid. Same as the range, ratios of SCVT grid and XU grid are non-convergence as resolution increases and have decreased to 0.554 and 0.555 at level 9, respectively. The ratios of SPRG grid and NOPT grid have converged slowly to 0.711 and 0.734, respectively. On the contrary, the ones of HR grid and OURS grid have increased as resolution increases and converged to 0.950 and 0.937, respectively at the same resolution.
To further emphasize the spherical distribution characteristic of grid area uniformity, taking the grid at level 9 as an example, the number of grids in different intervals of grid area relative deviation calculated by Equation (15) are counted in Table 3.
Refine these intervals to further demonstrate the area features of OURS grid (Table 4). Where, N is the number of grids in different intervals, and p is grid proportions.
The cumulative grid proportions of |Darea| are depicted in Table 5, and Figure 7 (|Darea| < 0.08%).
The Darea of all grids are in [−38.0%, 22.0%] (Table 3). The ones of HR grid and OURS grid are mainly in (−2.0%, 2.0%), which is smaller than other four grids. There are more than 90% grids with Darea in (−0.06%, 0.06%) in OURS grid. The cumulative proportions have increased logarithmically as |Darea| increases, in which the increasing rate of OURS grid is the fastest and that of NOPT grid is the slowest (Figure 7). The cumulative proportions of OURS grid with |Darea| of less than 0.044% is more than 90.00%, and the ones of other five grids are only 0.02% (NOPT grid), 0.62% (SPRG grid), 13.24% (HR grid), 0.37% (SCVT grid), and 1.25% (XU grid). Although the area ratio of HR grid is bigger than OURS grid, the proportion of grid with smaller is far less than that of OURS grid.
Meanwhile, the spherical distributions of the area relative deviations at level 9 are shown in Figure 8, where the lighter the color is, the closer the grid area deviation is to zero.
From Figure 8, the grids with maximum area deviation are the pentagons, and spherical distributions patterns of Darea are symmetry. The distribution of NOPT grid has fractal characteristics and is not continuous. Because of the larger number of grids with larger Darea, its distribution is shown by darker color. The distributions of the SPRG grid, SCVT grid and XU grid are similar. As the larger number of grids with smaller Darea (−0.06%, 0.06%) in OURS grid, its distribution is shown by the lightest color and some grids with larger Darea are mainly located along the triangular boundaries. Grids with larger Darea are also located around these boundaries, however, these regions are far larger than that of OURS grid.

3.1.2. The Grid Interval Uniformity

The geodesic distance, d, between grid points, the minimum distance, dmin and the maximum distance dmax (Table 6), the ratios between them (Table 7) and the grid length relative deviation (as Equation (16)) are calculated to describe the grid interval uniformity.
D length = L L avg L avg
The grid distance was normalized according to Equation (17) to compare the differences at various resolutions (Table 7 and Figure 9).
d ˜ = d × 2 L e v e l R
Ratios between the minimum and maximum distances, rd = dmin/dmax, of different grids at different resolutions are calculated in Table 7 and presented in Figure 10.
Same as the normalized area, d ˜ min , d ˜ max and ratios indicate the ranges of grid interval. The closer the ratio is to one, the smaller the interval range is, which indicates a more uniform grid interval. The d ˜ min , d ˜ max of the all grids except the SCVT grid and XU grid have gradually converged as the resolution increases (Figure 9 and Figure 10 and Table 6 and Table 7). The ratios of SCVT grid and XU grid are non-convergence as resolution increases, and ones of SPRG grid and HR grid have converged to 0.781 and 0.778 at level 9, respectively, the ratio of OURS grid is larger than them.
The number of grids in different intervals of grid length relative deviation have been counted (Table 8) to emphasize the spherical distribution characteristic of grid interval uniformity.
There are the similar characteristics of HR and OURS grid from Table 8, so the intervals is refined to further reveal their differences in Table 9. Where, N is the number of grids in different intervals, and p is grid proportions.
The cumulative proportions of grid with |Dlength| are depicted in Table 10 and Figure 11. Because of higher cumulative proportion in OURS grid in a smaller interval, only a part of results with |Dlength| < 1.00% are demonstrated.
The Dlength of all grids are in [−15%, 21%] (Table 8), and the ones of HR grid and OURS grid are mainly in (−1%, 1%). The proportion of grid in this range of OURS grid is about 7 times that of NOPT grid, although the interval ratio of NOPT grid is bigger than that of OURS grid. The cumulative proportions also increase logarithmically as |Dlength| increases, in which the increasing rate in OURS grid is the fastest and that of NOPT grid is the slowest. The cumulative proportions of OURS grid is more than 99.00% with |Dlength| less than 0.50%. Ones of other grids are only 3.86% (NOPT grid) 17.36% (SPRG grid), 75.17% (HR grid), 21.64% (SCVT grid), and 20.47% (XU grid) in the same range.
Meanwhile, the spherical distributions of Dlength are shown in Figure 12 at level 9.
Same as spherical distribution of grid area, the grids with maximum Dlength are these pentagons, and distributions of all grids are symmetry. Near these pentagons, the lengths of grids are less than the average length and shown by darker color in all grids. The distribution of NOPT grid has fractal characteristics and is not continuous, too. The Dlength in the remaining grids are reduced from the pentagons to triangular centers. The grids with larger Dlength are mainly located along the triangular boundaries in OURS grid, and those grids are spread like pentagon or star from these pentagons.

3.2. The Numerical Accuracy Evaluation

Spherical Laplacian operators of scalar field (as Equation (18)) is discretized to evaluate numerical accuracy on different grids.
f ( λ , φ ) = cos λ × cos 4 φ
where λ and φ are the longitude and latitude of point on the sphere, respectively.
The L2-norm error (as Equation (19)), L-norm error (as Equation (20)) of discretization have been calculated (Table 11 and Table 12 and Figure 13 and Figure 14).
L 2 = ( 1 i = 1 n v A i ( i = 1 n v A i ( f i num f i ana ) 2 ) )
L = max i = 1 , n v | f i num f i ana |
where fana and fnum are the analytical and numerical solutions of Laplacian operator, respectively.
The L2-norm errors and L-norm errors of all grids have been reduced as resolution increase, ones of OURS grid and HR grid are smaller than the others. The L-norm error (meaning the maximal error) of OURS grid has been reduced from 9.15 × 10−2 (NOPT grid), 2.57 × 10−2 (SPRG grid), 1.66 × 10−2 (HR grid), 3.87 × 10−2 (SCVT and XU grid) to 1.59 × 10−2 at level 8. The L2-norm error (meaning the RMS error) of OURS grid has been reduced from 4.52 × 10−4 (NOPT grid), 2.03 × 10−4 (SPRG grid), 6.05 × 10−5 (HR grid), 2.33 × 10−4 (SCVT and XU grid) to 5.86 × 10−4 at the same grid resolution, in which the enhancement of average accuracy is nearly 8 times that of NOPT grid and has more than 11.62% compared to the HR grid.
Meanwhile, the spherical distributions of the Laplacian operator error on different grids are shown in Figure 15 at level 5.
Grids in NOPT grid with larger errors are mainly located the boundaries of refined icosahedral triangular cells, and that of all optimization grids are mainly located on the icosahedral triangular boundaries. In OURS grid, there are some fluctuations near icosahedral vertices and the middle part of boundaries, grids’ error in these regions are smaller than that of HR grid. In addition, the error range of OURS grid is the smallest among all grids, and it has been narrowed to 84.81% of NOPT grid and 12.00% of HR grid.

3.3. Discussion

The raw grid (non-optimized grid) has been optimized by the proposed method. The area quasi-uniformity can be achieved by minimizing the grid area deviation cost, which improve the smoothness of grid area deviation to some extent. Although the grid area range of OURS grid is comparable to those of Heikes and Randall grid, the deformations of grid area and intervals of OURS grid are smoother, which is conducive to error control in simulation. In addition, the maximum error and RMS error of discritization of Laplacian operator have been decreased and converged as the resolution increases.
The grid quality is one of the aspects affecting the simulation accuracy. There is no high gradient of grid area and interval deformation in the optimized grid by our method, which can be helpful to improve the accuracy of discritization of Laplacian operator. A more extensive analysis of Laplace operator in a diffusion problem and some numerical experiments in terms of the accuracy and the numerical efficiency will be carried out in the future.

4. Conclusions

In this study, an overall uniformity and smoothness optimization method of the spherical icosahedral grid has been proposed based on the optimal transportation theory. The effectiveness of the proposed method was evaluated for grid uniformity and smoothness and the following conclusions may be drawn: (1) the area uniformity measured by the ratio between minimum and maximum grid area has been improved by 22.6% (SPRG grid), 38.3% (SCVT grid) and 38.2% (XU grid), and can be comparable to the HR grid. The interval uniformity has also been increased by 2.5% (SPRG grid), 2.8% (HR grid), 11.1% (SCVT grid) and 11.0% (XU grid). (2) the smoothness of grid area deformation measured by the number of grids with grid area deviation of less than 0.05% has been enhanced by 79.32% (HR grid) and more than 90% compared to the SPRG grid, SCVT grid and XU grid. The smoothness of grid length deformation has also been increased by 24.07% (HR grid) and more than 75% compared to other three grids. Although the uniformity of grid optimized by the proposed method is comparable to the HR grid, the smoothness of grid deformation is improved. In addition, the average accuracy of Laplacian operator is slightly better than the HR grid.
In this study, we only evaluated the geometric uniformity of the optimized grid, and the truncation error of Laplacian operator. Future work will focus on the accuracy of Laplacian operator of different grids in a diffusion equation problem and some numerical experiments in terms of the accuracy and the numerical efficiency.

Author Contributions

F.L. formed the research idea, designed and conducted the experiment, drafted the manuscript; X.Z. co-designed the research, analyzed the results and revised the manuscript; W.S. revised the manuscript; Y.L. and Y.D. contributed to the grammar modification. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2018YFB0505301) and the National Natural Science Foundation of China (No. 41671394, No. 41671383).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data set available on request to corresponding authors.

Acknowledgments

The authors would like to sincerely thank the two anonymous reviewers whose insightful comments have helped to substantially improve the manuscript. We would like to acknowledge and thank Peixoto for providing the Heikes and Randall grid and SCVT grid.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gassmann, A. A global hexagonal C-grid non-hydrostatic dynamical core (ICON-IAP) designed for energetic consistency. Q. J. R. Meteorol. Soc. 2013, 139, 152–175. [Google Scholar] [CrossRef]
  2. Schubert, J.J.; Stevens, B.; Crueger, T. Madden-Julian oscillation as simulated by the MPI Earth System Model: Over the last and into the next millennium. J. Adv. Model. Earth Syst. 2013, 5, 71–84. [Google Scholar] [CrossRef]
  3. Zängl, G.; Reinert, D.; Rípodas, P.; Baldauf, M. The ICON (ICOsahedral Non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core. Q. J. R. Meteorol. Soc. 2015, 141, 563–579. [Google Scholar] [CrossRef]
  4. Skamarock, W.C.; Klemp, J.B.; Duda, M.G.; Fowler, L.D.; Park, S.-H.; Ringler, T.D. A Multiscale Nonhydrostatic Atmospheric Model Using Centroidal Voronoi Tesselations and C-Grid Staggering. Mon. Weather. Rev. 2012, 140, 3090–3105. [Google Scholar] [CrossRef] [Green Version]
  5. Satoh, M. Atmospheric Circulation Dynamics and General Circulation Models; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  6. Dubos, T.; Dubey, S.; Tort, M.; Mittal, R.; Meurdesoif, Y.; Hourdin, F. DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility. Geosci. Model Dev. 2015, 8, 3131–3150. [Google Scholar] [CrossRef] [Green Version]
  7. Walko, R.L.; Avissar, R. The Ocean-Land-Atmosphere Model (OLAM). Part I: Shallow-Water Tests. Mon. Weather Rev. 2008, 136, 4033–4044. [Google Scholar] [CrossRef]
  8. Walko, R.L.; Avissar, R. The Ocean-Land-Atmosphere Model (OLAM). Part II: Formulation and Tests of the Nonhydrostatic Dynamic Core. Mon. Weather Rev. 2008, 136, 4045–4062. [Google Scholar] [CrossRef]
  9. Revokatova, A.; Nikitin, M.; Rivin, G.; Rozinkina, I.; Nikitin, A.; Tatarinovich, E. High-Resolution Simulation of Polar Lows over Norwegian and Barents Seas Using the COSMO-CLM and ICON Models for the 2019–2020 Cold Season. Atmosphere 2021, 12, 137. [Google Scholar] [CrossRef]
  10. Hsu, L.-H.; Chen, D.-R.; Chiang, C.-C.; Chu, J.-L.; Yu, Y.-C.; Wu, C.-C. Simulations of the East Asian Winter Monsoon on Subseasonal to Seasonal Time Scales Using the Model for Prediction Across Scales. Atmosphere 2021, 12, 865. [Google Scholar] [CrossRef]
  11. Hay, H.C.F.C.; Matsuyama, I. Nonlinear tidal dissipation in the subsurface oceans of Enceladus and other icy satellites. Icarus 2018, 319, 68–85. [Google Scholar] [CrossRef]
  12. Suzuki, K.; Nakajima, T.; Satoh, M.; Tomita, H.; Takemura, T.; Nakajima, T.; Stephens, G.L. Global cloud-system-resolving simulation of aerosol effect on warm clouds. Geophys. Res. Lett. 2008, 35, 610–616. [Google Scholar] [CrossRef]
  13. Cheng, Y.; Dai, T.; Zhang, H.; Xin, J.; Chen, S.; Shi, G.; Nakajima, T. Comparison and evaluation of the simulated annual aerosol characteristics over China with two global aerosol models. Sci. Total Environ. 2021, 763, 143003–143017. [Google Scholar] [CrossRef] [PubMed]
  14. Goto, D.; Kikuchi, M.; Suzuki, K.; Hayasaki, M.; Yoshida, M.; Nagao, T.M.; Choi, M.; Kim, J.; Sugimoto, N.; Shimizu, A.; et al. Aerosol model evaluation using two geostationary satellites over East Asia in May 2016. Atmos. Res. 2019, 217, 93–113. [Google Scholar] [CrossRef]
  15. Yamashita, Y.; Takigawa, M.; Goto, D.; Yashiro, H.; Satoh, M.; Kanaya, Y.; Taketani, F.; Miyakawa, T. Effect of Model Resolution on Black Carbon Transport from Siberia to the Arctic Associated with the Well-Developed Low-Pressure Systems in September. J. Meteorol. Soc. Jpn. 2021, 99, 287–308. [Google Scholar] [CrossRef]
  16. Cheng, Y.; Dai, T.; Goto, D.; Schutgens, N.A.J.; Shi, G.; Nakajima, T. Investigating the assimilation of CALIPSO global aerosol vertical observations using a four-dimensional ensemble Kalman filter. Atmos. Chem. Phys. 2019, 19, 13445–13467. [Google Scholar] [CrossRef] [Green Version]
  17. Korn, P. Formulation of an unstructured grid model for global ocean dynamics. J. Comput. Phys. 2017, 339, 525–552. [Google Scholar] [CrossRef]
  18. Peixoto, P.S. Accuracy analysis of mimetic finite volume operators on geodesic grids and a consistent alternative. J. Comput. Phys. 2016, 310, 127–160. [Google Scholar] [CrossRef]
  19. Peixoto, P.S.; Barros, S.R.M. Analysis of grid imprinting on geodesic spherical icosahedral grids. J. Comput. Phys. 2013, 237, 61–78. [Google Scholar] [CrossRef]
  20. Wang, N.; Lee, J.L. Geometric Properties of the Icosahedral-Hexagonal Grid on the Two-Sphere. Soc. Ind. Appl. Math. 2011, 33, 2536–2559. [Google Scholar] [CrossRef] [Green Version]
  21. Weller, H.; Thuburn, J.; Cotter, C.J. Computational modes and grid imprinting on five quasi-uniform spherical c-grids. Mon. Weather Rev. 2012, 140, 2734–2755. [Google Scholar] [CrossRef] [Green Version]
  22. Cheong, H.; Kang, H. Eigensolutions of the spherical Laplacian for the cubed-sphere and icosahedral-hexagonal grids. Q. J. R. Meteorol. Soc. 2015, 141, 3383–3398. [Google Scholar] [CrossRef]
  23. Miura, H.; Kimoto, M. A comparison of grid quality of optimized spherical hexagonal—pentagonal geodesic grids. Mon. Weather Rev. 2005, 133, 2817–2833. [Google Scholar] [CrossRef]
  24. Subich, C.J. Higher-order finite volume differential operators with selective upwinding on the icosahedral spherical grid. J. Comput. Phys. 2018, 368, 21–46. [Google Scholar] [CrossRef]
  25. Xu, G. Discrete Laplace-Beltrami Operator on Sphere and Optimal Spherical Triangulations. Int. J. Comput. Geom. Appl. 2006, 16, 75–93. [Google Scholar] [CrossRef]
  26. Korn, P.; Linardakis, L. A conservative discretization of the shallow-water equations on triangular grids. J. Comput. Phys. 2018, 375, 871–900. [Google Scholar] [CrossRef]
  27. Tomita, H.; Tsugawa, M.; Satoh, M.; Goto, K. Shallow Water Model on a Modified Icosahedral Geodesic Grid by Using Spring Dynamics. J. Comput. Phys. 2001, 174, 579–613. [Google Scholar] [CrossRef]
  28. Tomita, H.; Satoh, M.; Goto, K. An Optimization of the Icosahedral Grid Modified by Spring Dynamics. J. Comput. Phys. 2002, 183, 307–331. [Google Scholar] [CrossRef]
  29. Iga, S.; Tomita, H. Improved smoothness and homogeneity of icosahedral grids using the spring dynamics method. J. Comput. Phys. 2014, 258, 208–226. [Google Scholar] [CrossRef]
  30. Iga, S. An equatorially enhanced grid with smooth resolution distribution generated by a spring dynamics method. J. Comput. Phys. 2017, 330, 794–809. [Google Scholar] [CrossRef]
  31. Heikes, R.; Randall, D.A. Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid. Part I: Basic Design and Results of Tests. Mon. Weather Rev. 1995, 123, 1862–1880. [Google Scholar] [CrossRef]
  32. Heikes, R.; Randall, D.A. 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. Mon. Weather Rev. 1995, 123, 1881–1887. [Google Scholar] [CrossRef]
  33. Heikes, R.P.; Randall, D.A.; Konor, C.S. Optimized icosahedral grids: Performance of finite-difference operators and multigrid solver. Mon. Weather Rev. 2013, 141, 4450–4469. [Google Scholar] [CrossRef]
  34. Du, Q.; Gunzburger, M.D.; Ju, L. Constrained centroidal Voronoi tessellations for surfaces. Siam J. Sci. Comput. 2003, 24, 1488–1506. [Google Scholar] [CrossRef]
  35. Du, Q.; Gunzburger, M.D.; Ju, L. Voronoi-based finite volume methods, optimal Voronoi meshes, and PDEs on the sphere. Comput. Methods Appl. Mech. Eng. 2003, 192, 3933–3957. [Google Scholar] [CrossRef]
  36. Ju, L.; Ringler, T.; Gunzburger, M. Voronoi Tessellations and Their Application to Climate and Global Modeling; Springer: Berlin/Heidelberg, Germany, 2011; Volume 80, pp. 313–342. [Google Scholar]
  37. Du, Q.; Vance, F.; Max, G. Centroidal Voronoi Tessellations: Applications and Algorithms. SIAM Rev. 1999, 41, 637–676. [Google Scholar] [CrossRef] [Green Version]
  38. Miura, H. Application of the Synchronized B Grid Staggering for Solution of the Shallow-Water Equations on the Spherical Icosahedral Grid. Mon. Weather Rev. 2019, 147, 2485–2509. [Google Scholar] [CrossRef]
  39. Wang, N.; Bao, J.W.; Lee, J.L.; Moeng, F.; Matsumoto, C. Wavelet Compression Technique for High-Resolution Global Model Data on an Icosahedral Grid. J. Atmos. Ocean. Technol. 2015, 32, 1650–1667. [Google Scholar] [CrossRef]
  40. Jubair, M.; Alim, U.; Röber, N.; Clyne, J.; Mahdavi-Amiri, A. Icosahedral Maps for a Multiresolution Representation of Earth Data. In Proceedings of the VMV’16 Proceedings of the Conference on Vision, Modeling and Visualization, Bayreuth, Germany, 10–12 October 2016. [Google Scholar]
  41. Gu, X.; Luo, F.; Sun, J. Variational principles for Minkowski type problems, discrete optimal transport, and discrete Monge-Ampere equations. Math. Methods Solid State Superfluid Theory 2013, 20, 383–398. [Google Scholar] [CrossRef] [Green Version]
  42. Cui, L.; Qi, X.; Wen, C.; Lei, N.; Li, X.; Zhang, M.; Gu, X. Spherical optimal transportation. Comput. -Aided Des. 2019, 115, 181–193. [Google Scholar] [CrossRef]
  43. Tu, Y.; Wen, C.; Wen, Z.; Wu, J.F. Isometry Invariant Shape Descriptors for Abnormality Detection on Brain Surfaces Affected by Alzheimer’s Disease. In Proceedings of the 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Honolulu, HI, USA, 17–21 July 2018. [Google Scholar]
  44. Giri, A.; Choi, G.; Kumar, L. Open and closed anatomical surface description via hemispherical area-preserving map. Signal Process. 2020, 180, 107867–107880. [Google Scholar] [CrossRef]
  45. Su, Z.; Zeng, W.; Wang, Y.; Lu, L.Z.; Gu, X.F. Shape Classification Using Wasserstein Distance for Brain Morphometry Analysis. In Information Processing in Medical Imaging (IPMI); Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  46. Su, Z.; Wang, Y.; Shi, R. Optimal mass transport for shape matching and comparison. IEEE Trans. Pattern Anal. Mach. Intell. 2015, 37, 2246–2259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Lei, N.; Su, K.; Cui, L.; Yau, S.-T.; Gu, X.D. A Geometric View of Optimal Transportation and Generative Model. Comput. Aided Geom. Des. 2019, 68, 1–21. [Google Scholar] [CrossRef]
Figure 1. Area measure ψi of vi.
Figure 1. Area measure ψi of vi.
Atmosphere 12 01516 g001
Figure 2. Spherical power distance.
Figure 2. Spherical power distance.
Atmosphere 12 01516 g002
Figure 3. Spherical power diagram.
Figure 3. Spherical power diagram.
Atmosphere 12 01516 g003
Figure 4. Power cell and its dual triangle.
Figure 4. Power cell and its dual triangle.
Atmosphere 12 01516 g004
Figure 5. A ˜ min and A ˜ max of different grids.
Figure 5. A ˜ min and A ˜ max of different grids.
Atmosphere 12 01516 g005
Figure 6. Ratios between Amin and Amax of different grids.
Figure 6. Ratios between Amin and Amax of different grids.
Atmosphere 12 01516 g006
Figure 7. Curves of cumulative proportion for different grids at level 9 (|Darea| < 0.08%).
Figure 7. Curves of cumulative proportion for different grids at level 9 (|Darea| < 0.08%).
Atmosphere 12 01516 g007
Figure 8. Spherical distributions of Darea of different grids at level 9. (a) NOPT grid. (b) SPRG grid. (c) HR grid. (d) SCVT grid. (e) XU grid. (f) OURS grid.
Figure 8. Spherical distributions of Darea of different grids at level 9. (a) NOPT grid. (b) SPRG grid. (c) HR grid. (d) SCVT grid. (e) XU grid. (f) OURS grid.
Atmosphere 12 01516 g008aAtmosphere 12 01516 g008b
Figure 9. The normalized minimum and maximum distance of different grids.
Figure 9. The normalized minimum and maximum distance of different grids.
Atmosphere 12 01516 g009
Figure 10. Ratios between minimum and maximum distance of different grids.
Figure 10. Ratios between minimum and maximum distance of different grids.
Atmosphere 12 01516 g010
Figure 11. Curves of cumulative proportion with |Dlength| < 1.00% for different grids at level 9.
Figure 11. Curves of cumulative proportion with |Dlength| < 1.00% for different grids at level 9.
Atmosphere 12 01516 g011
Figure 12. Spherical distributions of Dlength at level 9. (a) NOPT grid. (b) SPRG grid. (c) HR grid. (d) SCVT grid. (e) XU grid. (f) OURS grid.
Figure 12. Spherical distributions of Dlength at level 9. (a) NOPT grid. (b) SPRG grid. (c) HR grid. (d) SCVT grid. (e) XU grid. (f) OURS grid.
Atmosphere 12 01516 g012
Figure 13. L2-norm error of Laplacian operator. (a) L2-norm error of all grids (b) L2-norm error of OURS and HR grid from level 5 to 8.
Figure 13. L2-norm error of Laplacian operator. (a) L2-norm error of all grids (b) L2-norm error of OURS and HR grid from level 5 to 8.
Atmosphere 12 01516 g013
Figure 14. L-norm error of Laplacian operator. (a) L-norm error of all grids (b) L-norm error of OURS and HR grid from level 5 to 8.
Figure 14. L-norm error of Laplacian operator. (a) L-norm error of all grids (b) L-norm error of OURS and HR grid from level 5 to 8.
Atmosphere 12 01516 g014
Figure 15. Spherical distributions of Laplacian operator error of different grids. (a) NOPT grid. (b) SPRG grid. (c) HR grid. (d) SCVT grid. (e) XU grid. (f) OURS grid.
Figure 15. Spherical distributions of Laplacian operator error of different grids. (a) NOPT grid. (b) SPRG grid. (c) HR grid. (d) SCVT grid. (e) XU grid. (f) OURS grid.
Atmosphere 12 01516 g015
Table 1. Amin, Amax of spherical icosahedral grids (unit: km2).
Table 1. Amin, Amax of spherical icosahedral grids (unit: km2).
LevelnvAavgNOPTOURS
AminAmaxAminAmax
01242, 505, 466.0642, 505, 466.0642, 505, 466.0642, 505, 466.0642, 505, 466.06
14212, 144, 418.8811, 115, 261.9112, 556, 071.1011, 115, 254.4812, 556, 081.85
21623, 148, 553.042, 812, 532.183, 339, 347.672, 885, 170.413, 216, 806.64
3642794, 494.69705, 296.93923, 852.78735, 897.42799, 720.75
42, 562199, 088.83176, 460.40237, 913.42185, 750.43199, 596.50
510, 24249, 801.3744, 123.6359, 942.4346, 638.4349, 880.62
640, 96212, 452.1711, 031.4415, 015.2811, 700.0012, 500.00
7163, 8423, 113.162, 757.893, 755.662, 919.713, 114.73
8655, 362778.30689.48939.35729.62778.67
92, 621, 442194.57172.37234.84182.34194.60
Table 2. The normalized area and the ratios between the Amin and Amax.
Table 2. The normalized area and the ratios between the Amin and Amax.
Level3456789
A ˜ avg 1.0471.0471.0471.0471.0471.0471.047
NOPT A ˜ min 0.9300.9280.9300.9280.9280.9280.928
A ˜ max 1.2181.2511.2601.2631.2631.2641.264
rA0.7640.7420.7400.7350.7350.7340.734
SPRG A ˜ min 0.8790.8480.8200.8030.7910.7810.776
A ˜ max 1.0801.0791.0801.0811.0831.0871.092
rA0.8140.7860.7600.7430.7300.7180.711
HR A ˜ min 0.9981.0081.0121.0141.0131.0161.016
A ˜ max 1.0621.0661.0681.0681.0691.0701.070
rA0.9400.9460.9500.9490.9480.9500.950
SCVT A ˜ min 0.8610.8120.7600.7190.6770.6390.599
A ˜ max 1.0821.0811.0801.0811.0811.0811.081
rA0.7960.7510.7100.6650.6260.5910.554
XU A ˜ min 0.8610.8120.7600.7190.6760.6380.600
A ˜ max 1.0821.0811.0801.0811.0821.0821.082
rA0.7960.7510.7100.6650.6250.5900.555
OURS A ˜ min 0.9700.9770.9810.9840.9820.9820.981
A ˜ max 1.0541.0501.0491.0511.0481.0481.047
rA0.9200.9300.9400.9360.9370.9370.937
Table 3. The number of grids in different Darea interval of different grids.
Table 3. The number of grids in different Darea interval of different grids.
Interval (%)NOPTSPRGHRSCVTXUOURS
[−100, −38.0]000000
(−38.0, −34.0)0001621620
(−34.0, −30.0)000000
(−30.0, −26.0)0009609600
(−26.0, −22.0)016201, 9201, 9200
(−22.0, −18.0)096002, 8803, 8400
(−18.0, −14.0)04, 800014, 40013, 4400
(−14.0, −10.0)19234, 560042, 24042, 2400
(−10.0, −6.0)436, 770180, 4800112, 320115, 20012
(−6.0, −2.0)1, 132, 800630, 72022, 786314, 880317, 7600
(−2.0, 2.0)366, 240887, 5202598, 4801, 403, 0401, 341, 5842, 621, 430
(2.0, 6.0)60, 480647, 040176728, 640784, 3360
(6.0, 10.0)0235, 2000000
(10.0, 14.0)380, 16000000
(14.0, 18.0)96, 00000000
(18.0, 22.0)148, 80000000
(22.0, 100.0)000000
total2, 621, 4422, 621, 4422, 621, 4422, 621, 4422, 621, 4422, 621, 442
Table 4. The number of grids in refined Darea interval of OURS grid.
Table 4. The number of grids in refined Darea interval of OURS grid.
Interval (%)Np (%)
[−100, −0.18)720.00
[−0.18, −0.14)600.00
[−0.14, −0.10)17, 6950.68
[−0.10, −0.06)44, 7331.71
[−0.06, −0.02)495, 83418.91
[−0.02, 0.02)1, 342, 12951.20
(−0.02, 0.60)720, 91927.50
(0.06, 100]00
total2, 621, 442100.00
Table 5. The cumulative proportion of |Darea| of different grids (%).
Table 5. The cumulative proportion of |Darea| of different grids (%).
Interval (%)NOPTSPRGHRSCVTXUOURS
[0, 0.004)001.0700.2612.57
[0.004, 0.012)0.020.073.270.070.5933.35
[0.012, 0.02)0.020.225.870.070.7351.20
[0.02, 0.028)0.020.228.630.220.9366.98
[0.028, 0.036)0.020.5111.060.221.2581.04
[0.036, 0.044)0.020.6213.240.371.2592.56
[0.044, 0.052)0.020.7715.100.661.3996.17
[0.052, 0.06)0.020.9916.730.921.3997.61
[0.06, 0.068)0.021.6522.952.381.7699.32
[0.068, 0.076)0.022.2328.113.042.1699.99
[0.076, 0.3)0.024.6943.546.235.3599.99
[0.3, 0.5)0.027.8756.6410.259.8799.99
[0.5, 0.7)0.0211.1766.6614.3914.2599.99
[0.7, 0.9)0.0214.3974.6418.8218.5399.99
[0.9, 2)13.9733.8699.1253.5251.1899.99
[6, 26)59.4982.610093.3393.32100
[26, 100]100100100100100100
Table 6. dmin, dmax, davg of different grids (unit: km).
Table 6. dmin, dmax, davg of different grids (unit: km).
LevelNOPTOURS
davgdmindmaxdavgdmindmax
07529.857529.857529.857529.857529.857529.85
13764.923526.834003.023764.923526.834003.02
21914.331763.412079.281915.301739.232087.62
3961.22881.711050.16961.85866.341059.12
4481.12440.85526.42481.47432.15532.56
5240.62220.43263.38240.81215.84267.18
6120.32110.21131.71120.42107.83134.40
760.1655.1165.8660.2054.1367.12
830.0827.5532.9330.1027.0633.56
915.0413.7816.4715.0513.5316.78
Table 7. Normalized distance and ratios between the dmin and dmax of different grids.
Table 7. Normalized distance and ratios between the dmin and dmax of different grids.
Level3456789
d ˜ avg 1.1821.1821.1821.1821.1821.1821.182
NOPT d ˜ min 1.1071.1071.1071.1071.1071.1071.107
d ˜ max 1.3191.3221.3231.3231.3231.3231.323
rd0.8390.8370.8370.8370.8370.8370.837
SPRG d ˜ min 1.0771.0591.0431.0171.0170.9990.983
d ˜ max 1.2761.2791.2791.2841.2851.2791.259
rd0.8440.8280.8150.8020.7910.7810.781
HR d ˜ min 1.0761.0741.0741.0741.0741.0741.075
d ˜ max 1.3511.3611.3641.3641.3651.3641.382
rd0.7960.7890.7870.7870.7870.7870.778
SCVT d ˜ min 1.0651.0351.0050.9750.9460.9170.890
d ˜ max 1.2751.2771.2771.2771.2771.2771.280
rd0.8350.8100.7870.7630.7410.7180.695
XU d ˜ min 1.0651.0351.0050.9750.9450.9160.888
d ˜ max 1.2751.2771.2771.2771.2771.2771.276
rd0.8350.8100.7870.7630.7400.7170.696
OURS d ˜ min 1.0881.0851.0841.0831.0871.0911.091
d ˜ max 1.3301.3371.3421.3501.3491.3541.353
rd0.8180.8120.8080.8020.8060.8060.806
Table 8. The number of grids with different Dlength interval of different grids.
Table 8. The number of grids with different Dlength interval of different grids.
NOPTSPRGHRSCVTXUOURS
[−100, −15]000000
(−15, −13]00060600
(−13, −11]0001201200
(−11, −9]0003604200
(−9, −7]05, 73007201, 0200
(−7, −5]028, 80002, 6402, 2800
(−5, −3]0151, 6800100, 800106, 0300
(−3, −1]1, 488, 960620, 1600301, 632301, 6320
(−1, 1)446, 880966, 2402, 447, 0561, 536, 1601, 471, 6802, 621, 430
[1, 3)60, 480679, 680174, 374678, 950738, 2000
[3, 5)0168, 9600000
[5, 7)446, 4001920000
[7, 9)144, 96000000
[9, 11)33, 75000000
[11, 13)000000
[13, 15)000000
[15, 17)1200000
[17, 19)0000012
[19, 21)0012000
[21, 100]000000
total2, 621, 4422, 621, 4422, 621, 4422, 621, 4422, 621, 4422, 621, 442
Table 9. The number of grids with refined Dlength interval of HR and OURS grid.
Table 9. The number of grids with refined Dlength interval of HR and OURS grid.
Interval (%)HROURS
Np (%)Np (%)
[−100, −1.1]160.0000
(−1.1, −0.9]46, 5921.7800
(−0.9, −0.7]137, 7765.2600
(−0.7, −0.5]835, 68031.8800
(−0.5, −0.3]468, 57617.87325, 32812.41
(−0.3, −0.1]288, 70411.01659, 53625.16
(−0.1,0.1)211, 2328.06676, 51225.81
[0.1, 0.3)166, 4326.35705, 72826.92
[0.3, 0.5)134, 1925.12234, 5288.95
[0.5, 0.7)110, 1184.2019, 7980.76
[0.7, 0.9)91, 3283.4800
[0.9, 1.0)73, 3282.8000
[1.1, 1.3)50, 6401.9300
[1.3, 1.5)6, 8160.2600
[1.5, 1.7)0000
[1.7, 2.3)00120
[2.3, 100]12000
total2, 621, 4421002, 621, 442100
Table 10. The cumulative proportion of |Dlength| of different grids (%).
Table 10. The cumulative proportion of |Dlength| of different grids (%).
Interval (%)NOPTSPRGHRSCVTXUOURS
[0, 0.02)00.512.200.731.325.02
[0.02, 0.10)0.023.4111.014.064.3225.81
[0.10, 0.18)0.795.9320.377.547.4745.63
[0.18, 0.26)1.789.0130.8710.7710.3366.58
[0.26, 0.34)3.8611.6544.4013.8413.4086.02
[0.34, 0.42)3.8614.1065.4917.8016.9697.28
[0.42, 0.50)3.8617.3675.1721.6420.4799.24
[0.50, 0.58)5.5120.4380.0325.2024.4599.96
[0.58, 0.66)5.5823.8083.9029.6628.6699.99
[0.66, 0.74)7.8226.9787.0034.4232.7599.99
[0.74, 0.82)12.5430.3489.5840.1737.3999.99
[0.82, 2.00)22.7266.5999.9990.0089.2399.99
[2.00, 5.00)76.1598.6799.9997.6297.6299.99
[5.00, 7.00)92.3099.6399.9998.9798.8399.99
[7.00, 100]100100100100100100
Table 11. L2-norm error of Laplacian operator of different grids.
Table 11. L2-norm error of Laplacian operator of different grids.
2345678
NOPT1.31 × 10−13.78 × 10−21.20 × 10−24.49 × 10−31.96 × 10−39.26 × 10−44.52 × 10−4
SPRG1.33 × 10−13.56 × 10−29.47 × 10−32.74 × 10−39.30 × 10−43.84 × 10−41.47 × 10−4
HR1.41 × 10−13.79 × 10−29.85 × 10−32.58 × 10−37.23 × 10−42.05 × 10−46.05 × 10−5
SCVT1.33 × 10−13.55 × 10−29.47 × 10−32.81 × 10−31.04 × 10−34.61 × 10−42.23 × 10−4
XU1.33 × 10−13.55 × 10−29.47 × 10−32.81 × 10−31.04 × 10−34.61 × 10−42.23 × 10−4
OURS1.38 × 10−13.75 × 10−29.78 × 10−32.57 × 10−36.98 × 10−41.99 × 10−45.42 × 10−5
Table 12. L-norm error of Laplacian operator of different grids.
Table 12. L-norm error of Laplacian operator of different grids.
2345678
NOPT3.52 × 10−11.28 × 10−18.08 × 10−28.89 × 10−29.10 × 10−29.10 × 10−29.15 × 10−2
SPRG2.81 × 10−17.62 × 10−23.56 × 10−23.23 × 10−23.00 × 10−22.70 × 10−22.57 × 10−2
HR2.84 × 10−18.67 × 10−23.23 × 10−21.41 × 10−26.75 × 10−33.31 × 10−31.66 × 10−3
SCVT2.81 × 10−17.61 × 10−23.92 × 10−23.88 × 10−23.87 × 10−23.87 × 10−23.87 × 10−2
XU2.81 × 10−17.61 × 10−23.92 × 10−23.88 × 10−23.87 × 10−23.87 × 10−23.87 × 10−2
OURS2.85 × 10−17.81 × 10−22.93 × 10−21.35 × 10−26.59 × 10−33.27 × 10−31.59 × 10−3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, F.; Zhao, X.; Sun, W.; Li, Y.; Duan, Y. An Overall Uniformity Optimization Method of the Spherical Icosahedral Grid Based on the Optimal Transformation Theory. Atmosphere 2021, 12, 1516. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12111516

AMA Style

Luo F, Zhao X, Sun W, Li Y, Duan Y. An Overall Uniformity Optimization Method of the Spherical Icosahedral Grid Based on the Optimal Transformation Theory. Atmosphere. 2021; 12(11):1516. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12111516

Chicago/Turabian Style

Luo, Fuli, Xuesheng Zhao, Wenbin Sun, Yalu Li, and Yuanzheng Duan. 2021. "An Overall Uniformity Optimization Method of the Spherical Icosahedral Grid Based on the Optimal Transformation Theory" Atmosphere 12, no. 11: 1516. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12111516

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop