Sunday, April 6, 2025

Using GEOET's Prospero model with minimal variations for simulating the non capacitive energy budget of snow and soil

This post is not self-explanatory and requires digging into other posts and some papers. 

Please review Section 2 of Concetta's paper (https://onlinelibrary.wiley.com/doi/10.1002/eco.70009?af=R) and verify the calculations presented there.

The model in D'Amato and Rigon (2025) uses a non-capacitive approach (it doesn't account for the thermal capacity of plants), and so will be if the same derivation is specialized for snow (or soil), which is a limitation. However, this approach is still more physically based than semi-empirical formulations or degree-day methods commonly used. In the literature, these are referred to as "stationary solutions" of the system. Despite the name, these solutions respond instantaneously to changing boundary conditions (radiation, latent and sensible heat fluxes), as evident in equation (10), which varies with radiation, wind velocity, and roughness.

Equation (10) and subsequent equations in Concetta's paper are essentially the Prospero solutions (though exact implementation should be verified in Concetta's code). The time interval of integration is, in principle instantaneous, but eventually you would like to integrate it over a finite time step (a hour, or a day, for instance). 

A non negligible aspect is that snow can melt into water and for any temperature you get from the energy budget, you need to partition the water in liquid water and ice. For this reason you probably need a partitioning function, like the one used for partitioning precipitation in rainfall and snowfall or you can simply use a melting law like  in simple models but now the temperature used should not be the air tempeature but the snow temperature.  See melting in simple models in  the links below

for further information. 

An important term is missing from the formulation: heat exchange by conduction with the ground, which should be represented as:

G = C_s T_Δg := C_s (T_g - T_s)

Where:

  • C_s is an appropriate exchange coefficient (can be taken as C_s = K/L, where K is the bulk thermal conductivity of the layer and L is its depth)
  • T_g is the ground temperature (which could be taken as the multi-annual air temperature average)
  • T_s is the snow temperature

Since this flux depends on the independent variable, it introduces additional terms that modify solution (10). Please derive these calculations independently.

Other terms that don't depend on the independent state variables can be included in the S_nk term. With these modifications, the Prospero code can effectively simulate the snow energy budget. Similar arguments apply to soil modeling.

A further consideration is the proper parameterization of the conductances C and C_E fluxes in equation (10), which differ from transpiration cases. For soil, according to Lehman-Or theory, evaporation should be modeled as potential until the water storage exceeds a threshold S_T, then decreasing proportionally with storage below this threshold when implementing an integrated model (Details ? I do not know).

I know that there are several missing aspects in this post. Who is interested, please ask. 

P.S. - These components have then to be carefully coupled to the other components. With respect to this, please consider the following: 

First review the presentation materials I've shared:

Getting new features to the linear systems (Vimeo2025)

The topic is that, based on my analysis, the second option is clearly the one should be applied in integrated distributed models (like GEOframe-NewAGE). However, this means we cannot simply subtract ET (or any other sink) from total rainfall - we need to incorporate this directly into the equation solver. While the example in the presentation uses a linear system with an analytical solution, the same principle applies to our non-linear fluxes where we use numerical integration. Therefore appropriate modifications could be necessary to the basic GEOframe-NewAGE codes. 

Thursday, April 3, 2025

A Ph.D. position on Advancing Physics-Informed Machine Learning for Environmental Modeling and Smart Irrigation Systems

Project Overview

This PhD grant, funded by Fondazione Bruno Kessler, aims to develop an advanced Physics-Informed Machine Learning (PIML) framework for modeling complex hydrodynamic and environmental systems. By integrating physical principles with data-driven methods, the research will focus on optimizing next-generation irrigation strategies. The project will harness the synergy between physical laws (as implemented in the GEOSPACE system) and machine learning to enable predictive, real-time, and scalable modeling tools for sustainable water resource management.

Key Objectives
  • Design hybrid PIML models that combine governing equations with data-driven predictive models (e.g., neural networks);
  • Improve predictive accuracy and generalizability across heterogeneous environmental conditions;
  • Incorporate real-time sensor data inputs to refine model states and parameters;
  • Benchmark PIML approaches against traditional numerical solvers and/or black-box machine learning models.

Methodological Approach

The core innovation of this project lies in the integration of physical constraints (as derived from GEOSPACE) into machine learning models. Building on recent advances in PIML, the goal is to design, develop, and validate models that enforce conservation laws and boundary conditions within neural architectures. This may involve:

  • Embedding partial differential equations (PDEs) directly into the loss functions of machine/deep learning models;
  • Developing adaptive training strategies to trade-oO data fidelity and physical consistency;
  • Utilizing sensor data to dynamically assimilate environmental variability into model predictions;
  • Leveraging high-performance computing to train and deploy models at scale across complex
  • domains.


Expected Outcomes

This integration will enable:
  • Accurate and efficient modeling of water distribution and use in precision irrigation systems;
  • Real-time monitoring and decision-support capabilities for agricultural and environmental applications;
  • Enhanced data efficiency and model robustness through physics-based regularization;
  • Improved understanding of system dynamics under data-scarce and/or non-stationary conditions.

Implementation Timeline

Months 1–6: Literature review on PIML methodologies;
Months 7–18: Design and develop a core PIML architecture, integrating IoT data sources;
Months 19–24: Validate models using lab-scale and field experimental datasets;
Months 25–36: Upscale models to real-world irrigation systems deployed within ongoing local and EU
projects (e.g., IRRITRE, AGRIF .OODTEF), and quantitatively assess their impact on water-saving strategies.

Possible Collaborations

Fabio Antonelli, Fondazione Bruno Kessler; Sara Bonetti and Concetta D'Amato, EPFL

Info: abouthydrology <at>  gmail.com

Sunday, March 30, 2025

A Ph.D. position ! Advanced Soil Biota-Hydraulics Interface for the WHETGEO-GEOSPACE system

Project Overview

This subproject, funded under the ICOSHELL project, aims to develop an integrated modeling
system that explicitly accounts for the dynamic interactions between soil biota activity and soil
hydraulic properties. Building upon the WHETGEO-1D and 2D frameworks, we will implement a
novel coupling between soil fauna population dynamics and plants root growth, evolving soil
hydraulic characteristics. The modelling system implemented will be eventually used for studying
the feedback between soil-vegetation hydrology.

Key Objectives


  • Extend the WHETGEO model architecture to incorporate time-varying soil hydraulic properties influenced by soil biota 
  • Implement the Kosugi soil water retention curve model with parameters that dynamically evolve based on biological activity
  • Develop and integrate a population dynamics module for key soil engineers (earthworms, ants, termites)
  • Create a comprehensive validation framework using laboratory and field experimental
  • data
Figure from Enrico Chiesa Master Thesis


Methodological Approach

The core innovation of this subproject is the implementation of a feedback loop between
 biological activity and soil physics. Following Meurer et al. (2020), we will start to model how earthworm populations modify soil structure, but significantly expand this approach by:

  • Replacing the van Genuchten model with the Kosugi water retention curve formulation, which provides a more direct physical interpretation of pore size distribution
  • Developing a differential equation system where the Kosugi parameters (median pore size and standard deviation) are directly modified by biological activity
  • Implementing these dynamics within the robust NCZ algorithm of WHETGEO, ensuring numerical stability across diverse conditions
  • The population dynamics will be modeled as a set of ordinary differential equations representing different functional groups of soil engineers, their reproduction, mortality, and activity rates as functions of environmental conditions (temperature, moisture, organic matter)

Expected Outcomes


This integration will allow to better capture:

• The temporal evolution of soil infiltration capacity following land-use changes

• The self-reinforcing positive feedback loops of ecosystem restoration, where initial

vegetation changes trigger soil biological activity that further enhances water retention

• The resilience of soil hydrological function under climate change scenarios


Implementation Timeline

Months 1-6: Preliminary studies, doctoral school activities

Months: 6-18 Implement Kosugi model in WHETGEO framework, develop and integrate

population dynamics module Months 18-24: Validate against experimental data Months 32-36:

Upscale to field applications and integration to estimate catchment scale effects. Study effects of

soil management

Possible collaborations

EPFL Lausannne, Prof. Sara Bonetti and Dr. Concetta D'Amato

Info: abouthydrology <at>  gmail.com

Wednesday, March 12, 2025

The Marvelous Physics of Plants: a personal Introduction

 "The Marvelous Physics of Plants" presents an exploration of the physics behind how plants function, particularly focusing on water transport mechanisms. The presentation begins with poetic descriptions of plant processes, then explores Erwin Schrödinger's fundamental question about how physics and chemistry can explain the events within living organisms. The authors examine various physics domains relevant to plants: quantum physics, thermodynamics, hydraulics, micrometeorology, stability, and light.

Among the other things,  the authors examine the physical limits of tree height, discussing how hydraulic restrictions ultimately limit how tall trees can grow. They also demonstrate synthetic tree models that scientists have created to replicate these natural mechanisms.

The slides combine mathematical formulations, anatomical diagrams, and experimental results to illustrate the physical principles governing plant function. A video of the talk is also available.


Friday, February 28, 2025

Three Batchelor Graduation Works

The first  Thesis-poster,  by Agnese Cavazzini, supervised by Gaia Roati and me, presents a hydrological study of the Secchia River basin using the GEOframe-NewAGE system. The research analyzes water balance and simulates river flow while generating soil moisture maps to identify drought-prone areas. Key elements include watershed division into sub-basins, mass balance equations, and calibration against measured data. Results show flow simulations at two monitoring stations and soil moisture anomaly maps. The successful implementation provides valuable insights into the basin's hydrological dynamics across Modena, Reggio Emilia, and Mantova provinces. You can get a high resolution poster by clicking on the Figure below..


The second Thesis-poster, by Lorenzo Dalsasso,  presents a statistical analysis of ground precipitation patterns by Lorenzo Dalsasso. Using hourly precipitation data from three weather stations, the study evaluates which probability distributions best represent precipitation duration, intensity, and intervals between events. A Python notebook with Kolmogorov-Smirnov tests determined that lognormal distributions best fit precipitation durations, Weibull distributions best represent precipitation intensities, and either Weibull (stations ID 40 and 1100) or lognormal (station ID 263) best characterize intervals between precipitation events. The results include detailed statistical parameters for each station. The high resolution poster can be found by clicking on the Figure below. 




Thr third thesis-poster presents Marco Feltrin's study on evapotranspirative fluxes in grapevines by integrating the GEOSPACE ecohydrological model with WiseConn sensor technology (dr. Marco Bezzi). The research compares two rainfall scenarios: a wet scenario (1147 mm of total precipitation) and a dry scenario (682.4 mm of total precipitation) to evaluate plant water stress. Using the one-dimensional GEOSPACE model with data from a vineyard near Verona, results show that the dry scenario led to half the plant transpiration during summer months. The model effectively demonstrates how water content throughout the soil column affects water stress in plants, with practical applications for irrigation management, water conservation, and predicting water availability for viticulture under changing climate conditions.  The high resoltion poster can be found by clicking on the Figure below.