scholarly journals On the use of Schwarz–Christoffel conformal mappings to the grid generation for global ocean models

2015 ◽  
Vol 8 (10) ◽  
pp. 3471-3485 ◽  
Author(s):  
S. Xu ◽  
B. Wang ◽  
J. Liu

Abstract. In this article we propose two grid generation methods for global ocean general circulation models. Contrary to conventional dipolar or tripolar grids, the proposed methods are based on Schwarz–Christoffel conformal mappings that map areas with user-prescribed, irregular boundaries to those with regular boundaries (i.e., disks, slits, etc.). The first method aims at improving existing dipolar grids. Compared with existing grids, the sample grid achieves a better trade-off between the enlargement of the latitudinal–longitudinal portion and the overall smooth grid cell size transition. The second method addresses more modern and advanced grid design requirements arising from high-resolution and multi-scale ocean modeling. The generated grids could potentially achieve the alignment of grid lines to the large-scale coastlines, enhanced spatial resolution in coastal regions, and easier computational load balance. Since the grids are orthogonal curvilinear, they can be easily utilized by the majority of ocean general circulation models that are based on finite difference and require grid orthogonality. The proposed grid generation algorithms can also be applied to the grid generation for regional ocean modeling where complex land–sea distribution is present.

2015 ◽  
Vol 8 (2) ◽  
pp. 1337-1373
Author(s):  
S. Xu ◽  
B. Wang ◽  
J. Liu

Abstract. In this article we propose two conformal mapping based grid generation algorithms for global ocean general circulation models (OGCMs). Contrary to conventional, analytical forms based dipolar or tripolar grids, the new algorithms are based on Schwarz–Christoffel (SC) conformal mapping with prescribed boundary information. While dealing with the basic grid design problem of pole relocation, these new algorithms also address more advanced issues such as smoothed scaling factor, or the new requirements on OGCM grids arisen from the recent trend of high-resolution and multi-scale modeling. The proposed grid generation algorithm could potentially achieve the alignment of grid lines to coastlines, enhanced spatial resolution in coastal regions, and easier computational load balance. Since the generated grids are still orthogonal curvilinear, they can be readily utilized in existing Bryan–Cox–Semtner type ocean models. The proposed methodology can also be applied to the grid generation task for regional ocean modeling where complex land–ocean distribution is present.


2007 ◽  
Vol 20 (17) ◽  
pp. 4293-4315 ◽  
Author(s):  
Jianjun Yin ◽  
Ronald J. Stouffer

Abstract Two coupled atmosphere–ocean general circulation models developed at GFDL show differing stability properties of the Atlantic thermohaline circulation (THC) in the Coupled Model Intercomparison Project/Paleoclimate Modeling Intercomparison Project (CMIP/PMIP) coordinated “water-hosing” experiment. In contrast to the R30 model in which the “off” state of the THC is stable, it is unstable in the CM2.1. This discrepancy has also been found among other climate models. Here a comprehensive analysis is performed to investigate the causes for the differing behaviors of the THC. In agreement with previous work, it is found that the different stability of the THC is closely related to the simulation of a reversed thermohaline circulation (RTHC) and the atmospheric feedback. After the shutdown of the THC, the RTHC is well developed and stable in R30. It transports freshwater into the subtropical North Atlantic, preventing the recovery of the salinity and stabilizing the off mode of the THC. The flux adjustment is a large term in the water budget of the Atlantic Ocean. In contrast, the RTHC is weak and unstable in CM2.1. The atmospheric feedback associated with the southward shift of the Atlantic ITCZ is much more significant. The oceanic freshwater convergence into the subtropical North Atlantic cannot completely compensate for the evaporation, leading to the recovery of the THC in CM2.1. The rapid salinity recovery in the subtropical North Atlantic excites large-scale baroclinic eddies, which propagate northward into the Nordic seas and Irminger Sea. As the large-scale eddies reach the high latitudes of the North Atlantic, the oceanic deep convection restarts. The differences in the southward propagation of the salinity and temperature anomalies from the hosing perturbation region in R30 and CM2.1, and associated different development of a reversed meridional density gradient in the upper South Atlantic, are the cause of the differences in the behavior of the RTHC. The present study sheds light on important physical and dynamical processes in simulating the dynamical behavior of the THC.


2000 ◽  
Vol 14 (4) ◽  
pp. 1267-1281 ◽  
Author(s):  
J. L. Sarmiento ◽  
P. Monfray ◽  
E. Maier-Reimer ◽  
O. Aumont ◽  
R. J. Murnane ◽  
...  

2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).


Sign in / Sign up

Export Citation Format

Share Document