Monday, December 31, 2012

Open Java resources for Hydrologists

Next question about Java, after having introduced it, and talked about options for doing numerics, is to see if there are effective resources on which hydrologists can build something useful, incrementally, and as a part of a community.
Also OpenMI uses Java, but, at the moment, Java there is just the second choice, well behind their C#. A second framework which we are looking to is:
  • Openda. It is a framework for data assimilation and calibration of models promoted by Deltares. 

There are out there probably other Java frameworks, and I will add them to the list, whenever they will be brought to my attention.
The field of more intense Java activity is actually the GIS one. There are at least three/four major open source efforts in Java:
In fact, it is not very clear to me, but also the
initiative seems having parts based on Java.
They are based on various type of resources among those in Java:
Both GeoTools and Sextante use 
  • JTS, the Java topology suite for several elementary operations.
An interesting,  old but not outdated initiative to give a third dimension to data visualisation is 
  • Visad, a library to visualise scientific datasets
from which sprout
  • the IDV  virtual globes family to visualize 3D geophysical data
IDV is not the unique Java Virtual globe. Another one is:
  • the Nasa WorldWind which gives (as open source) what Google Earth gives to the general public. Various prototype implementations of Nasa WorldWind in uDig (and in OMS3) were already made, and possibly sooner or later we will have a full mainstream implementation of it.
Last but not least an important resources is 
  • the NetCDF data format (of which IDV, but also ncBrowse,  is a viewer) has also an implementation in Java.
Even if my post is more concerned about the development of Hydrological resources and models, the above frameworks/programs/applications/models or part of them can just be used (and just not developed) with profit.  Starting from OMS3 there are at least three major modelling efforts that use it:
Not to forget is the whole Jgrasstools available for DEM manipulation (and further modelling) inside (and outside) uDig. All of them have more resources in a post dedicated to OMS 3 resources.
The Sextante library itself comes with more than 300 tools for manipulation of DEM and images. Also IDV and NasaWorldWind can be just used for the purposes of any research.
Another full-fledged Java Catchment model is

by Ricardo Mantilla. Its source code can be obtained following the instructions in the user manual.

A recent addition to the resources is


Why use Java, instead of  R ?

In one of my previous post, I talked about R software, and in fact, my group of people use both.  R is much more for some higher level operations, and is much less customisable than the resources developed in Java, and,  usually, less efficient.  My former students, using Java capabilities were able to built professional/industrial applications on the basis of Java (the GIS are a proof), that using just R would have been impossible.
However, the new incoming version of OMS3, should include the way to call R from the OMS console. Please contact the OMS version 3 developers for more information.

Why do not use Python ?

Python (or many other languages) could have been a good choice. Never say never, however I chose Java and to be successful one has to consistently invest his/her own limited resources in one direction.

There exist material/documentation on the above material ?

One introduction to programming GIS and essentially all I would require to a collaborator to know  can be found here by Andrea Antonello

Saturday, December 22, 2012

A Java Library (of books) for beginners, and a little path to self instruction for geophysicists


It is quite a few years that I insist for students moving to program in Java. Java is not anymore Java 1, it has evolved towards a mature language and now it is at its version 7 (Java 8, now and Java 9 soon) Java  13. Times ago I talked about doing matrix algebra in Java. Most of that post is still valid.

However, many would realize that for learning Java the right way s/he has to start from the beginning, and possibly, start to learn also the concepts of object oriented programming (the real reason to use Java).

My personal training began with the videos (in Italian) of one of my Trento university colleague Marco Ronchetti. For others, MIT Courseware could be a good starting, but also these less academical clips can help.

Programming is more than just knowing a language. This concept is well covered in

Introduction to programming using Java by David J.Eck
Java Precisely by P. Sestoft (suggested by Joshua Bloch, see below)

and nothing can substitute the reading of a good book. More specialised classics on Java are also:

Thinking in Java 4th edition

Thinking in Java annotated solutions book

the two by Bruce Eckel

The design pattern Java companion, by J. W. Cooper

These books sum up to a lot of pages to read and working them out is not a thing you do in a day.
People who teach Java, usually rely on command line. However, I would prefer to learn how to use an IDE. Our choice is, since many years, Eclipse. Learning to use it is easy but is made easier by reading Lars Vogel tutorials and material.

Obviously Oracle itself has a series of very good tutorials.

Going to more specialised book, for scientists, I also recommend

Mak, R. (2010). Prentice Hall - Java Number Cruncher. The Java Programmer's Guide to Numerical Computing (pp. 1–493).

Besset, D.H., Object-Oriented Implementation of Numerical Methods: An Introduction with Java & Smalltalk , 2000. (~70 Mb) (His code can be found here).

Nikishkov, G. (2010). Programming Finite Elements in Java (pp. 1–394). The preface of his book can give a good rational for using Java instead of other languages.

A general collection of Algorithms in Java is Algorithms by Kevin Wayne and Robert Sedgewick (also addressed here).

You can also give a look to old C Numerical Recipes (for algorithms) and use this free Java project for their application. (Do not support NR license policies).

Having got a clear idea of what Java is, you probably wants to write clean code. The right book is Clean Code by Robert Martin.

Now, you are probably ready to understand the subtleties of the Java Matrix Libraries I indicated in my older post.

If one wants to move to geosciences and GIS, the best choice is to give a look to the GEOtools site and to al the documentation s/he can find there. At that point s/he will be probably also ready for programming the Jgrasstools.  Programming OMS3 components would be easier, indeed, if you do not care of geographic features.  For GIS and working with us, one should consider the introduction by Andrea Antonello.

Eventually you can probably go back and make reacher your library (of books) by also reading and applying

Martin Fowler's UML Distilled: A Brief Guide to the Standard Object Modeling Language

Erich Gamma, Richard Helm, Ralph Johnson and John VlissidesDesign Patterns: Elements of Reusable Object-Oriented Software (but the easiest book, from the opinion of many is Head First Design Patterns by Eric Freeman and Elisabeth Robson.

and  the book that many consider the best one:

Effective Java Second Edition, by Joshua Bloch  by many considered the best book about Java.
Here it is the third Edition ;-).
Here there is also a video about.

For mastering Java generics and collections, I suggest:

Java Generics and collections by Maurcie Naftalin & Philip Wadler, by O'Reilly.

Obviously Java is one of the most used languages, and there is a lot of material on the web. Starting from the Oracles itself main site.

However, programming is about writing code, not just reading about it. Wishing I would have time to do it! 

Friday, December 21, 2012

The Boussinesq equation paper

Finally, after a long story (see the submission post) the paper on Boussinesq equation integration was accepted by Water Resources Research.  In the revision phase it was asked to add further comparison of the model with a model solving Richards equations.  We chose, for our comparison, and, a-posteriori , the choice was quite obvious, the model GEOtop.  Emanuele Cordano and I believe that the new added section could be quite interesting for understanding some characteristics of infiltration at hillslope scale.

To read the preprint, please follows the link. The code and documentation is available from the Boussinesq site. However, a more recent version, in Java and OMS compliant is also available at Github, after the joint work of Francesco Serafin (here) and Giuseppe Formetta.

Direct Solar Radiation Models by Formetta et al. 2012-2013

This paper presents two new modelling components based on the Object Modelling System v3 for the calculation of the shortwave incident radiation on complex topography settings, and the implementation of several ancillary tools. To understand it, it can be useful to give a look to the previous post on radiation of few weeks ago. The first component, NewAGE-SwRB, accounts for slope, aspect, shadow and the topographical information of the sites, and use suitable parametrisation for obtaining the cloudless irradiance. A second component, NewAGE-DEC-MOD's is implemented to estimate the irradiance reduction due to the presence of clouds, according to three parameterisations. To obtain a working modelling composition, suitable to be compared with ground data at measurement stations, the two components are connected to a Kriging component, and, with the use of a further component NewAGE-V (verification package), the performance of modeled is quantitatively evaluated. The two components (and the various parametrisations they contain) are tested using the data from three basins catchments, and some simple verification test is made to assess the goodness of the methods used. The components are part of a larger system, JGrass-NewAGE, their input and outputs are given as geometrical objects immediately visualisable in a GIS (for instance the companion uDig), and can be used seamlessly with the various modelling solutions available in JGrass-NewAGE for the estimation of long wave radiation, evapotranspiration, and snow melting, as well as stand-alone components to just estimate shortwave radiation for various uses. The modularity of the approach is shown to be extensible to more accurate physical-statistical studies aimed to assess in deep the components performances and to extend spatially their results, without the necessity of recoding any part of the component but just making use of connective scripts.

This is, obviously not the first effort in such direction, and some other good softwares were produced as 
SolarFlux (in ArcInfo GIS) (Dubayah and Paul, 1995; Hetrick et al., 1993), Solar Analyst (Fu and Rich, 2000), SRAD (by Moore, 1992, and documented in Wilson and Gallant, 2000- but you can see this), Solei (Miklanek, 1993 - see this) or r.sun (Hofierka and Suri, 2002), and often integrate the models in GIS. 
Our modelling, making treasure of these previous efforts, is also in line with those tools that try to respond to the increase demand of modularity and interchangeability in hydrological and biophysical models and have been developed in the last decades.

Executable, Data, and Documentation of the system can be found following this link (not yet operational - give me a little time for preparing it, with the objective that everyone can reproduce the results of the paper). Source code is partially available through the jgrasstool page

Friday, December 14, 2012

From Yesterday Just Open Data from the Italian Administrations ?

Yesterday Mr. Monti government Agenda Digitale decree-Law was transformed in state law. Hopefully this would make open all the environmental data produced by public administrations.
Let's see, and hope.
The text of the law can be found here.

Friday, December 7, 2012

Solar Radiation Physics and Geometry for hydrologists

I started to realize that a hydrologist has to learn about radiation, since the first moment I put my sight a little beyond rainfall-runoff moedelling, at least fifteen years ago. This was in building the first bricks of GEOtop: the first problem we faced was to avoid give radiation to point in shadows, i.e. calculating shadows. At the time this was not very much documented, even if the papers by Ralph Dubayah, Jeff Dozier, and others were already out since a few years. I started with giving a look to IPW and I tried to understand what and how it works, but eventually I finish to try my own modules.  To summarise some of the topic,  I renovated my slides on slideshare, which now better cover the notation, and can be considered as a preliminary reading before going to the Formetta et al. 2012 paper (see it on GMDD).
Radiation  treated in the slides mostly regards: geometry (where we finally followed the work of Javier Corripio), shortwave radiation treatment, and longwave radiation.
If, from many points of view working  with a real model of atmospheric absorption would be more exciting, hydrologists use mostly parameterizations, and that is what I mostly summarised in my slides together with a little of astronomical geometry.

Anyway, I found really useful the papers you see listed below.

References

(Mostly about short-wave radiation)

Bird, R. and Hulstrom, R.: Simplified clear sky model for direct and diffuse insolation on horizontal surfaces, Tech. rep., Solar Energy Research Inst., Golden, CO (USA), 1981.

Bird, R. E.,  Riordan, C, 1986: Simple Solar Spectral Model for Direct and Diffuse Irradiance on Horizontal and Tilted Planes at the Earth's Surface for Cloudless Atmospheres. J. Climate Appl. Meteor., 25, 87–97. doi: http://dx.doi.org/10.1175/1520-0450(1986)025<0087:SSSMFD>2.0.CO;2

Boland, J., Scott, L., and Luther, M.: Modelling the diffuse fraction of global solar radiation on a hori- zontal surface, Environmetrics, 12, 103–116, 2001.

Corripio, J.: Modelling the energy balance of high altitude glacierised basins in the Central Andes., PhD dissertation, University of Edinburgh, 2002.

Corripio, J.: Vectorial algebra algorithms for calculating terrain parameters from DEMs and solar radi- ation modelling in mountainous terrain, International Journal of Geographical Information Science, 17, 1–24, 2003.

Dozier, J. and Frew, J.: Rapid calculation of terrain parameters for radiation modeling from digital elevation data, Geoscience and Remote Sensing, IEEE Transactions on, 28, 963–969, 1990.

Dubayah, R.: Modeling a solar radiation topoclimatology for the Rio Grande River Basin., Journal of vegetation science : official organ of the International Association for Vegetation Science, 5, 627–640, http://ukpmc.ac.uk/abstract/AGR/IND20453099, 1994.

Dubayah, R. and Paul, M.: Topographic solar radiation models for GIS, International Journal of Geo- graphical Information Systems, 9, 405–419, 1995.

Duguay, C.: Radiation modeling in mountainous terrain review and status, Mountain Research and Development, pp. 339–357, 1993.

Erbs, D., Klein, S., and Duffie, J.: Estimation of the diffuse radiation fraction for hourly, daily and monthly-average global radiation, Solar Energy, 28, 293–302, 1982.

Gubler, S., Gruber, S., and Purves, R.: Uncertainties of parameterized surface downward clear-sky shortwave and all-sky longwave radiation., Atmos. Chem. Phys, 12, 5077–5098, 2012.

Helbig, N., Lowe, H., Mayer, B., and Lehning, M.: Explicit validation of a surface shortwave radiation
balance model over snow-covered complex terrain, Journal of Geophysical ResearchAtmospheres, 115, D18 113, 2010.

Long, C. and Ackerman, T.: Surface measurements of solar irradiance: A study of the spatial correlation between simultaneous measurements at separated sites, Journal of Applied Meteorology, 34, 1995. Mikla ́nek, P.: The estimation of energy income in grid points over the basin using simple digital elevation model, in: Annales Geophysicae, vol. 11, 11 European Geophysical Society, Springer, 1993.

Orgill, J. and Hollands, K.: Correlation equation for hourly diffuse radiation on a horizontal surface,
Solar energy, 19, 357–359, 1977.

Ranzi, R. and Rosso, R.: Distributed estimation of incoming direct solar radiation over a drainage basin, Journal of Hydrology, 166, 461–478, 1995.
Spencer, J.: Fourier series representation of the position of the sun, Search, 2, 172, 1971.

Tovar, J., Olmo, F., and Alados-Arboledas, L.: Local-Scale Variability of Solar Radiation in a Mountainous Region., Journal of Applied Meteorology, 34, 2316–2328, 1995.

Thursday, November 29, 2012

A simple and trivial consideration about rainfall return period in relation with a spatial analysis, and discharges

Statistics at station, i.e. evaluated at ground measurements stations do not reflect the statistics of the areal event. This is obvious. However, even today, I read a report where the guys:
  1. estimated the return period of rainfall in several point 
  2. used the estimated depths to infer the spatially varying rainfall intensities
  3. use these rainfalls as inputs of a rainfall-runoff model to obtain extreme discharges
In between there are a lot of technicalities, often useless. The point is, which is clear to the most, I hope, that, when moving from point 1 to point 2, one assumes that all the measured events are isochronous, which is not (otherwise we could not have let say 200 hundred return period events each year ia a space-wide area).  The above operation actually correspond to consider a precipitation with a higher return period than established (usually following some design criterion) and therefore maximise excessively the discharges.


What one should do is instead:
  1. studying the spatial statistics of precipitation to enable a stochastic weather generator^1^2
  2. run in continuos time her/his rainfall-runoff model for a long period, say 20 years if one wants to get some statistics a return period of 4 years,
  3. analyse the results and extracting the statistics of the discharges (i.e. their return period), eventually extracting the extreme events
The validity of each component of the modelling chain should have been tested against available data (of the basin) independently.  Traditionally engineers do not like to simulate events at continuos time, and prefer to model events. This latter approach, however, has several drawbacks, and especially:
  • one has to determine the initial conditions (which also introduce a bias in the return period) of the catchments (models that do not have this problem cannot be good models)
  • fall back into the issue of determining a spatially distributed rainfall with a certain return period
Engineers usually also neglect the role of snow in producing discharges. This cannot be neglect, except than in particular climatic conditions. Using continuos time simulations also implies the use of some parameterisation of evapotranspiration (and therefore requires a model like JGrass-NewAGE).

^1 - Remarkably using a weather generator can also allows the inclusion of foreseen trends (either in precipitation characteristics, as depth, interstorm inter-arrival time, or evapotranspiration or radiation).

^2 -  Usually these models are site depedent. Therefore, waiting for a spatial stochastic weather generator, one should run several copies of the weather generator, each one for each sites where there can be information to drive it, and subsequently, using the spatial data, spatially interpolate at each time-step the desired quantity.

Monday, November 19, 2012

Repertorio Nazionale dei dati



Prendo pari pari da JGrassTechTips una iniziativa che mi sembra utile promuovere.



Anche AboutHydrology inserisce oggi nella sua homepage il link del Repertorio Nazionale dei Dati Territoriali (RNDT). Attualmente questo strumento sta muovendo i primi passi, ma noi desideriamo che entri a fare parte della “cassetta degli attrezzi” di ogni geomatico e, in prospettiva, di ogni cittadino italiano. Non siamo da soli, ma insieme a:


Leggi l'articolo su TANTO.

European waters - assessment of status and pressures

This report's results present good and robust European overviews of the data reported by the first RBMPs, and of the ecological status and pressures affecting Europe's waters. Europe's waters are affected by several pressures, including water pollution, water scarcity and floods. Major modifications to water bodies also affect morphology and water flow. To maintain and improve the essential functions of our water ecosystems, we need to manage them well.
The report was produced by the European Environment Agency and can be downloaded following the link.

Saturday, November 17, 2012

uDig-JGrasstools Resources (in Italian)

Per usare uDig, la prima cosa è ovviamente installarlo seguendo le istruzioni. Fatto questo, il secondo passo è quello di istallare ed usare lo Spatial Toolbox.
Quindi non resta che usarlo!
Riepilogo qui alcune risorse alle quali gli utenti italiani di quel che fu JGrass possono fare ricorso

  • Il materiale sulla Horton Machine, ovvero per il trattamento dei dati digitali del terreno, è, ora, in un post separato
  • Il materiale di supporto alla progettazione di acquedotti e fognature è invece qui.
In questo link potete trovare anche il manuale di JGrass 3.0 ... che è sorpassato: ma, ad esempio, come costruire le legende si trova solo qui.

Per quelli che sono interessati alla capacità di scripting di uDig, allora è il caso di dare una occhiata:
Altre risorse sono certamente utili. Per esempio,  la serie di screen-cast che si puo' trovare su 


benche' siano in Inglese. Quanti volessero aiutare, li invito  a leggere il post di Andrea Antonello.


Friday, October 26, 2012

An introduction to Forecast, Errors, Uncertainty in Models, and Forecast Correction

Following the link below, you can find  the presentation I gave at the Pensa Trasversale initiative. I talked about modelling and its uncertainty, and I had, obviously, in mind my experience as a hydro-geomorpho-logical modeler.
The last topic in the title, forecast correction, is not as deep investigated as I would have. But in writing I realised that I had to speak to High School Students.
The presentation is in Italian, but I can translate it in English if someone is interested.  In the case, I should do a further investment in making some parts more technical and add some examples.

The presentation can be found, as usual, on slideshare. A companion of the presentation, also in Italian, of this presentation can be this paper that I wrote more than twenty years ago about modelling, entitled: Introduction to the Mathematical Models of Environment.  I dedicated this small contribution to Sandro Marani with whom I discussed a lot the topic.

Friday, October 12, 2012

Water Scarcity

Water scarcity is a reality for large populations in the World, and several projects started to understand it and to provide solutions.
Also regions where water is usually abundant suffer the threats of climate change and could be subject and conflicts about waters uses.
Today I was informed of the Glowasis blog   where many information are collected on the subject.

Thursday, October 4, 2012

Rivers in transitions

An editorial of Nature Geoscience to read.


"Rivers affect landscape structure and function to a much greater extent than might be expected from the fraction of the Earth's surface they cover. Rivers redistribute material as they flow, carving out canyons and building new land offshore. These morphological consequences of river flow are evident in any topographic map of the Earth's surface. ..."

Other links are also available at the page.

Browsing the same number also a commentary is available on the role of rivers, and a paper on river drainage patterns in New Zealand.

Monday, October 1, 2012

Guidelines for the Mapping of the Triggering of Landslides and Debris Flow


My studies on shallow landslides were not purely theoretical but directed to make safer the mountain environment in which I live. Therefore, since the beginning of my activities there was an effort to convert theoretical results into practical tools, which, in turn, helped research. These guidelines written for the Danube Flood Risk Project come with this attitude. The work was also supported by the IRASMOS EU Project and, more recently, from the Trento Province. While reading the guidelines themselves is probably the simplest way to approach the mapping of landslide triggering according to my perspective, I also make public the presentation that I gave last and this year on the subject.

The first presentation is an introduction to the subject of hydrological hazards in mountains areas and the topic of guidelines:



The second contains and comments some applications of the guidelines on catchments in Trentino.



All the operation seen (except for the very recent CI-SLAM model) can be reproduced using the tools in the Spatial Toolbox of uDig or using GEOtop.


References

Beven, K J and Kirkby, M J. 1979, A physically based variable contributing area model of basin hydrology Hydrol. Sci. Bull., 24(1),43-69

Beven, K, Rainfall-runoff modelling: the primer, Wiley, 2001

Borga, M., G. Dalla Fontana, F. Cazorzi, Analysis of topographic and climatic control on rainfall-triggered shallow landsliding using a quasi-dynamic wetness index, Jour. Hydrol., 268, 56-71, 2002
D’Odorico, P. and R. Rigon, Hillslope and channels contribution to the hydrologic response, Water Resour Res, 39(5) , 1-9, 2003

Lanni, C.; McDonnell, J. J.; Rigon, R., On the relative role of upslope and downslope 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 depth and 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 shallow landslide occurrence by means of a subsurface flow path connectivity index, Hydrol. Earth Syst. Sci. Discuss., 9, 4101-4134, www.hydrol-earth-syst-sci- discuss.net/9/4101/2012/ doi:10.5194/hessd-9-4101-2012, (in press at HESS)

Other papers and material about landslides can be found in this blog following the "Landslides" label.

Saturday, September 29, 2012

My Past Research on Cryopheric Hydrology


In [J22] it was demonstrated that a single-layer snowpack model can be sufficiently accurate in describing the evolution of the water equivalent of the snow, as long as the incident radiation is calculated accurately taking care of shadows and the complexity of mountain topography.


Subsequently, the single-layer model was replaced with a multilayer model in order to forecast the evolution of density and of metamorphism of the snow as well as the percolation phenomena within the snowpack, during the thesis of Stefano Endrizzi. Among the various studies carried out, one validates the snow model satellite data derived from MODIS [A41].  Furthermore, the same model was used to study the hydrological evolution of glaciers in Trentino (Alpine) and South America (Equatorial) [A39, A47].  Eventually, the modeling of the cryosphere moved towards considering evolutive processes of permafrost [thesis of Matteo Dall'Amico, and J30], that is the layer of soil subject to temperatures below zero centigrades for more than two consecutive years.  All of these research projects, as well as allowing the aforementioned studies, are necessary to modeling the entire yearly hydrological cycle in mountain environments such as Trentino.

[J30], drawing from an accurate work of reanalysis of process thermodynamics, implements a robust method for the integration of the freezing-soil equation.  The numeric algorithm used is globally convergent Newtonian method that is appropriate for the equations under study.  [J36] is a geomorphological survey of rock glaciers in Trentino, to be subsequently modelled with GEOtop.

References in English


[ J22] - Zanotti, F., Endrizzi, S, Bertoldi, G. e R. Rigon, The GEOTOP snow module, Hydrol. Proc., 18, 3667-3679 (2004), DOI 10.1002/hyp.5794

[j30]- M. Dall’Amico, S. Endrizzi, S. Gruber, and R. Rigon, An energy-conserving model of freezing variably-saturated soil, The Cryosphere, 5, 469-484, 2011, doi:10.5194/tc-5-469-2011

[J36] - R. Seppi, A. Carton, M. Zumiani, M. Dall’Amico, G. Zampedri, R. Rigon, "Inventory, distribution and topographic features of rock glaciers in the southern region of the Eastern Italian Alps (Trentino)" in Geografia Fisica e Dinamica Quaternaria, v. 2012, n. 35(2) (In press)

[A41] Endrizzi S., Bertoldi G., Neteler M., and Rigon R., Snow Cover Patterns and Evolution at Basin Scale: GEOtop Model Simulations and Remote Sensing Observations, Proceedings of the 63th Eastern Snow Conference,


References in Italian

[A47] Noldin I., Endrizzi S., Rigon R., Dall’Amico M, Sistema di drenaggio di un ghiacciaio alpino, Neve e Valanghe, n. 69, 48-52, 2010



My Past Research on Physico-Statistical Modelling of the Water Cycle at Basin Scale

While GEOtop [J24, J25] is for process-based modelling of the mass and energy budgets at a small scale, in order to model larger catchments, which include abstraction works or hydraulic structures, it was decided to implement a new modelling system JGrass-NewAGE [J34].  This system sacrifices process details in favour of  efficient calculations.  It is made of components apt at returning statistical hydrological quantities, opportunely averaged in time and space.  One of the goals of this implementation effort was to create the basis for a physico-statistical hydrology in which the hydrological spatially distributed dynamics is reduced into low dimensional components, when necessary surrogating the internal heterogeneities with "suitable noise" and a probabilistic description.


Unlike other efforts of synthesis, JGrass-NewAge wants to keep the spatial description explicit, at various degrees of simplicity.  This has been made possible by opportune processing of distributed information which, in this way, has become part of the model itself.
From the point of view of the information technology used to implement the modelling  [J41, A44, A49, A50], the system is based on the OMS v 3 system, which allows the use of modern, object-oriented strategies for the structuring of the deployment of the software and, at the same time, furnishing not a model, but various, interchangeable, modeling solutions (MS) that can be adapted to the problems in hand and the practical demands of the problem being solved.
The modeling system, as well as the components to model the physical processes themselves, also includes various tools for the processing of input data (for example, Kriging tools), including all the tools of the Horton Machine [eb3] for the processing of digital terrain data, and the tools for the treatment and interpretation of the output data, for the calibration of model parameters, and (in perspective) for continuous data assimilation.
With this in mind, an effort that is currently being made is that of creating an opportune digital watershed scheme that can accommodate the needs of the various modeling conceptualizations and the identification of areas that are hydrologically "similar" that can be treated conjointly during the calculation of flows and storage. At the moment, model solutions use standard implementations.  [J34, J41, A50] contains the description of the rainfall-runoff part of the modelling system; [J43] is a verification of the radiation budgets components; [J44] is an example of simplified snow modelling.  As a standard, any components is verified by itself against the data relative to the process that it describes, using various automatic calibration procedures, and quantitative objective functions. [J34, J41] using the infrastructure show how increased geomorphological (and processes) information affects the quality of reproduction of the hydrologic response. [j44] explains the watershed partition, based on a generalisation of the Pfafstetter numbering scheme, that guide the functioning of the JGrass-NewAGE system.


References 

In English:

[J24] - Rigon R., Bertoldi G e T. M. Over, GEOtop: A distributed hydrological model with coupled water and energy budgets, Vol. 7, No. 3, pages 371-388

[J25] Bertoldi G. R. Rigon e T. M. Over, Impact of watershed geomorphic char- acteristics on the energy and water budgets, Vol. 7, No. 3, pages 389-394, 2006

[J34] - Formetta, G.; Mantilla, R.; Franceschi, S., Antonello A., Rigon R., The JGrass- NewAge system for forecasting and managing the hydrological budgets at the basin scale: models of flow generation and propagation/routing, Geoscientific Model Development Volume: 4 Issue: 4 Pages: 943-955, DOI: 10.5194/gmd-4- 943-201, 2011

[A49] Formetta G., Antonello A., Franceschi S., David O. and Rigon R., The informatics of the hydrological modelling system JGrass-NewAge, 2012 International Congress on Environmental Modelling and Software Managing Resources of a Limited Planet, Sixth Biennial Meeting, Leipzig, Germany R. Seppelt, A.A. Voinov, S. Lange, D. Bankamp (Eds.) http://www.iemss.org/society/index.php/iemss- 2012-proceedings, 2012

[j36] - Formetta G., Rigon R., Chavez J.L., David O., The short wave radiation model in JGrass-NewAge System, Geosci. Model Dev., 6, 915-928, 2013, www.geosci-model-dev.net/6/915/2013/
doi:10.5194/gmd-6-915-2013

[J39] - Formetta G., Antonello A., Franceschi S., David O., and Rigon R., Hydrological modelling with components: A GIS-based open-source framework, Environmental Modelling & Software, 5 (2014), 190-200

[j42] - Formetta G., David O., Kampf S., Rigon R., The Cache la Poudre river basin snow water equivalent modeling with NewAge-JGrass, accepted GMD, 2014

[j44] Formetta G. , Antonello A. , Franceschi S. , David O., Rigon R.,  Digital watershed representation within the NewAge-JGrass system. Boletin Geologico y Minero, 125 (3): 371-381, 2014. ISSN: 0366-0176


In Italian:

[A44] Antonello A., Franceschi S., Formetta G., Rigon R., L’infrastruttura NewAGE per la previsione e la gestione dei bilanci idrici a scala di bacino: I - La struttura informatica, in Atti XXXII Convegno di Idraulica e Costruzioni Idrauliche, Palermo, 14-17 Settembre 2010

[A45] Formetta G., Franceschi S., Antonello A., Cordano E., Mantilla R., Rigon R., Il sistema NewAGE per la previsione e la gestione dei bilanci idrici a scala di bacino. II - I modelli di generazione, aggregazione e propagazione del deflusso. in Atti XXXII Convegno di Idraulica e Costruzioni Idrauliche, Palermo, 14-17 Settembre 2010

[A50] Formetta G., Rigon R, Le nuove componenti modellistiche di JGrass-NewAGE, Atti del XXXIII Convegno di Idraulica e costruzioni Idrauliche, Brescia, 10-15 settembre 2012

Friday, September 28, 2012

My Past Research on Shallow Landslide and Mass Flow Triggering


The role of hydrology in triggering mass movements was initially confronted with an implementation of the theories of Montgomery and Dietrich [1994] (MD), and the case of instability caused by surface runoff [A21, A26, A27].  The study then continued with the analysis of transient phenomena, that is the instabilities caused by the propagation of pressure waves in the unsaturated medium  [A31, A32], according to the theory by Iverson [2000] (I), and integrating  the two, MD1994 and I2000, views even in the case of rainfall of varying intensity [A21, J23].
Then, the simplified approach  (important in as so much as it highlighted some qualitative aspects of infiltration in the hillslopes) was supplanted by the use of the GEOtop model for the continual simulation of hydrological variables [A38, J26], and transient effects, within a minimal set of simplifications.  The use of GEOtop has allowed for the separation of the hydrological part, effectively modeled by GEOtop, and the geotechnical part, contained in the GEOtop-FS model [J26].  Particularly, the latter of these was the subject of a probabilistic treatment that introduced uncertainties into the main geotechnical parameters  [J26].
The paper [J26], and the thesis of Silvia Simoni introduced a systematic approach to the identification of areas of instability that made full use of the potential of on-site geophysical measurement campaigns and the a priori characterization of geotechnical properties of the soil in the laboratory, without using back analyses for the calibration of parameters as is generally done by simplified models.  The IRASMOS Reports [rep06, rep07 and rep08] represent a summary of the literature available on this subject which has been eventually refined in [rep09].


The most recent work  has been focused on trying to understand the dynamics of subsurface flow  in  by means of virtual experiments [A43] with GEOtop, and in more conceptualized terms to explicit the role of the variability of depth of soil [J33, J35, thesis of Cristiano Lanni]  with the model denominated CI-SLAM.  The result is the introduction of the concept of "hydrological connectivity" of the hillslopes, which is realized when a perched water table forms that covers the whole basin.  The connectivity concept bridged the gap between hillslope hydrology and basin hydrology, and has also consequences important for hillslopes' stability [J37]. In fact these concepts allows a better statistical identification of landslide areas, than previous similar models.  [J35] also contains a preliminary attempt to use the theories of self-organizing criticality in the context of instability propagation, which, evidently, heralds the actual landslide itself.

Paper [J46] faces the issues related to the choice of a certain parameterisation of the soil retention curves and analyses their relation to hillslope stability. It uses a new theory that uses double porosity, and estimates the stability with the use of the new theories by Lu, Likos and Godt.

References

In English:

[ J23] - D’Odorico, P., Fagherazzi G., Rigon R. Potential for landsliding: Dependenceon hyetograph characteristics J. Geophys. Res., Vol. 110, No. F1, F01007 10.1029/2004JF000127 10 February 2005

[J26] Simoni, S., F. Zanotti, G. Bertoldi and R. Rigon, Modelling the probability ofoccurrence of shallow landslides and channelized debris flows using GEOtop-FS, Hydrol. Process. 22, 532–545, 2008, DOI: 10.1002/hyp.6886

[J33] - 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

[J35] - 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

[J37] 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, (in press at HESS)

[A31] - E. Cordano, P., Bartolini, Rigon R. A flexible numerical approach to solving a generalized Richards’ equation problem and some applications, 2004

[rep06]- Rigon R., Rickenmann D., Catalogue of causes and triggering thresholds (Ed), IRASMOS EU Project Deliverable 1.1, 2007

[rep07] - Rigon R. (Ed), State-of-the-art models: their transferability and model application, IRASMOS EU ProjectDeliverable 1.2, 2007

[rep08] - R. Rigon, State of the art of prediction techniques, IRASMOS EU Project Deliverable 1.3, 2007

[rep09] - R.Rigon, Franceschi, S., Monacelli, G., and Formetta, G., The triggering of landslides and debris flows and their mapping, Danube Flood Risk EU Project, 2012

[J46] - Ciervo F. ,  Casini F. , Papa M.N. ,  Rigon R., Some remarks on bimodality effects of the hydraulic properties on shear strength of unsaturated soils, Vadose Zone Hydrology, published electronically, doi:10.2136/vzj2014.10.0152, 2015

In Italian:

[A21] - D’Odorico, P., Fagherazzi S., Rigon R. Frane superficiali e idrologia deiversanti: Un possibile metodo di indagine. Atti del XXVIII Convegno di Idraulica e Costruzioni Idrauliche, vol. V, pp.177-184, 2002

[A26] - Tiso, C., Bertoldi G. and R. Rigon. Il modello Geotop-SF per la determinazione dell’nnesco di fenomeni di franamento e di colata. Atti del Convegno Iterpraevent 2004, Riva del Garda, 24-28 Maggio 2004

[A27] - Rigon, R., A. Cozzini, S. Pisoni, G. Bertoldi e A. Armanini. A new simple method for the determination of the triggering of debris flows. Atti del Convegno Interpraevent 2004, Riva del Garda, 24-28 Maggio 2004

[A32] - Cordano, E., Panciera R., Rigon R., Bartolini P. Sulla soluzione diffusiva dell’equazione di Richards. Atti del XXIX Convegno di Idraulica e Costruzioni Idrauliche, Settembre 2004

[A43] Lanni C., Cordano E., Rigon R., Tarantino A., Analysis of the effect of normaland lateral subsurface water flow on the triggering of shallow landslides witha distributed hydrological model. in from geomorphology mapping to dynamic modelling, Strasbourg: CERG, 2009. Atti di: A Tribute to Prof. Dr. Theo van ASCH, Strasbourg, 6th-7th February 2009

Thursday, September 27, 2012

My Past Research on Hydroinformatics, GIS and Modelling by Components


Research in the aforementioned sectors was also carried out with the implementation of  open-source software, coded in C and Java and distributed with a GPL (v 3) license.  Involvement in this topic has been deemed necessary to easy cooperative research, and to improve reuse of codes among researchers and students,  and allow an incremental development of modelling solutions (avoiding to implement again and again the same algorithms at any new generation of students).
The software has been accurately documented [eb1 to eb13] so that it can be easily reused and modified for both research and didactic purposes.  The software originally included a series of C libraries for reading, writing and insertion of comments in the data files, dynamic allocation of memory, the statistical treatment of data aimed  especially at hydrology, hydraulics, and geomorphology, but not limited to these.  On the basis of these libraries, called "Fluid Turtles" and now obsolete, was implemented the initial version GEOtop model [j24, s3]  and an initial version of  the  Horton machine [e.g. eb-3, s3].


However, the traditional software architecture of the Fluid Turtles presented various limitations. Mainly: the lack of an interface for the processing of  input data and the treatment of output data and the difficulty of maintaining and testing the software and its parts each one independently from the others, with the growing number of processes being described [e.g. A44, A49]. In fact, with the increasing number of people working on the code, and with the success of the modelling ideas among users, it became necessary to be able to test and use groups of parts of the models separately (as in JGrass-NewAGE).  It was also envisioned necessary to predispose the models  to be linked (in the future) to external models, such as, for example, those simulating the evolution of the atmospheric boundary layer, or belonging to other domains than hydrology (for instance to build a Decision Support System).
These, and other reasons [J40], have brought through a decade of work, trials and errors, on the one hand, to the development of a new GIS,  JGrass, eventually embedded in uDig,  and, on the other, to the adoption of suitable informatics infrastructure in order to restructure the models in components according to the OMS standard.

In the latest version, JGrass has partially contributed to the uDig "core" [eb10, A49, J41], while the modeling part is migrating to the jgrasstools environment (based on OMS) called Spatial Toolbox.
The last version of the tools has been actually embedded in Hydrologis' S.T.A.G.E which is a stand-alone application connectable, in principle, to any Java GIS (thinking to future versions of uDig or GvSig). The tools for terrain analysis included in STAGE  (a.k.a "The Horton Machine") are well covered by [a57]

Along the years various prototypes where developed around the above infrastructures to connect models to SQL/Geographic databases (Postgresql/Postgis), to visualise results on the Nasa World Wind virtual globe, and to allows scripting to interact with models which were presented in various conferences, and on which we could discuss with those interested.

References

In English:

[ J24] - Rigon R., Bertoldi G e T. M. Over, GEOtop: A distributed hydrological model with coupled water and energy budgets, Vol. 7, No. 3, pages 371-388

[A49] Formetta G., Antonello A., Franceschi S., David O. and Rigon R., The informatics of the hydrological modelling system JGrass-NewAge, 2012 International Congress on Environmental Modelling and Software Managing Resources of a Limited Planet, Sixth Biennial Meeting, Leipzig, Germany R. Seppelt, A.A. Voinov, S. Lange, D. Bankamp (Eds.) http://www.iemss.org/society/index.php/iemss- 2012-proceedings, 2012

[J40] - Formetta G., Antonello A., Franceschi S., David O., and Rigon R., Hydrological modelling with components: A GIS-based open-source framework, Environmental Modelling & Software, 5 (2014), 190-200

[a57]- W. Abera, A. Antonello, S. Franceschi, G. Formetta, R Rigon , "The uDig Spatial Toolbox for hydro-geomorphic analysis" in Geomorphological Techniques, v. 4, n. 1 (2014), p. 1-19.

[eb2-b] Ghesla, E. and R. Rigon, A Tutorial for the Management of Digital Terrain Models, pg. 131, University of Trento, Department of Civil and Environmental Engineering, ISBN 10: 88-8443-155-7, 2006 (Now obsolete)

[eb3] - R.Rigon, E. Ghesla, C. Tiso and A. Cozzini, The Horton Machine, pg. viii, 136, ISBN 10:88-8443-147-6, University of Trento, 2006 (Now obsolete)

[eb4-a] Ghesla E and R. Rigon, A Tutorial for preparing GEOtop Input Files with JGrass, pg. vi, 62, ISBN 10:88-8443-153-0, University of Trento, 2006

* [eb05] Dall’Amico, A., Endrizzi, E., Gruber, S., Rigon R., The GEOtop Manual, Università di Trento, in press, 2013

In Italian:

[A44] Antonello A., Franceschi S., Formetta G., Rigon R., L’infrastruttura NewAGE per la previsione e la gestione dei bilanci idrici a scala di bacino: I - La struttura informatica in Atti XXXII Convegno di Idraulica e Costruzioni Idrauliche, Palermo, 14-17 Settembre 2010

[eb1-a] Antonello, A., S. Franceschi, A. Vitti e R. Rigon, Il Manuale JGRASS 2.0 (In Italiano), pg. 176, ISBN 10:88-8443-144-1, University of Trento, 2006 (Now obsolete)

[eb2-a] Ghesla, E. and R. Rigon, Un tutorial per il trattamento di modelli digitali del terreno con JGRASS - A Tutorial for the treatment of DEMs with JGRASS (in Italian), ISBN 10: 88-8443-146-8, 2006 (Now obsolete)

[eb4-b] Ghesla E and R. Rigon, Un tutorial per la generazione dei file di input per GEOtop utilizzando JGrass, pg. vi, 62, ISBN 10:88-8443-154-9, University of Trento, 2006 (Now obsolete)

* [eb6] Rigon R., Formetta G., Zini M., Franceschi S., Antonello A., La Horton Machine, Università di Trento, in press, 2013


* [eb7] Rigon R., Formetta G., Perathoner L., Iemma A., Franceschi S., Antonello A., Jiffle, una breve introduzione, Università di Trento, in press, 2013

* [eb8] Rigon R., Formetta G., Perathoner L., Franceschi S., Antonello A., Peakflow: teoria e pratica, Università di Trento, in press, 2013

* [eb9] Franceschi S., Rigon R., Formetta G., Perathoner L., Antonello A., Trentop, Manuale d’uso, Università di Trento, in press, 2013

* [eb10] - Antonello A., Franceschi S., Rigon R., Formetta G., Perathoner L., uDig: Installare lo Spatial Toolboox, Università di Trento, in press, 2013

[eb11] - Iemma, A., Antonello A., Franceschi S., Rigon R., Formetta G., Perathoner L., uDig walkthroughs, Lavorare con i formati di GRASS in uDig, Università di Trento, in press, 2013

My Past Research on Process Based Physical Modelling on the Hydrological Cycle

The first studies in this field took inspiration from analysis of moisture distribution in the soil, [J14, A8] where it was shown that, during relatively “dry periods, the moisture distribution in the soil can be understood and described with fractal analysis techniques. From these studies, and the desire to model evapotranspiration, eco-hydrological phenomena, hydrology and slope stability, and the evolution of the snowpack, there arose the need to develop an instrument capable of modelling the water cycle and soil moisture dynamics continually over time. These goals were the founding reasons for the implementation of the GEOtop Model [J24, A22]. GEOtop is “terrain-based (it is based on the use of digital terrain models and uses the knowledge of interac- tion between morphology and process) distributed (all the simulated variables are calculated for each pixel of the basin) model of “the water cycle (it simulates all the components of the water cycle, taking account of both the mass budget and the energy budget, the two budget equations being coupled through the temperature of the soil, which controls evaporation, hydraulic conductivity, and accumulation of the snowpack [J22]). A complete description of the model can be found in [J24, A22], articles that present the model system and a practical application to the Little Washita basin in Oklahoma, and, obviously in the manual [eb-05].


The GEOtop model was also applied during the study of the water cycle of Lake Serraia (Trentino, Italy) [A34]. [J25] demonstrates the effects of complex topography and morphology on the water cycle. In particular, one can observe that a more extensive channel network (as might arise in presence of greater slopes or more erodible soil) causes greater surface runoff and less evapotranspiration, which, in the energy budget, causes an increase in latent heat exchange with atmospheric boundary layer. The paper demonstrates, therefore, that topographic effects cannot be neglected in formulating the energy budget of the soil, as most global climate models normally do. Among the more theoretical studies, but essential to the distributed modelling of flows in unsaturated media, are [J27, A29]. In [J28] Richards Equation was perturbatively decomposed into a vertical component and a lateral one. The first dominates the initial phases of infiltration, the second the long-term redistribution of water volumes. With the work reported in [J30] the model was expanded with a soil freezing and thawing module, that allowed the analyses of the PermaNET project to be executed, and other studies performed by other researchers.

Recently, the GEOtop model has been used to estimate the impact of climate change on mountain catchments [rep05]. Ancillary studies have been dedicated to parameter calibration and uncertainties in hydrologic model forecasting [A37, A40].
The article [J38] envisages the restructuring of the GEOtop model with new numeric methods, developed together with Prof Vincenzo Casulli, and the adoption of a non-structured grid for the modeling. However, [J43] represents the state-of art of  GEOtop version 2.0, a milestone in the model history which contains Richards 3D integration, permafrost modelling, a multilayer snow model, renewed options for the treatment of meteo-data and radiation.
[J44] embrace the use of CLM and face the problem of data assimilation complemented by the use of Kriging techniques for filling the missing data.

A small community of users and developers has developed around GEOtop which is steadily growing.

The article [j38] envisages the restructuring of the GEOtop model with new numeric methods, developed together with Prof Vincenzo Casulli, and the adoption of a non-structured grid for the modeling.   [J43] represents the state-of art of  GEOtop version 2.0, a milestone in the model history which contains Richards 3D integration, permafrost modelling, a multilayer snow model, renewed options for the treatment of meteo-data and radiation.
[J44] embraces the use of CLM model, and faces the problem of data assimilation complemented by the use of Kriging techniques for filling the missing data.

References

In English:

[J14] - Rodriguez-Iturbe, I, Gregor K. Vogel, R. Rigon, D. Entekhabi, F. Castelli and A. Rinaldo, On the spatial organization of soil moisture fields, Geoph. Res. Letters, 22(20), 2757-2760, 1995.

[ J22] - Zanotti, F., Endrizzi, S, Bertoldi, G. e R. Rigon, The GEOTOP snow module, Hydrol. Proc., 18, 3667-3679 (2004), DOI 10.1002/hyp.5794

[ J24] - Rigon R., Bertoldi G e T. M. Over, GEOtop: A distributed hydrological model with coupled water and energy budgets, Vol. 7, No. 3, pages 371-388

[ J25] Bertoldi G. R. Rigon e T. M. Over, Impact of watershed geomorphic char- acteristics on the energy and water budgets, Vol. 7, No. 3, pages 389-394, 2006

[ J28] - Cordano E. R. Rigon, A perturbative view on the subsurface water pressure response at hillslope scale, Water Resour. Res., Vol. 44, No. 5, W05407- W05407, doi:10.1029/2006WR005740, 2008 

[j30]- M. Dall’Amico, S. Endrizzi, S. Gruber, and R. Rigon, An energy-conserving model of freezing variably-saturated soil, The Cryosphere Discussion, The Cryosphere Discuss., 4, 1243-1276, doi:10.5194/tcd-4-1243-2010, 2010

[J38] Cordano E., Rigon R., A mass-conservative method for the integration of the two-dimensional groundwater (Boussinesq) equation, submitted to Water Resour. Res., 2012  

[J43] - Endrizzi S., Gruber, S, Dall'Amico M., and Rigon R., 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, Geosci. Model Dev.,  7, 2831–2857, 2014, www.geosci-model-dev.net/7/2831/2014/ doi:10.5194/gmd-7-2831-2014

[J44] -  Han X., Lin, X, Rigon R., Jin R., Endrizzi S., Local Analysis of L-Band Microwave Brightness Temperature Assimilation With Geostatistics for Soil Moisture Estimation, PLOSONE, 2015

[A8] - Rodriguez-Iturbe, I., G. K. Vogel, R. Rigon, D. Entekhabi, F. Castelli and A. Rinaldo, Scaling properties of soil moisture, Proceeding of the workshop on climate change and hydrometereological hazards in the mediterranean area, Perugia, 1995 

[A34] - Bertola P., Bertoldi G.,Grisenti P., Piva G., Ragazzi M., Righetti M., Rigon R., Salvaterra M., Soppelsa G., Tomazzolli V., Integrated research on eutrophication processes on Caldonazzo lake (Trento, Italy). Atti del Convegno Simposio Internazionale di Ingegneria Sanitaria e Ambientale, Taormina, 23-26 Giugno, 2004

[A37] - Entezarolmahdi R., G. Bertoldi e R. Rigon, An Automatic “Watershed Model Calibration process, 7th International Congress of Civil Engineering, Teheran, 2006

[A31] - E. Cordano, P., Bartolini, Rigon R. A flexible numerical approach to solving a generalized Richards’ equation problem and some applications, 2004

[A40] Entezarolmahdi R., Rigon R., Bertoldi G., Assessment of parameter uncertainty for physically based hydrologic model, using automatic optimization approach, XXX Convegno di Idraulica e Costruzioni Idrauliche, Roma 2006 

[eb05] Dall’Amico, A., Endrizzi, E., Gruber, S., Rigon R., The GEOtop Manual, Università di Trento, A Draft here, in press, 2012

[rep05] - R. Rigon, A. Bellin, L. Forlin, H. Fowler, S. Blenkinsop, Testing of climate change scenarios on a case-study catchment using different methodologies, Deliverable C2.4 AQUATERRA EU Project, 2005

In Italian:

[A22] - Bertoldi, G., Rigon R., Overt T.M. Un indagine sugli effetti della topografia sul ciclo idrologico con il modello GEOtop. Atti del XXVIII Convegno di Idraulica e Costruzioni Idrauliche, Potenza, pp.313-324, 2002

[A32] - Cordano, E., Panciera R., Rigon R., Bartolini P. Sulla soluzione diffusiva dell’equazione di Richards. Atti del XXIX Convegno di Idraulica e Costruzioni Idrauliche, Settembre 2004

My Past Research on Hydro-geomorphology


The work on evolution of river networks, certainly enters also in this category. Here however, it is reported about those papers that deals directly with quantitative geomorphological analysis.
In most of the papers, the relationships between the various parts of a river basin are analyzed with fractal geometry techniques.  Such knowledge is useful not only for the evaluation of river basin evolution models, but also in  identifying the nature of their hydrological response to given events and their paleoclimate.  In [J2], using the Peano Basin, a mathematical reference structure, it is suggested that the amplitude function of natural networks can be reproduced with a multifractal multiplicative
process.   In [J11] this concept was rigorously formalized in the framework of random multifractal cascades theory.


The fractal properties, that is  power laws relating to  contributing areas and the length of stream reaches, were rigorously analyzed in [J15]. In [J17] more relationships between these characteristic quantities were found and an explanation of the nature of Hack's law is suggested.
In [J27] the structure of the river networks is further investigated by analysing the tributaries statistics. The relationships were verified experimentally by means of Digital Elevation Models (DEM).

Given programs for the extraction of digital terrain models (DEM) and their treatment, the natural development was the implementation of an open-source geographic information system: JGrass [s2].  JGrass, now part of uDig, contains within the package a large amount of GIS methods, jointly known as the Horton Machine [eb-3] to support the most common and some less common tools of analysis for river networks topology,  channel extraction, hillslope delineation. A review of these tools is in [a57].

Papers and work on landslide triggering are also concerned with geomorphology but referred in a different post.

References

[j2] - Marani, A., R. Rigon and A. Rinaldo, A Note on fractal channel networks,Water Resources Research, (27)5, 3041-3049, 1991.

[j11] - Marani, M., A. Rinaldo, R. Rigon and I. Rodriguez-Iturbe, Geomorphological width function and the random cascade, Geophysical Research Letters, 21(19), 2123-2126, 1994.

[j15] - Maritan, A., A. Rinaldo, R. Rigon, I. Rodriguez-iturbe and A. Giacometti, Scaling laws for river networks, Physical Review E , 53, 1510, 1996.

[j17] - Rigon, R., I. Rodriguez-Iturbe, A. Rinaldo, A. Maritan, A. Giacometti and D. Tarboton, On Hack’s law, Water Resources Research, 32(11), 3367, 1996

[s2] Jgrass (since 2002): it is a full featured GIS system built to support the hydro- geomorphological research of the group since the early 2000s. Recently it has been integrated in uDig. Current developments (which go far beyond my contribution) are available at: http://udig.refractions.net/

[eb3] - R.Rigon, E. Ghesla, C. Tiso and A. Cozzini, The Horton Machine, pg. viii, 136, ISBN 10:88-8443-147-6, University of Trento, 2006

[j27] - Convertino, M, Rigon, R; Maritan, A; I. Rodriguez-Iturbe and Rinaldo, A, Probabilistic structure of the distance between tributaries of given size in river networks, Water Resour. Res., Vol. 43, No. 11, W11418, doi:10.1029/2007WR006176, 2007

[a57]- W. Abera, A. Antonello, S. Franceschi, G. Formetta, R Rigon , "The uDig Spatial Toolbox for hydro-geomorphic analysis" in Geomorphological Techniques, v. 4, n. 1 (2014), p. 1-19

Monday, September 24, 2012

My Past Research on the Evolution of River Networks

Chronologically, one of my first interests was modeling the evolution of channel networks according to principles of minimal energy dissipation and self-organization by critical states. These two types of models proved to be capable of reproducing the two- and three-dimensional statistical characteristics of channel networks and natural basins, as well as the fractal and multifractal characteristics. This work, born with the intent of identifying a minimum set of characteristic dynamic elements in the evolution a hydrographic basin, has always been carried out in parallel to the refinement of measurement and analysis techniques of topographic data [J3]. 

The concept of optimality of a hydrographic basin was introduced in [J3, J4, J5]. In these works, the three postulates of optimal channel networks (OCNs) are stated and developed, proving how such principles can have quantitative effects on the morphology of river networks, particularly affecting the structure of slopes and of contributing areas, the geometry of the channels, and the characteristic velocity of the peak flow of a basin. All of these results explain numerous empirical laws and are still the basis of measurement campaigns.


In [J4], that which was postulated in [J3, J5] was verified by numeric simulation. It should be noted that the minimization of dissipated energy generates fractal forms that reproduce the quantitative characteristics of real basin. In [J6] the concept of optimality is further refined by introducing the hillslope contribution and presenting some case studies. In [J6], more tools are introduced for the qualitative comparison between the numeric models and the natural data. In [J8, J9, A6] a model of the evolution of river basins is presented that is based on the concepts of self-organization by critical states. This model proved to be equivalent to optimization model of [J3-J6]. In [J13] the impact of climatic variability on the morphology of the fluvial landscape is simulated, so offering an interpretative framework for some fluvial forms that can be found in nature.

Subsequently, the concept of optimality was refined observing that real basins do not have the configuration that would give an absolute minimum of dissipated energy, but rather that of states of local minimum that are dynamically accessible. From here the concept of feasible optimality was derived [A10, J18, J19]. It was also demonstrated that the states of absolute minimum, dynamically unreachable, have statistical properties that are not realistic, while accessible minimum states have the desired statistical characteristics. A relevant characteristic of the space-time dynamics of hydrographic networks is that they can be described by means of a parameter that can be linked to temperature [J16]. It is therefore possible to define the thermodynamics of the river networks. As with the thermodynamics of other physical systems, the relevant quantities are energy (dissipated in unit time), entropy, and the temperature.

It should be noted that the space-time evolution of river networks happens with an intermittent behavior similar to the concept of point equilibrium proposed for the evolu- tion of the biological species. It was also demonstrated that the temporal dynamics of river networks is coupled with the spatial activity at all scales and that natural networks, therefore, evolve according to conditions of minimum dissipation of energy but in the presence of a great variety of possible dynamic states.

In [J10] an accurate analysis of the fractal and multifractal properties of optimal river networks was carried out. The note [J18] is a review article, sent to the Annual Review of Earth and Planetary Sciences, that treats the aforementioned topics. In [J20] the results of a theorem on network topology that relates the sum of the contributing areas with the contributing areas themselves and hypothesizes that these quantities are analogous to the ratio of metabolic rhythm and mass of living beings. The two quantities are linked an exponential law with an exponent that was proved to be Hack’s exponent.

This work, born with the intent of identifying a minimum set of characteristic dynamic elements in the evolution of a hydrographic basin, has always been carried out in parallel to the refinement of measurement and analysis techniques of topographic data [s2,eb3]. Recently, this field of study has produced a work [J26] where the morphometric statistics of tributaries of natural rivers and OCNs are studied. These are related to the characteristics of peak flows and they have ecological implications such as, for example, the velocity of diffusion of waterborne diseases and the diffusion of species along the river network. In Rigon’s work, the morphological relations between the different parts of fluvial basins have been analyzed with ever more refined numeric instruments, to the point of creating a series of GIS methods known as the Horton Machine [eb3].

The paper [J41] is partially a review of old results, that were not collected before, and were overlooked by people because they did not appear in Rodriguez-Iturbe and Rinaldo 1997 book. It includes however some new set of simulation were injection of rainfall is assigned with certain distributions (with given correlation structure) producing differentiated power laws for discharge and contributing areas. Clearly a result to further explore.

References

In English:

[J3] - Rodriguez-Iturbe, I. , A. Rinaldo, R. Rigon, R. L. Bras, A. Marani and E.J. Ijjasz-Vasquez, Energy dissipation, runoff production, and the 3-dimensional structure of river basin, Water Resources Research, (28)4, 1095-1103, 1992.

[J4] - Rinaldo, A., I. Rodriguez-Iturbe, R. Rigon, R.L. Bras, E. J. Ijjasz-Vasquez e A. Marani, Minimum energy and fractal structures of drainage networks, Water Resources Research, (28), 2183, 1992.

[J5] - Rodriguez-Iturbe, I. , A. Rinaldo, R. Rigon, R. L. Bras, A. Marani and E.J. Ijjasz-Vasquez, Fractal structure as least energy patterns: The case of river networks, Geophysical Res. Letters, (19)9, 889-892, 1992.

[J6] - Rigon R., A. Rinaldo, I. Rodriguez-Iturbe, R. L. Bras and E. Ijjasz-Vasquez, Optimal channel networks: a framework for the study of river basin morphology, Water Resources Research, 29(6), 1635-1646, 1993.

[J7] - Ijjasz-Vasquez, E., R.L. Bras, I. Rodriguez-Iturbe, A. Rinaldo and R. Rigon, Are river networks OCN?, Advances in Water Resources, 16, 69-79, 1993.

[J8] - Rinaldo, A., I. Rodriguez-Iturbe, R. Rigon, E. Ijjasz-Vasquez, and R.L. Bras, Self organized fractal river networks, Physical Review Letters,70(6), 822-26, 1993.

[J9] - Rigon, R., A. Rinaldo and I. Rodriguez-Iturbe, On landscape self-organization, Journal of Geophysical Research, 99(B6), 11971-11993,1994.

[J10] - Rodriguez-Iturbe, I., M. Marani, R. Rigon and A. Rinaldo, Self-organized river basin landscapes: fractal and multifractal characteristics,Water Resources Research, 30(12), 3531-3539,1994.

[J16] - Rinaldo, A. Maritan, A. Flammini, F. Colaiori, R. Rigon, I. Rodriguez-Iturbe and J. R. Banavar, Thermodynamics of fractal networks, Physical Review Letters, 76(18), 3364-3367, 1996.

[J18] Rinaldo, A., I. Rodriguez-Iturbe and R. Rigon, Channel Networks, Annual Review of Earth and Planetary Sciences, 26, 289-327, 1998

[J27] - Convertino, M, Rigon, R; Maritan, A; I. Rodriguez-Iturbe and Rinaldo, A, Probabilistic structure of the distance between tributaries of given size in river networks, Water Resour. Res., Vol. 43, No. 11, W11418, doi:10.1029/2007WR006176, 2007

[J41] -Rinaldo,  A., Rigon R., Banavar, J., Maritan, A. and Rodriguez-Iturbe, I., Evolution and selection of river networks: Statics, dynamics, and complexity, PNAS 2014


In Italian:

[A04]- Rigon, R., Il clima è scritto nella forma del reticolo idrografico?, Rapporti e studi della commissione di studio dei provvedimenti per la conservazione e la difesa della cittá di Venezia, Tomo CLI, Classe di Scienze ff. mm. e nn., 1-21, 1992.

[A06] - Rigon, R. - Principi di auto-organizzazione nella dinamica evolutiva delle reti idrografiche, Tesi di Dottorato, Università degli Studi di Genova, Firenze, Padova, Trento, 1994

[A12] - Rigon, R., Che cosa guida i processi morfologici nei bacini fluviali? Reti ottime di canali e la legge di Hack, Atti XXVI Convegno di Idraulica e Costruzioni Idrauliche, Vol II, 121, 1998