Tuesday, July 16, 2013

Java for Hydrologists 101

There are a few postings on Java in this blog. Since I want to teach it to my students, I am quietly starting to populate this page with presentations which, eventually, will constitute the core of an informal class (;-)) the Java for Hydrologsts 101.
The focus of the JfH-101 is not simply to gain Java knowledge from the scratch (or so), but to address those topics and issues that have to do with my experience of hydrologist. So, I will try to cover Java as well as OMS3, and at the right time Geotools, and jgrasstools, not forgetting the tools of the tool (Ant, Maven, Git)^*, but in the meanwhile I will try to address the topics related to object oriented programming.  Programming is actually very much not talking about  that but doing it, so many of the slides will actually recall to actions.

Topics

0 - Getting Started (mostly things to read -or start to read- before the start) (YouTube 2018 video)
1 - Your first program (You Tube Video 2018)
2 - Solving a linear equation

3 -  A few diversions
4 - Reading  data from the system's console
5 - Reading data from a File

6 - Working with Git
-----Not yet implemented: ----

7 - Programming the heat diffusion equation
8 - Making the heat diffusion an OMS3 components
9 - Building Java projects with Ant, Maven and GRADLE
8 - GEOtools essentials
9 - Commenting the programming of the GEOframe-ET
10 - The Java REPL
11 - A little on Java Modules in Java 9
12 - Setting the continuous integration in GIT (using Travis)

The source code is available for download to from GitHub.

References

Please go to this blogpost.


^* - From the links you can quite understand the I rely very much on Lars Vogel site for the basic stuff. It is not obviously the only good resource available (stackoverflow is another one, for instance, and many others will be addressed).

Tuesday, July 9, 2013

Hillslope hydrology from the point of view of Richards equation

This is the lecture I gave at the second Summer school on Water Resources that I co-organised. The lecture was recorded in video and should be made available soon, together with the other lectures.
The lecture can be considered an expansion of the talk I gave at Ezio Todini's symposium with more details and, hopefully, more clear explanations. The presentation is, as usual, available on slideshare, and you can get it clicking on the image below.

The presentation mainly use a paper by Cordano and Rigon, 2008 to obtain various degree of simplification of the Richards equation. The scope (especially for what regards the Todini's symposium version) was to show that many approaches currently used derive, in fact, from a simplification of Richards, and there are not very much reasons to shoot to it as unphysical, at least if this is subsequently followed by the use of one of its simplifications.
However, the slides cover also a discussion over the time scales of flows in hillslopes, and the relative timing of vertical and lateral infiltration.

Here it is a short abstract:


The presentation covers Richards equation as applied to Hillslope Hydrology from its foundations. It is said that it assumes mass conservation and  the existence of the Darcy scale at which the soil medium can be treated as a continuum. Then it is specialised using some well known parameterisations (van-Genuchten Mualem), and subsequently is simplified to obtain other equations. In order: the 1-D Richards equation, the Boussinesq equation, the hillslope-storage Boussinesq equation, and finally its stationary approximations.  All of these are used in literature for various purposes, including soil moisture distribution and hillslope stability. The simplifications are based on the assumption that lateral (slope-parallel) flow is slower than slope-normal flow, which is subsequently shown not being necessarily true, true some simulation with a 3-D Richards equation solver. This eventuality is caused by hydraulic conductivity being (in some soils) high variable with water content. Eventually, a conceptual model is built on the knowledge acquired, in order to reduce the computational burden. Lastly some cases are discussed where Richards equation could fail using data from the Panola hillslope. It is shown that the fill-and-spill phenomena can be described properly, and that, on the contrary, the presence of macropores cannot.

The various schemes of simplification have a great effect on the identification of landslides' locations, and, in fact, many of the papers cited (and provided in the blog) are dealing with landslides hazard assessment. 

References cited in the presentation

Beven., K. J., M.J.. Kirkby, A physically based, variable contributing area model of basin hydrology, Hydrological Sciences bulletin-des Sciences Hydrologiques, 24, 1-3, 1979

Buckingham, E. 1907. Studies on the movement of soil moisture. Bulletin 38. USDA Bureau of Soils, Washington, DC.

Casulli, V,  Stelling GS, Semi‐implicit subgrid modelling of three‐dimensional free‐surface flows
International Journal for Numerical Methods in Fluids 67 (4), 441-449

Cordano, E., & Rigon, R. (2007). A perturbative view on the subsurface water pressure response at hillslope scale. Water Resources Research, 1–36.

Cordano, E., & Rigon, R. (2010). A mass-conservative method for the integration of two-dimensional groundwater (Boussinesq) equation. Water Resources Research, 1–24.
Dietrich, W. E., 1989, Slope morphology and erosion processes, in C. Wahrhaftig and D. Sloan (Eds.), Geology of San Francisco and Vicinity, Field Trip Guidebook T105, American Geophysical Union, p. 38-40.

D'Odorico, P., Fagherazzi, S., & Rigon, R. (2005). Potential for landsliding: Dependence on hyetograph characteristics. Journal of Geophysical Research, 110(F1), 1–10. doi:10.1029/2004JF000127

Iverson, R. M., Landslide triggering by rain infiltration, Water Resour. Res., Vol. 36, N0. 7,  1897-1910, 2000

 Lanni, C.; McDonnell, J. J.; Rigon, R., On the relative role of upslope anddownslope topography for describing water flow path and storage dynamics:a theoretical analysis, Hydrological Processes Volume: 25 Issue: 25 Pages: 3909-3923, DEC 15 2011, DOI: 10.1002/hyp.8263

 Lanni C., J. McDonnell JJ, Hopp L., Rigon R., "Simulated effect of soil depthand bedrock topography on near-surface hydrologic response and slope stability" in EARTH SURFACE PROCESSES AND LANDFORMS, v. 2012, (In press). - URL: http://onlinelibrary.wiley.com/doi/10.1002/esp.3267/abstract . - DOI: 10.1002/esp.3267

Lanni C., Borga M., Rigon R., and Tarolli P., Modelling catchment-scale shallowlandslide occurrence by means of a subsurface flow path connectivity index, Hydrol. Earth Syst. Sci. Discuss., 9, 4101-4134, 2013

Montgomery, D. R., and W. E. Dietrich, Where do channels begin?, 1988, Nature, v. 336, p. 232-234.

Mualem, Y., A New Model for predicting the hydraulic conductivity of unsaturated porous media, Water Resour. Res., vol 12, No 3, 1976

Narsilio, G. A., Buzzi, O., Fityus, S., Yun, T. S., & Smith, D. W. (2009). Upscaling Navier-Stokes Equations in porous media: Theoretical, numerical and experimental approach, Computers and Geotechnics, 36(7), 1200–1206. doi:10.1016/j.compgeo.2009.05.006

O'Loughlin, E.M., Prediction of Surface Saturation Zones in Natural Catchments by Topographic analysis, Water resour. Res., vol 22, no 5, 794-804, 1986

Orlandini, S., G. Moretti, M. Franchini, B. Aldighieri, and B. Testa (2003), Path-based methods for the determination of nondispersive drainage directions in grid-based digital elevation models, Water Resour. Res., 39(6), 1144, doi: 10.1029/2002WR001639.

Orlandini, S., P. Tarolli, G. Moretti, and G. Dalla Fontana (2011), On the prediction of channel heads in a complex alpine terrain using gridded elevation data, Water Resour. Res., 47(2), W02538, doi: 10.1029/2010WR009648.

Richards, L.A., Capillary conduction of liquids through porous mediums, Physics 1: 318-333, 1931

Troch P.A., Paniconi, C., van Loon E.E, Hillslope-storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes: 1. Formulation and characteristics response, Water Resour. Res., Vol 39, No 11, 1316, doi:10.1029/2002WR001728, 2003


Tromp-Van Meerveld, H. J., & Mcdonnell, J. J. (2006). Threshold relations in subsurface stormflow: 2. The fill and spill hypothesis. Water Resources Research, 42(2), W02411. doi:10.1029/2004WR003800

M. Th. van Genuchten, A Closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sc. Soc. of America, vol 44, no. 5, 1980

Whitaker, S., The Forcheimer equation: A theoretical development, Transport in Porous Media, October 1996, Volume 25, Issue 1, pp 27-61

Monday, July 1, 2013

2012 Water Resources Research Editors Award

Directly from the main page of Water Resources Research:

"With the approval of the AGU Hydrology Section Executive Committee, WRR has instituted an Editors’ Choice Award to be given to (about) the top one percent of published articles in any calendar year. Our goal is to provide professional recognition to scientists and students for their outstanding work. The selection is made by the WRR Editors based on technical significance, novelty, originality, presentation, and broader implications of the publication. Awards made in a given year are for articles published in the previous calendar year. We are pleased to announce that the following contributions are being recognized for 2012 award year:"



Grafton, R. Q., M. B. Ward, H. To, and T. Kompas (2011), Determinants of residential water consumption: Evidence and analysis from a 10-country household survey, Water Resour. Res.47, W08537, doi:10.1029/2010WR009685.

Sun, F., M. L. Roderick, W. H. Lim, and G. D. Farquhar (2011), Hydroclimatic projections for the Murray-Darling Basin based on an ensemble derived from Intergovernmental Panel on Climate Change AR4 climate models, Water Resour. Res.47, W00G02, doi:10.1029/2010WR009829.

Wigmosta, M. S., A. M. Coleman, R. J. Skaggs, M. H. Huesemann, and L. J. Lane (2011), National microalgae biofuel production potential and resource demand, Water Resour. Res.47, W00H04, doi:10.1029/2010WR009966, [printed 48(3), 2012].

Amiaz, Y., S. Sorek, Y. Enzel, and O. Dahan (2011), Solute transport in the vadose zone and groundwater during flash floods, Water Resour. Res.47, W10513, doi:10.1029/2011WR010747.

Gallant, J. C. and M. F. Hutchinson (2011), A differential equation for specific catchment area, Water Resour. Res.47, W05535, doi:10.1029/2009WR008540.


Huss, M. (2011), Present and future contribution of glacier storage change to runoff from macroscale drainage basins in Europe,Water Resour. Res.47, W07511, doi:10.1029/2010WR010299

Friday, June 14, 2013

A new mailing list for hydrological announcements

Dear All,

I opened a google group, a mailing list to communicate among hydrologists. This will act the very same way the Gilbert Club list does for geomorphologists in being used to send announces about events, Ph.D. and post-doc positions, new books, open source software, AGU and EGU sessions and so on.  



The list is:


abouthydrology@googlegroups.com

and below you can find instructions for subscribing to it (and eventually unsubscribe)

https://groups.google.com/forum/#!forum/abouthydrology

There you should see this window:


Click on the blue button "Apply to Join the group" to send an application.

I though it was more useful to have such a list than maintaining, each of us, a list of colleagues to send this kind of messages. I started the list and manage it, but I am looking for volunteers to manage it with me, and maybe without me, in the long run.

Friday, June 7, 2013

Two possible Ph.D. positions available

Dear reader,

in the new call of our doctoral school in Trento I proposed two topics which could be of interest for some of your friends, students or people you know. Here they are:





Distributed Modelling of the Hydrological Cycle at large scales (including human impacts)
This regards mainly the development of the  JGrass-NewAGE for the complete closure of the hydrological budgets, in medium to large scale modelling. This requires the implementation and testing of new physical-statistical model of the various terms of the hydrological cycle, and their application to case studies at the scale of hundreds to thousands of square kilometers. At the  moment the model has a first implementation of all the processes that is going to be thoroughly  tested, and the main interest in this research is to go beyond the simple forecasting of hydrological quantities (in space-time) to achieve  the estimation of error bounds in the predictions with the application of appropriate calibration methods, and data assimilation procedures.  The doctoral work is intended to achieve also the application of models and tools to a real case (which could be the river Adige) and will be pursued in coordination with ISTAT, the Italian Ministry of Environment, and the appropriate river Authorities.
The candidates are required to have good programming skills, especially in Java and R (or the will to pursue them), to improve the numerical modeling within the framework of the JGrass-NewAGE system. It is intended that all the tangible work in programming tools is produced as free software, and using free software.  Previous outcomes of the research line can be found at: 
Distributed Modelling and Data Assimilation of the Cryospheric processes
This study involves the modelling and forecasting of the evolution of the snow cover working with the model GEOtop.  Previous Ph.D researchers implemented a one dimensional energy budget of both the snowpack and  freezing soil. They also posed the bases for  further theoretical and numerical improvements of the model, to a 3D version, and eventually including also different constitutive relations, which could be pursued in this research.
The present proposal is especially dedicated to include (or embed) GEOtop modelling with a data assimilation system dedicated to real-time forecasting of the snow cover, depth, and status. The Ph.D. work could be oriented to assimilate either ground data than remote sensing data. 
The work will be made in coordination with Mountain-eering S.r.l, a spin-off of the University of Trento, and Stephan Gruber of University of Carleton (CA). The candidates are required to have good programming skills in Java/C++ (or the will to pursue them) to improve the numerical modeling within GEOtop (C++), and the system that is going to be built around it (Java). It is intended that all the tangible work in programming tools is produced as free software, and using free software.
Information about GEOtop can be found at:
Information about my research on Cryospheric processes at:

As general attitude in my research I believe that research must be reproducible:
and I require the same discipline to my Ph.D. students and collaborators. 
All my research topics, besides these two above are available at:

Programming skills and attitudes (at least to acquire programming skills) are mandatory. Knowledge of commercial softwares not necessary (and its use discouraged along the Ph.D. since, for the core research, I use softwares evolved by my group and my former students for almost everything). 


All the best,

riccardo rigon


Further information can be obtained by writing to claudia.fraizingher (at) unitn.it
She already wrote to some applicant:



" Dear ***,

with regard to your e-mail below, we inform you that the applications for our Doctorate programme are managed only via the online system. Please check our website: http://www.unitn.it/en/dricam/29815/current-call. The call for the 29th cycle of the Doctoral School in Civil, Environmental and Mechanical Engineering for the academic year 2013/2014 is open, application deadline is September 30, 2013. All the details relating to the application procedure are specified in the call (read it carefully). The list of research subjects proposed for this edition will be available in the next few days.

I also inform you that the assessment of qualifications and the oral examinations for candidates living abroad who wish to be interviewed in the country where they are currently living, can be arranged with a videoconference, subject to approval by the Head of the School. In this case you will have to complete the online section of the application “Offsite examination“ indicating the reason for such request (work, study, permanent residency, etc.), a Skype account and the address and phone number of a university or diplomatic representation as well as the name of a third person who will confirm your identity at the beginning of the examination. In any case you find all the information on-line. "

Wednesday, June 5, 2013

Ezio Todini 70th Symposium: my talk

Please you will find below the picture (click on it!) the presentation I will give the 6th of June in Bologna to mark the 70th birthday of Prof. Ezio Todini. I will say things that I already said in different contexts. However, reading them in this form can help to understand them more.
Ezio did a lot of interesting works on rainfall runoff modelling, starting from the Arno model, that can be considered a milestone in the field and Topkapi, his distributed hydrological model. He contributed also to data assimilation, models calibration, time series analysis, and I will dedicate a post in a few day to him and to his most important papers.

Below you can find also the link to our posters about:
The topics of the first poster where already covered in various posts. The second poster covers the same topics of the paper for the Hydrology days conference.

Monday, June 3, 2013

We need water

This time just a photo by the master Sebastiao Salgado. A photo that says more than thousand words.