Supplementary material to "Neodymium isotopes in the ocean model of the Community Earth System Model (CESM1.3)"

Author(s):  
Sifan Gu ◽  
Zhengyu Liu ◽  
Alexandra Jahn ◽  
Johannes Rempfer ◽  
Jiaxu Zhang ◽  
...  
2016 ◽  
Vol 9 (11) ◽  
pp. 4209-4225 ◽  
Author(s):  
Xiaomeng Huang ◽  
Qiang Tang ◽  
Yuheng Tseng ◽  
Yong Hu ◽  
Allison H. Baker ◽  
...  

Abstract. In the Community Earth System Model (CESM), the ocean model is computationally expensive for high-resolution grids and is often the least scalable component for high-resolution production experiments. The major bottleneck is that the barotropic solver scales poorly at high core counts. We design a new barotropic solver to accelerate the high-resolution ocean simulation. The novel solver adopts a Chebyshev-type iterative method to reduce the global communication cost in conjunction with an effective block preconditioner to further reduce the iterations. The algorithm and its computational complexity are theoretically analyzed and compared with other existing methods. We confirm the significant reduction of the global communication time with a competitive convergence rate using a series of idealized tests. Numerical experiments using the CESM 0.1° global ocean model show that the proposed approach results in a factor of 1.7 speed-up over the original method with no loss of accuracy, achieving 10.5 simulated years per wall-clock day on 16 875 cores.


2014 ◽  
Vol 7 (6) ◽  
pp. 7461-7503 ◽  
Author(s):  
A. Jahn ◽  
K. Lindsay ◽  
X. Giraud ◽  
N. Gruber ◽  
B. L. Otto-Bliesner ◽  
...  

Abstract. Carbon isotopes in the ocean are frequently used as paleo climate proxies and as present-day geochemical ocean tracers. In order to allow a more direct comparison of climate model results with this large and currently underutilized dataset, we added a carbon isotope module to the ocean model of the Community Earth System Model (CESM), containing the cycling of the stable isotope 13C and the radioactive isotope 14C. We implemented the 14C tracer in two ways: in the "abiotic" case, the 14C tracer is only subject to air–sea gas exchange, physical transport, and radioactive decay, while in the "biotic" version, the 14C additionally follows the 13C tracer through all biogeochemical and ecological processes. Thus, the abiotic 14C tracer can be run without the ecosystem module, requiring significantly less computational resources. The carbon isotope module calculates the carbon isotopic fractionation during gas exchange, photosynthesis, and calcium carbonate formation, while any subsequent biological process such as remineralization as well as any external inputs are assumed to occur without fractionation. Given the uncertainty associated with the biological fractionation during photosynthesis, we implemented and tested three parameterizations of different complexity. Compared to present-day observations, the model is able to simulate the oceanic 14C bomb uptake and the 13C Suess effect reasonably well compared to observations and other model studies. At the same time, the carbon isotopes reveal biases in the physical model, for example a too sluggish ventilation of the deep Pacific Ocean.


2017 ◽  
Vol 10 (12) ◽  
pp. 4723-4742 ◽  
Author(s):  
Sifan Gu ◽  
Zhengyu Liu

Abstract. The sediment 231Pa ∕ 230Th activity ratio is emerging as an important proxy for deep ocean circulation in the past. In order to allow for a direct model–data comparison and to improve our understanding of the sediment 231Pa ∕ 230Th activity ratio, we implement 231Pa and 230Th in the ocean component of the Community Earth System Model (CESM). In addition to the fully coupled implementation of the scavenging behavior of 231Pa and 230Th with the active marine ecosystem module (particle-coupled: hereafter p-coupled), another form of 231Pa and 230Th have also been implemented with prescribed particle flux fields of the present climate (particle-fixed: hereafter p-fixed). The comparison of the two forms of 231Pa and 230Th helps to isolate the influence of the particle fluxes from that of ocean circulation. Under present-day climate forcing, our model is able to simulate water column 231Pa and 230Th activity and the sediment 231Pa ∕ 230Th activity ratio in good agreement with available observations. In addition, in response to freshwater forcing, the p-coupled and p-fixed sediment 231Pa ∕ 230Th activity ratios behave similarly over large areas of low productivity on long timescales, but can differ substantially in some regions of high productivity and on short timescales, indicating the importance of biological productivity in addition to ocean transport. Therefore, our model provides a potentially powerful tool to help the interpretation of sediment 231Pa ∕ 230Th reconstructions and to improve our understanding of past ocean circulation and climate changes.


Sign in / Sign up

Export Citation Format

Share Document