Reposiroty for GOTM calibration of local lakes in ISIMIP3
The aim of this project is to compare the differences in results among simulations of lake data, including simulations with no calibration, simulations with global calibration, and simulations with specific calibration for each lake.
To do so, we started by gathering data from the ISIMIP3 Local Lakes repository at https://github.com/icra/ISIMIP_Local_Lakes/tree/main/LocalLakes, where we obtained temperature observations and lake hypsographs. Additionally, we retrieved data from a database where data providers included basic information for each lake, such as latitude, longitude, and elevation. This database is accessible on Google Drive at https://docs.google.com/spreadsheets/d/14f_co0zHNuYtWweguvnbdI_p5gxvTvqB9jvcjtRHz_E/edit#gid=0."
We then processed the data. To accomplish this, we began by converting all the temperature files into daily temperature files, all in the same format. Additionally, we converted the hypsograph files into a standard format to ensure compatibility with GOTM.
As the previous data consisted solely of lake observations, and we required additional parameters for GOTM, we also collected meteorological data from the ISMIMIP ************** repository, utilizing the ERAS5 dataset. This meteorological data included information such as temperature, relative humidity, and shortwave radiation. Additionally, we calculated cloud cover and dewpoint temperature using a Python function. After gathering all the meteorological data, we proceeded to compile distinct meteorological files for each lake, which would be utilized by GOTM during its execution.
For the simulations, we aimed to simulate both without calibration and by calibrating with global parameters. To achieve this, we maintained two folders containing identical ERAS5 meteorological data and hypsographs for each lake, but with different GOTM.yaml parameters for the 'no calibration' and 'global calibration' folders. All of this was automated using different scripts.
After we had gathered all the data (hypsographs, meteorological files and configuration yaml files) we automated the simulations for each lake, generating the corresponding "output.nc" files.
With the output files, we subsequently created various plots and computed specific statistics for each lake. In the graphs, we present a comparison between the simulation results and the actual temperature observations. The statistics we calculated include ***********. Subsequently, we compared the results obtained without calibration to those achieved with global calibration.
This repository consists of 4 directories:
- No calibration: It contains the results of the lakes that were simulated without calibration
- Global calibration: It contains the results of the lake simulations calibrated with global parameters
- Calibration: It will store the results of the lake simulations specificly calibrated for each lake
- Scripts: It contains various scripts used for generating files, creating graphics, calculating statistics, and automating the entire process.