Tuesday, April 16, 2024

Elementary Mathematics sheds light on plant Transpiration

By examining the derivation of Penn-Monteith-like equations for estimating evapotranspiration, one can uncover valuable insights into plant functionality. In essence, equations talk. For a more comprehensive and in-depth exploration of this topic, refer to this  erlier post
However, here you can find also the Jupyter  Notebooks and the data that were used to produce the figures in the presentation. The presentation itself can be found by clicking on the figure above.

Friday, April 5, 2024

A series of talks and material on Transit (Travel) time, Residence time and Response Time

Here below we started a little series of lectures about a statistical way of seeing water movements in catchments that, while having a long history (e.g. Niemi, 1977, Rigon et al, 2016) has been largely renewed recently starting from Botter et al., 2010 and Botter et al., 2011. The material is the same prepared for the Hydrological Modelling class however grouped here separately for the readers convenience. 

An alternative perspective is presented here regarding their concepts. While certain passages may pose some challenges, the enhanced comprehension of flux formation processes at the catchment scale is, in my opinion, immensely valuable and well worth the effort. The proposed approach involves the following line of thinking: a) the collective fluxes within catchments result from the cumulative movements of numerous small water volumes (water parcels); b) parcels can be understood through three key distributions: the travel time distribution, the residence time distribution, and the response time distribution; c) the interrelations among these distributions are elucidated; d) linking these distributions to catchment processes is achieved through the formulation of age-ranked distributions within ordinary differential equations; e) the theory developed here represents a generalization of the unit hydrograph theory.

Some References
  • Benettin, P., Soulsby, C., Birkel, C., Tetzlaff, D., , G. and Rinaldo, A. (2017) Using sas functions and high resolution isotope data to unravel travel time distributions in headwater catchments. Water Resources Research, 53, 1864–1878. URL: http: //doi.org/10.1002/2016WR020117. 
  • Benettin, Paolo, and Enrico Bertuzzo. 2018. “Tran-SAS v1.0: A Numerical Model to Compute Catchment-Scale Hydrologic Transport Using StorAge Selection Functions.” Geoscientific Model Development Discussions, January, 1–19.
  • Benettin, Paolo, Nicolas B. Rodriguez, Matthias Sprenger, Minseok Kim, Julian Klaus, Ciaran J. Harman, Ype van der Velde, et al. 2022. Transit Time Estimation in Catchments: Recent Developments and Future Directions.†Water Resources Research 58 (11). https://doi.org/10.1029/2022wr033096.
  • Botter, Gianluca, Enrico Bertuzzo, and Andrea Rinaldo. 2010. “Transport in the Hydrologic Response: Travel Time Distributions, Soil Moisture Dynamics, and the Old Water Paradox.” Water Resources Research 46 (3). http://doi.wiley.com/10.1029/2009WR008371.
  • Botter, Gianluca, Enrico Bertuzzo, and Andrea Rinaldo. 2011. “Catchment Residence and Travel Time Distributions: The Master Equation.” Geophysical Research Letters 38 (11). http://doi.wiley.com/10.1029/2011GL047666.
  • Drever, Mark C., and Markus Hrachowitz. 2017. “Migration as Flow: Using Hydrological Concepts to Estimate the Residence Time of Migrating Birds from the Daily Counts.” Methods in Ecology and Evolution / British Ecological Society 8 (9): 1146–57.
  • Harman, Ciaran J. 2015. “Time-Variable Transit Time Distributions and Transport: Theory and Application to Storage-Dependent Transport of Chloride in a Watershed.” Water Resources Research 51 (1): 1–30.
  • Harman, Ciaran J., and Esther Xu Fei. 2024. Mesas.py v1.0: A Flexible Python Package for Modeling Solute Transport and Transit Times Using StorAge Selection Functions.†Geoscientific Model Development 17 (2): 477–95. https://doi.org/10.5194/gmd-17-477-2024.
  • Hrachowitz, M., Benettin, P., van Breukelen, B. M., Fovet, O., Howden, N. J. K., Ruiz, L., van der Velde, Y. and Wade, A. (2016) Transit times-the link between hydrology and water quality at the catchment scale: Linking hydrology and transit times. Wiley Interdisciplinary Reviews: Water, 3, 629–657. 
  • McDonnell, Jeffrey J. 2014. The Two Water Worlds Hypothesis: Ecohydrological Separation of Water between Streams and Trees? Wiley Interdisciplinary Reviews: Water, April. http://doi.wiley.com/10.1002/wat2.1027.
  • Niemi, Antti J. 1977. “Residence Time Distributions of Variable Flow Processes.” The International Journal of Applied Radiation and Isotopes 28 (10): 855–60.
  • Rigon, Riccardo, Marialaura Bancheri, and Timothy R. Green. 2016. “Age-Ranked Hydrological Budgets and a Travel Time Description of Catchment Hydrology.” Hydrology and Earth System Sciences 20 (12): 4929–47.
  • Rigon, R., and M. Bancheri. “On the Relations between the Hydrological Dynamical Systems of Water Budget, Travel Time, Response Time and Tracer Concentrations.” http://abouthydrology.blogspot.com/2020/05/equivalences-and-differences-among.html.
  • Sprenger, M., Stumpp, C., Weiler, M., Aeschbach, W., ST, A., Benettin, P., Dubbert, M., Hartmann, A., Hrachowitz, M., Kirchner, J., McDonnel, J., Orlowski, N., Penna, D., Pfahl, S., Rinderer, M., Rodriguez, N., Schmidt, M. and Werner, C. (2019) The demographics of water: A review of water ages in the critical zone. Rev. Geophys., 2018RG000633. 
  • Schwemmle, Robin, and Markus Weiler. 2024. Consistent Modeling of Transport Processes and Travel Times: coupling Soil Hydrologic Processes with StorAge Selection Functions. Water Resources Research 60 (1). https://doi.org/10.1029/2023wr034441.
  • Velde, Y. van der, P. J. J. F. Torfs, S. E. A. T. M. Van der Zee, and R. Uijlenhoet. 2012. “Quantifying Catchment-Scale Mixing and Its Effect on Time-Varying Travel Time Distributions.” Water Resources Research 48 (6): W06536–13.
  • Velde, Ype van der, Ingo Heidbüchel, Steve W. Lyon, Lars Nyberg, Allan Rodhe, Kevin Bishop, and Peter A. Troch. 2014. “Consequences of Mixing Assumptions for Time-Variable Travel Time Distributions.” Hydrological Processes 29 (16): 3460–74.
  • Wilusz, Daniel C., Ciaran J. Harman, and William P. Ball. 2017. “Sensitivity of Catchment Transit Times to Rainfall Variability Under Present and Future Climates.” Water Resources Research 53 (12): 10231–56.

Friday, March 22, 2024

4DHydro website

 4DHydro is a project that came out from a call for tender by  ESA to which we had the pleasure to participate. All the making of the project is, since last week documented on the 4DHydro website that you can find following this link.


Not yet available, soon you'll see here a video explaining what the website is supposed to contain. 

Friday, March 8, 2024

Modelling and Hydrological Modelling

These lecture are actually part of the 2024 course in Hydrological Modelling. However because they can be of some more general interest, I am grouping them also here. They try to review the concepts of modelling in general and when applied to hydrology. In the series of lectures there is also a concise overview of catchment processes. The first lecture image, see below, it a Maurizo Cattelan artwork entitled "A donkey among doctors" which is my attitude when I approach the topic. 


The final idea about the practice to do model is expressed in the paper and in the various posts that regard DARTHs which can be find here.  Among people that more reflected on Hydrological Modelling there is certainly Keith Beven [GS]. To get a glimpse of his contributions, please see this other post which contains some of his relevant papers. 

Thursday, March 7, 2024

Stock and flow diagrams, a different way to represent dynamical systems

Stock and flow diagrams (see also here) are  a way to represent dynamical system which is the same area covered by EPN ((Extended Petri Nets). They were brought to my attention by the talk John Baez gave at  Edinburgh Mathematical Society last December. Fortunately the talk is available on Youtube.

Although I find that the visuals of EPN are more expressive and the accompanying infrastructure is easier for engineers to comprehend, I have come to realize that listening to the talk is incredibly instructive when it comes to realize that EPN falls in the objects of category theory. An intriguing aspect explored in the talk is the representation of open systems within stock-flow graphs. In EPN, it is assumed that a flow box not originating from a place indicates that the system is open. Additionally, when one EPN features an outgoing flow labeled A and another EPN has an input flow with the same label, they can be combined to create a composite graph. However, in this presentation, a new rectangular symbol is introduced for the same purpose.
Personally, I find the solution in EPN more intuitive, although it may be considered less abstract. Nevertheless, I have come to realize that a similar graphical approach could prove beneficial in extending EPN to represent not only ODE systems but also PDE systems.

Wednesday, March 6, 2024

On Hydrological Models and their choice (and a use of the AboutHydrology mailing list)

Initially, I was captivated by the visuals that I could incorporate into my presentations. To my pleasant surprise, I discovered that the AboutHydrology mailing list served as a valuable data source. Remarkably, this platform has been active for approximately a decade (I need to verify the exact date of its inception) and has amassed a wealth of information.

Reproduced from Melsen, 2022

Subsequently, I came across two intriguing papers authored by Melsen, delving into the "sociology of selecting a hydrological model." These papers proved to be quite engaging. Additionally, there are other noteworthy publications exploring similar themes. Notably, among the more recent works, Hamilton et al., 2022, and Horton et al., 2023, deserve special mention.  Please find their citation below. In the paper you can easily recover previous relevant literature. 

References

Hamilton, Serena H., Carmel A. Pollino, Danial S. Stratford, Baihua Fu, and Anthony J. Jakeman. 2022. “Fit-for-Purpose Environmental Modeling: Targeting the Intersection of Usability, Reliability and Feasibility.” Environmental Modelling & Software 148 (February): 105278. https://doi.org/10.1016/j.envsoft.2021.105278.

Horton, Pascal, Bettina Schaefli, and Martina Kauzlaric. 2022. “Why Do We Have so Many Different Hydrological Models? A Review Based on the Case of Switzerland.” WIREs. Water 9 (1). https://doi.org/10.1002/wat2.1574.

Melsen, Lieke A. 2023. “The Modeling Toolkit: How Recruitment Strategies for Modeling Positions Influence Model Progress.” Frontiers in Water 5 (May). https://doi.org/10.3389/frwa.2023.1149590.


Friday, February 16, 2024

Summarizing my (with a good company) cryospheric work

The hydrological cycle is significantly influenced by the presence of water in its condensed states in middle and extreme latitudes. Various hydrological parameters change below 0  Celsius, such as water viscosity, thermal capacity, and hydraulic conductivity. Consequently, mainstream hydrology treatments that neglect freezing provide incorrect results in winter, high elevations, and the far north and south for most of the year. In the current state of global warming that threatens the cryosphere which is progressively disappearing, it is even more crucial to address its dynamics


A little of-of-date itinerary can be found in a previous post here. To understand our progress, three milestone theses summarize the work done

  • Matteo's  Together, we worked out the Thermodynamics of non equilibrium for ice-systems and the theory of freezing soils. Matteo implemented also an integrator in GEOtop, not the perfect one, but acceptable. Matteo's 2011 paper is a benchmark paper in the topic. 
  • Stefano's brought GEOtop to some maturity and especially fine tuned the various tools related to snow and ice. Stefano's 2014 paper remains a landmark in our work. 
  • Niccolò's  pushes forward the previous work. Especially remarkable is his work on re-implementing the informatics according to new (for us) concepts in OO programming and using (finally) safe algorithms for the integration of the equations. His WHETGEO and FreeThaw papers are a must read for completeness and clarity.

Our work's focus was primarily on the critical zone, where we modified the Darcy-Buckingham law to account for freezing and thawing and their related hydrological and mechanical effects. We primarily focused on the hydrological effects neglecting the mechanical ones but not neglecting the energy budget, a common practice in hydrology, which is obviously not possible. Consequently, we faced the necessity to simultaneously solve both the mass budget and the energy budget.

The formulation of the equations can be found in the theses and papers cited above, and you will realize that establishing a correct relation between the Darcy scale energy content and the corresponding water (liquid or solid) is the main challenge. Proper physics requires the consideration of interfaces between the phases: air-water-soil-ice. While a complete understanding of this relation has not been yet achieved, some working approximations have been obtained. Looking at the two compartments, snow and ice in the soil, they differ in many aspects, with snow lacking soil and being affected by its aerial origin. Both snow and ice in the soil have their own complexities, which affect their evolution. They often interact and the fate of the soil with or without snow is quite different.

While determining the correct equations would be satisfactory goal for many, it remains unresolved how to numerically estimate these equations. It turns out that these mildly nonlinear equations pose problems when solved using the usual algorithms based on variations of the Newton method. Convergence of the numerical methods is not guaranteed, and many workarounds have been deployed to overcome these difficulties, often leading to issues with mass and energy conservation principles.

Fortunately, Casulli and Zanolli (2010) found a method to address these challenges. The fundamental paper can be found in bibliography, and a progressive approach to its formulation can be obtained by reading Casulli's lecture notes and completed by watching videos in this blog. Ideally, attending Casulli's annual school in Trento, held every second half of January after our GEOframe winter school, would provide the best understanding. Tubini's recent papers are the result of this approach, and FreeThaw and WHETGEO are concrete implementations of these algorithms in Java/OMS/GEOframe.

Another aspect to consider is the implementation of these algorithms in informatics. Concepts related to this can be found in parts of Tubini's thesis and the related papers, especially Tubini and Rigon, 2022.

Finally, as a source of information, all of these models' open-source codes can be found on GitHub, both for the older GEOtop and the more recent GEOframe model components.

References


Casulli, Vincenzo, and ZANOLLI. 2010. “A Nested Newton-Type Algorithm for Finite Colume Methods Solving Richards’ Equation in Mixed Form.SIAM Journal of Scientific Computing 32 (4): 2225–73.

Dall’Amico, Matteo. 2010. “Coupled Water and Heat Transfer in Permafrost Modeling.” Edited by Riccardo Rigon and Stephan Gruber. Phd, University of Trento. http://eprints-phd.biblio.unitn.it/335/.

Dall’Amico, M., S. Endrizzi, S. Gruber, and R. Rigon. 2011. “A Robust and Energy-Conserving Model of Freezing Variably-Saturated Soil.” The Cryosphere. https://tc.copernicus.org/articles/5/469/2011/.

Endrizzi, Stefano. 2007. “Snow Cover Modelling at a Local and Distributed Scale over Complex Terrain.” Ph.D. Thesis, January, 1–189.

Endrizzi, S., S. Gruber, M. Dall’Amico, and R. Rigon. 2014. “GEOtop 2.0: Simulating the Combined Energy and Water Balance at and below the Land Surface Accounting for Soil Freezing, Snow Cover and Terrain Effects.” Geoscientific Model Development 7 (6): 2831–57. https://doi.org/10.5194/gmd-7-2831-2014.

Tubini, N. 2021, June. “Theoretical and Numerical Tools for Studying the Critical Zone from Plots to Catchments.” Edited by R. Rigon and S. Gruber. Ph.D., Dipartimento di Ingegneria Civile, Ambientale e Meccanica, Università di Trento.

Tubini, Niccolò, Stephan Gruber, and Riccardo Rigon. 2021. “A Method for Solving Heat Transfer with Phase Change in Ice or Soil That Allows for Large Time Steps While Guaranteeing Energy Conservation.” The Cryosphere 15 (6): 2541–68. https://doi.org/10.5194/tc-15-2541-2021.

Tubini, Niccolò, and Riccardo Rigon. 2022. “Implementing the Water, HEat and Transport Model in GEOframe (WHETGEO-1D v.1.0): Algorithms, Informatics, Design Patterns, Open Science Features, and 1D Deployment.” Geoscientific Model Development 15 (1): 75–104. https://doi.org/10.5194/gmd-15-75-2022.