A Barotropic Solver for High-Resolution Ocean General Circulation Models

High-resolution global ocean general circulation models (OGCMs) play a key role in accurate ocean forecasting. However, the models of the operational forecasting systems are still not in high resolution due to the subsequent high demand for large computation, as well as the low parallel efficiency barrier. Good scalability is an important index of parallel efficiency and is still a challenge for OGCMs. We found that the communication cost in a barotropic solver, namely, the preconditioned conjugate gradient (PCG) method, is the key bottleneck for scalability due to the high frequency of the global reductions. In this work, we developed a new algorithm—a communication-avoiding Krylov subspace method with a PCG (CA-PCG)—to improve scalability and then applied it to the Nucleus for European Modelling of the Ocean (NEMO) as an example. For PCG, inner product operations with global communication were needed in every iteration, while for CA-PCG, inner product operations were only needed every eight iterations. Therefore, the global communication cost decreased from more than 94.5% of the total execution time with PCG to less than 63.4% with CA-PCG. As a result, the execution time of the barotropic modes decreased from more than 17,000 s with PCG to less than 6000 s with CA-PCG, and the total execution time decreased from more than 18,000 s with PCG to less than 6200 s with CA-PCG. Besides, the ratio of the speedup can also be increased from 3.7 to 4.6. In summary, the high process count scalability when using CA-PCG was effectively improved from that using the PCG method, providing a highly effective solution for accurate ocean simulation.

Simulation of global distribution of rare earth elements in the ocean using an ocean general circulation model

In this study, we report our ocean general circulation model simulations of the global distribution of rare earth elements (REEs) in the ocean. As previously reported (Oka et al. in Glob Biogeochem Cycles 23:1–16, 2009), the vertical profiles of REEs in the North Pacific Ocean are strongly controlled by the reversible scavenging process, and the systematic differences between REEs can be reproduced in the model by selecting an appropriate model parameter which controls affinity to particles. We here demonstrate that the external REE input from the coastal regions also plays a role in controlling the vertical profiles of dissolved REE and their inter-basin differences. The role of the external inputs is especially important for light REEs, such as neodymium (Nd). The linear increase in Nd concentration in the North Pacific Ocean cannot be sufficiently reproduced by the reversible scavenging alone; rather, a combination of the reversible scavenging and the external inputs is necessary. On the other hand, the distribution of heavy REEs, such as lutetium (Lu), can be broadly reproduced without the external inputs, suggesting that Lu has similarity with conservative nutrient-like tracer. When compared with REE observations compiled from both the recently obtained GEOTRACES dataset and pre-GEOTRACES reported data, our simulations successfully reproduced the overall features of these observations. Observational data suggested that the vertical profiles of REEs are not the same among the basins; our model simulations demonstrate that this feature can be clearly reproduced by considering both the reversible scavenging and the external REE inputs from the coastal regions.

Estimation of 137Cs inventories by a global ocean general circulation model for the global database interpolation

<p>Artificial radionuclide <sup>137</sup>Cs has been supplied into the ocean by global fallout due to atmospheric nuclear weapons tests since 1945, releases from reprocessing plants since 1952, and most recently by fallout and discharge due to the Fukushima Dai-ichi Nuclear Power Plant (1F NPP) accident since 2011.<sup>137</sup>Cs activities measured for scientific purposes as well as environmental health and safety monitoring have been summarized in a historical database by IAEA. The spatio-temporal density of the observations varies widely, therefore simulation by an ocean general circulation model (OGCM) can be helpful in the interpretation of these observations. We used the Parallel Ocean Program version 2 (POP2) of the Community Earth System Model version 2 (CESM2). The horizontal resolution is 1.125 degrees in longitude and 0.28 to 0.54 degrees in latitude. The simulation period was from 1945 to 2030, and the atmospheric conditions were forced to cycle through repeating normal years. The purposes of this study are to investigate the effect of the release from the reprocessing plants on the distribution of <sup>137</sup>Cs activity by global fallout in the Atlantic Ocean, and the effect of the release derived from the 1F NPP accident on the one by global fallout in the Pacific Ocean.</p><p>The simulated <sup>137</sup>Cs activities were in good agreement with the observed data in the database in the Atlantic Ocean and the Pacific Ocean. The simulated <sup>137</sup>Cs activity immediately after each release event in the North Pacific were inconsistent with the observed one because of the inadequate reproduction of the Kuroshio Current in this quasi-resolution ocean model. However, the influence of the dilution effect is expected to become smaller as the time after the release increases. The influence of the <sup>137</sup>Cs activity by release from the reprocessing plant on the one by global fallout in the Atlantic Ocean is limited to the northeast coast of the European continent and the Marginal Seas. It was also suggested that <sup>137</sup>Cs activity by global fallout has made detection difficult since the 1990s.The influence of the <sup>137</sup>Cs activity by the 1F NPP on the one by global fallout was found to be broadened by the Kuroshio extension area and extended to the California coast. This distribution was similar to that of the one by global fallout. However, there are few observed data off the California coast after 2011. It was also suggested that <sup>137</sup>Cs activity by global fallout has made detection difficult since the 2020 in the Pacific Ocean.</p><p>Even after 2020, it is still possible to detect <sup>137</sup>Cs activity by global fallout in the global ocean. The difference in the vertical distribution between the Pacific and Atlantic oceans reflects the ocean circulation, which is useful for the validation of ocean general circulation models. There is still room for improvement in setting the input conditions to the ocean for each event.</p>

Climate change effects Management with the approach of the uncertainty of Atmosphere-Ocean General Circulation Models in Hamadan Province, Iran

INTRODUCTION: Since Iran is located in the semi-arid belt, it has faced such issues as drought, dust crisis, and intensified migration. The assessment of the effects of climate change includes identifying some key aspects of uncertainties used to estimate its impacts, such as uncertainties in the context of Atmosphere-Ocean General Circulation Models (AOGCMs): in regional-scale climatology, in statistical or dynamic downscaling methods, and parametric and structural uncertainties in different models. One of the most important sources of uncertainty in climate change is the use of different AOGCMs that produce different outputs for climate variables. METHODS: In this study, to investigate the uncertainty of AOGCM models, the downscaled data of the NASA Earth Exchange Global Daily Downscaled Projections dataset obtained from 21 AOGCMs with medium Representative Concentration Pathway4.5 scenario were downloaded from the NASA site for 81 cells in Hamadan Province, Iran. After the validation of the models, they were evaluated against the criteria of the coefficient of determination and model efficiency coefficient in comparison with the data of the Hamedan synoptic station in the statistical period of 1976-2005. To reduce the uncertainty of AOGCMs, the ensemble performance (EP) of models was used in Climate Data Operators software. FINDINGS: It was revealed that MRI-CGCM3, MPI-ESM-LR, BNU-ESM, ACCESS1-0, MIROC-ESM, MIROC-ESM-CHEM, and MPI-ESM-MR models had better performance than similar models. It was also found that IPSL-CM5A-LR, CNRM-CM5, CSIRO-Mk3-6-0, CESM1-BGC, and GFDL-ESM2M had the lowest correlation between observational and simulation data of mean monthly precipitation. CONCLUSION: According to the results, this method could provide a good estimate in the base period (1976-2005), compared to the data of the Hamedan synoptic station, and was more accurate compared to the single implementation method of each AOGCM model. The results of EP of models in the future period (2020-2049) showed that precipitation will not change considerably in the future and will increase by 0.23 mm. In addition, the average, maximum, and minimum annual temperatures will increase by 1.54°C, 1.7°C, and 1.40°C, respectively.

FORTE 2.0: a fast, parallel and flexible coupled climate model

FORTE 2.0 is an intermediate-resolution coupled atmosphere–ocean general circulation model (AOGCM) consisting of the Intermediate General Circulation Model 4 (IGCM4), a T42 spectral atmosphere with 35σ layers, coupled to Modular Ocean Model – Array (MOMA), a 2∘ × 2∘ ocean with 15 z-layer depth levels. Sea ice is represented by a simple flux barrier. Both the atmosphere and ocean components are coded in Fortran. It is capable of producing a stable climate for long integrations without the need for flux adjustments. One flexibility afforded by the IGCM4 atmosphere is the ability to configure the atmosphere with either 35σ layers (troposphere and stratosphere) or 20σ layers (troposphere only). This enables experimental designs for exploring the roles of the troposphere and stratosphere, and the faster integration of the 20σ layer configuration enables longer duration studies on modest hardware. A description of FORTE 2.0 is given, followed by the analysis of two 2000-year control integrations, one using the 35σ configuration of IGCM4 and one using the 20σ configuration.