Improving the numerical solution of the energy equation in land models
Bibliographic record
Abstract
At its core, a hydrological model is comprised of the conservation of energy and mass for a myriad of modeled processes across a suite of spatio-temporal scales. Simulations over North America with the SUMMA hydrological model show that the form of the energy equation most commonly used in land models produces both large violations in energy conservation (especially in cold regions) as well as larger numerical errors in soil temperature and soil water content than is possible with more robust solvers. These numerical issues sabotage the success of efforts to improve process-representation. We present improved energy-conserving solutions for land models, testing five approaches over North America with the SUMMA model and evaluating tradeoffs between strict energy conservation and numerical errors in the energy equation. We include approaches that do not use time integration methods with rigorous error control (as is common in hydrological models) as well as approaches that do. The mixed form of the energy equation is discretized to conserve energy to within machine precision. Alternatively, the direct solution of the energy equation (i.e., using enthalpy as a primary variable) yields the smallest numerical errors because it allows error control to be placed on the inherent state variable. In the spirit of advancing process-representation, we illustrate the importance of accurate energy balance solutions for simulations of partially frozen soils, permafrost, and glaciers. In one prominent example, we demonstrate that debris-covered glaciers have substantially dissimilar runoff contributions when evolved using different solutions to the energy equation. The capability to accurately simulate the energy balance of terrestrial systems is essential to improve the theoretical underpinnings of process-based hydrologic models.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame distilled prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. Learned from the 10,348 direct Codex labels and 10,348 direct Gemma labels. Candidate is the union of thresholded teacher heads; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels or direct frontier model labels.
Codex and Gemma teacher scores by category
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.000 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one teacher head, not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".