In this post I am presenting another R exercise. Not that I believe very much on the physical content of the what Scottischsnow does here. (Actually this bold way of treating data and empirical formulas at large scale time and spatial is something that I usually blame).
However the information gathered in this way could be useful as a reference step for doing physically based things, at a more detailed scale (for instance by using our JGrass-NewAGE).
Moreover I like it because it uses R all the way for extracting the data and visualising them, and this can be equally likely instructive.
My reflections and notes about hydrology and being a hydrologist in academia. The daily evolution of my work. Especially for my students, but also for anyone with the patience to read them.
Wednesday, July 24, 2013
Monday, July 22, 2013
Rainfall statistics for the Netherlands
One of my most successful postings is the one of R resources for hydrologists. Sometimes contributions arrive from non -hydrologist. This is the case that of two posts from Wiekvoet whose authors download and analyse rainfall data from the Netherlands Royal Institute. The first post mainly performs an exploratory analysis of the data. The second one looks for rainfall changes in the past 100 years. Enjoy the reading.
As a byproduct one learn also how to use ggplot2 through some examples.
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
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).
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
- Short Intro to UML* class diagram (YouTube 2018)
- UML - stands fro Universal Modelling Language. It is a set of graphic tools that can be used to sketch OO programs. The classical reference is Martin Fowler's, UML 2.0 distilled, 2003. It is a short and nice and time-affordable reading.
- Naming Conventions (YouTube2018)
- The Structure of a Java Class and something about Javadoc (YouTube2018)
- Reading data from the console
- Storyboard of the classes
- Reading data from a File (no frills, no much design) (YouTube Live Video)
- Writing data to a File (no frills, no much design) (YouTube Live Video)
- Reading and writing data to a File. Designing the classes (YouTube Live Video)
- Reading and writing data to a File with Java 8. (YouTube Live Video)
- For getting more on Streams see the chapter 3 of Modern Java in Action (a book worth to buy).
- For Getting more on Lambdas, Chapter 6 to 9 of Modern Java in Action .
- Learning a little better what is a regular expression
- Reading and writing data to a File with Java 8 II and III (YouTube video)
- From command line
- from Eclipse (EGit)
-----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:
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
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.
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.
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.
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.
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
Whitaker, S., The Forcheimer equation: A theoretical development, Transport in Porous Media, October 1996, Volume 25, Issue 1, pp 27-61
Liu, Z, and Todini, E., Towards a comprehensive physically-based rainfall-runoff model,Hydrology and Earth System Sciences, 6(5), 859–881 (2002T)
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.006O'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
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,
She already wrote to some applicant:
" Dear ***,
" 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. "
Subscribe to:
Posts (Atom)






