An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme （TSPAS）, was generalized and implemented on a spherical icosahedral hexagonal grid （also referred to as a geodesic grid...An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme （TSPAS）, was generalized and implemented on a spherical icosahedral hexagonal grid （also referred to as a geodesic grid） to solve the transport equation. The C grid discretization was used for the spatial discretization. To implement TSPAS on an unstructured grid, the original finite-difference scheme was further generalized. The two-step integration utilizes a combination of two separate schemes （a low-order monotone scheme and a high-order scheme that typically cannot ensure monotonicity） to calculate the fluxes at the cell walls （one scheme corresponds to one cell wall）. The choice between these two schemes for each edge depends on a pre-updated scalar value using slightly increased fluxes. After the determination of an appropriate scheme, the final integration at a target cell is achieved by summing the fluxes that are computed by the different schemes. The conservative and shape-preserving properties of the generalized scheme are demonstrated. Numerical experiments are conducted at several horizontal resolutions. TSPAS is compared with the Flux Corrected Transport （FCT） approach to demonstrate the differences between the two methods, and several transport tests are performed to examine the accuracy, efficiency and robustness of the two schemes.展开更多
Global ocean surface currents estimated from satellite derived data based on a regular global grid are affected by the grid’s shape and placement.Due to different neighbourhood relationships,the rectangular lat/lon g...Global ocean surface currents estimated from satellite derived data based on a regular global grid are affected by the grid’s shape and placement.Due to different neighbourhood relationships,the rectangular lat/lon grids lose accuracy when interpolating andfitting elevation data.Hexagonal grids have shown to be advantageous due to their isotropic,uniform neighbourhood.Considering these merits,this paper aims to estimate global ocean surface current using a global isotropic hexagonal grid from satellite remote sensing data.First,gridded satellite altimeter data and wind data with different resolutions are interpolated into the centre of the global isotropic hexagonal grid.Then,geostrophic and Ekman currents components are estimated according to the Lagerlof Ocean currents theory.Finally,the inversion results are verified.By analyzing the results,we conclude that the ocean surface currents estimated based on the global isotropic hexagonal grid have considerable accuracy,with improvement over rectangular lat/lon grids.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.41505066)the Basic Scientific Research and Operation Foundation of Chinese Academy Meteorological Sciences(Grant Nos.2015Z002,2015Y005)the National Research and Development Key Program:Global Change and Mitigation Strategies(No.2016YFA0602101)
文摘An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme （TSPAS）, was generalized and implemented on a spherical icosahedral hexagonal grid （also referred to as a geodesic grid） to solve the transport equation. The C grid discretization was used for the spatial discretization. To implement TSPAS on an unstructured grid, the original finite-difference scheme was further generalized. The two-step integration utilizes a combination of two separate schemes （a low-order monotone scheme and a high-order scheme that typically cannot ensure monotonicity） to calculate the fluxes at the cell walls （one scheme corresponds to one cell wall）. The choice between these two schemes for each edge depends on a pre-updated scalar value using slightly increased fluxes. After the determination of an appropriate scheme, the final integration at a target cell is achieved by summing the fluxes that are computed by the different schemes. The conservative and shape-preserving properties of the generalized scheme are demonstrated. Numerical experiments are conducted at several horizontal resolutions. TSPAS is compared with the Flux Corrected Transport （FCT） approach to demonstrate the differences between the two methods, and several transport tests are performed to examine the accuracy, efficiency and robustness of the two schemes.
文摘Global ocean surface currents estimated from satellite derived data based on a regular global grid are affected by the grid’s shape and placement.Due to different neighbourhood relationships,the rectangular lat/lon grids lose accuracy when interpolating andfitting elevation data.Hexagonal grids have shown to be advantageous due to their isotropic,uniform neighbourhood.Considering these merits,this paper aims to estimate global ocean surface current using a global isotropic hexagonal grid from satellite remote sensing data.First,gridded satellite altimeter data and wind data with different resolutions are interpolated into the centre of the global isotropic hexagonal grid.Then,geostrophic and Ekman currents components are estimated according to the Lagerlof Ocean currents theory.Finally,the inversion results are verified.By analyzing the results,we conclude that the ocean surface currents estimated based on the global isotropic hexagonal grid have considerable accuracy,with improvement over rectangular lat/lon grids.
