Reliable low precision simulations in land surface models

Weather and climate models must continue to increase in both resolution and complexity in order that forecasts become more accurate and reliable. Moving to lower numerical precision may be an essential tool for coping with the demand for ever increasing model complexity in addition to increasing computing resources. However, there have been some concerns in the weather and climate modelling community over the suitability of lower precision for climate models, particularly for representing processes that change very slowly over long time-scales. These processes are difficult to represent using low precision due to time increments being systematically rounded to zero. Idealised simulations are used to demonstrate that a model of deep soil heat diffusion that fails when run in single precision can be modified to work correctly using low precision, by splitting up the model into a small higher precision part and a low precision part. This strategy retains the computational benefits of reduced precision whilst preserving accuracy. This same technique is also applied to a full complexity land surface model, resulting in rounding errors that are significantly smaller than initial condition and parameter uncertainties. Although lower precision will present some problems for the weather and climate modelling community, many of the problems can likely be overcome using a straightforward and physically motivated application of reduced precision.

[1]  Julien Langou,et al.  Accelerating scientific computations with mixed precision algorithms , 2008, Comput. Phys. Commun..

[2]  Kipton Barros,et al.  Solving lattice QCD systems of equations using mixed precision solvers on GPUs , 2009, Comput. Phys. Commun..

[3]  Peter D. Düben,et al.  Single Precision in Weather Forecasting Models: An Evaluation with the IFS , 2017 .

[4]  Florian Pappenberger,et al.  Improved seasonal prediction of the hot summer of 2003 over Europe through better representation of uncertainty in the land surface , 2016 .

[5]  Peter D. Düben,et al.  rpe v5: an emulator for reduced floating-point precision in large numerical simulations , 2016 .

[6]  Wanqiu Wang,et al.  A GCM assessment of atmospheric seasonal predictability associated with soil moisture anomalies over North America , 1998 .

[7]  Peter D. Düben,et al.  Benchmark Tests for Numerical Weather Forecasts on Inexact Hardware , 2014 .

[8]  Diana Verseghy,et al.  The reliability of single precision computations in the simulation of deep soil heat diffusion in a land surface model , 2016, Climate Dynamics.

[9]  R. Carsel,et al.  Developing joint probability distributions of soil water retention characteristics , 1988 .

[10]  B. Hurk,et al.  A Revised Hydrology for the ECMWF Model: Verification from Field Site to Terrestrial Water Storage and Impact in the Integrated Forecast System , 2009 .

[11]  Félix Francés,et al.  Spatial scale effect on the upper soil effective parameters of a distributed hydrological model , 2012 .

[12]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[13]  Florian Pappenberger,et al.  Evaluating uncertainty in estimates of soil moisture memory with a reverse ensemble approach , 2016 .

[14]  Roberto Buizza,et al.  The new ECMWF seasonal forecast system (system 4) , 2011 .

[15]  G. Balsamo,et al.  The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA‐Interim reanalysis data , 2014 .

[16]  Robert Strzodka,et al.  Performance and accuracy of hardware-oriented native-, emulated- and mixed-precision solvers in FEM simulations , 2007, Int. J. Parallel Emergent Distributed Syst..

[17]  Andreas W. Götz,et al.  SPFP: Speed without compromise - A mixed precision model for GPU accelerated molecular dynamics simulations , 2013, Comput. Phys. Commun..

[18]  Jonathan M. Borwein,et al.  High-precision computation: Mathematical physics and dynamics , 2010, Appl. Math. Comput..

[19]  D. Verseghy,et al.  The Canadian land surface scheme (CLASS): Its history and future , 2000 .

[20]  Wayne Luk,et al.  Architectures and Precision Analysis for Modelling Atmospheric Variables with Chaotic Behaviour , 2015, 2015 IEEE 23rd Annual International Symposium on Field-Programmable Custom Computing Machines.

[21]  Peter D. Düben,et al.  A study of reduced numerical precision to make superparameterization more competitive using a hardware emulator in the OpenIFS model , 2017 .

[22]  Peter D. Düben,et al.  On the use of scale‐dependent precision in Earth System modelling , 2017 .

[23]  F. Molteni,et al.  SEAS5: the new ECMWF seasonal forecast system , 2018, Geoscientific Model Development.

[24]  F. Pappenberger,et al.  ERA-Interim/Land: a global land surface reanalysis data set , 2015 .