Data

Tools

Indexes

Applications

Experiments

Open Sustainable Technology

A curated list of open technology projects to sustain a stable climate, energy supply, biodiversity and natural resources.

MATILDA

Tool for modeling water resources in glacierized catchments.
https://github.com/cryotools/matilda

Category: Natural Resources
Sub Category: Water Supply and Quality

Keywords

climate earth-science glaciers hbv hbv-model hydrology matilda melt modeling python spotpy

Keywords from Contributors

hydrological-modelling

Last synced: about 14 hours ago
JSON representation

Repository metadata

Tool for modeling water resources in glacierized catchments. Combines a temperature-index melt model with the conceptual catchment model HBV and a parameterized glacier area/volume re-scaling routine.

README.md

MATILDA: Modeling wATer resources In gLacierizeD cAtchments

DOI

MATILDA is a Python-based modeling framework for simulating water resources in glacierized catchments. This repository contains the core routines of the comprehensive MATILDA-Online workflow:

  1. The matilda.core module combines a temperature-index melt model with the HBV model for hydrological simulations.
  2. The matilda.mspot_glacier module provides a parameter optimization routine using the spotpy library.

While MATILDA can be used as a stand-alone model, the MATILDA-Online workflow provides comprehensive tools for data acquisition, pre- and post-processing with detailed documentation.

Repository Structure

.
├── Example
│   ├── example_workflow.py           # Example script for running MATILDA
│   ├── forcing_data.csv              # Input data for example workflow
│   └── runoff_data.csv               # Observed runoff data
├── matilda
│   ├── core.py                       # MATILDA core routines
│   └── mspot_glacier.py              # Parameter optimization routines
...

Installation

The tool's dependencies are set to integrate with MATILDA-Online. It requires Python 3.11 and the following libraries:

- xarray
- numpy
- pandas
- matplotlib
- scipy
- datetime
- hydroeval
- HydroErr
- plotly
- spotpy
- pyyaml

The MATILDA package and it's dependencies can be installed on your local machine using pip or a similar package manager. You can either install the package by using the link to this repository:

pip install git+https://github.com/cryotools/matilda.git

...or clone this repository to your local machine, navigate to the top-level directory and use it:

pip install .

Usage

A detailed walkthrough of the proposed modeling workflow and calibration strategy can be found at the MATILDA Online Webpage. For a quick start to the stand-alone model, see the application example and use the following guidelines.

Forcing Data

The minimum input is a CSV file containing time series of air temperature (°C or K), total precipitation (mm), and (if available) evapotranspiration (mm) data in the format shown below. If evapotranspiration is not provided, it is estimated from air temperature according to Oudin et.al. 2010. A series of discharge observations (mm) is used to calibrate the model. If no discharge data are provided, the model runs with default parameters. All datasets require daily resolution.

TIMESTAMP T2 RRR PE
2011-01-01 00:00:00 -18.2 0.00 0.00
2011-01-01 01:00:00 -18.3 0.1 0.00
2011-01-01 02:00:00 -18.2 0.1 0.00
... ... ... ...
Date Qobs
2011-01-01 0.17
2011-01-01 0.19
... ...

The forcing data are scaled to the mean elevations of the glacierized and ice-free subcatchments, respectively, using linear lapse rates. Reference elevations must be provided for the input data, the entire catchment, and the glacierized fraction. Automated routines for catchment delineation and public data download can be found in the MATILDA Online workflow.

Glacier Data

To apply the Δh parameterization of Huss and Hock 2015 within the DDM routine to calculate glacier evolution over the study period, you need to provide data on the initial glacier cover. The routine requires an initial glacier profile containing the spatial distribution of ice over elevation bands at the beginning of the study period in the form of a dataframe:

Elevation Area WE EleZone
3720 0.005 10786.061 3700
3730 0.001 13687.801 3700
3740 0.001 12571.253 3700
3750 0.002 12357.987 3800
.. ... ... ...
  • Elevation - elevation of each band (10 m intervals recommended)
  • Area - glacierized area in each band as a fraction of the total catchment area (column sum is the glacierized fraction of the total catchment)
  • WE - ice thickness in mm w.e.
  • EleZone - combined bands across 100-200 m.

Parameter List

MATILDA has 21 non-optional parameters, most of which are HBV standard parameters.

Parameter Description Unit Default Value
$\text{lr}_{\text{temp}}$ Temperature lapse rate °C m⁻¹ -0.006
$\text{lr}_{\text{prec}}$ Precipitation lapse rate mm m⁻¹ 0
$\text{PCORR}$ Precipitation correction factor - 1.0
$\text{TT}_{\text{snow}}$ Threshold temperature for snow °C 0
$\text{TT}_{\text{diff}}$ Temperature range for rain-snow transition °C 2
$\text{SFCF}$ Snowfall correction factor - 0.7
$\text{CFMAX}_{\text{snow}}$ Melt factor for snow mm °C⁻¹ day⁻¹ 5
$\text{CFMAX}_{\text{rel}}$ Melt factor for ice relative to snow - 2
$\text{CWH}$ Water holding capacity of snowpack - 0.1
$\text{CFR}$ Refreezing coefficient - 0.15
$\text{AG}$ Control parameter of the glacier storage-release scheme - 0.7
$\text{BETA}$ Shape coefficient for soil moisture routine - 1.0
$\text{CET}$ Correction factor for evapotranspiration - 0.15
$\text{FC}$ Field capacity of soil mm 250
$\text{LP}$ Fraction of field capacity for maximum evapotranspiration - 0.7
$\text{K}_0$ Recession coefficient for surface flow day⁻¹ 0.055
$\text{K}_1$ Recession coefficient for intermediate groundwater flow day⁻¹ 0.055
$\text{K}_2$ Recession coefficient for deep groundwater flow day⁻¹ 0.04
$\text{PERC}$ Percolation rate from upper to lower groundwater reservoir mm day⁻¹ 1.5
$\text{UZL}$ Threshold for quick flow from upper zone mm 120
$\text{MAXBAS}$ Length of triangular routing function day 3.0

Authors

See also the list of contributors who participated in the initial project funded by the GIZ.

License

This project is licensed under the MIT License.

References

The development of MATILDA integrated several well-established hydrological and glacier modeling tools. References for the primary methods and libraries used in the model are listed below:

PyPDD (Temperature-Index Model):

LHMP and HBV Models:

Δh (delta-h) Parametrization:

SPOTPY (Parameter Optimization):

Citation (CITATION.cff)

cff-version: 1.2.0
message: "If you use this software, please cite it using the following metadata."
title: "MATILDA: Modeling wATer resources In gLacierizeD cAtchments"
version: "1.0.0"
doi: "10.5281/zenodo.14267361"
date-released: "2024-12-03"
authors:
  - family-names: "Schuster"
    given-names: "Phillip"
    orcid: "0000-0002-0794-6689"
    affiliation: "Humboldt-Universität zu Berlin, Geography Department"
  - family-names: "Georgi"
    given-names: "Alexander"
    orcid: "0009-0000-9465-6761"
    affiliation: "Humboldt-Universität zu Berlin, Geography Department"
  - family-names: "Tappe"
    given-names: "Ana-Lena"
    affiliation: "Humboldt-Universität zu Berlin, Geography Department"
keywords:
  - glacier modeling
  - hydrology
  - water resources
  - temperature-index melt model
  - HBV model
  - parameter optimization
license: "MIT"
repository-code: "https://github.com/cryotools/matilda"
url: "https://matilda-online.github.io/jbook/"
abstract: >
  MATILDA is a Python-based modeling framework designed to simulate water 
  resources in glacierized catchments. It combines a temperature-index melt 
  model with the conceptual catchment model HBV and incorporates a 
  parameterized glacier area/volume re-scaling routine. 
  It also includes parameter optimization capabilities using the SPOTPY library.

Owner metadata

GitHub Events

Total
Last Year

Committers metadata

Last synced: 5 days ago

Total Commits: 665
Total Committers: 8
Avg Commits per committer: 83.125
Development Distribution Score (DDS): 0.484

Commits in past year: 53
Committers in past year: 2
Avg Commits per committer in past year: 26.5
Development Distribution Score (DDS) in past year: 0.094

Name Email Commits
Phillip Schuster p****r@p****e 343
Ana-Lena Tappe a****e@h****e 141
Correct Name a****e@g****m 138
geoalxx a****1@g****e 23
Phillip Schuster p****r@g****e 11
Phillip Schuster s****h@m****r 6
Phillip Schuster s****h@l****s 2
Georgi a****i@s****m 1

Committer domains:

Issue and Pull Request metadata

Last synced: 1 day ago

Total issues: 4
Total pull requests: 17
Average time to close issues: 12 months
Average time to close pull requests: 13 minutes
Total issue authors: 3
Total pull request authors: 2
Average comments per issue: 1.5
Average comments per pull request: 0.0
Merged pull request: 17
Bot issues: 0
Bot pull requests: 0

Past year issues: 1
Past year pull requests: 12
Past year average time to close issues: about 18 hours
Past year average time to close pull requests: 18 minutes
Past year issue authors: 1
Past year pull request authors: 2
Past year average comments per issue: 3.0
Past year average comments per pull request: 0.0
Past year merged pull request: 12
Past year bot issues: 0
Past year bot pull requests: 0

More stats: https://issues.ecosyste.ms/repositories/lookup?url=https://github.com/cryotools/matilda

Top Issue Authors

  • phiscu (2)
  • LKPatel1 (1)
  • KReeves092 (1)

Top Pull Request Authors

  • phiscu (15)
  • geoalxx (2)

Top Issue Labels

Top Pull Request Labels

Dependencies

matilda.egg-info/requires.txt pypi
  • HydroErr *
  • datetime *
  • hydroeval *
  • matplotlib *
  • numpy *
  • pandas *
  • plotly *
  • scipy *
  • xarray *
setup.py pypi
  • HydroErr *
  • datetime *
  • hydroeval *
  • matplotlib *
  • numpy *
  • pandas *
  • plotly *
  • scipy *
  • xarray *

Score: 4.969813299576