-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathSnowSeasonAnalysis.Rmd
50 lines (33 loc) · 5.43 KB
/
SnowSeasonAnalysis.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
---
title: "SnowSeasonAnalysis"
output: github_document
---
# Package: SnowSeasonAnalysis
## Introduction
The R package SnowSeasonAnalysis is a group of tools for the analysis of data during winter. In particular there are some scripts for quality check of precipitation data during winter, or there are some scripts for the analysis of snow presence using different methods, and combining results to improve the models. The third script is a procedure for analisys and filtering of snow height signal.
The goal of the project is to perform in R an algorithm to evaluate the precipitation recorded of a tipping bucket, based on paper ESOLIP, Mair E et al. [1]. Using other parameters measured on the station we classify 6 different events and give the opporunity of an user of data to decide if the precipitation recorded is realistic or not.
Another work done during the project was the implementation of an algorithm for snow cover duration. Classical method use daily amplitude and daily maximum of soil temperature for detection snow cover. But early snow fall wasn't detected because the soil freezing slowly, and expecially in early autumn, when temperatures are mild, and snow falls are weak, the soil doesn't react fastly. Using PAR we examine the part of radiation that hit soil and compare with a free sensor (2 m). In this way we improve the classical snow detection algorithm and can obtain some validation points for hydrological models or satellites products
The third part of this package is a snow filtering tool. This tool want to increase quality of snow height data excluding bad value (out of phisycal range and rate) and smoothing noise deriving by sampling condition. There is an open issue deriving to the behaviour of snow height signal during snow melting days. In warm and sunny days during spring the signal varies by a few centimeters during the day. We are not able to distinguish the cause of this strange behaviour. An hypotesis could be the insufficient signal correction with temperature, other analysis suggests that ultrasonic signal seems to penetrate on the "wet" snow more than in dry snow. A complete explanation need more investigation, and an appropriate filtering technique is not so easy to find.
In the package there are some scripts (folder **inst**):
* **Filtering_precipitation_readings.R** that run and save the ESOLIP quality check algorithm. [More details](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Filtering_precipitation_readings.Rmd)
* **Visualize_Filtering_precipitation_readings.R** that load results of Filtering_precipitation_readings.R and plot some graphs using *Shiny*.[More details](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Visualize_Filtering_precipitation_readings.Rmd)
* **Snow_detection_TS_PAR.R** that analyzes Soil Temperature (TS) and Photosynthetically active radiation (PAR) to detect the presence of snow at the station with a daily resolution.[More details](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Snow_detection_TS_PAR.Rmd)
* **Visualize_Snow_detection_TS_PAR.R** that load results of Snow_detection_TS_PAR.R and plot some dynamic graphs. [More details](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Visualize_Snow_detection_TS_PAR.Rmd)
* **Filtering_snow_height.R** that filter snow height signal using some thresholds and some filtering techniques.[More details](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Filtering_snow_height.Rmd)
* **Visualize_Filtering_snow_height.R** that load results of Filtering_snow_height.R and plot a dynamic graphs. [More details](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Visualize_Filtering_snow_height.Rmd)
## Required data
Put in folder *data* the following data:
1. In the folder *Climareport* add/update *Climareport.csv* [See on Section 1](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Filtering_precipitation_readings.Rmd)
2. In the folder *Input data* add/update *Station_xxx_yyy.csv*. Is has to be **hourly aggregated** and must have the following parameters:
+ **TIMESTAMP**: Format supported: **YYYY-MM-DD HH:MM**. Example: "2017-07-15 09:00"
+ **precipitation**: Units mm/h. Precipitation variable
+ **air_temperature**: Units deg C. Air temperature variable
+ **relative_humidity**: Units %. Relative humidity variable
+ **radiation**: Units W/m2. Solar Radiation variable
+ **SNOW_HEIGHT**: Units m. Snow height variable
+ **soil_temperature**: Units deg C . Soil temperature @ 0 cm
+ **phar_up**: Units W/m2. Photosynthetically active radiation on the station
+ **phar_down**: Units W/m2. Photosynthetically active radiation at soil level
3. In the folder *Support files* update (not necessary) **Range_min_max** to modify thresholds to detect ouliers using physical range of values
4. In the folder *Support files* update (not necessary) **Rate_min_max** to modify thresholds to detect ouliers using possible rate of increasing/decreasing
5. In the folder *Snow_Depth_Calibration* add files **Snow_Depth_Calibration_file.csv** with snow sureveys and virtual snow calibration point (see snow height time series and insert the date of end of winter season. Set the value of snow height at this date as 0 cm) [See on Section 3](https://github.com/EURAC-Ecohydro/SnowSeasonAnalysis/blob/master/vignettes/Filtering_snow_height.Rmd)