This project provides the basis in order to perform a reproducible clustering on water samples which can facilitate the detection of hydrochemical End-members . It has been created in the framework of the MedSal project.
It currently includes following steps:
-
Step 1: Data import
-
Step 2: Data quality assessment and filtering based on the Charge Balance Error
-
Step 3: Cluster analysis
-
Step 4: Data visualization
Input data are either excel/csv files that get imported.
| .xlsx | |
|---|---|
| Data quality assessment | |
| Charge balance error | ✔️ |
| Cluster analysis | |
| Data preparation | ✔️ |
| Clustering | ✔️ |
| Dendrogram | ✔️ |
| Cluster map | ✔️ |
| Data visualization | |
| Piper plot** | ✔️ |
| Shoeller plot | ✔️ |
| Ion ratios | ✔️ |
| Isotope data | ✔️ |
**Code from: https://github.com/markolipka/ggplot_Piper
- Create an .xlsx containing your water analyses. For the script to run the following variables must be provided with the exact column names shown below ( The column order is irrelevant). Longitude and Latitude in Decimal Degrees - WGS 84, ions in mg/L and isotopic data in per mil.
| Sample | Latitude | Longitude | Cl | HCO3 | SO4 | Ca | Mg | Na | K |
|---|
Other variables that can be included and are by default incorporated in the script are:
| d2H | d18O | D_Excess | NO3 | Br | CO3 |
|---|
Deuterium, d18O and Deuterium excess should be provided together!
- Find the r Provide the instructions for your run chunk
-
Provide the path to the .xlsx containing your data to the
water_analysis <- read_excel( ''C/.../your_data.xlsx'')
-
Choose the clustering parameters and the number of clusters
-
Specify EPSG Coordinate Reference System.
- Knit!