Monday, July 16, 2018

Towards a new SWMM, JSWMM

You want to design a storm water management system. What you usually get is a rainfall-runoff model. In this specific subfield, the model is SWMM. EPA SWMM contains many features that were implemented to simulate urban storm water depletion network. Meaning that there are specific model's parameters set for that, and that a community gathered around this tool.

However SWMM is not a system for designing sewers. Designing requires that you repeat the modeling actions several times. At the outlet of any pipe, you have to:
  • estimate the runoff under a "design rainfall" coming from some intensity-duration-frequency curves.
  • get the maximum discharge with an assigned return period (say 10 years)
  • Use simplified hydraulics for obtaining the size of the pipe apt to contain the maximum discharge
  • repeat the operation for the pipes downhill, without leaving out uphill branches.
Operation 2 above requires a search algorithm to find the rainfall duration that is responsible for the maximum discharge. The complete theory is in Rigon et al., 2011.
The point is that SWMM does not do the sequence of operation above. This is one of the reasons we implemented JSWMM. To see what it does, click on the figure.

Monday, June 11, 2018

On complex networks computation of mountain catchments

From June 12 to June 14 2018 in Trento there is the 5th European IAHR conference. I am convening (and also presenting a contribution) which derives from our modelling within the GEOframe system.
The scope of the presentation is to inform about functionalities of GEOframe and something of what it is hidden under the hood. Clicking on the figure above, please access the pdf of the presentation.

Wednesday, June 6, 2018

Snow related Ph.D. position (the Stradivari snow project)

The Stradivari project research aims to build better tools for analysing the processes of the hydrological cycle. The project is more focused on the tools (building the violin) but it does not forget the music that has to be played. It is conceived to account for hydrological processes interactions and feedbacks, and develop new mathematics (equations) for their description. Use of appropriate contemporary numerics is also part of the project. The overall project builds on the foundations given by the GEOtop ( model and the GEOframe-NewAGE ( infrastructure. 

It is time to move on the generation of snow models. The computation abilities and the physics of snow are now much better known than thirty years ago. However, most of the snow model are based on parameterisations which should be updated. We start from the experience made in GEOtop (versions 0.*,1.* and 2.*) which is capable of reliable snow height, temperature and densities, at operational level over all the Alps, but we look also to the experiences made by CROCUS, SnowPack and Alpine 3D. In particular GEOtop use highlighted various issues that we plan to overcome with a new version of the model, which requires both deepening the thermodynamics of snow and its numerical implementation.

In synthesis we identify the following aspects to be improved:
  • Actual GEOtop snow model is 1-d. It can improved including the vapor phase movements explicitly. Besides, we can add modules to have a better account for: density, viscosity and, adding the vapor phase, type of grain. 
  • The physics. Analysing the thermodynamics to keep out of the future formulation empirical parameterisations of the processes which revealed obsolete; 
  • Adding Richards equation for water percolation; 
  • Analysing the separation rainfall, snowfall through a more physical modelling than actual; 
  • the description of deposition on canopies, and subsequent effects of vegetation on sublimation; the effect of slope and topography characteristics on snow deposition and sublimation 
  • blowing snow in complex terrain 
  • assimilation of hydro-meteorological data and calibration 
  • integration with remote sensing data 
This research will be implemented in strict collaboration with MobyGIS, EURAC Research (Giacomo Bertoldi, Ph.D), Stephan Gruber,  Professor at Carleton University and with Niccolò Tubini, the Ph.D. student who is actually developing new theory codes for Richards equation and freezing soil. Collaborative and unselfish spirit is required in this research group.
The project has also some practical outcomes that are related to:
  • the avalanche triggering 
  • the water availability due to snow and glaciers melting (both in the short and long terms) 
  • the hydroelectric production 
All the code developed will be done in Github (or similar platform), inside the GEOframe community and will be Open Source according to the GPL v3 license.

The candidate will take care of implementing, besides the code, the appropriate procedures for continuous integration of the evolving source code, and s/he will be also asked to maintain a regular rate of commits to the common open platform. Despite these conditions, and being free and open source, the code will be intellectual property by the coder. This will be guaranteed also by the components-based infrastructure offered by OMS3, which allows to better define the contributions of anyone.
The implementation part will be followed, accompanied by testing activities, either for mathematical consistency, than for physical consistency with experiments and field measurements.
The Ph.D. student is intended to produce, besides working and tested codes, also at least three papers in major journals (VQR Class A), of which, at least one as first Author. Duration of the doctoral studies could be three or four years.

This project can enter either the curriculum C (Environmental Engineering) or the curriculum A (Modelling and Simulation) of our doctoral school.

Further information of the policies of the research group can be found:

Evaporation and Transpiration

This contains the video related to my interpretation of evapotranspiration.

The thermodynamics of evaporation

Transport of vapor in the atmospheric boundary layer
Evaporation from soils


The energy and mass budget 

The Penman-Monteith approach

Tuesday, June 5, 2018

From where do waters arrive

You cannot have an aqueduct if you do not have water to resupply it. So water comes from springs, wells (i.e. from groundwaters) or intakes (surface waters).


Tuesday, May 29, 2018

Do not do statistics if you do not have casual effects in mind

Statistics was believed after the master of the last century to be the science of correlation, not of causation. However it is clear to our contemporary researchers, at least some of them, that interpreting data without any guess about causation can bring to wrong conclusion. Here below, please find an example from: "The book of why: the new science of cause and effect" by Judea Pearl and Dana MacKenzie.
You should first look at the right figure. The scatterplot presents a roughly linear relation between Exercise and Cholesterol in blood. First observation, set this way, we probably have to reverse the axes. In a causal interpretation, it appears that exercise cannot cause cholesterol. On the contrary the cholesterol presence impose to the subjects to exercise more. Or there is something strange in data. More exercise cannot cause, by our normal belief, more cholesterol.
However, this is not actually even the main point. What the right figure suggests is that there is a positive correlation between the two variables: more cholesterol implies more exercise. However, as the left figure reveals, the real situation is not quite true. Because a cause of cholesterol is age, it appears that is reasonable to consider also this variable in the analysis. Then, when we separate the data among ages sets, we can see a further structure in the data and, in each class of age, in fact, the correlation between exercise and cholesterol is reversed. The less you exercise, the higher is your cholesterol. At the same time, the younger you are the less cholesterol you are expected to have in your blood. Now the picture is coherent with our causal expectations. I think there is something to learn.

Wednesday, May 23, 2018