How to Identify Good Coupling Methods With Error Analysis
The simulation of complex systems like the global atmosphere requires coupling of many physics processes. These include, for example, atmospheric motions, cloud and rain formation processes, as well as pollutant transport and transformation. Atmospheric models typically couple these processes by adopting methods found in earlier models or textbooks, but the implications of the methods and their impacts on the simulation results are often unclear. This study introduces a mathematically rigorous analysis framework that provides insight into the features and implications of coupling method choices. As an example, the framework is applied to the Energy Exascale Earth System Model version 1 (E3SMv1) to evaluate two coupling methods used for aerosols.
The framework enables scientists to evaluate the impacts of coupling method choices without implementing all possible methods or deriving lengthy mathematical expressions tailored to each method’s details. Furthermore, the framework provides an accuracy assessment for each interacting process. This can help avoid incidental error cancelation and increase confidence in estimating the overall accuracy of a numerical solution.
The need to choose, update, and evaluate coupling methods during the development of global atmospheric models for weather, climate, and Earth system prediction prompted researchers to launch this study. Many of the existing mathematical analyses focused on details of specific physics problems and numerical models; some analyses were more general but produced mathematical expressions that were hard to interpret by physical scientists. This study introduces an analysis framework that is general and intuitive to use. The framework ignores some details of the mathematical formulation and numerical algorithms used for each interacting process and focuses on coupling issues. The study starts by describing two fundamental sources of numerical error as two processes are coupled. It then demonstrates how these insights can be applied through simple arithmetic to derive numerical errors in multi-process problems. Two coupling methods used for the aerosol life cycles in E3SMv1 are analyzed as a concrete example of applying the tool.
The computational resources used in this study were provided by the National Energy Research Scientific Computing Center, a U.S. Department of Energy (DOE) Office of Science user facility located at Lawrence Berkeley National Laboratory, and by the Compy supercomputer operated for DOE by Pacific Northwest National Laboratory.