Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS
Abstract Numerical simulation of land models without error control can be highly inaccurate. We present the incorporation of the Suite of Nonlinear and Differential‐Algebraic Equation Solvers (SUNDIALS) package to solve the equations that simulate thermodynamics and hydrologic processes in the Struc...
Saved in:
| Main Authors: | , , , , , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
American Geophysical Union (AGU)
2024-12-01
|
| Series: | Journal of Advances in Modeling Earth Systems |
| Online Access: | https://doi.org/10.1029/2024MS004256 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1850079743125225472 |
|---|---|
| author | Raymond J. Spiteri Ashley E. VanBeusekom Kyle Klenk Reza Zolfaghari Sean J. Trim Wouter J. M. Knoben Andrew M. Ireson Martyn P. Clark |
| author_facet | Raymond J. Spiteri Ashley E. VanBeusekom Kyle Klenk Reza Zolfaghari Sean J. Trim Wouter J. M. Knoben Andrew M. Ireson Martyn P. Clark |
| author_sort | Raymond J. Spiteri |
| collection | DOAJ |
| description | Abstract Numerical simulation of land models without error control can be highly inaccurate. We present the incorporation of the Suite of Nonlinear and Differential‐Algebraic Equation Solvers (SUNDIALS) package to solve the equations that simulate thermodynamics and hydrologic processes in the Structure for Unifying Multiple Modeling Alternatives (SUMMA) land model. The algorithmic features of SUNDIALS, such as error estimation and adaptive order and step‐size control, result in a SUMMA‐SUNDIALS model that delivers substantially improved accuracy and relative computational efficiency compared to integration with the previous SUMMA model, which uses the low‐order backward Euler method with no rigorous error control. The results are demonstrated through simulations over the North American continent with more than 500,000 spatial elements. Compared to the previous SUMMA model, we find that the simulations produced by the SUMMA‐SUNDIALS model are orders of magnitude closer to converged solutions for the same computational cost. Being able to efficiently perform more reliable simulations makes the SUMMA‐SUNDIALS model a powerful tool for improving our understanding of the terrestrial component of the Earth System. |
| format | Article |
| id | doaj-art-9d8642606d724366b8e23d7f9d1780bc |
| institution | DOAJ |
| issn | 1942-2466 |
| language | English |
| publishDate | 2024-12-01 |
| publisher | American Geophysical Union (AGU) |
| record_format | Article |
| series | Journal of Advances in Modeling Earth Systems |
| spelling | doaj-art-9d8642606d724366b8e23d7f9d1780bc2025-08-20T02:45:07ZengAmerican Geophysical Union (AGU)Journal of Advances in Modeling Earth Systems1942-24662024-12-011612n/an/a10.1029/2024MS004256Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALSRaymond J. Spiteri0Ashley E. VanBeusekom1Kyle Klenk2Reza Zolfaghari3Sean J. Trim4Wouter J. M. Knoben5Andrew M. Ireson6Martyn P. Clark7Department of Computer Science University of Saskatchewan Saskatoon SK CanadaCentre for Hydrology University of Saskatchewan Canmore AB CanadaDepartment of Computer Science University of Saskatchewan Saskatoon SK CanadaDepartment of Computer Science University of Saskatchewan Saskatoon SK CanadaDepartment of Computer Science University of Saskatchewan Saskatoon SK CanadaDepartment of Civil Engineering University of Calgary Calgary AB CanadaCentre for Hydrology University of Saskatchewan Canmore AB CanadaDepartment of Civil Engineering University of Calgary Calgary AB CanadaAbstract Numerical simulation of land models without error control can be highly inaccurate. We present the incorporation of the Suite of Nonlinear and Differential‐Algebraic Equation Solvers (SUNDIALS) package to solve the equations that simulate thermodynamics and hydrologic processes in the Structure for Unifying Multiple Modeling Alternatives (SUMMA) land model. The algorithmic features of SUNDIALS, such as error estimation and adaptive order and step‐size control, result in a SUMMA‐SUNDIALS model that delivers substantially improved accuracy and relative computational efficiency compared to integration with the previous SUMMA model, which uses the low‐order backward Euler method with no rigorous error control. The results are demonstrated through simulations over the North American continent with more than 500,000 spatial elements. Compared to the previous SUMMA model, we find that the simulations produced by the SUMMA‐SUNDIALS model are orders of magnitude closer to converged solutions for the same computational cost. Being able to efficiently perform more reliable simulations makes the SUMMA‐SUNDIALS model a powerful tool for improving our understanding of the terrestrial component of the Earth System.https://doi.org/10.1029/2024MS004256 |
| spellingShingle | Raymond J. Spiteri Ashley E. VanBeusekom Kyle Klenk Reza Zolfaghari Sean J. Trim Wouter J. M. Knoben Andrew M. Ireson Martyn P. Clark Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS Journal of Advances in Modeling Earth Systems |
| title | Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS |
| title_full | Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS |
| title_fullStr | Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS |
| title_full_unstemmed | Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS |
| title_short | Accurate and Efficient Numerical Simulation of Land Models Using SUMMA With SUNDIALS |
| title_sort | accurate and efficient numerical simulation of land models using summa with sundials |
| url | https://doi.org/10.1029/2024MS004256 |
| work_keys_str_mv | AT raymondjspiteri accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT ashleyevanbeusekom accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT kyleklenk accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT rezazolfaghari accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT seanjtrim accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT wouterjmknoben accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT andrewmireson accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials AT martynpclark accurateandefficientnumericalsimulationoflandmodelsusingsummawithsundials |